Procedures

Return the solver state variable for a given set of elements

Coefficient alpha from the recursive formulation of Legendre polynomials, for the Legendre mode 'mode'.

Prodcut of alpha(numerator) * beta(denominator) / alpha(denominator) as needed by the Clenshaw algorithm in ply_split_legendre_matrix.

Quotient of two alpha values.

append a value to the dynamic array and return its position.

append a value to the dynamic array and return its position.

append a value to the dynamic array and return its position.

append a value to the dynamic array and return its position.

append a value to the dynamic array and return its position.

append a value to the dynamic array and return its position.

append a value to the dynamic array and return its position.

append a value to the dynamic array and return its position.

append a value to the dynamic array and return its position.

This routine convert plane to triangle and add triangle to growing array and add position of triangle into spatialObj%primitive_position

This routine convert box to planes and then planes to triangles and add triangle to growing array and add position of triangle into spatialObj%primitive_position

This routine convert canonical line to seeder line array

This routine convert plane to triangle and add triangle to growing array and add position of triangle into spatialObj%primitive_position

This routine convert canonical point to seeder point array

This routine single cylinder from object table

appending a value to the dynamic array

appending a value to the dynamic array

appending a sorted list of values to the dynamic array

appending a sorted list of values to the dynamic array

Append a new node to the protoTree.

This routine append periodic type to growing array of periodic

This routine converts sdr_spacer_type to spacerInterwoven_type for as each filament and add it to the growing array of spacerInterwoven_type

This routine converts sdr_spacer_type to spacerInterwoven_type for as each filament and add it to the growing array of spacerInterwoven_type

Coefficient beta from the recursive formulation of Legendre polynomials, for the Legendre mode 'mode'.

This function checks for intersection of a axis aligned box and a parallelepiped.

This function checks for intersection of box and cube

This routine checks if neighbor node with intersected boundary is level higher than current node level.

Copy the FPT header information.

This routine append 8 children to protoTree and inherit property bits from parent. leaf bit is removed from parent.

This routine creates children for each parent if children intersect with boundary object.

Routine to create an element with subelement resolution.

destroy the dynamic array

destroy the dynamic array

destroy the dynamic array

destroy the dynamic array

destroy the dynamic array

destroy the dynamic array

destroy the dynamic array

destroy the dynamic array

destroy the dynamic array

destruction of a dynamic array

destruction of a dynamic array

empty the entries without changing arrays

empty the entries without changing arrays

empty the entries without changing arrays

empty the entries without changing arrays

empty the entries without changing arrays

empty the entries without changing arrays

empty the array, reset nvals to be 0

empty the entries without changing arrays

empty the array, reset nvals to be 0

empty all contents of the array without changing the size or status of any array

empty all contents of the array without changing the size or status of any array

increase the size of the container for the array.

increase the size of the container for the array.

increase the size of the container for the array.

increase the size of the container for the array.

increase the size of the container for the array.

increase the size of the container for the array.

increase the size of the container for the array.

increase the size of the container for the array.

increase the size of the container for the array.

expanding the dynamic array

expanding the dynamic array

Mark all virtual nodes, which contain a flooded child as flooded starting from the second finest level moving up to the root. This allows to easily avoid non-flooded domains later on.

This routine loops over all intersected with geoemtry nodes and fluidify some node according to the following rule:\n 1. one of its link does noe intersect with any geometry that requires qVal 2. it has fluid neighbor on that direction. i.e. it is wet in that side. Jiaxing Qi

This routine checks for qVal of the periphery and floods if qVal < 0 and the node in that direction is fluid and not intersected by boundary

This routine loop over all nodes are flood non-interesting leaf node with wet face and inherit the wetness of the virtual node to the eligble childrens

Get sampled data.

Get sampled data.

This routine gets minBCID of the given node position in the protoTree. If the minBcid is periodic then it bcID is set to treeID of fluid node on the opposite side of periodic plane. It also computes the qVal if calc_dist = true. If qVal = -1 then there is no intersection and if qVal > 1 then the geometry is intersected after the link distance.

This function returns the position of treeID of given coordReal in the the given mesh Start from minLevel which is the level of neighbor and find the treeID which is a leaf in protoTree

This routine inherit the intersected boundary objects from parent to childrens

initialize the dynamic array

initialize the dynamic array

initialize the dynamic array

initialize the dynamic array

