Geometry package: Difference between revisions

From Octave
Jump to navigation Jump to search
(using templates and octave guidelines for code formatting)
(using the Forge and Codeline templates)
Line 13: Line 13:


This tutorial shows the workflow to generate a triangular mesh inside an arbitrary region.
This tutorial shows the workflow to generate a triangular mesh inside an arbitrary region.
This tutorial requires that you install the package [http://octave.sourceforge.net/fpl/index.html fpl] and [http://octave.sourceforge.net/msh/index.html msh] (which requires [http://geuz.org/gmsh/ Gmsh] installed in your system). Alternatively, the core function <tt>delaunay</tt> could be used (the tutorial explains how) but the result aren't so pretty.
This tutorial requires that you install the package {{Forge|fpl}} and {{Forge|msh}} (which requires [http://geuz.org/gmsh/ Gmsh] installed in your system). Alternatively, the core function {{Codeline|delaunay}} could be used (the tutorial explains how) but the result aren't so pretty.


The first part of the tutorial requires an interesting shape. If you have Inkscape you can use the previous tutorial to load it into octave. Here I will be using [http://ubuntuone.com/5pNS12ZChUXeGNBniWNa3J this SVG].  
The first part of the tutorial requires an interesting shape. If you have Inkscape you can use the previous tutorial to load it into octave. Here I will be using [http://ubuntuone.com/5pNS12ZChUXeGNBniWNa3J this SVG].  
Line 25: Line 25:
P        = bsxfun (@minus, P, centroid (P));
P        = bsxfun (@minus, P, centroid (P));
</pre>}}
</pre>}}
Now we have our SVG as a polygon compatible with the Geometry package format. You can plot the polygon using the function <tt>drawPolygon</tt>. the next step is to mesh the interior of the polygon. To do this we could just call <tt>delaunay</tt> on the polygon and be done with it, but in general such mesh wont be so nice (you will need to add interior points). A very effective way of generating a good mesh is to use the package [http://octave.sourceforge.net/msh/index.html msh], which requires [http://geuz.org/gmsh/ Gmsh] installed in your system. The function data2geo in the Geometry package makes our work very easy:
Now we have our SVG as a polygon compatible with the Geometry package format. You can plot the polygon using the function {{Codeline|drawPolygon}}. the next step is to mesh the interior of the polygon. To do this we could just call {{Codeline|delaunay}} on the polygon and be done with it, but in general such mesh wont be so nice (you will need to add interior points). A very effective way of generating a good mesh is to use the package {{Forge|msh}}, which requires [http://geuz.org/gmsh/ Gmsh] installed in your system. The function {{Codeline|data2geo}} in the Geometry package makes our work very easy:


{{Code|Generating mesh for plot with msh package|<pre>
{{Code|Generating mesh for plot with msh package|<pre>
Line 36: Line 36:


<!-- </syntaxhighlight> -->
<!-- </syntaxhighlight> -->
To plot the generated mesh we use the function <tt>pdemesh</tt> from the [http://octave.sourceforge.net/fpl/index.html fpl] package.
To plot the generated mesh we use the function {{Codeline|pdemesh}} from the {{Forge|fpl}} package.


{{Code|Plotting mesh with fpl package|<pre>
{{Code|Plotting mesh with fpl package|<pre>

Revision as of 12:17, 8 March 2012

The geometry package is multipackage providing functions to manipulate geometrical entities in 2D and 3D.

Relation to matGeom

octCLIP

Piece-wise 2D polynomial polygons

Tutorials

Loading SVG files

Meshing Octave

Under construction

This tutorial shows the workflow to generate a triangular mesh inside an arbitrary region. This tutorial requires that you install the package fpl and msh (which requires Gmsh installed in your system). Alternatively, the core function delaunay could be used (the tutorial explains how) but the result aren't so pretty.

The first part of the tutorial requires an interesting shape. If you have Inkscape you can use the previous tutorial to load it into octave. Here I will be using this SVG.

Octave.png
Code: Loading the file as polygon compatible with geometry package
octavesvg = svg ("octave.svg");
ids       = octavesvg.pathid();
P         = octavesvg.path2polygon (ids{1}, 12)(1:end-1,:);
P         = bsxfun (@minus, P, centroid (P));

Now we have our SVG as a polygon compatible with the Geometry package format. You can plot the polygon using the function drawPolygon. the next step is to mesh the interior of the polygon. To do this we could just call delaunay on the polygon and be done with it, but in general such mesh wont be so nice (you will need to add interior points). A very effective way of generating a good mesh is to use the package msh, which requires Gmsh installed in your system. The function data2geo in the Geometry package makes our work very easy:

Code: Generating mesh for plot with msh package
pkg load msh
filename = tmpnam ();
meshsize = sqrt (mean (sumsq (diff (P, 1, 1), 2)))/2;
data2geo (P, meshsize, "output", [filename ".geo"]);
T = msh2m_gmsh (filename);

To plot the generated mesh we use the function pdemesh from the fpl package.

Code: Plotting mesh with fpl package
pkg load fpl 
pdemesh (T.p, T.e, T.t)
Octave meshed.png

From piece-wise polynomial shapes to polygons

See also