def __init__(self, problem, basisrule, nquadpoints, usecache = False):
self.problem = problem
self.basis = psp.constructBasis(problem, basisrule)
fquad, equad = puq.quadrules(problem.mesh.dim, nquadpoints)
self.facequads = pmmu.MeshQuadratures(problem.mesh, fquad)
self.elementquads = pmmu.MeshElementQuadratures(problem.mesh, equad)
self.facevandermondes = self.faceVandermondes(pcbu.FaceToBasis(problem.mesh, self.basis), usecache=usecache)
self.elementvandermondes = pcv.ElementVandermondes(problem.mesh, self.basis, self.elementquads)
def __init__(self, mesh, quadpoints, etods, k, nfb):
self.k = k
self.mesh = mesh
self.etods = etods
self.nfb= nfb
_, equad = puq.quadrules(mesh.dim, quadpoints)
self.mqs = pmmu.MeshElementQuadratures(mesh, equad)
self.fbrule = pcb.FourierBesselBasisRule(np.arange(0,self.nfb) - self.nfb/2)
def __init__(self, mesh, quadpoints, etob, ibc):
self.mesh = mesh
fquad, _ = puq.quadrules(mesh.dim, quadpoints)
self.mqs = pmmu.MeshQuadratures(mesh, fquad)
self.etobc = dict([(e, ibc(e)) for e in mesh.partition])
self.etonbcs = {}
self.etob = etob
self.populate()
def __init__(self, problem, quadpoints, basis):
_, self.equad = puq.quadrules(problem.mesh.dim, quadpoints)
self.basis = basis
self.problem = problem
