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}$',
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.pde.poisson_model_2d import CrackData, LShapeRSinData, CosCosData, KelloggData, SinSinData, ffData
from fealpy.vem import PoissonVEMModel
from fealpy.mesh.adaptive_tools import AdaptiveMarker
from fealpy.tools.show import showmultirate
from fealpy.quadrature import TriangleQuadrature
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
m = int(sys.argv[1])
maxit = int(sys.argv[2])
p = int(sys.argv[3])
if m == 1:
pde = KelloggData()
quadtree = pde.init_mesh(n=4)
elif m == 2:
pde = LShapeRSinData()
quadtree = pde.init_mesh(n=4)
elif m == 3:
pde = CrackData()
quadtree = pde.init_mesh(n=4)
elif m == 4:
pde = CosCosData()
quadtree = pde.init_mesh(n=2)
elif m == 5:
pde = SinSinData()
quadtree = pde.init_mesh(n=3)
elif m == 6:
pde = ffData()