/jigsaw-matlab

JIGSAW is a Delaunay-based unstructured mesh generator for two- and three-dimensional geometries.

Primary LanguageC++OtherNOASSERTION

JIGSAW: An unstructured mesh generator

     

JIGSAW is a computational library for unstructured mesh generation and tessellation; designed to generate high-quality triangulations and polyhedral decompositions of general planar, surface and volumetric domains. JIGSAW includes refinement-based algorithms for the construction of new meshes, optimisation-driven techniques for the improvement of existing grids, as well as routines to assemble (restricted) Delaunay tessellations and Voronoi complexes.

This package provides a MATLAB / OCTAVE based scripting interface to the underlying JIGSAW mesh generator, including a range of additional facilities for file I/O, mesh visualisation and post-processing operations.

JIGSAW has been compiled and tested on various 64-bit Linux , Windows and Mac based platforms.

Code Structure

JIGSAW is a multi-part library, consisting of a MATLAB / OCTAVE front-end, and a core c++ back-end. All of the heavy-lifting is done in the c++ layer - the interface contains additional scripts for file I/O, visualisation and general data processing:

├── JIGSAW  :: MATLAB/OCTAVE top-level functions
├── script  -- MATLAB/OCTAVE utilities
└── jigsaw
    ├── src -- JIGSAW source files
    ├── inc -- JIGSAW header files (for libjigsaw)
    ├── bin -- put JIGSAW exe binaries here
    ├── lib -- put JIGSAW lib binaries here
    ├── geo -- default folder for JIGSAW inputs
    ├── out -- default folder for JIGSAW output
    └── uni -- unit tests and libjigsaw programs

The MATLAB / OCTAVE interface is provided for convenience - you're not forced to use it, but it's perhaps the easiest way to get started!

It's also possible to interact with the JIGSAW back-end directly, either through (i) scripting: building text file inputs and calling the JIGSAW executable from the command-line, or (ii) programmatically: using JIGSAW data-structures within your own applications and calling the library via its API.

Getting Started

The first step is to compile and configure the code! JIGSAW can either be built directly from src, or installed using the conda package manager.

Building from src

The full JIGSAW src can be found in ../jigsaw/src/.

JIGSAW is a header-only package - the single main jigsaw.cpp file simply #include's the rest of the library directly. JIGSAW does not currently dependent on any external packages or libraries.

JIGSAW consists of several pieces: (a) a set of command-line utilities that read and write mesh data from/to file, and (b) a shared library, accessible via a C-format API.

Using cmake

JIGSAW can be built using the cmake utility. To build, follow the steps below:

* Ensure you have the cmake utility installed.
* Clone or download this repository.
* Navigate to the main `../jigsaw/` directory.
* Create a new temporary directory BUILD (to store the cmake build files).
* Navigate into the temporary directory.
* Execute: cmake -D CMAKE_BUILD_TYPE=BUILD_MODE ..
* Execute: make
* Execute: make install
* Delete the temporary directory.

This process will build a series of executables and the shared library: jigsaw itself - the main command-line meshing utility, tripod - JIGSAW's tessellation infrastructure, as well as libjigsaw - JIGSAW's shared API. BUILD_MODE can be used to select different compiler configurations and should be either RELEASE or DEBUG.

See example.jig for documentation on calling the command-line executables, and the headers in ../jigsaw/inc/ for details on the API.

Using g++ / llvm

JIGSAW has been successfully built using various versions of the g++ and llvm compilers. The build process is a simple one-liner (from ../jigsaw/src/):

g++ -std=c++11 -pedantic -Wall -O3 -flto -D NDEBUG
-D __cmd_jigsaw jigsaw.cpp -o ../bin/jigsaw

will build the main jigsaw command-line executable,

g++ -std=c++11 -pedantic -Wall -O3 -flto -D NDEBUG
-D __cmd_tripod jigsaw.cpp -o ../bin/tripod