initialize the dynamic array

initialize the dynamic array

initialize the dynamic array

initialize the dynamic array

initialize the dynamic array

initialization of a dynamic array

initialization of a dynamic array

Initialize the node type.

Test the intersection between the given cube and the object specified by obj_pos.

This function provides the test for equality of two attributes.

This function provides the test for equality of two projections.

This function provides the test for equality of two projections.

This function provides the test for equality of two nodes descriptions.

This function provides the test for equality of two projections.

This function provides the test for equality of the header for two projections.

This function provides the test for equality of two projections.

This function provides a comparison of two attributes.

This function provides a > comparison of two projections.

This function provides a > comparison of two projections.

This function provides a > comparison of nodes descriptions.

This function provides a > comparison of two projections.

This function provides a > comparison of the header of two projections.

This function provides a > comparison of two projections.

This function provides a >= comparison of two projections.

This function provides a >= comparison of two projections.

This function provides a >= comparison of two nodes descriptions.

This function provides a >= comparison of two projections.

This function provides a >= comparison of the header of two projections.

This function provides a >= comparison of two projections.

This function provides a comparison of two attributes.

This function provides a < comparison of two projections.

This function provides a < comparison of two projections.

This function provides a < comparison of two nodes descriptions.

This function provides a < comparison of two projections.

This function provides a < comparison of the header of two projections.

This function provides a < comparison of two projections.

This function provides a comparison of two attributes.

This function provides a <= comparison of two projections.

This function provides a <= comparison of two projections.

This function provides a <= comparison of two nodes descriptions.

This function provides a <= comparison of two projections.

This function provides a <= comparison of the header of two projections.

This function provides a <= comparison of two projections.

This function provides the test for unequality of two attributes.

This function provides the test for unequality of two projections.

This function provides the test for unequality of two projections.

This function provides the test for unequality of two nodes descriptions.

This function provides the test for unequality of two projections.

This function provides the test for unequality of the header of two projections.

This function provides the test for unequality of two projections.

This routine loads the distance refine table from attribute table and appends the loaded distance in radius and level to refine nodes within radius the given radius.

load single distance refine table

Load the value fill definition for a single color.

This routine load the spacer filament type information for both length and width

Small helping routine to keep track of leaf nodes.

This routine checks if a boundary need calc qVal for a given BCID It is used in identify_boundary routine

This routine checks if a boundary need flood periphery for diagonal directions for a given BCID. It is used in identify_boundary routine

insert an element at a given position

insert an element at a given position

insert an element at a given position

insert an element at a given position

insert an element at a given position

insert an element at a given position

insert an element at a given position

adds the value to a given position inside the growing array.

adds the values starting from a given position inside the growing array.

adds the value to a given position inside the growing array.

adds the values starting from a given position inside the growing array.

adds the value to a given position inside the growing array.

adds the values starting from a given position inside the growing array.

adds the value to a given position inside the growing array.

adds the values starting from a given position inside the growing array.

adds the value to a given position inside the growing array.

adds the values starting from a given position inside the growing array.

adds the value to a given position inside the growing array.

adds the values starting from a given position inside the growing array.

adds the value to a given position inside the growing array.

adds the values starting from a given position inside the growing array.

Coefficients from the recursive formulation of legendre polynomials. L_n = alpha * x * L_n-1 + beta * L_n-2

Prodcut of alpha(numerator) * beta(denominator) / alpha(denominator)

Quotient of two alpha values.

Coefficients from the recursive formulation of legendre polynomials. L_n = alpha * x * L_n-1 + beta * L_n-2

do IDeg1 = 1, mPd+1 do IDeg2 = 1, mPd=1 !! iDeg2 = mod(iDeg-1,mpd+1)+1 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

Compute the derivative in X direction for 3D Legendre polynomial.

Compute the derivative in Y direction for 3D Legendre polynomial.

Compute the derivative in Y direction for 3D Legendre polynomial.

Subroutine to change the polynomial space (Q or P) of an atl_statedata_type from Q-space to P-space and vice versa.

Copy a single element state into a larger array and pad it with zeroes.

Copy a single 1D element state into a larger array and pad it with zeroes.

Copy a single 2D element state into a larger array and pad it with zeroes.

Copy a single element state into a larger array and pad it with zeroes.

Truncating an oversampled polynomial representation back to the original representation.

