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))