#include "Basic/Utilities.hpp"
#include "Db/Db.hpp"
#include "Db/DbGrid.hpp"
#include "Mesh/AMesh.hpp"
#include "Mesh/MeshEStandard.hpp"
#include "Mesh/Delaunay.hpp"
#include "Core/io.hpp"
#include <math.h>
#include <stdio.h>

Functions
int	MSS (int ndim, int ipol, int icas, int icorn, int idim)

VectorDouble	extend_grid (DbGrid db, const VectorDouble &gext, int nout)

VectorDouble	extend_point (Db db, const VectorDouble &gext, int nout)

VectorDouble	get_db_extension (Db dbin, Db dbout, int *nout)

static int	st_load_segment (DbGrid *dbgrid, VectorInt &mesh, VectorInt &order, int ipos, int ix1, int ix2)

static int	st_load_triangle (DbGrid *dbgrid, VectorInt &mesh, VectorInt &order, int ipos, int ix1, int iy1, int ix2, int iy2, int ix3, int iy3)

static int	st_load_tetra (DbGrid *dbgrid, VectorInt &mesh, VectorInt &order, int ipos, int ix1, int iy1, int iz1, int ix2, int iy2, int iz2, int ix3, int iy3, int iz3, int ix4, int iy4, int iz4)

static MeshEStandard *	st_ultimate_regular_grid (Db *dbgrid, int ndim, int nmesh, int ncorner, VectorInt &meshes, VectorInt &order)

AMesh *	meshes_turbo_2D_grid_build (DbGrid *dbgrid)

int	meshes_2D_write (const char file_name, const char obj_name, int verbose, int ndim, int ncode, int ntri, int npoints, const VectorInt &ntcode, const VectorInt &triangles, const VectorDouble &points)

void	mesh_stats (int ndim, int ncorner, int nmesh, const int meshes, const double points)

AMesh *	meshes_turbo_3D_grid_build (DbGrid *dbgrid)

AMesh *	meshes_turbo_1D_grid_build (DbGrid *dbgrid)

Function Documentation

◆ extend_grid()

VectorDouble extend_grid	(	DbGrid *	db,
		const VectorDouble &	gext,
		int *	nout
	)

Extend a Grid by gext

Returns: The coordinates of the extension points (Dimension: number * ndim)

Parameters

[in]	db	Output Db grid structure
[in]	gext	Array of domain dilation
[out]	nout	Number of extension points

Remarks: The returned array 'ext' must be freed by the calling function

◆ extend_point()

VectorDouble extend_point	(	Db *	db,
		const VectorDouble &	gext,
		int *	nout
	)

Extend a Point Domain by gext

Returns: The coordinates of the extension points (Dimension: number * ndim)

Parameters

[in]	db	Output Db grid structure
[in]	gext	Array of domain dilation
[out]	nout	Number of extension points

Remarks: The returned array 'ext' must be freed by the calling function

◆ get_db_extension()

VectorDouble get_db_extension	(	Db *	dbin,
		Db *	dbout,
		int *	nout
	)

Define the maximum extension between dbin and/or dbout

Returns: Array of extensions

Parameters

[in]	dbin	Db input (optional)
[in]	dbout	Db output (optional)
[out]	nout	Number of extension points

Remarks: The calling function must free the returned array 'ext'; with its dimension: ndim * 2^ndim

◆ mesh_stats()

void mesh_stats	(	int	ndim,
		int	ncorner,
		int	nmesh,
		const int *	meshes,
		const double *	points
	)

Calculate the statistics on a set of meshes

Parameters

[in]	ndim	Space dimension
[in]	ncorner	Number of corners for each mesh
[in]	nmesh	Number of meshes
[in]	meshes	Array of vertex indices for each mesh
[in]	points	Array of 'ndim' coordinates for mesh vertex

◆ meshes_2D_write()

int meshes_2D_write	(	const char *	file_name,
		const char *	obj_name,
		int	verbose,
		int	ndim,
		int	ncode,
		int	ntri,
		int	npoints,
		const VectorInt &	ntcode,
		const VectorInt &	triangles,
		const VectorDouble &	points
	)

Dump the contents of a triangulation in an ASCII file according to the STL format

Parameters

[in]	file_name	Nmae of the created ASCII file
[in]	obj_name	Name assigned to the object
[in]	verbose	Verbose flag
[in]	ndim	Space dimension
[in]	ncode	Number of different codes
[in]	ntri	Number of triangles (expected)
[in]	npoints	Number of poins (expected)
[in]	ntcode	Array of number of triangles per code
[in]	triangles	Array of vertex indices for each triangle
[in]	points	Array of 3-D coordinates for triangle vertices

◆ meshes_turbo_1D_grid_build()

AMesh* meshes_turbo_1D_grid_build ( DbGrid * dbgrid )

Build the regular meshing from a 1-D grid

Returns: The newly created AMesh structure

Parameters

[in] dbgrid Db structure

◆ meshes_turbo_2D_grid_build()

AMesh* meshes_turbo_2D_grid_build ( DbGrid * dbgrid )

Build the regular 2-D grid meshing

Returns: Pointer to the newly created AMesh structure

Parameters

[in] dbgrid Db structure

◆ meshes_turbo_3D_grid_build()

AMesh* meshes_turbo_3D_grid_build ( DbGrid * dbgrid )

