Releases: llnl/axom
Axom-v0.12.0
Please download the Axom-v0.12.0.tar.gz tarball below, which includes all of the Axom submodules as well. The AxomData-v0.12.0.tar.gz tarball contains extra files to put in the data directory only if you want to run certain Axom tests.
Added
- Added a new Python interface for sidre, using nanobind to generate Python bindings.
- Added a new "BUMP" (Blueprint Utilities for Mesh Processing) component in Axom, which includes utilities that were formerly included in the Axom MIR component. BUMP is useful for writing algorithms that process Blueprint meshes.
- New
axom::MALLOC_ALLOCATOR_IDis for using malloc and free even when axom is configured with Umpire support. - The
axom::mir::ElviraAlgorithmclass, which performs material interface reconstruction using the ELVIRA algorithm, was enhanced so it supports 3D structured mesh inputs. The output mesh is a Blueprint mesh with a 3D unstructured polyhedral topology. - The
axom::mir::ElviraAlgorithmclass, was enhanced to accept a "plane" option that causes it to return clipping plane origin and normal as fields on the mesh. - The
axom::mir::ElviraAlgorithmclass, was enhanced to accept a "pointmesh" option that causes it to return a mesh consisting of points located at clipping plane origins for each clipped material fragment, instead of returning polygonal or polyhedral meshes. This option is off by default. - Exposed primal clip operations for clipping various shapes with a plane.
- Adds constructs in the
axomnamespace that wrap RAJA atomics, reductions, scans, and sorts. When RAJA is not available, serial-only substitutes are provided, allowing algorithms to still compile and run. These constructs are templated on theExecSpace(execution space) so it is not necessary to query RAJA policies via theexecution_spacetype traits classes. - 2D and 3D implementations for
axom::for_allwere added. - Adds
axom::FlatMapView, a helper class associated withaxom::FlatMapto support queries from within a GPU kernel. - Adds
axom::FlatMap::create<ExecSpace>()andaxom::FlatMap::insert<ExecSpace>()to support constructing or inserting a hash map over a batch of keys and values on the GPU or with OpenMP. - Adds support for custom allocators to
axom::FlatMap. - Primal: Adds ability to perform sample-based shaping on tetrahedral shapes.
- Improves efficiency of volume fraction computation from quadrature samples during sample-based shaping.
- Adds a
axom::DeviceHashtype as a GPU-enabled version of thestd::hashinterface. - Added a new
quest::STLWriterclass that writes mint meshes to STL format. - Adds
assignmethods toaxom::Array. - Adds
assign,fill,setmethods toaxom::ArrayView. - Core: Adds a
TempFileclass to Axom's FileUtilities to help with generating temp files with unique file names that can be automatically removed when the instance goes out of scope. - Klee: We now support optional specification of a per-shape
dimensionsfield for the geometry of a shape. These can be used to override the globaldimensionsof a Klee input file. - Lumberjack: Added sorting of log messages by creation time. This is now the default.
- Slic: Added check for pending messages which is done by default when calling flushStreams() and pushStreams().
- Sina: Records can now be provided with a curve ordering to use when writing to file. By default, all records will now use oldest-first ordering (ULTRA-like)
- Sina: Documents can optionally be written as HDF5 instead of JSON. HDF5 should provide better performance for large, curve-rich sets
- Sina: Written documents can now be appended to. This is a flexible system with a few different uses, ex: continuous writing of timeseries, overwriting values that change over the course of a simulation, and snapshot handling. See documentation for details.
- Adds
quest::MFEMReaderfor reading 1D MFEM contours in 2D space. - Adds an option to
quest::SamplingShaperto allow in/out tests based on winding numbers for MFEM contours. - The
shaping_driverexample program can select--sampling inoutto do the default In/Out sampling and--sampling windingnumberto select winding number in/out tests for MFEM data.
Changed
- Updates blt submodule to [BLT version 0.7.1][https://github.com/LLNL/blt/releases/tag/v0.7.1]
- Updates to [Conduit version 0.9.5][https://github.com/LLNL/conduit/releases/tag/v0.9.5]
- Updates to [RAJA version 2025.09.0][https://github.com/LLNL/RAJA/releases/tag/v2025.09.0]
- Updates to [camp version 2025.09.2][https://github.com/LLNL/camp/releases/tag/v2025.09.2]
- Updates to [Umpire version 2025.09.0][https://github.com/LLNL/Umpire/releases/tag/v2025.09.0]
- Axom now requires
C++17and will default to that if not specified viaBLT_CXX_STD. - Fixed
Timer::elapsed*()methods so they properly report the sum of all start/stop cycles since the lastreset(). - Adds support for allocations using
mallocandfreeeven when Axom is configured with Umpire support. - Adds a new utility tool,
mesh_converter, which converts between mesh formats. The first conversion is from a Pro-E tetrahedral mesh to an STL mesh of its boundary triangles. - Primal: Adds a method to determine if a point is contained within a Tetrahedron.
- The
primal::BoundingBoxclass'expand()andshift()methods were modified so they do nothing when called on invalid bounding boxes. - Updates to [MFEM version 4.8.0][https://github.com/mfem/mfem/releases/tag/v4.8]
- Readers in Quest were moved from a
quest/readersdirectory toquest/io. - Sina: Renames a Fortran module to
sina_hdf5_config(fromhdf5_config) - Spin: Uses
axom::FlatMapinSparseOctreeLevelimplementation. We have observed a performance regression of about 20% during InOutOctree construction and queries over STL surface meshes relative to the previous sparsehash implementation. Please reach out to Axom developers if this affects you while we work on fixes for these. - Klee: Moves source files related to IO into a new
iosubdirectory in the Klee component - Primal: Consolidates construction logic for
BezierCurve,BezierPatch,KnotVector,NURBSCurveandNURBSPatchclasses and add overloads fromaxom::ArrayView - Primal: 2D and 3D winding number methods are now accelerated via memoization (dynamic caching + reuse) when supplied arrays of query points
- Core: Updates behavior of
FlatMap::reserve()to only trigger a rehash if maximum load factor would be exceeded. - Quest: Moves curve linearization from the
quest::C2CReaderintoquest::LinearizeCurvesso the logic can be used with other curve data.
Fixed
- Core: prevent incorrect instantiations of
axom::Arrayfrom a host-only compile, when Axom is compiled with GPU support. Instances where this occurs will now trigger a static assertion during compile time. - Fixes build with
ninjagenerator - Primal: Fixes a
BoundingBoxconstructor with zero (or fewer) points - Sina: Fixes configuration variables related to inclusion of
AdiakWriter.hppand to hdf5 support insina_fortran_interface.f - Spin: Fixes undefined behavior in BVH tree construction associated with using signed indexes
- Spin: Fixes undefined behavior in UniformGrid construction associated with invalid geometry bounding boxes
- Core: Fixes undefined behavior in MapCollection when searching empty collections
- Core: Fixes some edge cases in the
joinPathfile utility - Core: Fixes
FlatMap::erase()to update reported size. - Core: Fixes
FlatMap::rehash()when allocated in device-only memory.
Deprecated
- Primal: Deprecates
Triangle::checkInTriangle(pt). UseTriangle::contains(pt)instead.
Axom-v0.11.0
Please download the Axom-v0.11.0.tar.gz tarball below, which includes all of the Axom submodules as well. The AxomData-v0.11.0.tar.gz tarball contains extra files to put in the data directory only if you want to run certain Axom tests.
Added
- Added a new "Mir" Axom component for accelerated Material Interface Reconstruction (MIR) algorithms.
MIR algorithms take Blueprint meshes with a matset and they use the matset's material information
to split any input zones that contain multiple materials into zones that contain a single material.
The Mir component provides an implementation of the Equi-Z MIR algorithm, which is a visualization-
oriented algorithm that produces smooth interfaces between zones and their neighbors. The Mir
component also provides a 2D ELVIRA algorithm, which reconstructs polygonal shapes and preserves
volume fractions. - Support in
quest::IntersectionShaperfor Blueprint mesh stored in aconduit::Node
orsidre::Group. - Adds new CMake configuration options,
AXOM_ENABLE_ASANandAXOM_ENABLE_UBSAN, to enable/disable AddressSanitizer and UndefinedBehaviorSanitizer respectively in Axom. Default is OFF for both. - A number of new
klee::Geometryconstructors are added, for the different shapes now supported.
This is a temporary change. The class will be subclassed in the future to support a diversity of geometries. - Support some analytical shapes in
IntersectionShaper. - Support generating shapes in memory (not requiring input files).
sidre::Viewholding array data may now be re-shaped. Seesidre::View::reshapeArray.- Sina C++ library is now a component of Axom
- Adds optional dependency on Open Cascade. The initial intention is
to use Open Cascade's file I/O capabilities in support of Quest applications. - Adds
primal::NURBSCurveandprimal::NURBSPatchclasses, supported byprimal::KnotVector. - Adds a Quest example that reads in a STEP file using Open Cascade and processes its geometry
- Adds a piecewise method to load external data using
sidre::IOManager. This adds new overloaded methods
ofloadExternalDatainsidre::IOManagerandsidre::Group. - Adds intersection routines between
primal::Rayobjects andprimal::NURBSCurve/primal::NURBSPatchobjects. - Adds LineFileTagCombiner to Lumberjack to allow combining messages if line number, file, and tag are equal.
- Adds some support for 2D shaping in
quest::IntersectionShaper, using STL meshes with zero for z-coordinates or in-memory triangles as input. - Adds ability in Lumberjack to own and set communicators.
- Adds
NonCollectiveRootCommunicatorto Lumberjack to provide an MPI-based communicator for logging messages non-collectively.
Changed
- Updates blt submodule to [BLT version 0.7.0][https://github.com/LLNL/blt/releases/tag/v0.7.0]
- Updates to [MFEM version 4.7.0][https://github.com/mfem/mfem/releases/tag/v4.7]
- Updates to [Caliper version 2.12.1][https://github.com/LLNL/Caliper/releases/tag/v2.12.1]
- Updates to [Conduit version 0.9.3][https://github.com/LLNL/conduit/releases/tag/v0.9.3]
- Updates to [RAJA version 2025.03.0][https://github.com/LLNL/RAJA/releases/tag/v2025.03.0]
- Updates to [camp version 2025.03.0][https://github.com/LLNL/camp/releases/tag/v2025.03.0]
- Updates to [Umpire version 2025.03.0][https://github.com/LLNL/Umpire/releases/tag/v2025.03.0]
primal::NumericArrayhas been moved tocore. The header iscore/NumericArray.hpp.quest::Shaperandquest::IntersectionShaperconstructors require a runtime policy.
Changing the policy after construction is no longer supported.- Importing Conduit array data into
sidre::Viewnow allocates destination
data using theView's parent's allocator ID, instead of always using
host memory. This is consistent with the behavior of deep-copying data
from Sidre. - ItemCollection and its child classes MapCollection, ListCollection, and IndexedCollection were moved from Sidre
to core. The namespace prefix for these classes is nowaxom::instead ofaxom::sidre. The internal usage of
these types within Sidre Datastore and Group is unchanged. MFEMSidreDataCollection::LoadExternalDatanow takes two optional string parameters, one that is a
filename (defaults to thenamemember variable) and the other is aGrouppath relative to the base of
the Data Collection itself (defaults to the root of theDataStore).SLIC_ASSERT,SLIC_ASSERT_MSG,SLIC_CHECK, andSLIC_CHECK_MSGmacros delegate to assert() within HIP device kernels.
Fixed
- Fixes compilation issue with [email protected] on 32-bit Windows configurations.
This required a RAJA fix to avoid 64-bit intrinsics,
as well as support for 32-bitWords in Slam'sBitSetclass. - Minor bugfix to
primal::intersect(segment, ray)to better handle cases when segment and ray overlap. - Fixes a memory leak in
axom::Arraycopy constructor.
Axom-v0.10.1
Please download the Axom-v0.10.1.tar.gz tarball below, which includes all of the Axom submodules. The AxomData-v0.10.1.tar.gz tarball contains extra files to put in the data directory only if you want to run certain Axom tests.
Added
- Constructor to Axom BVH that avoids an unnecessary copy of each bounding box.
Fixed
- Issue with uninitialized state in axom::Array class was causing host-initialization of device-allocated memory in certain situations. This could cause warnings about uninitialized memory or crashes in the axom::Array constructor.
- Added a guard for sidre-related mint API usage in a quest example.
- Removed
std::endsusage fromSLIC_ASSERT,SLIC_ASSERT_MSG,SLIC_CHECK, andSLIC_CHECK_MSGmacros that prevented Lumberjack from combining messages. - Some line numbers linking file contents into the Axom Quickstart Guide were incorrect causing the docs to appear incomplete.
Axom-v0.10.0
Please download the Axom-v0.10.0.tar.gz tarball below, which includes all of the Axom submodules. The AxomData-v0.10.0.tar.gz tarball contains extra files to put in the data directory only if you want to run certain Axom tests.
Added
- Added SLIC constructors that take in a
std::stringfor the stream. If string is interpreted as a file name, the file is not opened until SLIC flushes and the stream has at least one message logged. - Primal: Adds a
clip()operator overload for clipping a 2D polygon against another 2D polygon. - Primal: Adds
Polygon::reverseOrientation()to reverse orientation of a polygon in-place. - Adds
StaticArray, a wrapper forStackArraywith a size member variable. - Multidimenional
core::Arraysupports column-major and arbitrary stride ordering, in addition to the default row-major ordering. - Adds new
PolygonArrayandMAX_VERTStemplate parameters toprimal::Polygonfor dynamic or static allocation. - Adds support for the optional
caliperandadiakdependencies to axom. These dependencies are added through axom'sspackpackage via the new+profilingvariant, and are enabled in axom's build system via theCALIPER_DIRandADIAK_DIRconfiguration paths. - Adds new annotation macros to axom:
AXOM_ANNOTATE_{BEGIN,END,SCOPE,METADATA}. These replace the previous annotation macrosAXOM_PERF_MARK_{FUNCTION,SECTION}. - Adds a RAII-based
MPIWrapperutility class to axom's core component. This can help setup/teardown MPI in axom's examples. It can also be used in configurations with MPI. - Primal: Adds a
closest_pointoperator for finding the closest point on aSegment - Primal: Adds a
reflectPointmethod to thePlaneprimitive - Primal: Makes several primitive methods available in device code
- Improves support for
axom::Arrayallocated in unified and pinned memory on GPU platforms. Use of GPU-based operations for Arrays allocated in a unified memory space is controlled with a new method,Array::setDevicePreference(). - Adds
svg2contoursscript to convert paths in an SVG file to an MFEM NURBS mesh. - Quest: Adds an example to query winding numbers on an MFEM NURBS mesh.
Changed
- Updates to [Conduit version 0.9.2][https://github.com/LLNL/conduit/releases/tag/v0.9.2]
- Updates to [RAJA version 2024.07.0][https://github.com/LLNL/RAJA/releases/tag/v2024.07.0]
- Updates to [camp version 2024.07.0][https://github.com/LLNL/camp/releases/tag/v2024.07.0]
- Updates to [Umpire version 2024.07.0][https://github.com/LLNL/Umpire/releases/tag/v2024.07.0]
axom::CLI::ExitCodes::Successhas been changed toaxom::CLI::ExitCodes::CLI11_Successto avoid conflict when X11#definesSuccess.MarchingCubesmasking now uses the mask field's integer values instead of converting them to booleans. The new behavior lets you select a value to mask for. If you want to continue the boolean behavior, use only 0 or 1 in your mask field.- Primal:
Polyhedron::centroid()function changed to return center of mass of the polyhedron.Polyhedron::vertexMean()added to return average of polyhedron's vertices.Polyhedron::moments()returns the volume and centroid of the polyhedron, the 0th and 1st moments respectively. quest::ArrayIndexeris nowaxom::MDMapping, adopting conventional terminology and moving out ofquest.mint::structured_execis nowaxom::nested_for_exec, to support nested for loops for all of Axom. Seesrc/axom/core/execution/nested_for_exec.hpp.- Set default Umpire allocator id to device instead of unified for CUDA and HIP execution policies.
- Upgrades
vcpkgusage for axom's automated Windows builds to its 2024.03.19 release. Also updates vcpkg port versions for axom dependencies. Temporarily removesumpirefrom axom's default dependencies on Windows due to incompatibility between umpire's externalfmtand axom's vendored copy. - Turn off CMake finding dependencies on system paths.
axom::Array: trivially-copyable types with a non-trivial constructor are now initialized on the GPU.- SLIC no longer outputs the rank count in the
RANKformat string in parallel loggers. You can access the rank count via new format optionRANK_COUNT.
Removed
- Removes config option
AXOM_ENABLE_ANNOTATIONS. Annotations are now provided bycaliper(andadiakfor metadata) and are available when axom is configured withCALIPER_DIRandADIAK_DIRconfig variables. - Removes caching of
{PACKAGE}_FOUNDvariables inSetupAxomThirdParty.cmake - We no longer test Axom with the XL compiler. So users should consider XL unsupported.
Fixed
numerics::eigen_solve()has been corrected to avoid an early return with error state.
Axom-v0.9.0
Please download the Axom-v0.9.0.tar.gz tarball below, which includes all of the Axom submodules as well. The AxomData-v0.9.0.tar.gz tarball contains extra files to put in the data directory only if you want to run certain Axom tests.
Added
- Primal: Adds a
Quadrilateralprimitive - Primal: Adds a
compute_bounding_box()operator for computing the bounding
box of aQuadrilateral - Primal: Adds a
clip()operator for clipping a tetrahedron against the
half-space defined by a plane - Primal: Adds a
checkAndFixOrientation()function toprimal::Tetrahedron
that swaps the order of vertices if the signed volume of the Tetrahedron is
negative, resulting in the signed volume becoming positive. - Adds
FlatMap, a generic key-value store which aims for drop-in compatibility
withstd::unordered_map, but utilizes an open-addressing design. - Adds support for device-side use of
Array::push_back()andArray::emplace_back(). - Adds initial support for using Slic streams with tags
- Adds an example that finds intersection candidate pairs between two Silo
hexahedral meshes using either a BVH or Implicit Grid spatial index - Quest: Adds
setTetPredFromBoundingBox()andsetTetPred()functions to
quest::ProEReaderandPProEReaderthat set a tet predicate, allowing
user code to read in a subset of a Pro/E ASCII tetrahedron mesh file.
Changed
MarchingCubeshas optimizations to improve GPU performance, particularly for
repeated computations. The constructor has changed and a newsetMeshmethod
is added to set (or change) the mesh. New accessors present output data
without moving them from device to host. These accessors are an interim
solution and likely to be updated in the future.DistributedClosestPointoutputs are now controlled by thesetOutputmethod.MarchingCubesallows user to select the underlying data-parallel implementationfullParallelworks best on GPUs.hybridParallelreduces the amount of data processed and works best with
MarchingCubesRuntimePolicy::seq.byPolicy(the default) selects the implementation based on the runtime policy.
MarchingCubesandDistributedClosestPointclasses identify domains by their
state/domain_idparameters if provided, or the local iteration index if not.MarchingCubesandDistributedClosestPointclasses changed from requiring the Blueprint
coordset name to requiring the Blueprint topology name. The changed interface methods are:DistributedClosestPoint::setObjectMeshDistributedClosestPoint::computeClosestPointsMarchingCubes::MarchingCubesMarchingCubesSingleDomain::MarchingCubesSingleDomain
- Primal:
Polyhedron::volume()function changed from returning a signed
volume to an unsigned volume. The addedPolyhedron::signedVolume()function
returns the signed volume. - Primal:
intersection_volume()operators changed from returning a signed
volume to an unsigned volume. - Primal's
BoundingBox::contains(BoundingBox)now returnstruewhen the input is empty - Renamed axom's bit utility functions to conform to
C++20standard:popCount() -> popcount(),
trailingZeros() -> countr_zero()andleadingZeros() -> countl_zero() - Renamed
axom::utilities::swapEndian() -> byteswap()to conform toC++23standard
Fixed
- quest's
SamplingShapernow properly handles material names containing underscores - quest's
SamplingShapercan now be used with an mfem that is configured for (GPU) devices - primal's
Polygonarea computation in 3D previously only worked when the polygon was aligned with the XY-plane. It now works for arbitrary polygons. - Upgrades our
vcpkgusage for automated Windows builds of our TPLs to its 2023.12.12 release - Fixed a bug in the bounds checks for
primal::clip(Triangle, BoundingBox) - Fixed a bug when loading Sidre groups with attributes that already exist
- Fixed
std::localeerror when when compilingsrc/axom/core/utilities/System.cppusing nvcc - Include
cstdintfor higher gcc version support (e.g. gcc-13) - Fixed several memory leaks in
axom::Array,quest::Shapingandsidre::MFEMSidreDataCollection
Axom-v0.8.1
Please download the Axom-v0.8.1.tar.gz tarball below, which includes all of the Axom submodules as well. The AxomData-v0.8.1.tar.gz tarball contains extra files to put in the data directory only if you want to run certain Axom tests.
Changed
- Updates to [RAJA version 2023.06.0][https://github.com/LLNL/RAJA/releases/tag/v2023.06.0]
- Updates to [camp version 2023.06.0][https://github.com/LLNL/camp/releases/tag/v2023.06.0]
- Updates to [Umpire version 2023.06.0][https://github.com/LLNL/Umpire/releases/tag/v2023.06.0]
Fixed
- Fixed MFEMSidreDataCollection finite element space bug
- Various fixes to CMake machinery
Axom-v0.8.0
Please download the Axom-v0.8.0.tar.gz tarball below, which includes all of the Axom submodules as well. The AxomData-v0.8.0.tar.gz tarball contains extra files to put in the data directory only if you want to run certain Axom tests.
Added
- Adds MarchingCubes class implementing the marching cubes algorithm for surface detection.
- Adds the following methods to
axom::Arrayto conform more closely with thestd::vectorinterface:Array::front(): returns a reference to the first elementArray::back(): returns a reference to the last elementArray::resize(size, T value): resizes the array, and sets any new elements tovalue.
- Adds an
ArrayView::empty()method to return whether the view is empty or not. - Adds an
area()function toprimal::Polygon - Adds initial support for using Slam types on the GPU
- Adds support for using
ArrayViewIndirectionindirection policy withslam::Mapandslam::BivariateMap - Adds
const_iteratorsupport toslam::BivariateMapandslam::SubMap - Primal: Adds a
Hexahedronprimitive - Primal: Adds a
clip()operator for computing the intersection of aHexahedronand anotherTetrahedronas aPolyhedron - Primal: Adds an
intersection_volume()operator for computing the volume of intersection between a primitive and aTetrahedron - Primal: Adds a
primal::Polyhedron::from_primitive()operator that returns aPolyhedronobject from a given primitive. - Adds
DataStore::getBufferInfo()andGroup::getDataInfo()methods that insert information into a ConduitNodeabout buffers in aDataStoreobject or data in aGroupsubtree. The information can be accessed from theNodeby the caller from specifically named fields in theNode. - Quest: Adds a
quest::ProEReaderfor reading in Pro/E tetrahedral meshes - Quest: The
quest::IntersectionShaperclass can now use a percent error to determine whether the revolved volume for a shape is sufficiently accurate or whether the shape must be further refined. This new dynamic method of shaping complements the existing segment-based curve refinement method and it is activated usingShaper::setRefinementType()and by callingShaper::setPercentError()to set a refinement error percentage. - Multimat: adds initial support for fields stored on the GPU.
MultiMat::setAllocatorID()orMultiMat::setFieldAllocatorID()can be called to change the memory space in which a field is allocated. - Multimat: adds an
MultiMat::addExternalField()function to support fields where the memory is managed externally. Fields which are externally-managed cannot be transposed between sparse and dense layouts, or moved between allocator spaces. - Multimat: adds a
MultiMat::removeField()function to remove fields from the Multimat instance. - Multimat: adds an overload of
MultiMat::setCellMatRel()that supports setting a multi-material relation in a compressed sparse-row (CSR) representation. - Quest: Adds ability to import volume fractions into
SamplingShaperbefore processingKleeinput - Slam: adds a
slam::MappedVariableCardinalitypolicy to accelerate mapping flat indices back to first-set indices when used in aStaticRelation - Adds an
ArrayView(data, shape, strides)constructor to support column-major and custom striding layouts. - Adds an
ArrayView::subspan()overload for multi-dimensional subspans - Adds an
axom::utilities::insertionSort()method. - Quest: Adds Pro/E tetrahedral meshes as input to the
IntersectionShaper - Quest: For sample-based shaping, users can register callback functions to modify the input points before querying the spatial index. This allows, e.g. querying 3D points against 2D surfaces of revolution provided as c2c contour files.
- Adds an
axom::utilities::localeutility function to guard against platforms that do not have the requested locales via thestd::localefunction. If the system does not have the requested locale (e.g.en_US.UTF8), it returns the user's default locale. - Sidre: Add new protocols
sidre_layout_jsonandconduit_layout_jsonto provide output of DataStore layout in a user-readable format that excludes the numerical arrays held by Views and Buffers. - Sidre: Add methods to methods for destroying Groups that will also destroy Buffers if the destruction of a Group and the Views in its subtree cause the Buffer to become detached from all Views.
- Sidre: Add two Group methods -- one to return a vector of the valid I/O protocols (based on compilation options), and one that returns a default protocol.
- Klee: Add support in shaping driver and MFEMSidreDataCollection to write Blueprint datasets that contain matset metadata needed for VisIt to treat volume fraction arrays as a material. This enables VisIt plots such as FilledBoundary.
Changed
- Fixed bug in
mint::mesh::UnstructuredMeshconstructors, affecting capacity. A missing factor was added. If you worked around this by adding the factor yourself, you may want to undo that work-around. - Updates blt submodule to [email protected]
- Updates uberenv submodule to HEAD of main on 12May2023
- Updates to conduit version 0.8.6
- Updates to mfem version 4.5
- Updates to fmt version 9.1.0
- Updates to
c2cversion 1.8.0 - The Axom library has been broken down into its component libraries (prefixed with
axom_). This change requires no change to downstream CMake users who import our targets. The exported CMake targetaxomincludes all components, but users who do not import our targets will need to create the link line themselves. The following replacement can be used:-laxom->-laxom_quest -laxom_multimat -laxom_slam -laxom_mint -laxom_klee -laxom_inlet -laxom_sidre -laxom_slic -laxom_lumberjack -laxom_core. If you only need a subset of the components, you can now use those targets directly, ie.axom::inlet. IntersectionShapernow implements material replacement rules.axom::Arraymove constructors are nownoexcept.- Exported CMake targets,
cli11,fmt,sol, andsparsehash, have been prefixed withaxom::to guard against conflicts. DistributedClosestPointquery now supports any blueprint-valid mesh format, including multidomain. Domain underloading and overloading can be expressed using multidomain format. Closest points are identified by cp_rank, cp_domain_index, and cp_index. The new cp_domain_index specifies the domain containing the closest point.DistributedClosestPointinterfacing variable namesclosest_pointandmin_distancehave been changed tocp_coordsandcp_distance, respectively, to match the naming convention of other interfacing variables.- Adds
vcpkgports forRAJA,Umpirewith optionalOpenMPfeature for automated Windows build. - Reduce size of
ArrayView::subspanto prevent accessing invalid memory. - Adds
vcpkgport forluaas optional dependency on Windows. - Adds additional parameters to quest's
PointInCellquery to control the Newton solve from physical to reference space for a given element. - Remove function pointer call in IteratorBase::advance().
- Slam:
IndirectionPolicy::data()now returns a reference to the underlying buffer. Rebinding an indirection to a new buffer is now achieved throughIndirectionPolicy::ptr(), which returns a mutable pointer to the buffer. - Quest:
Shaper::applyTransforms()is no longer a public method. - Multimat: fields are now returned as shallow, device-copyable views of a field instead of full copies of field data.
- Multimat:
MultiMat::addField()andMultiMat::setVolfracField()API now useaxom::ArrayViewto accept data. - Multimat: Ported field data/sparsity layout conversion methods to GPU.
- Multimat:
MultiMat::makeOtherRelation()now runs on the GPU with an appropriately-set allocator ID. - Multimat:
MultiMat::setCellMatRel(counts, indices)now runs on the GPU, and accepts GPU-side data. - Renames
ArrayView::spacing()toArrayView::minStride(). - Klee: A shape's geometry no longer needs a
pathfield when itsformatis "none" - Quest: Shapes without geometry can participate in replacement rules for sample-based shaping. Volume fractions for the associated materials must be supplied before shaping.
- Primal: Expose Polyhedron::getFaces() as public method.
- Removed custom Spack package recipes in favor of using the radiuss-spack-configs repository as a submodule.
Fixed
- Fixed issues with CUDA build in CMake versions 3.14.5 and above. Now require CMake 3.18+ for CUDA and non-gpu builds.
- Fix to allow Axom to build when using RAJA, but not Umpire.
- Checks validity of bounding boxes in
primalintersection operators against planes and triangles before using the geometry. - Improves import logic for
luadependency - Improves import logic for
mfemdependency in device builds whenmfemis configured withcaliper - Fixes ambiguity when calling
Array::resize(size, value)forArray<bool> - Sidre: Group methods
hasChildView()andhasChildGroup()changed to return false when Group holds items in list format. This fixes an IOManager issue causing failures for multi-domain meshes when writing Blueprint index file. - Sidre: Changed Group::copyToConduitNode() method to properly handle Groups using the list format. Previously, it was assumed that all child items in a Group have a name, which is not the case for list format.
- Sidre: Fixed Group method
importConduitTreeExternal()to handle lists with unnamed members the same asimportConduitTree(). - Slic and Lumberjack: Add missing calls to
flush()for parallel output log streams.
Deprecated
- Integer types in
src/axom/core/Types.hppare deprecated becausec++11supports their equivalents.
Axom-v0.7.0
Please download the Axom-v0.7.0.tar.gz tarball below, which includes all of the Axom submodules as well.
Added
- Adds a
view()method toaxom::Arrayclass to simplify creation of a correspondingaxom::ArrayView - Adds GPU/OpenMP support to
spin::ImplicitGrid.
The following functions run with the user-specified execution space (specified as a template argument
onImplicitGrid):ImplicitGrid::insert(nelems, bboxes): insert a batch of bounding boxes into the implicit gridImplicitGrid::getCandidatesAsArray(nquery, queryObjs, ...): query the implicit grid for a batch of
query objects, and generate a CSR-format array for the candidates.
In addition,ImplicitGrid::getQueryObject()returns an object that may be used within a GPU kernel
to query the implicit grid.
- Added initial implementation of GPU/OpenMP-accelerated point-in-cell queries
- Added an alternative surface mesh tester function to Quest, based on
ImplicitGrid - Add
constversions ofbegin()andend()forArrayandArrayView - Add support for passing compatible custom allocator IDs to
axom::Arraywith explicitly specified
memory space - Adds constructor overloads to
axom::Arraythat uses both uninitialized data (ArrayOptions::Uninitialized)
and custom allocators - Adds a comparison operator (
operator<()) toStackArray. This allows it to be used as a key forstd::map - Use compiler intrinsics for axom's bit utility functions:
popCount(),trailingZeros()andleadingZeros() - Adds a random-access iterator to
slam::DynamicSet - Adds an overload to
ImpicitGrid::getCandidates()from a grid cell of the lattice.
This makes it easier to iterate over the bins of the spatial index. - Defines iterator traits on
axom::Array<T>/ArrayView<T>iterators, to allow passing iterator
pairs to standard library functions - Adds full support for calling methods on
axom::Array<T>allocated in device-only memory. - Adds ability to index into subarrays of a multidimensional
axom::Array<T>using
operator[]andoperator() - Adds ability to build axom using the
hipcompiler. Support for running device
kernels with hip will be added in the future. - Adds new host-configs for HIP on LLNL platforms
- Adds GPU/OpenMP support to
spin::UniformGrid.
The following functions run with a user-specified execution space (specified as a template argument
onUniformGrid):UniformGrid::initialize(): creates/re-creates a uniform grid with an array of objects and their
corresponding bounding boxesUniformGrid::getCandidatesAsArray(): query the uniform grid for objects that share a grid
cell with the query bounding box
In addition,UniformGrid::getQueryObject()returns an object that may be used within a GPU kernel
to query the uniform grid.
- Adds ability to specify a storage policy in
UniformGrid. Two policies are provided:DynamicGridStoragestores the bins as an array of arrays (default)FlatGridStoragestores the bins as a flat array of elements, where each bin is a slice of the
array
- Adds a templated uniform grid-based surface mesh tester function to Quest
- Adds an initializer list constructor and assignment operator to
axom::Array - Adds an overload of
axom::Array::resize(ArrayOptions::Uninitialized, dims)to support resizes
without constructing or initializing new elements - Adds examples and tests for using Slic interface in Fortran
- Adds examples for using the BVH device traversal API
- Adds a
ScatteredInterpolationquery to quest, which enables interpolating scalar fields at arbitrary
points from a given 2D or 3D point mesh in the Mesh Blueprint format. The current implementation
generates a Deluanay complex over the points and performs linear interpolation over the
triangle/tetrahedron at the query points. - Adds a HIP execution policy for device kernels to run on AMD GPU hardware
- Adds new Slic macros that allow you to selectively print messages only on root ranks. For example,
SLIC_ERROR_ROOT(msg)andSLIC_ERROR_ROOT_IF(EXP, msg). This can be set via
slic::initialize(bool is_root = true)orslic::setIsRoot(). - Adds forward iterators to the
Views andGroups of asidre::Group.
These can be accessed via the range-for syntax asfor(auto& view: grp.views()){...},
Or using the iterator syntax as
for(auto& it = grp.views().begin(), itEnd = grp.views().end(); it ! itEnd ; ++it) {...}, and similarly for the groups of a group. - Adds forward iterators to the
Attributes andBufferss of asidre::DataStore,
with a similar syntax, e.g.for(auto& buf : datastore.buffers()){...}. - Adds an overload of
ImplicitGrid::getCandidatesAsArray()to accept query points/bounding boxes
as anaxom::ArrayView. - Adds a
primal::closest_point(point,sphere)overload to find the closest point on a sphere to a given point - Adds an overload to quest's
SignedDistancequery to return the closest point on the surface
to the query point and the surface normal at that point. Also exposes this functionality
in quest's signed_distance C API. - Adds utility function for linear interpolation (
lerp) of two numbers - Adds utility function to compute binomial coefficients
- Adds a
CurvedPolygonclass to primal representing a polygon withBezierCurves as edges - Adds functions to compute the moments (area, centroid) of a
CurvedPolygon - Adds functions to evaluate integrals over
BezierCurveandCurvedPolygonobjects - Adds a
ArrayViewIndirectionstorage policy to Slam - Adds set accessor methods to
slam::DynamicVariableRelation - Adds a new component to Axom,
multimat, to simplify the handing of multi-material meshes and
fields. - Adds functions to compute winding numbers and in/out queries for
PolygonandCurvedPolygonobjects. - Adds
constants.hppto primal to track geometric constants. Initially includes
a value forPRIMAL_TINY, a small constant that can be added to
denominators to avoid division by zero. DistributedClosestPointquery now supports "domain underloading" -- ranks that are passed in can
have empty object meshes and/or empty query meshes- 'BezierCurve' objects now support Rational Bezier curve functionality
- Primal: Adds a
clip()operator for computing the intersection of aTetrahedronand anotherTetrahedronas aPolyhedron - Added
slic::outputLocalMessages()to output messages from the current rank to the console for MPI-enabled LogStreams.
Changed
- Axom now requires C++14 and will default to that if not specified via
BLT_CXX_STD. - Moved bit-twiddling functions to core component
axom::Arraynow default-initializes its data by default. To create an Array with uninitialized
elements, pass anaxom::ArrayOptions::Uninitializedas the first constructor argument.axom::ArrayView<const T>can now be created from aconst Array<T>- Added new
ExecSpacetemplate parameter tospin::ImplicitGrid.
ExecSpaceis now the second template parameter (out of three) and defaults toaxom::SEQ_EXEC. - Instead of saving the entire
DataStore,MFEMSidreDataCollectionwill now save only
its domain and global groups - When an
inlet::Fieldfails a range or valid value constraint, the provided value and
corresponding range/set of valid values are now included in the error message - IOManager::write now allows the calling code to pass in the full name of
the root file that it will produce - Improved consistency of orientation operations in
primal::Plane,primal::Sphere,primal::orientation()
andprimal::in_sphere() - Improved efficiency for
primal::in_sphere()-- the computations are now based on(D+1)x(D+1)determinants
instead of(D+2)x(D+2)determinants forD-dimensional spheres - Improved efficiency for (signed) area/volume functions of
primal::Segment,primal::Triangleandprimal::Tetrahedron - Updates interface for
primal::Sphereto use more ofprimal, e.g. usesprimal::Pointinstead ofT*to represent points - Adds
circumsphere()functions toprimal::Triangleandprimal::Tetrahedronto return theSphere
that circumscribes the triangle/tetrahedron's vertices - Adds operator overloads to subtract a
primal::Vectorfrom aprimal::Pointto yield a newprimal::Point - Consolidates
quest::findTriMeshIntersections*()implementations forBVHandImplicitGrid BVH::find*()batch functions now return the total number of candidate intersections found- Enables empty
axom::Array<T>to be iterated over withbegin()/end() - Removed
AXOM_VERSION_EXTRAin favor ofaxom::gitSHA()and adding the SHA toaxom::getVersion()and
axom::about() - Use more specific type trait checks in
ArrayOps, to avoid generating unnecessary copies in
fill/destroy operations on otherwise trivially-copyable/destructible types. axom::Arraynow consistently propagates the allocator ID on copy, move, and swap operations when possible.
This is a breaking change; copy-construction of a dynamic array from a device array will no longer automatically
move the array to host memory, and will instead maintain the same allocator ID as the source array.- The device traversal method
BVH::TraverserType::traverse_tree()now supports passing in arbitrary query objects
for BVH traversal. - Moved
inlet::LuaReader::solState()to be a protected function that now returns astd::shared_ptr<axom::sol::state>.
This is an advanced feature that could cause users to break an input file state after verification. This also allows us
to not exposeaxom/sol.hppto all users of Inlet. This greatly reduces compile times. Using this feature requires
both a derived class and includingaxom/sol.hppin the user code. - Renamed some overloads of function
createViewof...
Axom-v0.6.1
Please download the Axom-v0.6.1.tar.gz tarball below, which includes all of the Axom submodules as well.
Added
- Added a config variable,
AXOM_DEBUG_DEFINEto control whether theAXOM_DEBUGcompiler define is enabled.
ByDEFAULT, it is enabled forDebugandRelWithDebInfoconfigurations, but this can be overriden
by settingAXOM_DEBUG_DEFINEtoONorOFF. axom::Arrayis now GPU-compatible, in particular via a memory space template parameter and via
extensions toaxom::ArrayViewthat allow for copying into kernels and transfers between memory spaces.- Adds some utility arithmetic operators for adding and subracting
primal::Points andprimal::Vectors
Changed
- Renamed
AXOM_NOT_USEDmacro toAXOM_UNUSED_PARAMfor better consistency with other Axom macros - Added
explicittoaxom::Inlet::InletVectorconstructors and added a constructor that accepts adouble* AXOM_ENABLE_MFEM_SIDRE_DATACOLLECTIONconfiguration option is nowONby default (rather thanOFF).
This option should be disabled ifmfemwas configured withMFEM_USE_SIDRE.
Fixed
- The
AXOM_DEBUGcompiler define is now properly exported via theaxomCMake target when it is enabled - Added tolerance parameter
EPStoprimal::closest_point()operator. This effectively snaps
closest points to the triangle boundaries vertices and edges when they are withinEPS,
improving consistency when, e.g., querying multiple triangles from the same mesh. - Fixed regression in
SignedDistancequeries for query points closest to edges or vertices
of the input triangle mesh
Axom-v0.6.0
Please download the Axom-v0.6.0.tar.gz tarball below, which includes all of the Axom submodules.
Added
- Added new CMake option to allow users to turn off Axom created tools:
AXOM_ENABLE_TOOLS - Inlet can now log verification errors to a user-processable list instead of using Slic
- Sidre parallel I/O: Added new mapping arrays to the automatically-generated Blueprint
index to support new schema for multi-domain parallel meshes. - Added support for optional third-party
c2c("contours to codes") library for parsing 2D spline data.
c2cis currently only available for Axom configurations on LLNL platforms. - Primal's
intersect(Ray, Segment)can now return the parametric coordinates of the intersection
along both the ray and the segment (when the intersection exists) - Primal's
intersect(Segment, BoundingBox)can now return the parametric coordinates bounding the
portion of the segment contained within the BoundingBox (when the intersection exists) - Generalizes Quest's
InOutOctreeclass to work with 2D line segment meshes. Previously,
it only worked with 3D triangle meshes - Added support for reading
c2c".contour" files in Quest. Contours that enclose a 2D region can be linearized
into line segment meshes and loaded into Quest'sInOutOctreefor in/out queries. - Updated the Quest
inoutC API to support 2D queries using thec2clibrary, when Axom is configured withc2c - Updated the C++ Quest "containment" example to support 2D in/out queries
(in addition to the already supported 3D queries) - Added
axom::Arraymodeled afterstd::vector. Previousaxom::Arrayrenamed toaxom::MCArray. Future changes to both arrays are expected. - Added a
data_collection_utiltool to generate Mesh Blueprint compliant high order distributed meshes from
an mfem mesh or over a Cartesian domain - Added utility functions
axom::utilities::getHostName()andaxom::utilities::getUserName(). - Added new
axom::primal::ZipIterable<T>type to convert structure-of-arrays data to a given
Primal geometric primitive. - Quest: Added a
computeDistances()function toSignedDistanceclass for batched
signed-distance queries. - Spin: Added a
getTraverser()function toBVH, enabling the customized traversal of a
BVH from within a device kernel. - Primal: Adds an
Octahedronprimitive - Primal: Adds a
Polyhedronprimitive for representing convex polyhedra bounded by planar polygons in 3D - Primal: Adds a
clip()operator for computing the intersection of aTetrahedronand anOctahedronas aPolyhedron - Klee: Adds a new component,
klee, for specifying non-conformal shape overlays for materials onto simulation meshes.
This component defines a schema for defining, transforming and overlaying 2D and 3D shapes
and validates klee input files. See the klee documentation for more information. - Quest: Adds a new query for sampling-based "shaping" onto low- or high-order computational meshes
- Quest: Adds a new query for intersection-based "shaping" of revolved contours onto 3D hexahedral meshes.
This capability uses a RAJA policy operate on various execution spaces (host, openmp, device). - Quest: Adds a "shaping" example for embedding a klee specification onto an MFEM mesh
- Added Sidre function
View::clear(). - Core now provides an
axom::ArrayViewthat provides view/indexing semantics over a raw pointer.
This replaces the external buffer logic previously provided byaxom::Array.
Changed
-
MFEMSidreDataCollectionnow reuses FESpace/QSpace objects with the same basis -
Harden configuration options for BLT tools (style, code quality, etc.) against accidentally being enabled for users. Developers will
always give a full path (e.g.CLANGFORMAT_EXECUTABLE) -
Inlet:
Writers are passed directly toInlet::writeinstead of being registered -
Inletobjects can now be constructed without a user-providedsidre::DataStore -
Conduit version changed to v. 0.7.2
-
Renames
AXOM_DEBUG_VARmacro toAXOM_UNUSED_VARsince there were many cases where the latter
was not the appropriate name. This macro elides warnings about unused variables -
Inlet's
isUserProvidedcan now be used to query the status of subobjects of aContainervia a name parameter -
Upgrades our
vcpkgusage for automated Windows builds of our TPLs to its 2021.05.12 release -
Upgrades built-in
cli11library to its v1.9.1 release -
Quest's
inoutC API has two new functions:inout_set_dimension()andinout_set_segments_per_knot_span().
The latter is only applicable for 2D queries onc2ccontours -
Spin: Refactored
BVHpublic API based on user suggestions
BVHconstructor only handles setting up default values, while the actual building of the BVH is
now done in aBVH::initialize(primal::BoundingBox*, int)method.
Alternate Umpire allocator IDs are supplied viaBVH::setAllocatorID(int).
OtherBVHmethods have been modified to accept or return Primal primitives. -
Spin: Removed hard dependency on RAJA and Umpire from
BVH. -
Moved
slam::IteratorBasetoaxom::IteratorBase. -
sidre::Arraynow derives fromaxom::MCArray. -
axom::Arrayis now multidimensional; it intends to behave likestd::vectorin the 1D case
andnumpy.ndarrayin the multidimensional case -
Quest:
SignedDistancehas been modified to usespin::BVHinstead ofBVHTree. This
enables signed-distance queries to run on the GPU, as specified via a new template
parameter. -
Spin: Removed
BVHTreeclass in favor ofBVH. -
Quest's
signed_distanceC API: Removed functions related to oldBVHTreeclass
and added functions related toBVHclass- Removed:
void signed_distance_set_max_levels( int maxLevels ) - Removed:
void signed_distance_set_max_occupancy( int maxOccupancy ) - Added:
void signed_distance_set_allocator( int allocatorID ) - Added:
void signed_distance_set_execution_space( SignedDistExec execSpace )
- Removed:
-
All built-in third-party libraries (
fmt,cli11,sol, andsparsehash) have been guarded to allow downstream users to
have their own versions. This includes moving their headers underinclude/axominstead ofinclude/and
moving their C++ namespace underaxom(eg.fmttoaxom::fmt). If you don't use our built-n TPLs this has no
affect on you, but if you do these are some the changes you will need to make:Library Namespace changes Header include changes fmt fmt::→axom::fmt::#include "fmt/fmt.hpp"→#include "axom/fmt.hpp"sol sol::→axom::sol::#include "sol/sol.hpp"→#include "axom/sol.hpp"sparsehash google::→axom::google::#include "sparsehash→#include "axom/sparsehashcli11 CLI::→axom::CLI::#include "CLI11/CLI11.hpp"→#include "axom/CLI11.hpp" -
Moved
axom::MCArrayand thesidre::Arrayit was based on intomint
asaxom::deprecated::MCArrayandsidre::deprecated::MCArray, respectively.
sidre::Arrayis now based onaxom::Array. -
utilities::string::splitnow returns a vector instead of using an out-parameter,
Inlet's string utilities were moved to Core, andsplitLastNTokenswas renamed torsplitN -
axom::Array-related classes have been moved into individual files. -
RAJA dependency updated to 0.14.0
-
Umpire dependency updated to 0.6.0. Support for versions prior to v2.1.0 was removed.
-
Conduit dependency updated to 0.7.2+ (develop as of Sept 13, 2021). This was required because Spack
is now usingHDF5's CMake build system. -
Internal BLT dependency updated to 0.4.1
Fixed
- Fixed Primal's
intersect(Ray, Segment)calculation for Segments that do not have unit length - Fixed problem with Cray Fortran compiler not recognizing MSVC pragmas in
axom/config.hpp.
The latter are now only added in MSVC configurations. - Fixed bug in
Mint's VTK output for fields of typeint64andfloat - Improved loading of data collections in
MFEMSidreDataCollection - Added workaround to
MFEMSidreDataCollectionforC++14standard library feature that was not available in[email protected] - Delayed finalizing reloaded mesh in
MFEMSidreDataCollectionuntil after setting
the nodalGridFunction(when applicable) - Transposed
RandZcoordinates when linearizing NURBS curves inc2creader - Fixed user-reported in/out ambiguity within some
InOutOctreecases with grazing triangles