will build the tripod command-line utility (JIGSAW's tessellation infrastructure) and,

g++ -std=c++11 -pedantic -Wall -O3 -flto -fPIC -D NDEBUG
-D __lib_jigsaw jigsaw.cpp -shared -o ../lib/libjigsaw.so

will build JIGSAW as a shared library (libjigsaw).

Install via conda

JIGSAW is also available as a conda environment. To install and use, follow the steps below:

* Ensure you have conda installed. If not, consider miniconda as a lightweight option.
* Add conda-forge as a channel: conda config --add channels conda-forge
* Create a jigsaw environment: conda create -n jigsaw jigsaw

Each time you want to use JIGSAW simply activate the environment using: conda activate jigsaw

Once activated, the various JIGSAW command-line utilities will be available in your run path, JIGSAW's shared library (libjigsaw) will be available in your library path and its include files in your include path.

Example Problems

After compiling and configuring the code, navigate to the JIGSAW installation directory in your MATLAB / OCTAVE environment and run the following set of example problems:

meshdemo(0); % simple 2-dim. examples to get started
meshdemo(1); % simple 3-dim. examples to get started
meshdemo(2); % frontal-delaunay methods for surfaces
meshdemo(3); % frontal-delaunay methods for volumes
meshdemo(4); % dealing with sharp-features in piecewise smooth domains
meshdemo(5); % dealing with user mesh-size controls
meshdemo(6); % dealing with topological constraints
meshdemo(7); % mesh iso-surface geometry -- (case 1)
meshdemo(8); % mesh iso-surface geometry -- (case 2)

Additionally, the ../jigsaw/example.jig file provides a description of JIGSAW's configuration options, and can be used as a command-line example. A set of unit-tests and libjigsaw example programs are contained in ../jigsaw/uni/. The JIGSAW-API is documented via the header files in ../jigsaw/inc/. A repository of 3D surface models generated using JIGSAW can be found here.

License

This program may be freely redistributed under the condition that the copyright notices (including this entire header) are not removed, and no compensation is received through use of the software. Private, research, and institutional use is free. You may distribute modified versions of this code UNDER THE CONDITION THAT THIS CODE AND ANY MODIFICATIONS MADE TO IT IN THE SAME FILE REMAIN UNDER COPYRIGHT OF THE ORIGINAL AUTHOR, BOTH SOURCE AND OBJECT CODE ARE MADE FREELY AVAILABLE WITHOUT CHARGE, AND CLEAR NOTICE IS GIVEN OF THE MODIFICATIONS. Distribution of this code as part of a commercial system is permissible ONLY BY DIRECT ARRANGEMENT WITH THE AUTHOR. (If you are not directly supplying this code to a customer, and you are instead telling them how they can obtain it for free, then you are not required to make any arrangement with me.)

DISCLAIMER: Neither I nor: Columbia University, the Massachusetts Institute of Technology, the University of Sydney, nor the National Aeronautics and Space Administration warrant this code in any way whatsoever. This code is provided "as-is" to be used at your own risk.

References

There are a number of publications that describe the algorithms used in JIGSAW in detail. If you make use of JIGSAW in your work, please consider including a reference to the following:

[1] - Darren Engwirda: Generalised primal-dual grids for unstructured co-volume schemes, J. Comp. Phys., 375, pp. 155-176, https://doi.org/10.1016/j.jcp.2018.07.025, 2018.

[2] - Darren Engwirda, Conforming Restricted Delaunay Mesh Generation for Piecewise Smooth Complexes, Procedia Engineering, 163, pp. 84-96, https://doi.org/10.1016/j.proeng.2016.11.024, 2016.

[3] - Darren Engwirda, Voronoi-based Point-placement for Three-dimensional Delaunay-refinement, Procedia Engineering, 124, pp. 330-342, http://dx.doi.org/10.1016/j.proeng.2015.10.143, 2015.

[4] - Darren Engwirda, David Ivers, Off-centre Steiner points for Delaunay-refinement on curved surfaces, Computer-Aided Design, 72, pp. 157-171, http://dx.doi.org/10.1016/j.cad.2015.10.007, 2016.

[5] - Darren Engwirda, Locally-optimal Delaunay-refinement and optimisation-based mesh generation, Ph.D. Thesis, School of Mathematics and Statistics, The University of Sydney, http://hdl.handle.net/2123/13148, 2014.