from fealpy.recovery import FEMFunctionRecoveryAlg
from fealpy.mesh.adaptive_tools import mark
from fealpy.quadrature import TriangleQuadrature
from mpl_toolkits.mplot3d import Axes3D
from fealpy.tools.show import showmultirate
m = int(sys.argv[1])
maxit = int(sys.argv[2])
p = int(sys.argv[3])
if m == 1:
pde = LShapeRSinData()
mesh = pde.init_mesh(n=4, meshtype='tri')
elif m == 2:
pde = KelloggData() #TODO
mesh = pde.init_mesh(n=4, meshtype='tri')
elif m == 3:
pde = CrackData()
mesh = pde.init_mesh(n=4, meshtype='tri')
elif m == 4:
pde = SinSinData()
mesh = pde.init_mesh(n=1, meshtype='tri')
elif m == 5:
pde = CosCosData()
mesh = pde.init_mesh(n=2, meshtype='tri')
theta = 0.1
k = maxit - 15
errorType = [
'$|| u_I - u_h ||_{l_2}$',
from fealpy.recovery import FEMFunctionRecoveryAlg
from fealpy.mesh.adaptive_tools import mark
from fealpy.quadrature import TriangleQuadrature
from mpl_toolkits.mplot3d import Axes3D
from fealpy.tools.show import showmultirate
m = int(sys.argv[1])
maxit = int(sys.argv[2])
p = int(sys.argv[3])
if m == 1:
pde = LShapeRSinData()
mesh = pde.init_mesh(n=4, meshtype='tri')
elif m == 2:
pde = KelloggData()
mesh = pde.init_mesh(n=4, meshtype='tri')
elif m == 3:
pde = CrackData()
mesh = pde.init_mesh(n=4, meshtype='tri')
elif m == 4:
pde = SinSinData()
mesh = pde.init_mesh(n=1, meshtype='tri')
elif m == 5:
pde = CosCosData()
mesh = pde.init_mesh(n=2, meshtype='tri')
theta = 0.1
k = maxit - 15
errorType = [
'$|| u_I - u_h ||_{l_2}$', '$|| u - u_h||_{0}$',