Variational Polygonal/Polyhedral Shape Functions (Version 1.00)

links executables usage compilation changes

This software supports finite-elements-type calculations over polygonal and polyhedral meshes. Supported applications include:

Simulation of deformable solids (with linear elasticity) in 2D and 3D,
Solution of the Franke test in 2D and 3D,
Calculation of Geodesics in Heat on polygonal meshes, and
Gradient domain processing of signals on polygonal meshes.

LINKS

Papers:

[Bunge, Herholz, Sorkine-Hornung, Botsch, and Kazhdan, 2022]

[Crane, Weischedel, and Wardetzky, 2013]

Executables:

Source Code:

Data:

Older Versions:

EXECUTABLES

--in <input polygonal/polyhedral mesh>: This string specifies the the name of the mesh.
For 2D simulations, the input polygonal mesh is assumed to be in PLY format, giving the set of vertices with the x-, y-, and z-coordinates of the positions encoded by the properties x, y, and z and the set of polygons encoded by a list of vertex indices.
For 3D simulations, the input polyhedral mesh is assumed to be in OVM format.
[--xForm <linear transform>]: This 2x2 (resp. 3x3) set of floating point values describes the entries of the linear transformation initially applied to the solid.
The default values spcify the identity transformation.
[--lock]: If enabled, this flag specifies that the values on the y-axis (resp. yz-plane) should be locked during the course of the animation.
[--gravity <gravitational force>]: This floating point value describes the force of gravity acting on the solid. (Note that without the --lock parameter, using a non-zero value for gravity will have the solid fall off the screen.)
The default value for this parameter is -500,000,000.
[--mg]: If enabled, this flag specifies that a multigrid solver should be used (instead of the default sparse Cholesky solver).
[--vCycles <number of v-cycles per animation step>]: If a multigrid solver is used, ths integer value specifies the number of v-cycles to be performed at each step of the animation.
The default value for this parameter is 1.
[--gsIters <number of Gauss-Seidel iterations per level>]: If a multigrid solver is used, ths integer value specifies the number of Gauss-Seidel iterations to be done within each level of the v-cycle.
The default value for this parameter is 5.

--in <input polygonal/polyhedral mesh>: This string specifies the the name of the mesh.
For 2D simulations, the input polygonal mesh is assumed to be in PLY format, giving the set of vertices with the x-, y-, and z-coordinates of the positions encoded by the properties x, y, and z and the set of polygons encoded by a list of vertex indices.
For 3D simulations, the input polyhedral mesh is assumed to be in OVM format.
[--mg]: If enabled, this flag specifies that a multigrid solver should be used (instead of the default sparse Cholesky solver).
[--vCycles <number of v-cycles per animation step>]: If a multigrid solver is used, ths integer value specifies the number of v-cycles to be performed at each step of the animation.
The default value for this parameter is 3.
[--gsIters <number of Gauss-Seidel iterations per level>]: If a multigrid solver is used, ths integer value specifies the number of Gauss-Seidel iterations to be done within each level of the v-cycle.
The default value for this parameter is 5.

--in <input polygonal mesh>: This string specifies the the name of the mesh.
The input polygonal mesh is assumed to be in PLY format, giving the set of vertices with the x-, y-, and z-coordinates of the positions encoded by the properties x, y, and z and the set of polygons encoded by a list of vertex indices.
[--time <diffusion time>]: This floating point values specifies the time for diffusing the source delta function .
The default value for this parameter is 1e-3.
[--mg]: If enabled, this flag specifies that a multigrid solver should be used (instead of the default sparse Cholesky solver).
[--vCycles <number of v-cycles per animation step>]: If a multigrid solver is used, ths integer value specifies the number of v-cycles to be performed at each step of the animation.
The default value for this parameter is 1.
[--gsIters <number of Gauss-Seidel iterations per level>]: If a multigrid solver is used, ths integer value specifies the number of Gauss-Seidel iterations to be done within each level of the v-cycle.
The default value for this parameter is 5.

--in <input polygonal mesh>: This string specifies the the name of the input polygonal mesh.
The polygonal mesh is assumed to be in PLY format, giving the set of vertices with the x-, y-, and z-coordinates of the positions encoded by the properties x, y, and z and the set of polygons encoded by a list of vertex indices.
[--out <output polygonal mesh>]: This string specifies the the name of the output (processed) polygonal mesh.
The polygonal mesh is assumed to be in PLY format, giving the set of vertices with the x-, y-, and z-coordinates of the positions encoded by the properties x, y, and z and the set of polygons encoded by a list of vertex indices.
[--gWeight <gradient interpolation weight>]: This floating point value specifies the weight that should be given to gradient interpolation.
The default value for this parameter is 1e-5.
[--gScale <gradient dampening/amplification factor>]: This floating point value specifies the scale that is to be appled to the gradients.
The default value for this parameter is 1, specifying unmodified output.
[--mg]: If enabled, this flag specifies that a multigrid solver should be used (instead of the default sparse Cholesky solver).
[--vCycles <number of v-cycles per animation step>]: If a multigrid solver is used, ths integer value specifies the number of v-cycles to be performed at each step of the animation.
The default value for this parameter is 3.
[--gsIters <number of Gauss-Seidel iterations per level>]: If a multigrid solver is used, ths integer value specifies the number of Gauss-Seidel iterations to be done within each level of the v-cycle.
The default value for this parameter is 5.
[--color]: If enabled and the mesh contains per-vertex color information, this flag specifies that the color at the vertices, not the position, should be processed.

