Re: Distmesh, mesh refinement

Andreas Stahel <[email protected]>
Newsgroups gmane.comp.gnu.octave.general
Message-ID <[email protected]>
Mishal0488 wrote
> Hi Guys
> 
> Does anyone have any experience with the use of Distmesh? I have attached
> a
> link below.
> 
> http://persson.berkeley.edu/distmesh/
> 
> Below is my code for a mesh. I am not sure to refine the mesh to attain
> better results for my finite element code. The 0.01 which is in bold
> allows
> the mesh to get finer when compared to a value of 0.1, however anything
> smaller than 0.01 does not create a better mesh.
> 
> If anyone has any experience in this package please help.
> 
> P1 = [0 0];
> P2 = [0 60];
> P3 = [30 30]; #load point
> P4 = [60 0];
> pv = [P1(1,1) P1(1,2);P2(1,1) P2(1,2);P3(1,1) P3(1,2);P4(1,1)
> P4(1,2);P1(1,1) P1(1,2)];
> fd = { 'l_dpolygon', [], pv };
> fh = @(p) (min(p(:,3)));#ones(size(p,1),1);
> [p,T] = distmesh( fd, fh, *0.01*, [-2,-2;2,2],pv);
> patch( 'vertices', p, 'faces', T, 'facecolor', [.9, .9, .9] )
> 
> 
> 
> --
> Sent from:
> https://octave.1599824.n4.nabble.com/Octave-General-f1599825.html

Dear Mishal

With distmesh I have no experience, but will have a look soon.
For my FEM code I used triangle, which woked without major problems.
If you like give FEMoctave a try. Install by
pkg install https://github.com/AndreasStahel/FEMoctave/archive/v2.0.3.tar.gz
or use github at https://github.com/AndreasStahel/FEMoctave

I hope it helps

Andreas






--
Sent from: https://octave.1599824.n4.nabble.com/Octave-General-f1599825.html
lmpx.com only provides a reader for public news (NNTP) servers. It is not affiliated with the servers or forums shown here and is not responsible for the content of articles, which is written by their respective authors.