657
edits
Geometry package: Difference between revisions
Jump to navigation
Jump to search
→Meshing Octave
| Line 35: | Line 35: | ||
n = 1e6; | n = 1e6; | ||
thrsd = 0.38; | thrsd = 0.38; | ||
while n>size (Ps,1) | while n > size (Ps,1) | ||
n = size (Ps,1); | n = size (Ps,1); | ||
ds = sqrt (sumsq (Ps(1:end,:) - Ps([2:end 1],:),2)); | ds = sqrt (sumsq (Ps(1:end,:) - Ps([2:end 1],:),2)); | ||
Ps = Ps(ds>mean (ds)*(1-thrsd),:); | Ps = Ps(ds > mean (ds)*(1-thrsd),:); | ||
end | end | ||
whos P Ps | whos P Ps | ||
| Line 54: | Line 54: | ||
<!-- </syntaxhighlight> --> | <!-- </syntaxhighlight> --> | ||
To plot the generated mesh we use the function {{Codeline|pdemesh}} from the {{Forge|fpl}} package. | After this code finishes, the structure T contains our mesh. To plot the generated mesh we use the function {{Codeline|pdemesh}} from the {{Forge|fpl}} package. In general is a good idea to use the openGL render (called ''fltk'') to plot meshes. | ||
{{Code|Plotting mesh with fpl package|<pre> | {{Code|Plotting mesh with fpl package|<pre> | ||
pkg load fpl | pkg load fpl | ||
pdemesh (T.p, T.e, T.t) | graphics_toolkit('fltk') | ||
pdemesh (T.p, T.e, T.t); | |||
view (2) | |||
axis tight | |||
axis equal | |||
</pre>}} | </pre>}} | ||
The output should look something like this | |||
[[File:octave_meshed.png|600x400px|center]] | [[File:octave_meshed.png|600x400px|center]] | ||
Questions? Ideas? Join us in the mailing list of in the IRC channel. | |||
[[File:octave_meshed.png|200px|center]] | |||
=== From piece-wise polynomial shapes to polygons === | === From piece-wise polynomial shapes to polygons === | ||