USAGE EXAMPLES (WITH SAMPLE DATA)
For testing purposes, a number of polgonal/polyhedral models are provided (using the .ply and .ovm extensions respectively).

Linear Elasticity in 2D

To run this executable you must specify the input polygonal mesh. For example, to see the deformation of the unit square, tessellated by a Voronoi diagram, deforming under the action of gravity, with the vertices on the left side locked, and using a direct solver to advance time-steps, execute:

% Bin/*/DeformableSolids2D --in ../VPSF.Data/square.voronoi.3.ply --lock

To see the deformation of the unit square, tessellated using concave polygons, evolving to its rest state after an initial anisotropic scaling is applied, using a hierarchical solver to advance time-steps, execute:

% Bin/*/DeformableSolids2D --in ../VPSF.Data/square.concave.3.ply --gravity 0 --xForm 2 0 0 0.5 --mg

You can toggle the animtation by hitting [SPACE] and you can step through the animation by hitting "+".
You can also interact with the animation by left-clicking to drag a vertex.

Linear Elasticity in 3D

To run this executable you must specify the input polyhedral mesh. For example, to see the deformation of a unit cube, tessellated by a Voronoi diagram, deforming under the action of gravity, with the vertices on the left side locked, and using a direct solver to advance time-steps, execute:

% Bin/*/DeformableSolids3D --in ../VPSF.Data/cube.voronoi.3.ovm --lock

To see the deformation of the unit cube, tessellated using truncated cells, evolving to its rest state after an initial anisotropic scaling is applied, using a hierarchical solver to advance time-steps, execute:

% Bin/*/DeformableSolids3D --in ../VPSF.Data/cube.truncated.3.ovm --gravity 0 --xForm 2 0 0 0 1 0 0 0 0.5 --mg

You can toggle the animtation by hitting [SPACE] and you can step through the animation by hitting "+".
You can pan by by dragging with the left mouse button while holding down the [CTRL] key.
You can rotate by dragging with the left mouse button.
You can also rotate by using the "q", "w" , "a", "z", "s", and "x" keys.

Franke Test in 2D

To run this executable you must specify the input polygonal mesh. For example, to run the test on the unit square tessellated by a Voronoi diagram and using a direct solver, execute:

% Bin/*/FrankeTest2D --in ../VPSF.Data/square.voronoi.3.ply

To run the test on the unit square tessellated using concave polygons and using a hierarchical solver, execute:

% Bin/*/FrankeTest2D --in ../VPSF.Data/square.concave.3.ply --mg

Franke Test in 3D

To run this executable you must specify the input polyhedral mesh. For example, to run the test on the unit cube, tessellated by a Voronoi diagram, and using a direct solver, execute:

% Bin/*/FrankeTest3D --in ../VPSF.Data/cube.voronoi.3.ovm

To run the test on the unit cube, tessellated using truncated cells, and using a hierarchical solver, execute:

% Bin/*/FrankeTest3D --in ../VPSF.Data/cube.truncated.3.ovm --mg

Geodesics in Heat

To run this executable you must specify the input polygonal mesh. For example, to visualize single-source geodesics on the model of the Armadillo Man, using a direct solver, execute:

% Bin/*/GeodesicsInHeat --in ../VPSF.Data/armadillo.ply

To visualize single-source geodesics on the model of the Fanblade, using a hierarchical solver, execute:

% Bin/*/GeodesicsInHeat --in ../VPSF.Data/fanblade.ply --mg

You can specify the geodesic source by left-clicking while holding down the [SHIFT] key.
You can pan by by dragging with the left mouse button while holding down the [CTRL] key.
You can rotate by dragging with the left mouse button.
You can also rotate by using the "q", "w" , "a", "z", "s", and "x" keys.

Gradient Domain Processing

To run this executable you must specify the input and output polygonal meshes as well as the gradient interpolation weight and the gradient dampening/amplification scale. For example, to smooth the Bunny model using a direct solver, execute:

% Bin/*/GradientDomainProcessing --in ../VPSF.Data/bunny.ply --gScale 0 --out bunny.smooth.ply

For more aggressive smoothing, you can increase the gradient interpolation weight:

% Bin/*/GradientDomainProcessing --in ../VPSF.Data/bunny.ply --gScale 0 --out bunny.smooth.ply --gWeight 1e-3

To sharpen the Armadillo Man model using a hierarchical solver, execute:

% Bin/*/GradientDomainProcessing --in ../VPSF.Data/armadillo.ply --gScale 2 --out armadillo.sharp.ply --mg

COMPILATION AND EXECUTION

The Windows executables require both the glew and glut dynamically linked libraries to run. These can be found here and should be included either in the directory with the executables, or in the directory from which the executables are run.
Compiling under Windows requires both the glew and glut libraries. These can be found here and should be placed in the output directory for linkage.

HISTORY OF CHANGES

HOME