v0.11.0

• CLI: Add a switch --check_mesh_orientation to check if the mesh orientation is consistent (i.e. vertex normals have the same orientation as all adjacent face normals)
• CLI: Support for reading point/particle attributes (integers, floats and vectors of floats) from BGEO files
• Lib: Remove octree-based domain decomposition (superseded by regular grid subdivision approach --subdomain-grid)
• CLI: Set --subdomain-grid=on by default
• CLI: Remove all arguments for octree-based domain decomposition
• CLI: Remove options to output some debug files (octree grid, density map, etc.)
• Lib: Fixed a rare bug where the "density map" did not always contain a full zero one-ring around the function values (See #215)

v0.10.0

This release implements "Weighted Laplacian Smoothing for Surface Reconstruction of Particle-based Fluids" (Löschner, Böttcher, Jeske, Bender; 2023), mesh cleanup based on "Mesh Displacement: An Improved Contouring Method for Trivariate Data" (Moore, Warren; 1991) and a new, more efficient domain decomposition (see README.md for more details).

• Lib: Implement new spatial decomposition based on a regular grid of subdomains (subdomains are dense marching cubes grids)
• CLI: Make new spatial decomposition available in CLI with --subdomain-grid=on
• Lib: Implement weighted Laplacian smoothing to remove bumps from surfaces according to paper "Weighted Laplacian Smoothing for Surface Reconstruction of Particle-based Fluids" (Löschner, Böttcher, Jeske, Bender 2023)
• CLI: Add arguments to enable and control weighted Laplacian smoothing --mesh-smoothing-iters=..., --mesh-smoothing-weights=on etc.
• Lib: Implement marching_cubes_cleanup function: a marching cubes "mesh cleanup" decimation inspired by "Mesh Displacement: An Improved Contouring Method for Trivariate Data" (Moore, Warren 1991)
• CLI: Add argument to enable mesh cleanup: --mesh-cleanup=on
• Lib: Add functions to TriMesh3d to find non-manifold edges and vertices
• CLI: Add arguments to check if output meshes are manifold (no non-manifold edges and vertices): --mesh-check-manifold=on, --mesh-check-closed=on
• Lib: Support for mixed triangle and quad meshes
• Lib: Implement convert_tris_to_quads function: greedily merge triangles to quads if they fulfill certain criteria (maximum angle in quad, "squareness" of the quad, angle between triangle normals)
• CLI: Add arguments to enable and control triangle to quad conversion with --generate-quads=on etc.
• Lib: Support for reading and writing PLY meshes (MixedTriQuadMesh3d)
• CLI: Support for filtering input particles using an AABB with --particle-aabb-min/--particle-aabb-max
• CLI: Support for clamping the triangle mesh using an AABB with --mesh-aabb-min/--mesh-aabb-max

v0.9.3

Version 0.9.3

• CLI: Make input filename/input sequence pattern a positional argument instead of separate arguments, arguments --input-file and --input-sequence are removed
• CLI: Introduce some short arguments: -r for --particle-radius, -l for --smoothing-length, -c for --cube-size, -t for --surface-threshold
• CLI: Add arguments -s/--start-index and -e/--end-index to specify the first and last index to process for a sequence of files
• CLI: Process sequences of files in natural sort order instead of alphabetical order
• CLI: Show progress bar when reconstructing a sequence of files

Version 0.9.2

• Lib: In the UniformCartesianCubeGrid3d::from_aabb constructor, re-align the min-coordinate of the AABB to multiples of the cube size by default. This way multiple frames of an animation will be automatically consistent in terms of marching cubes grid alignment without having to manually specify an AABB for the entire simulation.

v0.9.1

• CLI: Fix bug where no input sequence files were found if the sequence pattern path did not contain any "parent" path (e.g. file_{}.vtk instead of ./file_{}.vtk)

v0.9.0

• Upgrade to Rust edition 2021
• Lib: Update to nalgebra 0.32.1
• Lib: Support for writing particles as JSON
• Lib: Support for writing particles as BGEO
• CLI: Support for new writers in convert subcommand
• Lib: Rename AxisAlignedBoundingBox*d typedefs to Aabb3d and Aabb2d
• Lib: Support reading VTU (VTK XML) files
• CLI: Support reconstruction of particles from VTU (VTK XML) files, including attributes
• CLI: Support for specifying an output pattern (e.g. "surface_{}.obj") using --output-file/-o when processing sequences of files, this allows to specify a different output format than VTK for sequences (which was the only available output format for sequences before)
• CLI: Files belonging to a sequence are now detected by treating the placeholder "{}" as a "(\d+)" regex (i.e. any number of digits) and sorting the matching files lexicographically (previously, the tool just counted up from 1)

v0.8.0

• Lib: Update to nalgebra 0.29
• Lib: All the particle/mesh file IO code was moved from the splashsurf binary crate to the splashsurf_lib crate. They are behind the optional io feature flag. This allows using the IO code in tests of the library itself (previously the code was duplicated). But these functions may also be helpful to library users building their own pipeline if they wan to to quickly test code by e.g. loading particle data from files.
• Lib: Replace some unsafe blocks by using bytemuck instead
• Lib: Redesign the API for SPH kernel functions
• Lib: Add module sph_interpolation for the interpolation of normals, scalar fields and vector fields to arbitrary points using SPH based interpolation
• CLI: The --output-normals=on/off argument was renamed to --normals=on/off
• CLI: Add command line switch --sph-normals=on/off (enabled by default) to toggle whether to use SPH interpolation (instead of area weighted triangle normals) to compute normals when --normals is on
• CLI: Support for writing surfaces meshes (including normals) to OBJ files directly from the reconstruct subcommand (previously this was only supported for the convert subcommand) (See #29)
• CLI: Support for specifying names of point attribute fields via the argument --interpolate-attributes=... that should be interpolated to the reconstructed surface. Currently it is only possible to load these fields from VTK files and the interpolated values can also only be written if the output file format is VTK.
• CLI: All on/off arguments now require using an equal sign to set them. Primarily, this was done to prevent the CLI from accepting empty values which would be replaced by the respective default arguments. E.g. previously --normals was equivalent to the default --normals=off which is quite confusing. Now the default is only applied when the argument is not provided at all while --normals results in an error and --normals=on can be used to override the default.

v0.7.0

• Lib: Fix a bug that caused the neighborhood search results to be incomplete
• Lib: Update to nalgebra 0.27 with const generics and removed usage of nalgebra types that are deprecated
• Lib: Add a "naive" O(n^2) neighborhood search for testing purposes
• Lib: Rename other neighborhood search functions to neighborhood_search_spatial_hashing and neighborhood_search_spatial_hashing_parallel
• CLI: Add a command line option --output-normals to the reconstruct sub-command to write mesh normals to the output file. Note that currently the normals are only computed using an area weighted average of triangle normals.
• CLI: Add JSON files as supported particle input file
• CLI: Change the convert sub-command to allow both conversion of particle files and conversion of mesh formats. Instead of using the --input option you now have to either use the --particles or --mesh option to specify the input file depending on whether you want to convert a particle or mesh file. For particles VTK, BGEO, PLY, XYZ, JSON -> VTK, PLY is supported. For meshes only VTK, PLY -> OBJ is supported.
• CLI: Rename the command line option to specify a sequence of input files from --input-sequence-pattern to --input-sequence

v0.6.1

• Lib: Add docs.rs metadata to manifest to document all features and to build only for a single target

v0.6.0

This release fixes a couple of bugs that may lead to inconsistent surface reconstructions when using domain decomposition (i.e. reconstructions with artificial bumps exactly at the subdomain boundaries, especially on flat surfaces). Currently there are no other known bugs and the domain decomposed approach appears to be really fast and robust.

In addition the CLI now reports more detailed timing statistics for multi-threaded reconstructions.

Otherwise this release contains just some small changes to command line parameters.

• Lib: Add a ParticleDensityComputationStrategy enum to the SpatialDecompositionParameters struct. In order for domain decomposition to work consistently, the per particle densities have to be evaluated to a consistent value between domains. This is especially important for the ghost particles. Previously, this resulted inconsistent density values on boundaries if the ghost particle margin was not at least 2x the compact support radius (as this ensures that the inner ghost particles actually have the correct density). This option is now still available as the IndependentSubdomains strategy. The preferred way, that avoids the 2x ghost particle margin is the SynchronizeSubdomains where the density values of the particles in the subdomains are first collected into a global storage. This can be faster as the previous method as this avoids having to collect a double-width ghost particle layer. In addition there is the "playing it safe" option, the Global strategy, where the particle densities are computed in a completely global step before any domain decomposition. This approach however is really slow for large quantities of particles. For more information, read the documentation on the ParticleDensityComputationStrategy enum.
• Lib: Fix bug where the workspace storage was not cleared correctly leading to inconsistent results depending on the sequence of processed subdomains
• Lib: Fix bug in the octree splitting where ghost particles of a parent node were actually classified as non-ghost particles in the child node, leading to many subsequent splits
• Lib: Change AABB's contains_point method such that it now considers the AABB as "open" to its max value, i.e. it checks if the point is in the half-open set [min, max[
• Lib: Implemented an own profile! macro that also works in multi-threaded code, i.e. with together with rayon
• CLI: The CLI now also prints detailed profiling/timing output when running in parallel with domain decomposition thanks to the new profile macro
• CLI: Add --domain-min and --domain-max flags to the convert subcommand that allows to filter out particles
• CLI: Remove the --splash-detection-radius flag as it did not work for a couple of releases
• Lib: Update to vtkio 0.6.0

v0.5.1

• Lib: Fixed bug in triangulation that resulted in at most 4 out 5 possible triangles being emitted. This resulted in holes in surface meshes in some cases.