proteus.DiagUtils module

Module for diagnostic utilities

proteus.DiagUtils.L2errorFEMvsAF(analyticalFunction, quadraturePointArray, quadratureWeightArray, functionValueArray, T=None)

supposed to be L2 norm of error in vector quantity I think just using dot would cover both scalar and vector case

proteus.DiagUtils.getQuadraturePhysPointsAndWeights(mesh, femSpace, quadrature, verbose=0)

for convenience, hide steps for generating quadrature points and weights, with Jacobians, on physical mesh based on points and weights on reference element

returns points array that’s nelem x nquadloc x 3

weight array that’s nelem x nquadloc

proteus.DiagUtils.getFEMvals(u, xiArray, verbose=0)

for convenience, hide steps for generating finite element solution at physical points corresponding to reference points held in xiArray

returns array that’s nelem x npointloc

proteus.DiagUtils.testCrRavNodalBasis(nd, verbose=0)

test local Crouzeix-Raviart element space

proteus.DiagUtils.testQuadNodalBasis(nd, verbose=0)

test local P2 nodal finite element space

proteus.DiagUtils.testEdgeDOFMap(mesh, nd)

test edge dof map to see what its doing

proteus.DiagUtils.testQuadRefMats(nd, verbose=0)

test quad reference matrices to see what its doing

proteus.DiagUtils.testQuadDOFMap(mesh, nd, verbose=0)

test quad dof map to see what its doing