def functionSpace(space="Pch", mesh=mesh.mesh(2, 1), order=1):
"""create a function space
"""
key = space + '(' + str(mesh.dimension()) + ',' + str(order) + ')'
if key not in _spaces:
raise RuntimeError('FunctionSpace ' + key +
' not existing in dictionary')
return _spaces[key](mesh)
import sys
from pyfeelpp import core
e=core.Environment(sys.argv)
from pyfeelpp import mesh,discr,filters,vf
geo={
'2':core.download( "github:{repo:feelpp,path:feelpp/quickstart/laplacian/feelpp2d/feelpp2d.geo}", worldComm=core.Environment.worldCommPtr() )[0],
'3':core.download( "github:{repo:feelpp,path:feelpp/quickstart/laplacian/feelpp3d/feelpp3d.geo}", worldComm=core.Environment.worldCommPtr() )[0]
}
def run( m, geo ):
m2d = mesh.load(m,geo,0.1)
Xh=discr.functionSpace( space="Pch", mesh=m2d, order=1 )
P0h = discr.functionSpace( space="Pdh", mesh=m2d, order=0 )
#u=Xh.elementFromExpr("{sin(2*pi*x)*cos(pi*y)}:x:y")
u=Xh.element()
u.on(range=mesh.elements(m2d),expr=vf.expr("x*x:x"))
e = filters.exporter(mesh=m2d,name="feelpp"+str(m.dimension())+"d")
e.addScalar("un", 1.)
e.addP1c("u",u);
e.addP0d("pid",discr.pid( P0h ));
e.save()
run( mesh.mesh( dim=2 ), geo['2'] )
run( mesh.mesh( dim=3, realdim=3 ), geo['3'] )