External libraries
Curvature
extractborder
Extract a border of a mesh.
1 | int Res = extractborder(Th, Label, Points);
|
Parameters:
Th
(mesh
ormesh3
)Label
(int
) Label of the border to extractPoints
(real[int, int]
) Extracted points Must be allocated asreal[int, int] Points(3, 1);
Output:
Res
(real
) Length of the extracted border
dfft
Refer to the FFTW documentation for more informations.
distance
Need
1 | load "distance"
|
ff_gsl_awk
Refer to the GSL documentation for more informations
fflapack
Refer to the LAPACK documentation for more informations.
ff-NLopt
Refer to the NLOPT documentation for more informations.
gmsh
Need
1 | load "gsmh"
|
The gmsh
software is available here
gmshload
Load a 2D mesh build with Gmsh.
1 | mesh Th = gmshload(MeshFile, [reftri=RefTri], [renum=Renum]);
|
Parameters:
MeshFile
(string
) Mesh file namereftri=
(.. todo:: todo)renum=
(.. todo:: todo)
Output:
Th
(mesh
)
iovtk
savevtk
Save mesh or solution in vtk/vtu format.
1 | savevtk(FileName, Th, [Ux, Uy, Uz], p, [dataname=DataName], [withsurfacemesh=WithSurfaceMesh], [order=Order], [floatmesh=FloatMesh], [floatsol=FloatSol], [bin=Bin], [swap=Swap]);
|
Parameters:
FileName
(string
) File name:*.vtk
or*.vtu
Th
(mesh
ormesh3
)[Ux, Uy, Uz], p
(fespace
function of vector offespace
functions) Solutions to save, as much as wanteddataname=
(string
) Name of solutions, seprated by a spacewithsurfacemesh=
(bool
) .. todo:: todoorder=
(int[int]
) Order of solutions.Available: 0 or 1
floatmesh=
(bool
) .. todo:: todofloatsol=
(bool
) .. todo:: todobin=
(bool
) If true, save file in binary formatswap
(bool
) .. todo:: todo
Output:
- None
isoline
Need
1 | load "isoline"
|
msh3
cube
Construct a cubic mesh.
1 | mesh3 Th = cube(nnX, nnY, nnZ, [X(x), Y(y), Z(z)], [label=Label], [flags=Flags], [region=Region]);
|
Parameters:
nnX
(int
) Number of discretization point along \(x\)nnY
(int
) Number of discretization point along \(y\)nnZ
(int
) Number of discretization point along \(z\)X(x)
(func
) [Optional] Affine function of \(x\) to define the length Default:x
Y(y)
(func
) [Optional] Affine function of \(y\) to define the width Default:y
Z(z)
(func
) [Optional] Affine function of \(z\) to define the height Default:z
label=
(int[int]
) [Optional]List of surface labels Default:
[1, 2, 3, 4, 5, 6]
flags=
(int
) [Optional]Refer to square
region=
(int
) [Optional]Region number of the cube volume Default:
0
Output:
Th
(mesh3
) Cube mesh
tetgen
Refer to the Tetgen documentation for more informations.