def testDerivatives(self):
k = 10
# meshes = [tum.regularsquaremesh(2)]
meshes = tum.meshes2d()
structuredpoints = pug.StructuredPoints([[0.01,0.01],[0.99,0.99]], [20,30])
rules = [pcbu.planeWaveBases(2, k),
pcbu.FourierBesselBasisRule(range(-4,5)),
pcbu.FourierHankelBasisRule([[-1,-1]], range(-4,5)),
pcbu.ProductBasisRule(pcbu.planeWaveBases(2, k, 3),pcbr.ReferenceBasisRule(pcbr.Dubiner(2)))]
basisDerivatives(rules, meshes, structuredpoints, k)
def testPToE2D(self):
self.compareStructuredAndUnstructured(tum.meshes2d(), 2)
def test2DQuadratures(self):
""" Ensure quadratures on neighbouring faces of various 2D meshes match """
self.compareQuadratures(tum.meshes2d(), puq.legendrequadrature(6))
