Releases: LLNL/axom
Axom-v0.10.1
Please download the Axom-v0.10.1.tar.gz tarball below, which includes all of the Axom submodules.
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::ends
usage fromSLIC_ASSERT
,SLIC_ASSERT_MSG
,SLIC_CHECK
, andSLIC_CHECK_MSG
macros 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.
Added
- Added SLIC constructors that take in a
std::string
for 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 forStackArray
with a size member variable. - Multidimenional
core::Array
supports column-major and arbitrary stride ordering, in addition to the default row-major ordering. - Adds new
PolygonArray
andMAX_VERTS
template parameters toprimal::Polygon
for dynamic or static allocation. - Adds support for the optional
caliper
andadiak
dependencies to axom. These dependencies are added through axom'sspack
package via the new+profiling
variant, and are enabled in axom's build system via theCALIPER_DIR
andADIAK_DIR
configuration 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
MPIWrapper
utility 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_point
operator for finding the closest point on aSegment
- Primal: Adds a
reflectPoint
method to thePlane
primitive - Primal: Makes several primitive methods available in device code
- Improves support for
axom::Array
allocated 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
svg2contours
script 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::Success
has been changed toaxom::CLI::ExitCodes::CLI11_Success
to avoid conflict when X11#define
sSuccess
.MarchingCubes
masking 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::ArrayIndexer
is nowaxom::MDMapping
, adopting conventional terminology and moving out ofquest
.mint::structured_exec
is 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
vcpkg
usage for axom's automated Windows builds to its 2024.03.19 release. Also updates vcpkg port versions for axom dependencies. Temporarily removesumpire
from axom's default dependencies on Windows due to incompatibility between umpire's externalfmt
and 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
RANK
format 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
(andadiak
for metadata) and are available when axom is configured withCALIPER_DIR
andADIAK_DIR
config variables. - Removes caching of
{PACKAGE}_FOUND
variables 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.
Added
- Primal: Adds a
Quadrilateral
primitive - 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::ProEReader
andPProEReader
that set a tet predicate, allowing
user code to read in a subset of a Pro/E ASCII tetrahedron mesh file.
Changed
MarchingCubes
has optimizations to improve GPU performance, particularly for
repeated computations. The constructor has changed and a newsetMesh
method
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.DistributedClosestPoint
outputs are now controlled by thesetOutput
method.MarchingCubes
allows user to select the underlying data-parallel implementationfullParallel
works best on GPUs.hybridParallel
reduces the amount of data processed and works best with
MarchingCubesRuntimePolicy::seq
.byPolicy
(the default) selects the implementation based on the runtime policy.
MarchingCubes
andDistributedClosestPoint
classes identify domains by their
state/domain_id
parameters if provided, or the local iteration index if not.MarchingCubes
andDistributedClosestPoint
classes changed from requiring the Blueprint
coordset name to requiring the Blueprint topology name. The changed interface methods are:DistributedClosestPoint::setObjectMesh
DistributedClosestPoint::computeClosestPoints
MarchingCubes::MarchingCubes
MarchingCubesSingleDomain::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 returnstrue
when the input is empty - Renamed axom's bit utility functions to conform to
C++20
standard:popCount() -> popcount()
,
trailingZeros() -> countr_zero()
andleadingZeros() -> countl_zero()
- Renamed
axom::utilities::swapEndian() -> byteswap()
to conform toC++23
standard
Fixed
- quest's
SamplingShaper
now properly handles material names containing underscores - quest's
SamplingShaper
can now be used with an mfem that is configured for (GPU) devices - primal's
Polygon
area computation in 3D previously only worked when the polygon was aligned with the XY-plane. It now works for arbitrary polygons. - Upgrades our
vcpkg
usage 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::locale
error when when compilingsrc/axom/core/utilities/System.cpp
using nvcc - Include
cstdint
for higher gcc version support (e.g. gcc-13) - Fixed several memory leaks in
axom::Array
,quest::Shaping
andsidre::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.
[v0.8.1] - Release date 2023-08-16
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.
Added
- Adds MarchingCubes class implementing the marching cubes algorithm for surface detection.
- Adds the following methods to
axom::Array
to conform more closely with thestd::vector
interface: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
ArrayViewIndirection
indirection policy withslam::Map
andslam::BivariateMap
- Adds
const_iterator
support toslam::BivariateMap
andslam::SubMap
- Primal: Adds a
Hexahedron
primitive - Primal: Adds a
clip()
operator for computing the intersection of aHexahedron
and anotherTetrahedron
as 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 aPolyhedron
object from a given primitive. - Adds
DataStore::getBufferInfo()
andGroup::getDataInfo()
methods that insert information into a ConduitNode
about buffers in aDataStore
object or data in aGroup
subtree. The information can be accessed from theNode
by the caller from specifically named fields in theNode
. - Quest: Adds a
quest::ProEReader
for reading in Pro/E tetrahedral meshes - Quest: The
quest::IntersectionShaper
class 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
SamplingShaper
before processingKlee
input - Slam: adds a
slam::MappedVariableCardinality
policy 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::locale
utility function to guard against platforms that do not have the requested locales via thestd::locale
function. 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_json
andconduit_layout_json
to 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::UnstructuredMesh
constructors, 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
c2c
version 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 targetaxom
includes 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
. IntersectionShaper
now implements material replacement rules.axom::Array
move constructors are nownoexcept
.- Exported CMake targets,
cli11
,fmt
,sol
, andsparsehash
, have been prefixed withaxom::
to guard against conflicts. DistributedClosestPoint
query 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.DistributedClosestPoint
interfacing variable namesclosest_point
andmin_distance
have been changed tocp_coords
andcp_distance
, respectively, to match the naming convention of other interfacing variables.- Adds
vcpkg
ports forRAJA
,Umpire
with optionalOpenMP
feature for automated Windows build. - Reduce size of
ArrayView::subspan
to prevent accessing invalid memory. - Adds
vcpkg
port forlua
as optional dependency on Windows. - Adds additional parameters to quest's
PointInCell
query 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::ArrayView
to 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
path
field when itsformat
is "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
primal
intersection operators against planes and triangles before using the geometry. - Improves import logic for
lua
dependency - Improves import logic for
mfem
dependency in device builds whenmfem
is 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.hpp
are deprecated becausec++11
supports 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::Array
class 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
const
versions ofbegin()
andend()
forArray
andArrayView
- Add support for passing compatible custom allocator IDs to
axom::Array
with explicitly specified
memory space - Adds constructor overloads to
axom::Array
that 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
hip
compiler. 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:DynamicGridStorage
stores the bins as an array of arrays (default)FlatGridStorage
stores 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
ScatteredInterpolation
query 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
View
s andGroup
s 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
Attribute
s andBuffers
s 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
SignedDistance
query 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
CurvedPolygon
class to primal representing a polygon withBezierCurve
s as edges - Adds functions to compute the moments (area, centroid) of a
CurvedPolygon
- Adds functions to evaluate integrals over
BezierCurve
andCurvedPolygon
objects - Adds a
ArrayViewIndirection
storage 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
Polygon
andCurvedPolygon
objects. - Adds
constants.hpp
to 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. DistributedClosestPoint
query 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 aTetrahedron
and anotherTetrahedron
as 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::Array
now default-initializes its data by default. To create an Array with uninitialized
elements, pass anaxom::ArrayOptions::Uninitialized
as the first constructor argument.axom::ArrayView<const T>
can now be created from aconst Array<T>
- Added new
ExecSpace
template parameter tospin::ImplicitGrid
.
ExecSpace
is now the second template parameter (out of three) and defaults toaxom::SEQ_EXEC
. - Instead of saving the entire
DataStore
,MFEMSidreDataCollection
will now save only
its domain and global groups - When an
inlet::Field
fails 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::Triangle
andprimal::Tetrahedron
- Updates interface for
primal::Sphere
to use more ofprimal
, e.g. usesprimal::Point
instead ofT*
to represent points - Adds
circumsphere()
functions toprimal::Triangle
andprimal::Tetrahedron
to return theSphere
that circumscribes the triangle/tetrahedron's vertices - Adds operator overloads to subtract a
primal::Vector
from aprimal::Point
to yield a newprimal::Point
- Consolidates
quest::findTriMeshIntersections*()
implementations forBVH
andImplicitGrid
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_EXTRA
in 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::Array
now 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.hpp
to all users of Inlet. This greatly reduces compile times. Using this feature requires
both a derived class and includingaxom/sol.hpp
in the user code. - Renamed some overloads of function
createView
of...
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_DEFINE
to control whether theAXOM_DEBUG
compiler define is enabled.
ByDEFAULT
, it is enabled forDebug
andRelWithDebInfo
configurations, but this can be overriden
by settingAXOM_DEBUG_DEFINE
toON
orOFF
. axom::Array
is now GPU-compatible, in particular via a memory space template parameter and via
extensions toaxom::ArrayView
that allow for copying into kernels and transfers between memory spaces.- Adds some utility arithmetic operators for adding and subracting
primal::Point
s andprimal::Vector
s
Changed
- Renamed
AXOM_NOT_USED
macro toAXOM_UNUSED_PARAM
for better consistency with other Axom macros - Added
explicit
toaxom::Inlet::InletVector
constructors and added a constructor that accepts adouble*
AXOM_ENABLE_MFEM_SIDRE_DATACOLLECTION
configuration option is nowON
by default (rather thanOFF
).
This option should be disabled ifmfem
was configured withMFEM_USE_SIDRE
.
Fixed
- The
AXOM_DEBUG
compiler define is now properly exported via theaxom
CMake target when it is enabled - Added tolerance parameter
EPS
toprimal::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
SignedDistance
queries 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.
c2c
is 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
InOutOctree
class 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'sInOutOctree
for in/out queries. - Updated the Quest
inout
C API to support 2D queries using thec2c
library, 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::Array
modeled afterstd::vector
. Previousaxom::Array
renamed toaxom::MCArray
. Future changes to both arrays are expected. - Added a
data_collection_util
tool 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 toSignedDistance
class 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
Octahedron
primitive - Primal: Adds a
Polyhedron
primitive for representing convex polyhedra bounded by planar polygons in 3D - Primal: Adds a
clip()
operator for computing the intersection of aTetrahedron
and anOctahedron
as 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::ArrayView
that provides view/indexing semantics over a raw pointer.
This replaces the external buffer logic previously provided byaxom::Array
.
Changed
-
MFEMSidreDataCollection
now 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:
Writer
s are passed directly toInlet::write
instead of being registered -
Inlet
objects can now be constructed without a user-providedsidre::DataStore
-
Conduit version changed to v. 0.7.2
-
Renames
AXOM_DEBUG_VAR
macro toAXOM_UNUSED_VAR
since there were many cases where the latter
was not the appropriate name. This macro elides warnings about unused variables -
Inlet's
isUserProvided
can now be used to query the status of subobjects of aContainer
via a name parameter -
Upgrades our
vcpkg
usage for automated Windows builds of our TPLs to its 2021.05.12 release -
Upgrades built-in
cli11
library to its v1.9.1 release -
Quest's
inout
C API has two new functions:inout_set_dimension()
andinout_set_segments_per_knot_span()
.
The latter is only applicable for 2D queries onc2c
contours -
Spin: Refactored
BVH
public API based on user suggestions
BVH
constructor 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)
.
OtherBVH
methods have been modified to accept or return Primal primitives. -
Spin: Removed hard dependency on RAJA and Umpire from
BVH
. -
Moved
slam::IteratorBase
toaxom::IteratorBase
. -
sidre::Array
now derives fromaxom::MCArray
. -
axom::Array
is now multidimensional; it intends to behave likestd::vector
in the 1D case
andnumpy.ndarray
in the multidimensional case -
Quest:
SignedDistance
has been modified to usespin::BVH
instead ofBVHTree
. This
enables signed-distance queries to run on the GPU, as specified via a new template
parameter. -
Spin: Removed
BVHTree
class in favor ofBVH
. -
Quest's
signed_distance
C API: Removed functions related to oldBVHTree
class
and added functions related toBVH
class- 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/axom
instead ofinclude/
and
moving their C++ namespace underaxom
(eg.fmt
toaxom::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/sparsehash
cli11 CLI::
→axom::CLI::
#include "CLI11/CLI11.hpp"
→#include "axom/CLI11.hpp"
-
Moved
axom::MCArray
and thesidre::Array
it was based on intomint
asaxom::deprecated::MCArray
andsidre::deprecated::MCArray
, respectively.
sidre::Array
is now based onaxom::Array
. -
utilities::string::split
now returns a vector instead of using an out-parameter,
Inlet's string utilities were moved to Core, andsplitLastNTokens
was 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 typeint64
andfloat
- Improved loading of data collections in
MFEMSidreDataCollection
- Added workaround to
MFEMSidreDataCollection
forC++14
standard library feature that was not available in[email protected]
- Delayed finalizing reloaded mesh in
MFEMSidreDataCollection
until after setting
the nodalGridFunction
(when applicable) - Transposed
R
andZ
coordinates when linearizing NURBS curves inc2c
reader - Fixed user-reported in/out ambiguity within some
InOutOctree
cases with grazing triangles
Axom-v0.5.0
Please download the Axom-v0.5.0.tar.gz tarball below, which includes all of the Axom submodules as well.
Added
- Added the MFEMSidreDataCollection class for describing [MFEM] meshes and associated fields. This
class was adapted from MFEM's SidreDataCollection and is enabled when Axom is built with MFEM
and theAXOM_ENABLE_MFEM_SIDRE_DATACOLLECTION
CMake option is enabled. - Added
slic::setAbortFunction
to configure a custom callback when SLIC aborts. - Added a
batched
option to quest'sInOutOctree
containment query example application.
This uses a kernel to test for containment on an array of points.
The query uses OpenMP threading, when available. - Inlet: Added support for user-defined conversions from Inlet tables to user-defined
types, and support for arrays of user-defined types - Added compiler define
NOMINMAX
toaxom/config.hpp.in
to avoid problems with
the Windowsmin
andmax
macros. - Added
cpp14
variant to Spack package to allowInlet::LuaReader
to be used more easily. - Inlet: Added support for string-keyed associative arrays (dictionaries)
- Inlet: Added support for defining and retrieving functions from the input file
- Inlet: Added support for YAML and JSON input files
- Inlet: Added support for mixed-key (integer and string) associative arrays
- Inlet: Added support for deeply nested containers of structs
- Inlet: Added support for
void
and strings in Lua-defined functions - Inlet: Added
get<std::vector<T>>
for retrieving arrays without index information - Inlet: Added a new
Writer
for generating JSON schemas which can be used by text editors
for autocompletion - Inlet: SphinxWriter will now document the signature of function callbacks added to a schema
- Axom::Path - New class for performing basic path operations with user-selectable
delimiter characters - Inlet: Added a method to
inlet::Inlet
that retrieves the set of unexpected names
in the input file - Inlet: Added an option to mark
Container
s as strict, which fail verification when unexpected
entries are present - Added support in
MFEMSidreDataCollection
for registeringQFunctions
(data associated with quadrature points on a mesh) - Inlet: The internal hierarchy of an
Inlet
object can be reconstructed from a Sidre group,
excluding callback functions - Added new overloaded version of method
sidre::DataStore::generateBlueprintIndex to incorporate new MPI
features in conduit and allow for generation of a blueprint index on
an under-decomposed parallel mesh - Added new method sidre::View::importArrayNode to import a
conduit::Node holding array data directly into a sidre::View - Added support for registering material and species sets in
MFEMSidreDataCollection
. These correspond to
matset
s
and
specset
s
in the Mesh Blueprint
Changed
- Converted [Uberenv] to a git submodule. We previously vendored a copy of this script.
- The Sidre Datastore no longer rewires Conduit's error handlers to SLIC by default.
It can be explicitly rewired using the static
DataStore::setConduitSLICMessageHandlers()
method. - Inlet: Changed
SchemaCreator
to an abstract class and added missing functions - Inlet: Added ability to access the
Reader
class fromInlet
and Sol Lua state
from theLuaReader
class - Inlet: Switched accessor interface to match that of the STL with
operator[]
and
T get<T>()
- Inlet:
std::shared_ptr<T>
has been replaced withT&
in non-owning contexts
andstd::unique_ptr<T>
in owning contexts - specifically, within Inlet's internal
tree structure - Unified core and SPIO unit tests into fewer executables to limit size of build directory
- Renamed
axom::slic::UnitTestLogger
toaxom::slic:SimpleLogger
because it's used in
more than just unit tests. - Inlet: Input file functions can now be of arbitrary signature subject to type and arity
restrictions - Exported all symbols on Windows by default when compiling a dynamic library
- Updated TPL
conduit
to version 0.6.0 released Nov 2, 2020. - Updated built-in TPL
sparsehash
to version 2.0.4 released Aug 11, 2020. - Inlet: Exposed
primal::Vector
in Lua for use in input-file-defined functions - The
MFEMSidreDataCollection
will now reconstruct fields and the mesh when a
datastore isLoad
ed in - Inlet: Cleaned up
Table
interface to eliminate ambiguity and duplicated functionality - Inlet: Renamed
DocWriter
toWriter
and refactored its interface - Inlet: Renamed
Table
toContainer
- Inlet collections of mixed or incorrect type will now fail verification, even if they're
not marked as required - Required Inlet collections no longer fail Inlet verification if they are empty in the input file
- Inlet:
operator bool
forField
andContainer
has been replaced with more preciseisUserProvided
andexists
, which also returnstrue
if a default value was specified. - Updated built-in TPL
fmt
to master branch snapshot, March 26, 2021. - Inlet:
SphinxWriter
will now print only one element schema per container instead of
printing the same schema for each element in the container - Updated BLT to version 0.4.0 released 9 Apr 2021
- Updated MFEM to version 4.2 released 30 Oct 2020. Axom no longer requires MFEM to be built serially
- The macro for exporting symbols is now
AXOM_EXPORT
instead ofAXOM_API
- Updated Conduit to v0.6.0
- Updated SCR to compatibility with v3.0rc1
Fixed
- Updated to new BLT version that does not fail when ClangFormat returns an empty
version string. BLT/Axom now issues a warning and disables thestyle
build
target if version is unknown or wrong. - Inlet: Apply lambda verifiers on generic containers to individual elements
for consistency - Inlet: Fixed a bug relating to nested table lookups of primitive arrays and functions
- Fixed a bug relating to deeply nested callback functions in Inlet
- Inlet: Always ignore primitive array elements that do not match the requested type
- Inlet: Empty structs/collections of structs with required sub-elements no longer fail
verification - Quest: Fixed a bug with InOutOctree for triangles that lie on faces of octree blocks
- Updated to use newer Conduit config directory
- Add support for legacy hdf5 cmake build system
Axom-v0.4.0
[Version 0.4.0] - Release date 2020-09-23
Please download the Axom-v0.4.0.tar.gz tarball below, which includes all of the Axom submodules as well.
Added
- Exposed the tolerance parameter
EPS
that is used to determine intersections between
triangles inprimal:intersect()
as an optional final parameter. - Added BVH spatial index option to the
mesh_tester
utility for calculating
triangle-triangle intersection. - Added
axom::execution_space< ExecSpace >::onDevice()
to check if execution
space is on device. - Added Axom macro
AXOM_SUPPRESS_HD_WARN
to silence host device compiler
warnings. - Added option to quest's
SignedDistance
class and C API to toggle whether
the distance query computes the sign. - Added a
batched
option to quest's signed distance query example application.
This computes all distance queries on an array of points using a single call tocomputeDistance
.
The query uses OpenMP threading, when available. - Added new component, Inlet, to assist in retrieving and storing data from
an input deck. - Added the ability to specify an [Umpire] allocator ID to use with the
BVH. This allows the application to use a device allocator for the BVH and
avoid use of Unified Memory (UM) on the GPU, which can hinder perfomrmance,
or use a pool allocator to mitigate the latencies associated with allocation/deallocation.
The allocator ID is specified as an optional argument to the BVH constructor. - Added new CMake option,
AXOM_ENABLE_ANNOTATIONS
, to enable/disable code
annotations in Axom. Default is OFF. - Added Axom annotation macros. The macros can be used to annotate functions,
using theAXOM_PERF_MARK_FUNCTION
macro, or at a more fine grain level,
different sections of code can be annotated by wrapping them within an
AXOM_PERF_MARK_SECTION
block. As a first cut, this works with NVTX tools.
However, the hooks are in place to add support for Caliper in the future. - Added a simple interface to NVTX that allows an application to set the color
and category for NVTX ranges corresponding to annotated code in Axom. The
application can now callaxom::nvtx:set_color()
and
axom::nvtx::set_category()
to set the corresponding parameters respectively.
This facilitates in the performance evaluation by allowing developers to easily
filter out calls by category or visually by setting a different color to use
in GUI tools, such as, NVVP and NSight. - Added a portable floating_point_limits traits class, to return min(), max(), lowest()
and epsilon() of afloat
ordouble
type. The functionality is equivalent to that provided by
std::numeric_limits, but, the code is host/device decorated accordingly such that it
can also be called on the device. - Added initial support for ray queries using the BVH. The caller may now supply a set of rays to
a BVH and the BVH will return a set of candidate BVH bins that intersect each ray. - Added initial support for bounding box queries using the BVH. The caller may
now supply a set of bounding boxes to a BVH and the BVH will return a set of
candidate BVH bins that intersect each bounding box. - Added an
axom-config.cmake
file to axom's installation to streamline incorporating axom
into user applications. See<axom-install>/examples/axom
for example usages. - Added [Sol] as a built-in TPL for fast and simple
C++
andLua
binding.
Sol is automatically enabled whenLUA_DIR
is found.
The version of Sol used in this release isv2.20.6
, which requiresC++14
.
Removed
- Removed the
AXOM_ENABLE_CUB
option, sinceCUB
is no lonher used directly in
Axom code. Instead, we useRAJA::stable_sort
with RAJA-v0.12.1 and fallback
tostd::stable_sort
with older versions of RAJA and when the code is built
without RAJA.
Changed
- Updated Axom to support RAJA-v0.12.1 and Umpire-v4.01, but the code remains
backwards compatible with previous versions of RAJA and Umpire. - Transitioned Axom's code formatting tool from
Uncrustify
to [clang-format].
Axom's clang-format rules depend on clang 10. - Modified the command line interface for
mesh_tester
utility. Interface
now uses a -m, --method option to select the spatial index, and -p, policy
option now accepts a string or integer value. - Renamed the
AXOM_USE_MPI3
option toAXOM_ENABLE_MPI3
for consistency. - Modified the API for the BVH to accomodate different query types. The queries are now
more explicitly calledBVH::findPoints()
andBVH::findRays()
. - Modified the API of Axom's memory management routines to not leak usage of Umpire. Instead of
passing anumpire::Allocator
object to specify an allocator, we now use the corresponding
integer ID associated with the allocator. - All names in the C API now preserve the case of the C++ function.
ex.SIDRE_datastore_new
is nowSIDRE_DataStore_new
. - Fortran API in slic module.
axom::slic::message
Level enums are changed
from enum-name_enumerator to namespace_enumerator.
ex.level_error
is nowmessage_error
. - Fortran derived-type constructors are now generic functions named afer the derived type.
datastore_new
is nowSidreDataStore
iomanager_new
is nowIOManager
Fixed
- Fixed a bug in
primal::intersect(Segment, BoundingBox)
and added regression tests. - Spin's octrees can now be used with 64-bit indexes. This allows octrees
with up to 64 levels of resolution when using a 64-bit index type. - Resolved issue with
AXOM_USE_64BIT_INDEXTYPE
configurations. Axom can once again
be configured with 64-bit index types. - Fixed a triangle-triangle intersection case in primal that produced inconsistent results
depending on the order of the triangle's vertices. - Fixed issue in the parallel construction of the BVH on GPUs, due to incoherent
L1 cache that could result in some data corruption in the BVH. The code now
calls__threadfence_system()
after the parent is computed and stored back
to global memory to ensure that the write is visible to all threads. - Fixed issue in Mint that would cause the [email protected] compiler to segfault. The
mint_cell_types.cpp
test was causing a segfault in the compiler. The main
issue triggering this compiler bug was the use ofconstexpr
when defining the
staticcell_info
array of structs. The code has been modified to useconst
instead. - Fixed issue in Quest's Signed Distance query that would prevent consecutive
calls to Quest when MPI-3 shared memory is enabled due to not properly
nullifying internal pointers when finalize is called. - Fixed issue where the BVH would dispatch to the CPU sort() routine when the
specified execution policy was CUDA_EXEC async. Now, when the execution policy
is CUDA_EXEC the code would correctly dispatch to the GPU sort, using CUB
(when CUB is enabled), regardless of whether it's synchronous or asynchronous. - Fixed issue with missing the bvh_traverse.hpp from the install prefix, which was preventing
applications from using the BVH when pointing to an Axom install prefix. - Fixed usage of cuda kernel policies in Mint. Raja v0.11.0 changed the way max threads
launch bounds is calculated. Consequently, a large number of threads was being launched
leading to max registry count violation when linking. We are now using fixed kernel size
of 256 threads (16x16 in 2D and 8x8x4 in 3D). - Third-party libraries can now build on the Windows platform through uberenv using vcpkg
("zero-to-axom support on Windows")
Known Bugs
- Encountered a compiler bug on IBM LC platforms when using the IBM XL C/C++
compiler. The issue is manifested in thegenerate_aabbs_and_centroids
method
in thespin_bvh.cpp
unit test. It seems that the compiler does not handle
the lambda capture of the arrays correctly which leads to a segfault. A
workaround for the IBM XL compiler is provided. - There is a known bug in MVAPICH that prevents consecutive creation/deletion
of MPI windows. This was encountered on LC platforms when enabling shared
memory in the Signed Distance Query. See the corresponding
Github Issue for details.