OFELI's Logo

An Object Oriented Finite Element Library

To adapt mesh in function of given solution. More...

#include <MeshAdapt.h>

Public Member Functions

 MeshAdapt ()
 Default constructor.
 
 MeshAdapt (Mesh &ms)
 Constructor using initial mesh. More...
 
 MeshAdapt (Domain &dom)
 Constructor using a reference to class Domain. More...
 
 ~MeshAdapt ()
 Destructor.
 
DomaingetDomain () const
 Get reference to Domain instance.
 
MeshgetMesh () const
 Get reference to current mesh.
 
void set (Domain &dom)
 Set reference to Domain instance.
 
void set (Mesh &ms)
 Set reference to Mesh instance.
 
void setSolution (const Vect< real_t > &u)
 Define label of node.
 
void setJacobi (int n)
 Set number of Jacobi iterations for smoothing.
 
void setSmooth (int n)
 Set number of smoothing iterations.
 
void AbsoluteError ()
 Metric is constructed with absolute error.
 
void RelativeError ()
 Metric is constructed with relative error.
 
void setError (real_t err)
 Set error threshold for adaption.
 
void setHMin (real_t h)
 Set minimal mesh size.
 
void setHMax (real_t h)
 Set maximal mesh size.
 
void setHMinAnisotropy (real_t h)
 Set minimal mesh size and set anisotropy.
 
void setRelaxation (real_t omega)
 Set relaxation parameter for smoothing. More...
 
void setAnisotropic ()
 Set that adapted mesh construction is anisotropic.
 
void MaxAnisotropy (real_t a)
 Set maximum ratio of anisotropy.
 
void setMaxSubdiv (real_t s)
 Change the metric such that the maximal subdivision of a background's edge is bounded by the given number (always limited by 10)
 
void setMaxNbVertices (size_t n)
 Set maximum number of vertices. More...
 
void setRatio (real_t r)
 Set ratio for a smoothing of the metric. More...
 
void setNoScaling ()
 Do not scale solution before metric computation. More...
 
void setNoKeep ()
 Do not keep old vertices. More...
 
void setHessian ()
 set computation of the Hessian
 
void setOutputMesh (string file)
 Create mesh output file.
 
void setGeoFile (string file)
 Set Geometry file.
 
void setGeoError (real_t e)
 Set error on geometry.
 
void setBackgroundMesh (string bgm)
 Set background mesh.
 
void SplitBoundaryEdges ()
 Split edges with two vertices on boundary.
 
void CreateMetricFile (string mf)
 Create a metric file.
 
void setMetricFile (string mf)
 Set Metric file.
 
void getSolutionMbb (string mbb)
 Set solution defined on background mesh for metric construction.
 
void getSolutionMBB (string mBB)
 Set solution defined on background mesh for metric construction.
 
void getSolutionbb (string rbb)
 Read solution defined on the background mesh in bb file. More...
 
void getSolutionBB (string rBB)
 Read solution defined on the background mesh in BB file. More...
 
void getSolution (Vect< real_t > &u, int is=1)
 Get the interpolated solution on the new mesh. More...
 
void getInterpolatedSolutionbb ()
 Write the file of interpolation of the solutions in bb file.
 
void getInterpolatedSolutionBB ()
 Write the file of interpolation of the solutions in BB file.
 
void setTheta (real_t theta)
 Set angular limit for a corner (in degrees) More...
 
void Split ()
 Split triangles into 4 triangles.
 
void saveMbb (string file, const Vect< real_t > &u)
 Save a solution in metric file. More...
 
int run ()
 Run adaptation process. More...
 
int run (const Vect< real_t > &u)
 Run adaptation process using a solution vector. More...
 
int run (const Vect< real_t > &u, Vect< real_t > &v)
 Run adaptation process using a solution vector and interpolates solution on the adapted mesh. More...
 

Detailed Description

To adapt mesh in function of given solution.

Class MeshAdapt enables modifying mesh according to a solution vector defining at nodes. It concerns 2-D triangular meshes only.

Remarks
Class MeshAdapt is mainly based on the software 'Bamg' developed by F. Hecht, Universite Pierre et Marie Curie, Paris. We warmly thank him for accepting incoporation of Bamg in the OFELI package
Author
Rachid Touzani

Constructor & Destructor Documentation

◆ MeshAdapt() [1/2]

MeshAdapt ( Mesh ms)

Constructor using initial mesh.

Parameters
[in]msReference to initial mesh

◆ MeshAdapt() [2/2]

MeshAdapt ( Domain dom)

Constructor using a reference to class Domain.

Parameters
[in]domReference to Domain class

Member Function Documentation

◆ getSolution()

void getSolution ( Vect< real_t > &  u,
int  is = 1 
)

Get the interpolated solution on the new mesh.

The solution must have been saved on an output bb file

Parameters
[out]uVector that contains on output the obtained solutions. This vector is resized before being initialized
[in]is[Default: 1]

◆ getSolutionbb()

void getSolutionbb ( string  rbb)

Read solution defined on the background mesh in bb file.

Solution is interpolated on created mesh

◆ getSolutionBB()

void getSolutionBB ( string  rBB)

Read solution defined on the background mesh in BB file.

Solution is interpolated on created mesh

◆ run() [1/3]

int run ( )

Run adaptation process.

Returns
Return code:
  • = 0: Adaptation has been normally completed
  • = 1: An error occured

◆ run() [2/3]

int run ( const Vect< real_t > &  u)

Run adaptation process using a solution vector.

Parameters
[in]uSolution vector defined on the input mesh
Returns
Return code:
  • = 0: Adaptation has been normally completed
  • = 1: An error occured

◆ run() [3/3]

int run ( const Vect< real_t > &  u,
Vect< real_t > &  v 
)

Run adaptation process using a solution vector and interpolates solution on the adapted mesh.

Parameters
[in]uSolution vector defined on the input mesh
[in]vSolution vector defined on the (adapted) output mesh
Returns
Return code:
  • = 0: Adaptation has been normally completed
  • = 1: An error occured

◆ saveMbb()

void saveMbb ( string  file,
const Vect< real_t > &  u 
)

Save a solution in metric file.

Parameters
[in]fileFile name where the metric is stored
[in]uSolution vector to store

◆ setMaxNbVertices()

void setMaxNbVertices ( size_t  n)

Set maximum number of vertices.

Default value is 500000

◆ setNoKeep()

void setNoKeep ( )

Do not keep old vertices.

By default, old vertices are kept

◆ setNoScaling()

void setNoScaling ( )

Do not scale solution before metric computation.

By default, solution is scaled (between 0 and 1)

◆ setRatio()

void setRatio ( real_t  r)

Set ratio for a smoothing of the metric.

Parameters
[in]rRatio value.
Note
If r is 0 then no smoothing is performed, if r lies in [1.1,10] then the smoothing changes the metric such that the largest geometrical progression (speed of mesh size variation in mesh is bounded by r) (by default no smoothing)

◆ setRelaxation()

void setRelaxation ( real_t  omega)

Set relaxation parameter for smoothing.

Default value for relaxation parameter is 1.8

◆ setTheta()

void setTheta ( real_t  theta)

Set angular limit for a corner (in degrees)

The angle is defined from 2 normals of 2 consecutive edges






Copyright © 1998-2022 Rachid Touzani