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()