Truncating an oversampled 1D polynomial representation back to the original representation.

Truncating an oversampled 2D polynomial representation back to the original representation.

Truncating an oversampled polynomial representation back to the original representation.

Subroutine to convert linearized dof index to ansatz function number for Q-Polynomials.

Evaluate three-dimensional tensor product Legendre polynomials (not-normalized) at a given set of coordinates.

Returns the value of the non-normalized Legendre polynomial at the left boundary of the reference element, i.e. at -1.

Returns the value of the non-normalized Legendre polynomial at the left boundary of the reference element, i.e. at -1.

Returns the value of the non-normalized differentiated Legendre polynomial at the leftboundary of the reference element, i.e. at -1.

Returns the value of the gradient of the dual Legendre polynomial at the left boundary of the reference element, i.e. at -1.

Returns the value of the gradient of the dual Legendre polynomial at the left boundary of the reference element, i.e. at -1.

Returns the value of the dual Legendre polynomial at the left boundary of the reference element, i.e. at -1.

Returns the value of the dual Legendre polynomial at the left boundary of the reference element, i.e. at -1.Vectorized version.

Returns the value of the derivaitve of the dual Legendre polynomial at the left boundary of the reference element, i.e. at -1.

Returns the value of the derivaitve of the dual Legendre polynomial at the left boundary of the reference element, i.e. at -1.Vectorized version.

Returns the value of the non-normalized differentiated Legendre polynomial at the right boundary of the reference element, i.e. at +1.

Returns the value of the gradient of dual Legendre polynomial at the right boundary of the reference element, i.e. at +1.

Returns the value of the gradient of dual Legendre polynomial at the right boundary of the reference element, i.e. at +1. Vectorized version.

Returns the value of the dual Legendre polynomial at the right boundary of the reference element, i.e. at +1.

Returns the value of the dual Legendre polynomial at the right boundary of the reference element, i.e. at +1. Vectorized Version.

Returns the value of the derivaitve of the dual Legendre polynomial at the right boundary of the reference element, i.e. at +1.

Returns the value of the derivaitve of the dual Legendre polynomial at the right boundary of the reference element, i.e. at +1.vectoized version.

Load settings to describe a projection method from a Lua table.

Fill ups the bodys accroding to the DA.

Loading parameters for the filtering from the configuration script. This needs to be performed before any call of the actual transformation ply_split_element_1D.

Filter a polynomial representation in elements in one dimension according to its odd mode fraction.

Filter one-dimensional elements of degree element_degree.

Filter two-dimensional elements of degree element_degree.

Filter three-dimensional elements of degree element_degree.

Convert strip of coefficients of a modal representation in terms of Legendre polynomials to modal coefficients in terms of Chebyshev polynomials.

Convert coefficients of a modal representation in terms of Legendre polynomials to modal coefficients in terms of Chebyshev polynomials.

Define settings for the Fast Polynomial Transformation.

Print the FPT settings to the log output.

Read the FPT configuration options from the provided Lua script in conf.

Write FPT settings into a Lua table.

Convert strip of coefficients of a modal representation in terms of Legendre polynomials to modal coefficients in terms of Chebyshev polynomials.

Load settings to describe a projection method from a Lua table.

Write FXT settings into a Lua table.

Convert modal data to nodal data in 1D using flpt.

Convert modal data to nodal data in 2D using flpt.

Convert modal data to nodal data in 3D using flpt.

Convert nodal data to modal data using flpt.

subroutine to create gauss points and weights for one-dimensional integration on the interval [x1,x2].

Create Gauss-Legendre integration points and weights for one-dimensional integration on the interval [x1,x2].

function to provide the coordinates from the quadrature points on the faces

Initialize the flpt data structure for fast legendre polynomial transformation via the fxtpack.

Initialize the transformations via L2 projections.

Subroutine to initialize the fast polynomial transformation for Legendre expansion.

Integral of combination of all anzatz functions for projection onto finer element

Integral of combination of all anzatz functions for projection onto finer element

Routine to initialize the projection coefficients for a usage in the subsampling routine to project degrees of freedoms of a parent cell to the degrees of freedoms of a child cell if the degrees of freedoms are Q-Legendre polynomials.

Integrate the integrand function in Legendre basis, and represent the integral again in the Legendre basis up to the maximal degree.

Actual implementation of the matrix operation to change between nodal and modal representations.