Build the regular 3-D grid meshing

Returns: Pointer to the newly created AMesh structure

Parameters

[in] dbgrid Db structure

◆ MSS()

int MSS	(	int	ndim,
		int	ipol,
		int	icas,
		int	icorn,
		int	idim
	)

Returns the shift value for the apex('icorn')

Returns: Shift value

Parameters

[in]	ndim	Space dimension (1, 2 or 3)
[in]	ipol	Rank of the polarization (starting from 0)
[in]	icas	Rank of the case (starting from 0)
[in]	icorn	Rank of the corner (starting from 0)
[in]	idim	Rank of the coordinate (starting from 0)

Remarks: This function returns the shift value to be applied to each; coordinate index, for each apex for each mesh which constitutes; the partition of a cell into meshes (using Turbo facility); When 'ndim' is provided as negative, a special case is programmed

◆ st_load_segment()

static int st_load_segment	(	DbGrid *	dbgrid,
		VectorInt &	mesh,
		VectorInt &	order,
		int	ipos,
		int	ix1,
		int	ix2
	)

static

Load the vertices in a segment and check if the segment is masked

Returns: 1 If the segment is valid because at least one vertex is active

Parameters

[in]	dbgrid	Db structure
[in]	ipos	Position of newly created mesh information
[in]	ix1	Grid index along X for the vertex #1
[in]	ix2	Grid index along X for the vertex #2
[out]	mesh	Array of triangle ranks (dimension = 3)
[out]	order	Array of relative ranks

Remarks: The values in 'order' are the absolute indices (starting from 1),; negative if the grid node is masked off

◆ st_load_tetra()

static int st_load_tetra	(	DbGrid *	dbgrid,
		VectorInt &	mesh,
		VectorInt &	order,
		int	ipos,
		int	ix1,
		int	iy1,
		int	iz1,
		int	ix2,
		int	iy2,
		int	iz2,
		int	ix3,
		int	iy3,
		int	iz3,
		int	ix4,
		int	iy4,
		int	iz4
	)

static

Load the vertices in a tetrahedron and check if the tetrahedron is masked

Returns: 1 If the tetrahedron is valid because at least one vertex is active

Parameters

[in]	dbgrid	Db structure
[in]	mesh	Array of triangle ranks (dimension = 4)
[in]	order	Array of relative ranks
[in]	ipos	Position of newly created mesh information
[in]	ix1	Grid index along X for the vertex #1
[in]	iy1	Grid index along Y for the vertex #1
[in]	iz1	Grid index along Z for the vertex #1
[in]	ix2	Grid index along X for the vertex #2
[in]	iy2	Grid index along Y for the vertex #2
[in]	iz2	Grid index along Z for the vertex #2
[in]	ix3	Grid index along X for the vertex #3
[in]	iy3	Grid index along Y for the vertex #3
[in]	iz3	Grid index along Z for the vertex #3
[in]	ix4	Grid index along X for the vertex #4
[in]	iy4	Grid index along Y for the vertex #4
[in]	iz4	Grid index along Z for the vertex #4

Remarks: The values in 'order' are the absolute indices (starting from 1),; negative if the grid node is masked off

◆ st_load_triangle()

static int st_load_triangle	(	DbGrid *	dbgrid,
		VectorInt &	mesh,
		VectorInt &	order,
		int	ipos,
		int	ix1,
		int	iy1,
		int	ix2,
		int	iy2,
		int	ix3,
		int	iy3
	)

static

Load the vertices in a triangle and check if the triangle is masked

Returns: 1 If the triangle is valid because at least one vertex is active

Parameters

[in]	dbgrid	Db structure
[in]	ipos	Position of newly created mesh information
[in]	ix1	Grid index along X for the vertex #1
[in]	iy1	Grid index along Y for the vertex #1
[in]	ix2	Grid index along X for the vertex #2
[in]	iy2	Grid index along Y for the vertex #2
[in]	ix3	Grid index along X for the vertex #3
[in]	iy3	Grid index along Y for the vertex #3
[out]	mesh	Array of triangle ranks (dimension = 3)
[out]	order	Array of relative ranks

Remarks: The values in 'order' are the absolute indices (starting from 1),; negative if the grid node is masked off

◆ st_ultimate_regular_grid()

static MeshEStandard* st_ultimate_regular_grid	(	Db *	dbgrid,
		int	ndim,
		int	nmesh,
		int	ncorner,
		VectorInt &	meshes,
		VectorInt &	order
	)

static

Perform the ultimate task for the regular grid

Returns: Pointer to the newly created MeshEStandard structure

Parameters

[in]	dbgrid	Db structure
[in]	ndim	Space dimension
[in]	nmesh	Current number of meshes
[in]	ncorner	Number of corner per mesh
[in]	meshes	Initial array of meshes
[in]	order	Array of relative ranks

Functions

Function Documentation

◆ extend_grid()

◆ extend_point()

◆ get_db_extension()

◆ mesh_stats()

◆ meshes_2D_write()

◆ meshes_turbo_1D_grid_build()

◆ meshes_turbo_2D_grid_build()

◆ meshes_turbo_3D_grid_build()

◆ MSS()

◆ st_load_segment()

◆ st_load_tetra()

◆ st_load_triangle()

◆ st_ultimate_regular_grid()