Skip to content

Commit

Permalink
update mesh doc
Browse files Browse the repository at this point in the history
  • Loading branch information
lyc102 committed Aug 27, 2021
1 parent 842ba40 commit e8ce3cf
Show file tree
Hide file tree
Showing 29 changed files with 149 additions and 183 deletions.
15 changes: 6 additions & 9 deletions docs/_mesh/basic.md
Original file line number Diff line number Diff line change
Expand Up @@ -35,9 +35,8 @@ findnode(node); % plot indices of all vertices
```


![png]({{ site.baseurl }}/assets/images/mesh/meshbasicdoc_2_0.png)
<img src="mesh_figures/meshbasicdoc_2_0.png" alt="meshbasicdoc_2_0" style="zoom:75%;" />



Apply the uniform refinement sevreal times to obtain a fine mesh.
Expand All @@ -51,9 +50,8 @@ showmesh(node,elem);
```


![png]({{ site.baseurl }}/assets/images/mesh/meshbasicdoc_4_0.png)
<img src="mesh_figures/meshbasicdoc_4_0.png" alt="meshbasicdoc_4_0" style="zoom:75%;" />



## Example: A Cube in 3-D
Expand All @@ -70,9 +68,8 @@ findnode3(node);
```


![png]({{ site.baseurl }}/assets/images/mesh/meshbasicdoc_6_0.png)
<img src="mesh_figures/meshbasicdoc_6_0.png" alt="meshbasicdoc_6_0" style="zoom:75%;" />



Unlike the 2-D case, to apply uniform refinement to obtain a fine mesh with good mesh quality, a different ordering of the initial mesh, which may violate the positive ordering, should be used. See [3-D Red Refinement]({{ site.baseurl }}{% link _mesh/uniformrefine3.md %}).
File renamed without changes
File renamed without changes
File renamed without changes
File renamed without changes
File renamed without changes
File renamed without changes
File renamed without changes
File renamed without changes
File renamed without changes
File renamed without changes
File renamed without changes
File renamed without changes
File renamed without changes
File renamed without changes
File renamed without changes
311 changes: 140 additions & 171 deletions docs/_mesh/sc.md

Large diffs are not rendered by default.

2 changes: 1 addition & 1 deletion example/solver/compareMaxwellsolvers3.m
Original file line number Diff line number Diff line change
Expand Up @@ -23,7 +23,7 @@
%% Lshape adaptive grids
mesh.shape = 'Lshape';
mesh.type = 'adaptive';
mesh.size = 4e4;
mesh.size = 3e4;
pde = 'Maxwell';
fem = 'ND1';
% get the matrix
Expand Down
4 changes: 2 additions & 2 deletions solver/amgMaxwell.m
Original file line number Diff line number Diff line change
Expand Up @@ -134,8 +134,8 @@
bdidx(isBdNode) = 1;
end
Tbd = spdiags(bdidx,0,N,N);
% BP = gradt*A(1:NE,1:NE)*grad + Tbd;
BP = gradt*spdiags(edgeLength.*beta,0,NE,NE)*grad + Tbd;
BP = gradt*A(1:NE,1:NE)*grad + Tbd;
% BP = gradt*spdiags(edgeLength.*beta,0,NE,NE)*grad + Tbd;

%% Transfer operators between multilevel meshes
setupOption.solver = 'NO';
Expand Down

0 comments on commit e8ce3cf

Please sign in to comment.