Load settings to describe a projection method from a Lua table.

Write L2P settings into a Lua table.

Transformation between modal and nodal values in 1D via L2 projection.

Transformation between modal and nodal values in 2D via L2 projection.

Transformation between modal and nodal values in 3D via L2 projection.

Define a new polynomial in the Lagrange basis.

Evaluate a polynomial in the Lagrange basis at some point x.

Evaluate the given Lagrangian mode (which is 1 at coord(mode) and 0 in all other points) at a given point x.

\todo: as we use a relation of gamma, it might be better to use the gammln function provided by the numerical recipes, and just use the difference in an exponential function.

Evaluate all 1D Legendre polynomials at a given set of points up to the given degree.

Subroutine to convert Chebyshev (A) to Legendre (B) coefficients.

Subroutine to transform Legendre expansion to point values at Chebyshev nodes.

Subroutine to transform Legendre expansion to point values at Chebyshev nodes.

Subroutine to transform Legendre expansion to point values at Chebyshev nodes. VK: no multivar fashion of this routine is used anymore

Subroutine to transform Legendre expansion to point values at Chebyshev-Lobatto nodes.

Vectorizing subroutine to transform Legendre expansion to point values at Chebyshev-Lobatto nodes.

Subroutine to transform Legendre expansion to point values at Chebyshev nodes.

Vectorizing subroutine to transform Legendre expansion to point values at Chebyshev nodes.

Evaluate a given set of Legendre polynomials a given set of 1D points.

Initialize points with the Chebyshev quadrature points, 3D

Create the integration points on the surface of (cubical) elements.

Create multidimensional points from given 1D set of nodes in the cubic reference element.

Generates a given number of Chebyshev-Lobatto points on the unit interval [-1;+1].

Generates a given number of Chebyshev points on the unit interval [-1;+1].

Compute Gauss-Legendre integration points on the interval [-1,1].

Subroutine to transform Legendre expansion to point values at Chebyshev nodes.

Subroutine to transform Legendre expansion to point values at Chebyshev nodes.

Subroutine to transform Legendre expansion to point values at Chebyshev nodes.

Subroutine to transform Legendre expansion to point values at Chebyshev nodes.

Subroutine to transform point values at Chebyshev-Lobatto nodes to a Legendre expansion.

Vectorizing subroutine to transform point values at Chebyshev-Lobatto nodes to a Legendre expansion.

Subroutine to transform point values at Chebyshev nodes to a Legendre expansion.

Vectorizing subroutine to transform point values at Chebyshev nodes to a Legendre expansion.

Compute a multi-dimensional tensor for the given set of nodes.

Fill the body of the projection with all required data, ply_poly_project_define has to be used beforehand to set necessary header information.

Convert nDoF modes to nodal values.

Convert nodal values to nDoFs modes.

Projection of polynomial data from parent elements to child elements. The projection is done by a direct transformation of the modal coeffiecients to another coordinate system with z=ax+b.

Load settings to describe a projection method from a Lua table.

Load settings to describe a projection method from a Lua table.

Define a projection, without filling its body.

Subroutine to project elemental data from a parent cell to one of its children.

Subroutine to project element data from a parent cell to its children.

Routine to create one-dimensional projection coefficient for a coarse element to a fine element.

Function to calculate the squared L2-Norm of a given Legendre polynomial on the reference element [-1,+1].

Subsampling by L2-Projection of the Q-Tensorproduct Legendre polynomials.

Sample data described by varsys in orig_mesh according to the tracking object trackInst with adaptive refinements.

Sampling polynomial data from a given array and mesh to a new mesh with a new data array, where just a single degree of freedom per element is used.

Initialize the sampled tracking entities.

Output sampled tracking data.

Load the configuration of sampled tracking objects.

Load the configuration for adaptive subsampling.

Free previously allocated methodData of variable.

This subroutine reads the sampling configuration from the Lua script provided in conf and fills the sampling data in 'me' accordingly.

Allocate memory for a sampled variable.

This routine computes for each element whether the solution in it is considered to be deviating from the mean above the given threshold or not. The logical result is stored in var%deviates for each element.

Move the variable data from source to destination.

Create a variable system for the given tracking instance.

Function to calculate the scalar product between a Legendre polynomial (ansatz function) and a dual Legendre polynomial (test function) on the reference element [-1;+1].

