test convergence on non-uniform meshes.

lyc102 · lyc102 · commit 5747bd7a5bf5 · 2021-12-03T17:29:12.000-08:00
diff --git a/example/fem/Poisson/cubePoisson.m b/example/fem/Poisson/cubePoisson.m
@@ -44,7 +44,7 @@
 
 %% Plot convergence rates
 figure;
-showrateh3(h,errH1,2,'-*', '$|| Du - Du_h ||', ...
+showrateh3(h,errH1,2,'-*', '|| Du - Du_h ||', ...
            h,errL2,2,'k-+', '|| u - u_h ||', ...
            h,erruIuh,2,'m-+','|| D u_I - D u_h ||');
 
diff --git a/mesh/cubegradmesh.m b/mesh/cubegradmesh.m
@@ -9,6 +9,11 @@
 %     [node,elem] = cubegradmesh(xh,yh,zh);
 %     showmesh3(node,elem);
 %
+%   xh = (0:0.125:1).^2;
+%   yh = (0:0.125:1).^2;
+%   zh = 0:0.125:0.5;
+%   [node,elem] = cubegradmesh(xh,yh,zh);
+%   showmesh3(node,elem);
 %
 % Copyright (C) Long Chen. See COPYRIGHT.txt for details. 
 

Original file line number	Diff line number	Diff line change
`@@ -9,6 +9,11 @@`
`9`	`9`	`% [node,elem] = cubegradmesh(xh,yh,zh);`
`10`	`10`	`% showmesh3(node,elem);`
`11`	`11`	`%`
	`12`	`+% xh = (0:0.125:1).^2;`
	`13`	`+% yh = (0:0.125:1).^2;`
	`14`	`+% zh = 0:0.125:0.5;`
	`15`	`+% [node,elem] = cubegradmesh(xh,yh,zh);`
	`16`	`+% showmesh3(node,elem);`
`12`	`17`	`%`
`13`	`18`	`% Copyright (C) Long Chen. See COPYRIGHT.txt for details.`
`14`	`19`