search

LutzGross / esys-escript.github.io

Other
29 stars 13 forks source link

finley: For higher order elements interpolation frm ReducedFunction to Function Fails. #35

Open aellery opened 2 years ago

aellery commented 2 years ago

from esys.escript import * from esys.finley import Rectangle

domain1=Rectangle(2,2,l0=50., l1=50.)

c=Scalar(1., ReducedFunction(domain1)) c.expand() # why is this needed? c=interpolate(c, Function(domain1)) print("c=" ,c)

domain2=Rectangle(2,2,l0=50., l1=50., order=2) c=Scalar(1., ReducedFunction(domain2)) c.expand() # why is this needed? c=interpolate(c, Function(domain2)) # Fails! print("c=" ,c)

domain3=Rectangle(2,2,l0=50., l1=50., order=2, fullOrder=True) c=Scalar(1., ReducedFunction(domain3)) c.expand() # why is this needed? c=interpolate(c, Function(domain3)) # Fails! print("c=" ,c)

LutzGross commented 2 years ago

This is also the case for FunctionOnBoundary!