Vectorized Function to calculate the scalar product between a Legendre polynomial (ansatz function) and a dual Legendre polynomial (test function) on the reference element [-1;+1] and to calculate the scalar product between a Legendre polynomial (ansatz function) and a differentiated dual Legendre polynomial (test function) on the reference element [-1;+1].

Function to calculate the scalar product between a Legendre polynomial (ansatz function) and a differentiated dual Legendre polynomial (test function) on the reference element [-1;+1].

Function to calculate the L2 scalar product of a Legendre polynomial with itself on the reference element [-1,+1].

Split one-dimensional elements of degree parent_degree into two elements with polynomials of degree child_degree.

Testing the 1D splitting.

Split two-dimensional elements of degree parent_degree into four elements with polynomials of degree child_degree.

Testing the 2D splitting.

Split three-dimensional elements of degree parent_degree into eight elements with polynomials of degree child_degree.

Testing the 3D splitting.

Initialization of the module. This needs to be performed before any call of the actual transformation ply_split_element_1D.

Project a polynomial representation in elements in one dimension to its two halves in that direction.

Testing routine for the functions of this module.

Compute the transformation matrix for a projection to the left and right half-interval of Legendre polynomials for the given maximal number of modes.

A small testing routine to check the functions of this module.

Get the color of an element.

Get the subresolution data for all elements for a given color and in the requested format.

Subroutine to load subresolution information for a given tree.

Routine to subsample mesh information for one refinement level.

Small helping routine to wrap transfers in all allowed dimensions.

Transfer of degrees of freedom from one polynomial to another in 1D.

Transfer of degrees of freedom from one polynomial to another in 2D.

Transfer of degrees of freedom from one polynomial to another in 3D.

Transfer the polynomial in P representation from on dimension to another one.

Compute the transformation matrix for a projection to the left and right half-interval of Legendre polynomials for the given maximal number of modes.

return the position of a given value in the array val, which is what you usually want to know. it is the index of a given value

return the position of a given value in the array val, which is what you usually want to know. it is the index of a given value

the actual position of a given value in the dynamic array

the actual position of a given value in the dynamic array

This function project given vector on an given plane

Routine to convert protoTree to Treelm data format. append all leaves to the temData%treeID

Small helping routine to get the variable data from a leaf.

Setup timers to assess the runtime of various parts of Seeder

Returns if any bc attribute has the distance refine option with reach_level>0.

Returns if any bc attribute has the subresolution option set.

This routine appends temporary child intersected object to actual growing array of intersected objects. To reduce memory usuage, the child with maximum number of intersected objects from parent is appended at same position as its parent. Rest of the childrens intersected objects are appended to the end of growing array

This routine appends temporary child intersected object to actual growing array of intersected objects. To reduce memory usuage, the child with maximum number of intersected objects from parent is appended at same position as its parent. Rest of the childrens intersected objects are appended to the end of growing array.

This routine created sphere objects and new attribute and extend a list of spatial objects if node intersected boundary has distance refine.

Get all attributes of the given object kind and (unique) id from the list of attributes in "attribute".

Return all colors of the given bitfield encoded in an array of characters.

This routine builds the preliminary tree with geometry intersection and neighbor identification

Clear a bit in the last byte of the node properties in a given node.

Return all colors encoded in an array of characters based on an array of logicals indicating for each color wether it should be set or not.

Get the color at all given points.

This routine dumps the final fluid tree leaves in the disk

Performance results are written to a file for statistical review

This routine find the treeID on the opposite side neighbor of the periodic plane for current leaf node

This routine identifies the nodes, which are supposed to be part of the computational domain, as defined by the seed objects.

Read the configuration for the Seeder harvesting from a Lua script.

Clean up the properties data type.

Import subresolved color data.

Load the properties from the mesh.

Identify all unique boundary colors, after all boundary attributes and seeds are known.

This routine checks for boundary neighbors and level of the boundary node

Look up colors that should be inverted and set their inversion flag accordingly.

This routines inherit distance refine sphere object from root node down to leaf node. Only the object with level greater than node level are inHerited

If parent has hasBoundary_bit then this function will inherit this property to eligible childrens

Subroutine to initialize the list for all attributes

Prominently let the user now, what he actually is running right now.

Returns if periodic bc attribute is defined

This routine loads the attribute information from the config file.

This routine loads canonical geometrial objects like point, line, plane and box and add them to the growing array of each primitive geometries and the position of this geometries are stored in the growing array of spatialObject

Load the configuration from the Lua script provided on the command line or from seeder.lua by default, if no file name is given as program argument.

This routine loads the boundCube table from config file

\brief Loading cylinder information from config file \n

This routine loads the deformation table from transformation table

Load ellipsoid information from config file.

This routine loads the geometry table from the config file i.e loading different geometry kinds like canoND, cube, periodic, STL etc.

This routine reads a single geometry table from the config file i.e loading different geometry kinds like canoND, cube, periodic, STL etc.

Routine to load spatial object defined in config file and store in geometry type

\brief load periodic table from config file.\n

\brief Loading spacer information from config file \n

This routine single spacer from object table

Routine to load single spatial object table defined in config file

Load sphere information from config file.

This routine loads STL files from config and reads the triangles from the files into the dynamic array of triangles.

This routine loads the transformation table for each spatial object table in config file

This routine loads the translation table from transformation table

Load triangle information from config file.

This routine floods the node with the given color and increases nFloodedLeaves.

This routine marks 26 direct neighbors as has boundary bit

Find the neighbor position in protoTree for iDir on the same level or on any one above.

This routine identifies the 6 direct neighbors of a node in the prototree

Return all colors of the given node encoded in an array of characters.

Set a bit in the last byte of the node properties in a given node.

Function compute intersection of plane with cube by checking intersection of two triangle of a plane with cube

This subroutine creates the treelmesh from the flooded prototree.

This routine computes the minimum distance of a given link and all the geometries in a given node:\n the link is given by a vector and a origin point.\n the node is given by the node position in the protoTree.\n If there is no intersection, qVal returns -1.0

Read the restart file into the state vectors

This routine extends the protoTree with max of minlevel or level of refinement object.

Set a bit in the last byte of the node properties in a given node.

This routine smoothens fluid domain with maximum level jumps of 1.

This function checks intesection of solid cube and spacer interwoven

Add a value definition for a color to the list of fillings.

Load the filling definition for subresolved colors.

Find the value definitions for all unique colors.

Loading the settings for the subresolution in the mesh.

This routine apply transformations to canonical objects

This routine gives special treatment when qVal > 1.0 or qVal == -1.0 for flooded neighbor, treat it as normal fluid: clean BCID, set qVal to -1 (no itersection). for non-flooded neighbor, treat it as high order wall: set qVal to 1

This routine will wet the neighbors sides of the neighbor node in all colors, that are flooded.

This routine will wet the neighbors sides of the neighbor node

Write current leaves of the prototree as treelm restart.

return the position of a given value in the list 'sorted'. this is mainly for internal usage. the sorted list is only a pointer list to the actual values thus, in order to get the index of a given value, you need to look up the entry in the sorted list. this is done by the positionofval routine

return the position of a given value in the list 'sorted'. this is mainly for internal usage. the sorted list is only a pointer list to the actual values thus, in order to get the index of a given value, you need to look up the entry in the sorted list. this is done by the positionofval routine

return the sorted position of a value in the given dynamic array

return the sorted position of a value in the given dynamic array

fix the dynamic array, meaning: store the array in the sorted order and cut off the trailing empty entries

fix the dynamic array, meaning: store the array in the sorted order and cut off the trailing empty entries

fixing the dynamic array

fixing the dynamic array

Internal subroutine to load the list of colors, which by default should apply subelement resolution to its boundaries.

This routine apply transformation to canonical objects.

This routine apply transformation to canonical objects.

Recursively traverse the tree in a depth first manner to obtain the the ordering required by treelm.

truncate the array, meaning cut off the trailing empty entries

truncate the array, meaning cut off the trailing empty entries

truncate the array, meaning cut off the trailing empty entries

truncate the array, meaning cut off the trailing empty entries

truncate the array, meaning cut off the trailing empty entries

truncate the array, meaning cut off the trailing empty entries

truncate the array, meaning cut off the trailing empty entries

truncate the array, meaning cut off the trailing empty entries

truncate the array, meaning cut off the trailing empty entries

truncate the array after the last valid entry and hence cut off the empty trailing empty entries

truncate the array after the last valid entry and hence cut off the empty trailing empty entries

Truncate the growing arrays in the node list to their actual size.

Small helping routine to write leaves in order into a treelmesh formatted file.