xref: /libCEED/doc/sphinx/source/releasenotes.md (revision 43e1b16f0c7c5dbfae04543c0675e255f265ff67)
1bcb2dfaeSJed Brown# Changes/Release Notes
2bcb2dfaeSJed Brown
3bcb2dfaeSJed BrownOn this page we provide a summary of the main API changes, new features and examples
4bcb2dfaeSJed Brownfor each release of libCEED.
5bcb2dfaeSJed Brown
6bcb2dfaeSJed Brown(main)=
7bcb2dfaeSJed Brown
8bcb2dfaeSJed Brown## Current `main` branch
9bcb2dfaeSJed Brown
107e7773b5SJeremy L Thompson### Interface changes
117e7773b5SJeremy L Thompson
127e7773b5SJeremy L Thompson- Update {c:func} `CeedQFunctionGetFields` and {c:func} `CeedOperatorGetFields` to include number of fields.
13ce4822f6SJeremy L Thompson- Promote to the public API: QFunction and Operator field objects, `CeedQFunctionField` and `CeedOperatorField`, and associated getters, {c:func}`CeedQFunctionGetFields`; {c:func}`CeedQFunctionFieldGetName`; {c:func}`CeedQFunctionFieldGetSize`; {c:func}`CeedQFunctionFieldGetEvalMode`; {c:func}`CeedOperatorGetFields`; {c:func}`CeedOperatorFieldGetElemRestriction`; {c:func}`CeedOperatorFieldGetBasis`; and {c:func}`CeedOperatorFieldGetVector`.
14f04ea552SJeremy L Thompson- Clarify and document conditions where `CeedQFunction` and `CeedOperator` become immutable and no further fields or suboperators can be added.
1570a7ffb3SJeremy L Thompson- Add {c:func} `CeedOperatorLinearAssembleQFunctionBuildOrUpdate` to reduce object creation overhead in assembly of CeedOperator preconditioning ingredients.
164db537f9SJeremy L Thompson- Promote {c:func} `CeedOperatorCheckReady`to the public API to facilitate interactive interfaces.
17dcc1e3ecSJeremy L Thompson- Warning added when compiling OCCA backend to alert users that this backend is experimental.
189a1d3511SJeremy L Thompson- `ceed-backend.h`, `ceed-hash.h`, and `ceed-khash.h` removed. Users should use `ceed/backend.h`, `ceed/hash.h`, and `ceed/khash.h`.
19*43e1b16fSJeremy L Thompson- Added {c:func} `CeedQFunctionGetKernelName`; refactored {c:func} `CeedQFunctionGetSourcePath` to exclude function kernel name.
207e7773b5SJeremy L Thompson
21f479eb23SJeremy L Thompson### New features
22f479eb23SJeremy L Thompson
23f479eb23SJeremy L Thompson- `CeedScalar` can now be set as `float` or `double` at compile time.
24f479eb23SJeremy L Thompson
25bcb2dfaeSJed Brown### Maintainability
26bcb2dfaeSJed Brown
27bcb2dfaeSJed Brown- Refactored preconditioner support internally to facilitate future development and improve GPU completeness/test coverage.
28db52d626SJeremy L Thompson- `Include-what-you-use` makefile target added as `make iwyu`.
29bcb2dfaeSJed Brown
30bcb2dfaeSJed Brown(v0-9)=
31bcb2dfaeSJed Brown
32bcb2dfaeSJed Brown## v0.9 (Jul 6, 2021)
33bcb2dfaeSJed Brown
34bcb2dfaeSJed Brown### Interface changes
35bcb2dfaeSJed Brown
36bcb2dfaeSJed Brown- Minor modification in error handling macro to silence pedantic warnings when compiling with Clang, but no functional impact.
37bcb2dfaeSJed Brown
38bcb2dfaeSJed Brown### New features
39bcb2dfaeSJed Brown
40bcb2dfaeSJed Brown- Add {c:func}`CeedVectorAXPY` and {c:func}`CeedVectorPointwiseMult` as a convenience for stand-alone testing and internal use.
41bcb2dfaeSJed Brown- Add `CEED_QFUNCTION_HELPER` macro to properly annotate QFunction helper functions for code generation backends.
42bcb2dfaeSJed Brown- Add `CeedPragmaOptimizeOff` macro for code that is sensitive to floating point errors from fast math optimizations.
43bcb2dfaeSJed Brown- Rust support: split `libceed-sys` crate out of `libceed` and [publish both on crates.io](https://crates.io/crates/libceed).
44bcb2dfaeSJed Brown
45bcb2dfaeSJed Brown### Performance improvements
46bcb2dfaeSJed Brown
47bcb2dfaeSJed Brown### Examples
48bcb2dfaeSJed Brown
49bcb2dfaeSJed Brown- Solid mechanics mini-app updated to explore the performance impacts of various formulations in the initial and current configurations.
50bcb2dfaeSJed Brown- Fluid mechanics example adds GPU support and improves modularity.
51bcb2dfaeSJed Brown
52bcb2dfaeSJed Brown### Deprecated backends
53bcb2dfaeSJed Brown
54bcb2dfaeSJed Brown- The `/cpu/self/tmpl` and `/cpu/self/tmpl/sub` backends have been removed. These backends were intially added to test the backend inheritance mechanism, but this mechanism is now widely used and tested in multiple backends.
55bcb2dfaeSJed Brown
56bcb2dfaeSJed Brown(v0-8)=
57bcb2dfaeSJed Brown
58bcb2dfaeSJed Brown## v0.8 (Mar 31, 2021)
59bcb2dfaeSJed Brown
60bcb2dfaeSJed Brown### Interface changes
61bcb2dfaeSJed Brown
62bcb2dfaeSJed Brown- Error handling improved to include enumerated error codes for C interface return values.
63bcb2dfaeSJed Brown- Installed headers that will follow semantic versioning were moved to {code}`include/ceed` directory. These headers have been renamed from {code}`ceed-*.h` to {code}`ceed/*.h`. Placeholder headers with the old naming schema are currently provided, but these headers will be removed in the libCEED v0.9 release.
64bcb2dfaeSJed Brown
65bcb2dfaeSJed Brown### New features
66bcb2dfaeSJed Brown
67bcb2dfaeSJed Brown- Julia and Rust interfaces added, providing a nearly 1-1 correspondence with the C interface, plus some convenience features.
68bcb2dfaeSJed Brown- Static libraries can be built with `make STATIC=1` and the pkg-config file is installed accordingly.
69bcb2dfaeSJed Brown- Add {c:func}`CeedOperatorLinearAssembleSymbolic` and {c:func}`CeedOperatorLinearAssemble` to support full assembly of libCEED operators.
70bcb2dfaeSJed Brown
71bcb2dfaeSJed Brown### Performance improvements
72bcb2dfaeSJed Brown
73bcb2dfaeSJed Brown- New HIP MAGMA backends for hipMAGMA library users: `/gpu/hip/magma` and `/gpu/hip/magma/det`.
74bcb2dfaeSJed Brown- New HIP backends for improved tensor basis performance: `/gpu/hip/shared` and `/gpu/hip/gen`.
75bcb2dfaeSJed Brown
76bcb2dfaeSJed Brown### Examples
77bcb2dfaeSJed Brown
78bcb2dfaeSJed Brown- {ref}`example-petsc-elasticity` example updated with traction boundary conditions and improved Dirichlet boundary conditions.
79bcb2dfaeSJed Brown- {ref}`example-petsc-elasticity` example updated with Neo-Hookean hyperelasticity in current configuration as well as improved Neo-Hookean hyperelasticity exploring storage vs computation tradeoffs.
80bcb2dfaeSJed Brown- {ref}`example-petsc-navier-stokes` example updated with isentropic traveling vortex test case, an analytical solution to the Euler equations that is useful for testing boundary conditions, discretization stability, and order of accuracy.
81bcb2dfaeSJed Brown- {ref}`example-petsc-navier-stokes` example updated with support for performing convergence study and plotting order of convergence by polynomial degree.
82bcb2dfaeSJed Brown
83bcb2dfaeSJed Brown(v0-7)=
84bcb2dfaeSJed Brown
85bcb2dfaeSJed Brown## v0.7 (Sep 29, 2020)
86bcb2dfaeSJed Brown
87bcb2dfaeSJed Brown### Interface changes
88bcb2dfaeSJed Brown
89bcb2dfaeSJed Brown- Replace limited {code}`CeedInterlaceMode` with more flexible component stride {code}`compstride` in {code}`CeedElemRestriction` constructors.
90bcb2dfaeSJed Brown  As a result, the {code}`indices` parameter has been replaced with {code}`offsets` and the {code}`nnodes` parameter has been replaced with {code}`lsize`.
91bcb2dfaeSJed Brown  These changes improve support for mixed finite element methods.
92bcb2dfaeSJed Brown- Replace various uses of {code}`Ceed*Get*Status` with {code}`Ceed*Is*` in the backend API to match common nomenclature.
93bcb2dfaeSJed Brown- Replace {code}`CeedOperatorAssembleLinearDiagonal` with {c:func}`CeedOperatorLinearAssembleDiagonal` for clarity.
94bcb2dfaeSJed Brown- Linear Operators can be assembled as point-block diagonal matrices with {c:func}`CeedOperatorLinearAssemblePointBlockDiagonal`, provided in row-major form in a {code}`ncomp` by {code}`ncomp` block per node.
95bcb2dfaeSJed Brown- Diagonal assemble interface changed to accept a {ref}`CeedVector` instead of a pointer to a {ref}`CeedVector` to reduce memory movement when interfacing with calling code.
96bcb2dfaeSJed Brown- Added {c:func}`CeedOperatorLinearAssembleAddDiagonal` and {c:func}`CeedOperatorLinearAssembleAddPointBlockDiagonal` for improved future integration with codes such as MFEM that compose the action of {ref}`CeedOperator`s external to libCEED.
97bcb2dfaeSJed Brown- Added {c:func}`CeedVectorTakeAray` to sync and remove libCEED read/write access to an allocated array and pass ownership of the array to the caller.
98bcb2dfaeSJed Brown  This function is recommended over {c:func}`CeedVectorSyncArray` when the {code}`CeedVector` has an array owned by the caller that was set by {c:func}`CeedVectorSetArray`.
99bcb2dfaeSJed Brown- Added {code}`CeedQFunctionContext` object to manage user QFunction context data and reduce copies between device and host memory.
100bcb2dfaeSJed Brown- Added {c:func}`CeedOperatorMultigridLevelCreate`, {c:func}`CeedOperatorMultigridLevelCreateTensorH1`, and {c:func}`CeedOperatorMultigridLevelCreateH1` to facilitate creation of multigrid prolongation, restriction, and coarse grid operators using a common quadrature space.
101bcb2dfaeSJed Brown
102bcb2dfaeSJed Brown### New features
103bcb2dfaeSJed Brown
104bcb2dfaeSJed Brown- New HIP backend: `/gpu/hip/ref`.
105bcb2dfaeSJed Brown- CeedQFunction support for user `CUfunction`s in some backends
106bcb2dfaeSJed Brown
107bcb2dfaeSJed Brown### Performance improvements
108bcb2dfaeSJed Brown
109bcb2dfaeSJed Brown- OCCA backend rebuilt to facilitate future performance enhancements.
110bcb2dfaeSJed Brown- Petsc BPs suite improved to reduce noise due to multiple calls to {code}`mpiexec`.
111bcb2dfaeSJed Brown
112bcb2dfaeSJed Brown### Examples
113bcb2dfaeSJed Brown
114bcb2dfaeSJed Brown- {ref}`example-petsc-elasticity` example updated with strain energy computation and more flexible boundary conditions.
115bcb2dfaeSJed Brown
116bcb2dfaeSJed Brown### Deprecated backends
117bcb2dfaeSJed Brown
118bcb2dfaeSJed Brown- The `/gpu/cuda/reg` backend has been removed, with its core features moved into `/gpu/cuda/ref` and `/gpu/cuda/shared`.
119bcb2dfaeSJed Brown
120bcb2dfaeSJed Brown(v0-6)=
121bcb2dfaeSJed Brown
122bcb2dfaeSJed Brown## v0.6 (Mar 29, 2020)
123bcb2dfaeSJed Brown
124bcb2dfaeSJed BrownlibCEED v0.6 contains numerous new features and examples, as well as expanded
125bcb2dfaeSJed Browndocumentation in [this new website](https://libceed.readthedocs.io).
126bcb2dfaeSJed Brown
127bcb2dfaeSJed Brown### New features
128bcb2dfaeSJed Brown
129bcb2dfaeSJed Brown- New Python interface using [CFFI](https://cffi.readthedocs.io/) provides a nearly
130bcb2dfaeSJed Brown  1-1 correspondence with the C interface, plus some convenience features.  For instance,
131bcb2dfaeSJed Brown  data stored in the {cpp:type}`CeedVector` structure are available without copy as
132bcb2dfaeSJed Brown  {py:class}`numpy.ndarray`.  Short tutorials are provided in
133bcb2dfaeSJed Brown  [Binder](https://mybinder.org/v2/gh/CEED/libCEED/main?urlpath=lab/tree/examples/tutorials/).
134bcb2dfaeSJed Brown- Linear QFunctions can be assembled as block-diagonal matrices (per quadrature point,
135bcb2dfaeSJed Brown  {c:func}`CeedOperatorAssembleLinearQFunction`) or to evaluate the diagonal
136bcb2dfaeSJed Brown  ({c:func}`CeedOperatorAssembleLinearDiagonal`).  These operations are useful for
137bcb2dfaeSJed Brown  preconditioning ingredients and are used in the libCEED's multigrid examples.
138bcb2dfaeSJed Brown- The inverse of separable operators can be obtained using
139bcb2dfaeSJed Brown  {c:func}`CeedOperatorCreateFDMElementInverse` and applied with
140bcb2dfaeSJed Brown  {c:func}`CeedOperatorApply`.  This is a useful preconditioning ingredient,
141bcb2dfaeSJed Brown  especially for Laplacians and related operators.
142bcb2dfaeSJed Brown- New functions: {c:func}`CeedVectorNorm`, {c:func}`CeedOperatorApplyAdd`,
143bcb2dfaeSJed Brown  {c:func}`CeedQFunctionView`, {c:func}`CeedOperatorView`.
144bcb2dfaeSJed Brown- Make public accessors for various attributes to facilitate writing composable code.
145bcb2dfaeSJed Brown- New backend: `/cpu/self/memcheck/serial`.
146bcb2dfaeSJed Brown- QFunctions using variable-length array (VLA) pointer constructs can be used with CUDA
147bcb2dfaeSJed Brown  backends.  (Single source is coming soon for OCCA backends.)
148bcb2dfaeSJed Brown- Fix some missing edge cases in CUDA backend.
149bcb2dfaeSJed Brown
150bcb2dfaeSJed Brown### Performance Improvements
151bcb2dfaeSJed Brown
152bcb2dfaeSJed Brown- MAGMA backend performance optimization and non-tensor bases.
153bcb2dfaeSJed Brown- No-copy optimization in {c:func}`CeedOperatorApply`.
154bcb2dfaeSJed Brown
155bcb2dfaeSJed Brown### Interface changes
156bcb2dfaeSJed Brown
157bcb2dfaeSJed Brown- Replace {code}`CeedElemRestrictionCreateIdentity` and
158bcb2dfaeSJed Brown  {code}`CeedElemRestrictionCreateBlocked` with more flexible
159bcb2dfaeSJed Brown  {c:func}`CeedElemRestrictionCreateStrided` and
160bcb2dfaeSJed Brown  {c:func}`CeedElemRestrictionCreateBlockedStrided`.
161bcb2dfaeSJed Brown- Add arguments to {c:func}`CeedQFunctionCreateIdentity`.
162bcb2dfaeSJed Brown- Replace ambiguous uses of {cpp:enum}`CeedTransposeMode` for L-vector identification
163bcb2dfaeSJed Brown  with {cpp:enum}`CeedInterlaceMode`.  This is now an attribute of the
164bcb2dfaeSJed Brown  {cpp:type}`CeedElemRestriction` (see {c:func}`CeedElemRestrictionCreate`) and no
165bcb2dfaeSJed Brown  longer passed as `lmode` arguments to {c:func}`CeedOperatorSetField` and
166bcb2dfaeSJed Brown  {c:func}`CeedElemRestrictionApply`.
167bcb2dfaeSJed Brown
168bcb2dfaeSJed Brown### Examples
169bcb2dfaeSJed Brown
170bcb2dfaeSJed BrownlibCEED-0.6 contains greatly expanded examples with {ref}`new documentation <Examples>`.
171bcb2dfaeSJed BrownNotable additions include:
172bcb2dfaeSJed Brown
173bcb2dfaeSJed Brown- Standalone {ref}`ex2-surface` ({file}`examples/ceed/ex2-surface`): compute the area of
174bcb2dfaeSJed Brown  a domain in 1, 2, and 3 dimensions by applying a Laplacian.
175bcb2dfaeSJed Brown
176bcb2dfaeSJed Brown- PETSc {ref}`example-petsc-area` ({file}`examples/petsc/area.c`): computes surface area
177bcb2dfaeSJed Brown  of domains (like the cube and sphere) by direct integration on a surface mesh;
178bcb2dfaeSJed Brown  demonstrates geometric dimension different from topological dimension.
179bcb2dfaeSJed Brown
180bcb2dfaeSJed Brown- PETSc {ref}`example-petsc-bps`:
181bcb2dfaeSJed Brown
182bcb2dfaeSJed Brown  - {file}`examples/petsc/bpsraw.c` (formerly `bps.c`): transparent CUDA support.
183bcb2dfaeSJed Brown  - {file}`examples/petsc/bps.c` (formerly `bpsdmplex.c`): performance improvements
184bcb2dfaeSJed Brown    and transparent CUDA support.
185bcb2dfaeSJed Brown  - {ref}`example-petsc-bps-sphere` ({file}`examples/petsc/bpssphere.c`):
186bcb2dfaeSJed Brown    generalizations of all CEED BPs to the surface of the sphere; demonstrates geometric
187bcb2dfaeSJed Brown    dimension different from topological dimension.
188bcb2dfaeSJed Brown
189bcb2dfaeSJed Brown- {ref}`example-petsc-multigrid` ({file}`examples/petsc/multigrid.c`): new p-multigrid
190bcb2dfaeSJed Brown  solver with algebraic multigrid coarse solve.
191bcb2dfaeSJed Brown
192bcb2dfaeSJed Brown- {ref}`example-petsc-navier-stokes` ({file}`examples/fluids/navierstokes.c`; formerly
193bcb2dfaeSJed Brown  `examples/navier-stokes`): unstructured grid support (using PETSc's `DMPlex`),
194bcb2dfaeSJed Brown  implicit time integration, SU/SUPG stabilization, free-slip boundary conditions, and
195bcb2dfaeSJed Brown  quasi-2D computational domain support.
196bcb2dfaeSJed Brown
197bcb2dfaeSJed Brown- {ref}`example-petsc-elasticity` ({file}`examples/solids/elasticity.c`): new solver for
198bcb2dfaeSJed Brown  linear elasticity, small-strain hyperelasticity, and globalized finite-strain
199bcb2dfaeSJed Brown  hyperelasticity using p-multigrid with algebraic multigrid coarse solve.
200bcb2dfaeSJed Brown
201bcb2dfaeSJed Brown(v0-5)=
202bcb2dfaeSJed Brown
203bcb2dfaeSJed Brown## v0.5 (Sep 18, 2019)
204bcb2dfaeSJed Brown
205bcb2dfaeSJed BrownFor this release, several improvements were made. Two new CUDA backends were added to
206bcb2dfaeSJed Brownthe family of backends, of which, the new `cuda-gen` backend achieves state-of-the-art
207bcb2dfaeSJed Brownperformance using single-source {ref}`CeedQFunction`. From this release, users
208bcb2dfaeSJed Browncan define Q-Functions in a single source code independently of the targeted backend
209bcb2dfaeSJed Brownwith the aid of a new macro `CEED QFUNCTION` to support JIT (Just-In-Time) and CPU
210bcb2dfaeSJed Browncompilation of the user provided {ref}`CeedQFunction` code. To allow a unified
211bcb2dfaeSJed Browndeclaration, the {ref}`CeedQFunction` API has undergone a slight change:
212bcb2dfaeSJed Brownthe `QFunctionField` parameter `ncomp` has been changed to `size`. This change
213bcb2dfaeSJed Brownrequires setting the previous value of `ncomp` to `ncomp*dim` when adding a
214bcb2dfaeSJed Brown`QFunctionField` with eval mode `CEED EVAL GRAD`.
215bcb2dfaeSJed Brown
216bcb2dfaeSJed BrownAdditionally, new CPU backends
217bcb2dfaeSJed Brownwere included in this release, such as the `/cpu/self/opt/*` backends (which are
218bcb2dfaeSJed Brownwritten in pure C and use partial **E-vectors** to improve performance) and the
219bcb2dfaeSJed Brown`/cpu/self/ref/memcheck` backend (which relies upon the
220bcb2dfaeSJed Brown[Valgrind](http://valgrind.org/) Memcheck tool to help verify that user
221bcb2dfaeSJed Brown{ref}`CeedQFunction` have no undefined values).
222bcb2dfaeSJed BrownThis release also included various performance improvements, bug fixes, new examples,
223bcb2dfaeSJed Brownand improved tests. Among these improvements, vectorized instructions for
224bcb2dfaeSJed Brown{ref}`CeedQFunction` code compiled for CPU were enhanced by using `CeedPragmaSIMD`
225bcb2dfaeSJed Browninstead of `CeedPragmaOMP`, implementation of a {ref}`CeedQFunction` gallery and
226bcb2dfaeSJed Brownidentity Q-Functions were introduced, and the PETSc benchmark problems were expanded
227bcb2dfaeSJed Brownto include unstructured meshes handling were. For this expansion, the prior version of
228bcb2dfaeSJed Brownthe PETSc BPs, which only included data associated with structured geometries, were
229bcb2dfaeSJed Brownrenamed `bpsraw`, and the new version of the BPs, which can handle data associated
230bcb2dfaeSJed Brownwith any unstructured geometry, were called `bps`. Additionally, other benchmark
231bcb2dfaeSJed Brownproblems, namely BP2 and BP4 (the vector-valued versions of BP1 and BP3, respectively),
232bcb2dfaeSJed Brownand BP5 and BP6 (the collocated versions---for which the quadrature points are the same
233bcb2dfaeSJed Brownas the Gauss Lobatto nodes---of BP3 and BP4 respectively) were added to the PETSc
234bcb2dfaeSJed Brownexamples. Furthermoew, another standalone libCEED example, called `ex2`, which
235bcb2dfaeSJed Browncomputes the surface area of a given mesh was added to this release.
236bcb2dfaeSJed Brown
237bcb2dfaeSJed BrownBackends available in this release:
238bcb2dfaeSJed Brown
23968e843eeSJed Brown| CEED resource (`-ceed`)  | Backend                                             |
24068e843eeSJed Brown|--------------------------|-----------------------------------------------------|
24168e843eeSJed Brown| `/cpu/self/ref/serial`   | Serial reference implementation                     |
24268e843eeSJed Brown| `/cpu/self/ref/blocked`  | Blocked reference implementation                    |
24368e843eeSJed Brown| `/cpu/self/ref/memcheck` | Memcheck backend, undefined value checks            |
24468e843eeSJed Brown| `/cpu/self/opt/serial`   | Serial optimized C implementation                   |
24568e843eeSJed Brown| `/cpu/self/opt/blocked`  | Blocked optimized C implementation                  |
24668e843eeSJed Brown| `/cpu/self/avx/serial`   | Serial AVX implementation                           |
24768e843eeSJed Brown| `/cpu/self/avx/blocked`  | Blocked AVX implementation                          |
24868e843eeSJed Brown| `/cpu/self/xsmm/serial`  | Serial LIBXSMM implementation                       |
24968e843eeSJed Brown| `/cpu/self/xsmm/blocked` | Blocked LIBXSMM implementation                      |
25068e843eeSJed Brown| `/cpu/occa`              | Serial OCCA kernels                                 |
25168e843eeSJed Brown| `/gpu/occa`              | CUDA OCCA kernels                                   |
25268e843eeSJed Brown| `/omp/occa`              | OpenMP OCCA kernels                                 |
25368e843eeSJed Brown| `/ocl/occa`              | OpenCL OCCA kernels                                 |
25468e843eeSJed Brown| `/gpu/cuda/ref`          | Reference pure CUDA kernels                         |
25568e843eeSJed Brown| `/gpu/cuda/reg`          | Pure CUDA kernels using one thread per element      |
25668e843eeSJed Brown| `/gpu/cuda/shared`       | Optimized pure CUDA kernels using shared memory     |
25768e843eeSJed Brown| `/gpu/cuda/gen`          | Optimized pure CUDA kernels using code generation   |
25868e843eeSJed Brown| `/gpu/magma`             | CUDA MAGMA kernels                                  |
259bcb2dfaeSJed Brown
260bcb2dfaeSJed BrownExamples available in this release:
261bcb2dfaeSJed Brown
26268e843eeSJed Brown:::{list-table}
26368e843eeSJed Brown:header-rows: 1
26468e843eeSJed Brown:widths: auto
26568e843eeSJed Brown* - User code
26668e843eeSJed Brown  - Example
26768e843eeSJed Brown* - `ceed`
26868e843eeSJed Brown  - * ex1 (volume)
26968e843eeSJed Brown    * ex2 (surface)
27068e843eeSJed Brown* - `mfem`
27168e843eeSJed Brown  - * BP1 (scalar mass operator)
27268e843eeSJed Brown    * BP3 (scalar Laplace operator)
27368e843eeSJed Brown* - `petsc`
27468e843eeSJed Brown  - * BP1 (scalar mass operator)
27568e843eeSJed Brown    * BP2 (vector mass operator)
27668e843eeSJed Brown    * BP3 (scalar Laplace operator)
27768e843eeSJed Brown    * BP4 (vector Laplace operator)
27868e843eeSJed Brown    * BP5 (collocated scalar Laplace operator)
27968e843eeSJed Brown    * BP6 (collocated vector Laplace operator)
28068e843eeSJed Brown    * Navier-Stokes
28168e843eeSJed Brown* - `nek5000`
28268e843eeSJed Brown  - * BP1 (scalar mass operator)
28368e843eeSJed Brown    * BP3 (scalar Laplace operator)
28468e843eeSJed Brown:::
285bcb2dfaeSJed Brown
286bcb2dfaeSJed Brown(v0-4)=
287bcb2dfaeSJed Brown
288bcb2dfaeSJed Brown## v0.4 (Apr 1, 2019)
289bcb2dfaeSJed Brown
290bcb2dfaeSJed BrownlibCEED v0.4 was made again publicly available in the second full CEED software
291bcb2dfaeSJed Browndistribution, release CEED 2.0. This release contained notable features, such as
292bcb2dfaeSJed Brownfour new CPU backends, two new GPU backends, CPU backend optimizations, initial
293bcb2dfaeSJed Brownsupport for operator composition, performance benchmarking, and a Navier-Stokes demo.
294bcb2dfaeSJed BrownThe new CPU backends in this release came in two families. The `/cpu/self/*/serial`
295bcb2dfaeSJed Brownbackends process one element at a time and are intended for meshes with a smaller number
296bcb2dfaeSJed Brownof high order elements. The `/cpu/self/*/blocked` backends process blocked batches of
297bcb2dfaeSJed Browneight interlaced elements and are intended for meshes with higher numbers of elements.
298bcb2dfaeSJed BrownThe `/cpu/self/avx/*` backends rely upon AVX instructions to provide vectorized CPU
299bcb2dfaeSJed Brownperformance. The `/cpu/self/xsmm/*` backends rely upon the
300bcb2dfaeSJed Brown[LIBXSMM](http://github.com/hfp/libxsmm) package to provide vectorized CPU
301bcb2dfaeSJed Brownperformance. The `/gpu/cuda/*` backends provide GPU performance strictly using CUDA.
302bcb2dfaeSJed BrownThe `/gpu/cuda/ref` backend is a reference CUDA backend, providing reasonable
303bcb2dfaeSJed Brownperformance for most problem configurations. The `/gpu/cuda/reg` backend uses a simple
304bcb2dfaeSJed Brownparallelization approach, where each thread treats a finite element. Using just in time
305bcb2dfaeSJed Browncompilation, provided by nvrtc (NVidia Runtime Compiler), and runtime parameters, this
306bcb2dfaeSJed Brownbackend unroll loops and map memory address to registers. The `/gpu/cuda/reg` backend
307bcb2dfaeSJed Brownachieve good peak performance for 1D, 2D, and low order 3D problems, but performance
308bcb2dfaeSJed Browndeteriorates very quickly when threads run out of registers.
309bcb2dfaeSJed Brown
310bcb2dfaeSJed BrownA new explicit time-stepping Navier-Stokes solver was added to the family of libCEED
311bcb2dfaeSJed Brownexamples in the `examples/petsc` directory (see {ref}`example-petsc-navier-stokes`).
312bcb2dfaeSJed BrownThis example solves the time-dependent Navier-Stokes equations of compressible gas
313bcb2dfaeSJed Browndynamics in a static Eulerian three-dimensional frame, using structured high-order
314bcb2dfaeSJed Brownfinite/spectral element spatial discretizations and explicit high-order time-stepping
315bcb2dfaeSJed Brown(available in PETSc). Moreover, the Navier-Stokes example was developed using PETSc,
316bcb2dfaeSJed Brownso that the pointwise physics (defined at quadrature points) is separated from the
317bcb2dfaeSJed Brownparallelization and meshing concerns.
318bcb2dfaeSJed Brown
319bcb2dfaeSJed BrownBackends available in this release:
320bcb2dfaeSJed Brown
32168e843eeSJed Brown| CEED resource (`-ceed`)  | Backend                                             |
32268e843eeSJed Brown|--------------------------|-----------------------------------------------------|
32368e843eeSJed Brown| `/cpu/self/ref/serial`   | Serial reference implementation                     |
32468e843eeSJed Brown| `/cpu/self/ref/blocked`  | Blocked reference implementation                    |
32568e843eeSJed Brown| `/cpu/self/tmpl`         | Backend template, defaults to `/cpu/self/blocked`   |
32668e843eeSJed Brown| `/cpu/self/avx/serial`   | Serial AVX implementation                           |
32768e843eeSJed Brown| `/cpu/self/avx/blocked`  | Blocked AVX implementation                          |
32868e843eeSJed Brown| `/cpu/self/xsmm/serial`  | Serial LIBXSMM implementation                       |
32968e843eeSJed Brown| `/cpu/self/xsmm/blocked` | Blocked LIBXSMM implementation                      |
33068e843eeSJed Brown| `/cpu/occa`              | Serial OCCA kernels                                 |
33168e843eeSJed Brown| `/gpu/occa`              | CUDA OCCA kernels                                   |
33268e843eeSJed Brown| `/omp/occa`              | OpenMP OCCA kernels                                 |
33368e843eeSJed Brown| `/ocl/occa`              | OpenCL OCCA kernels                                 |
33468e843eeSJed Brown| `/gpu/cuda/ref`          | Reference pure CUDA kernels                         |
33568e843eeSJed Brown| `/gpu/cuda/reg`          | Pure CUDA kernels using one thread per element      |
33668e843eeSJed Brown| `/gpu/magma`             | CUDA MAGMA kernels                                  |
337bcb2dfaeSJed Brown
338bcb2dfaeSJed BrownExamples available in this release:
339bcb2dfaeSJed Brown
34068e843eeSJed Brown:::{list-table}
34168e843eeSJed Brown:header-rows: 1
34268e843eeSJed Brown:widths: auto
34368e843eeSJed Brown* - User code
34468e843eeSJed Brown  - Example
34568e843eeSJed Brown* - `ceed`
34668e843eeSJed Brown  - * ex1 (volume)
34768e843eeSJed Brown* - `mfem`
34868e843eeSJed Brown  - * BP1 (scalar mass operator)
34968e843eeSJed Brown    * BP3 (scalar Laplace operator)
35068e843eeSJed Brown* - `petsc`
35168e843eeSJed Brown  - * BP1 (scalar mass operator)
35268e843eeSJed Brown    * BP3 (scalar Laplace operator)
35368e843eeSJed Brown    * Navier-Stokes
35468e843eeSJed Brown* - `nek5000`
35568e843eeSJed Brown  - * BP1 (scalar mass operator)
35668e843eeSJed Brown    * BP3 (scalar Laplace operator)
35768e843eeSJed Brown:::
358bcb2dfaeSJed Brown
359bcb2dfaeSJed Brown(v0-3)=
360bcb2dfaeSJed Brown
361bcb2dfaeSJed Brown## v0.3 (Sep 30, 2018)
362bcb2dfaeSJed Brown
363bcb2dfaeSJed BrownNotable features in this release include active/passive field interface, support for
364bcb2dfaeSJed Brownnon-tensor bases, backend optimization, and improved Fortran interface. This release
365bcb2dfaeSJed Brownalso focused on providing improved continuous integration, and many new tests with code
366bcb2dfaeSJed Browncoverage reports of about 90%. This release also provided a significant change to the
367bcb2dfaeSJed Brownpublic interface: a {ref}`CeedQFunction` can take any number of named input and output
368bcb2dfaeSJed Brownarguments while {ref}`CeedOperator` connects them to the actual data, which may be
369bcb2dfaeSJed Brownsupplied explicitly to `CeedOperatorApply()` (active) or separately via
370bcb2dfaeSJed Brown`CeedOperatorSetField()` (passive). This interface change enables reusable libraries
371bcb2dfaeSJed Brownof CeedQFunctions and composition of block solvers constructed using
372bcb2dfaeSJed Brown{ref}`CeedOperator`. A concept of blocked restriction was added to this release and
373bcb2dfaeSJed Brownused in an optimized CPU backend. Although this is typically not visible to the user,
374bcb2dfaeSJed Brownit enables effective use of arbitrary-length SIMD while maintaining cache locality.
375bcb2dfaeSJed BrownThis CPU backend also implements an algebraic factorization of tensor product gradients
376bcb2dfaeSJed Brownto perform fewer operations than standard application of interpolation and
377bcb2dfaeSJed Browndifferentiation from nodes to quadrature points. This algebraic formulation
378bcb2dfaeSJed Brownautomatically supports non-polynomial and non-interpolatory bases, thus is more general
379bcb2dfaeSJed Brownthan the more common derivation in terms of Lagrange polynomials on the quadrature points.
380bcb2dfaeSJed Brown
381bcb2dfaeSJed BrownBackends available in this release:
382bcb2dfaeSJed Brown
38368e843eeSJed Brown| CEED resource (`-ceed`) | Backend                                             |
38468e843eeSJed Brown|-------------------------|-----------------------------------------------------|
38568e843eeSJed Brown| `/cpu/self/blocked`     | Blocked reference implementation                    |
38668e843eeSJed Brown| `/cpu/self/ref`         | Serial reference implementation                     |
38768e843eeSJed Brown| `/cpu/self/tmpl`        | Backend template, defaults to `/cpu/self/blocked`   |
38868e843eeSJed Brown| `/cpu/occa`             | Serial OCCA kernels                                 |
38968e843eeSJed Brown| `/gpu/occa`             | CUDA OCCA kernels                                   |
39068e843eeSJed Brown| `/omp/occa`             | OpenMP OCCA kernels                                 |
39168e843eeSJed Brown| `/ocl/occa`             | OpenCL OCCA kernels                                 |
39268e843eeSJed Brown| `/gpu/magma`            | CUDA MAGMA kernels                                  |
393bcb2dfaeSJed Brown
394bcb2dfaeSJed BrownExamples available in this release:
395bcb2dfaeSJed Brown
39668e843eeSJed Brown:::{list-table}
39768e843eeSJed Brown:header-rows: 1
39868e843eeSJed Brown:widths: auto
39968e843eeSJed Brown* - User code
40068e843eeSJed Brown  - Example
40168e843eeSJed Brown* - `ceed`
40268e843eeSJed Brown  - * ex1 (volume)
40368e843eeSJed Brown* - `mfem`
40468e843eeSJed Brown  - * BP1 (scalar mass operator)
40568e843eeSJed Brown    * BP3 (scalar Laplace operator)
40668e843eeSJed Brown* - `petsc`
40768e843eeSJed Brown  - * BP1 (scalar mass operator)
40868e843eeSJed Brown    * BP3 (scalar Laplace operator)
40968e843eeSJed Brown* - `nek5000`
41068e843eeSJed Brown  - * BP1 (scalar mass operator)
41168e843eeSJed Brown    * BP3 (scalar Laplace operator)
41268e843eeSJed Brown:::
413bcb2dfaeSJed Brown
414bcb2dfaeSJed Brown(v0-21)=
415bcb2dfaeSJed Brown
416bcb2dfaeSJed Brown## v0.21 (Sep 30, 2018)
417bcb2dfaeSJed Brown
418bcb2dfaeSJed BrownA MAGMA backend (which relies upon the
419bcb2dfaeSJed Brown[MAGMA](https://bitbucket.org/icl/magma) package) was integrated in libCEED for this
420bcb2dfaeSJed Brownrelease. This initial integration set up the framework of using MAGMA and provided the
421bcb2dfaeSJed BrownlibCEED functionality through MAGMA kernels as one of libCEED’s computational backends.
422bcb2dfaeSJed BrownAs any other backend, the MAGMA backend provides extended basic data structures for
423bcb2dfaeSJed Brown{ref}`CeedVector`, {ref}`CeedElemRestriction`, and {ref}`CeedOperator`, and implements
424bcb2dfaeSJed Brownthe fundamental CEED building blocks to work with the new data structures.
425bcb2dfaeSJed BrownIn general, the MAGMA-specific data structures keep the libCEED pointers to CPU data
426bcb2dfaeSJed Brownbut also add corresponding device (e.g., GPU) pointers to the data. Coherency is handled
427bcb2dfaeSJed Browninternally, and thus seamlessly to the user, through the functions/methods that are
428bcb2dfaeSJed Brownprovided to support them.
429bcb2dfaeSJed Brown
430bcb2dfaeSJed BrownBackends available in this release:
431bcb2dfaeSJed Brown
43268e843eeSJed Brown| CEED resource (`-ceed`) | Backend                         |
43368e843eeSJed Brown|-------------------------|---------------------------------|
43468e843eeSJed Brown| `/cpu/self`             | Serial reference implementation |
43568e843eeSJed Brown| `/cpu/occa`             | Serial OCCA kernels             |
43668e843eeSJed Brown| `/gpu/occa`             | CUDA OCCA kernels               |
43768e843eeSJed Brown| `/omp/occa`             | OpenMP OCCA kernels             |
43868e843eeSJed Brown| `/ocl/occa`             | OpenCL OCCA kernels             |
43968e843eeSJed Brown| `/gpu/magma`            | CUDA MAGMA kernels              |
440bcb2dfaeSJed Brown
441bcb2dfaeSJed BrownExamples available in this release:
442bcb2dfaeSJed Brown
44368e843eeSJed Brown:::{list-table}
44468e843eeSJed Brown:header-rows: 1
44568e843eeSJed Brown:widths: auto
44668e843eeSJed Brown* - User code
44768e843eeSJed Brown  - Example
44868e843eeSJed Brown* - `ceed`
44968e843eeSJed Brown  - * ex1 (volume)
45068e843eeSJed Brown* - `mfem`
45168e843eeSJed Brown  - * BP1 (scalar mass operator)
45268e843eeSJed Brown    * BP3 (scalar Laplace operator)
45368e843eeSJed Brown* - `petsc`
45468e843eeSJed Brown  - * BP1 (scalar mass operator)
45568e843eeSJed Brown* - `nek5000`
45668e843eeSJed Brown  - * BP1 (scalar mass operator)
45768e843eeSJed Brown:::
458bcb2dfaeSJed Brown
459bcb2dfaeSJed Brown(v0-2)=
460bcb2dfaeSJed Brown
461bcb2dfaeSJed Brown## v0.2 (Mar 30, 2018)
462bcb2dfaeSJed Brown
463bcb2dfaeSJed BrownlibCEED was made publicly available the first full CEED software distribution, release
464bcb2dfaeSJed BrownCEED 1.0. The distribution was made available using the Spack package manager to provide
465bcb2dfaeSJed Browna common, easy-to-use build environment, where the user can build the CEED distribution
466bcb2dfaeSJed Brownwith all dependencies. This release included a new Fortran interface for the library.
467bcb2dfaeSJed BrownThis release also contained major improvements in the OCCA backend (including a new
468bcb2dfaeSJed Brown`/ocl/occa` backend) and new examples. The standalone libCEED example was modified to
469bcb2dfaeSJed Browncompute the volume volume of a given mesh (in 1D, 2D, or 3D) and placed in an
470bcb2dfaeSJed Brown`examples/ceed` subfolder. A new `mfem` example to perform BP3 (with the application
471bcb2dfaeSJed Brownof the Laplace operator) was also added to this release.
472bcb2dfaeSJed Brown
473bcb2dfaeSJed BrownBackends available in this release:
474bcb2dfaeSJed Brown
47568e843eeSJed Brown| CEED resource (`-ceed`) | Backend                         |
47668e843eeSJed Brown|-------------------------|---------------------------------|
47768e843eeSJed Brown| `/cpu/self`             | Serial reference implementation |
47868e843eeSJed Brown| `/cpu/occa`             | Serial OCCA kernels             |
47968e843eeSJed Brown| `/gpu/occa`             | CUDA OCCA kernels               |
48068e843eeSJed Brown| `/omp/occa`             | OpenMP OCCA kernels             |
48168e843eeSJed Brown| `/ocl/occa`             | OpenCL OCCA kernels             |
482bcb2dfaeSJed Brown
483bcb2dfaeSJed BrownExamples available in this release:
484bcb2dfaeSJed Brown
48568e843eeSJed Brown:::{list-table}
48668e843eeSJed Brown:header-rows: 1
48768e843eeSJed Brown:widths: auto
48868e843eeSJed Brown* - User code
48968e843eeSJed Brown  - Example
49068e843eeSJed Brown* - `ceed`
49168e843eeSJed Brown  - * ex1 (volume)
49268e843eeSJed Brown* - `mfem`
49368e843eeSJed Brown  - * BP1 (scalar mass operator)
49468e843eeSJed Brown    * BP3 (scalar Laplace operator)
49568e843eeSJed Brown* - `petsc`
49668e843eeSJed Brown  - * BP1 (scalar mass operator)
49768e843eeSJed Brown* - `nek5000`
49868e843eeSJed Brown  - * BP1 (scalar mass operator)
49968e843eeSJed Brown:::
500bcb2dfaeSJed Brown
501bcb2dfaeSJed Brown(v0-1)=
502bcb2dfaeSJed Brown
503bcb2dfaeSJed Brown## v0.1 (Jan 3, 2018)
504bcb2dfaeSJed Brown
505bcb2dfaeSJed BrownInitial low-level API of the CEED project. The low-level API provides a set of Finite
506bcb2dfaeSJed BrownElements kernels and components for writing new low-level kernels. Examples include:
507bcb2dfaeSJed Brownvector and sparse linear algebra, element matrix assembly over a batch of elements,
508bcb2dfaeSJed Brownpartial assembly and action for efficient high-order operators like mass, diffusion,
509bcb2dfaeSJed Brownadvection, etc. The main goal of the low-level API is to establish the basis for the
510bcb2dfaeSJed Brownhigh-level API. Also, identifying such low-level kernels and providing a reference
511bcb2dfaeSJed Brownimplementation for them serves as the basis for specialized backend implementations.
512bcb2dfaeSJed BrownThis release contained several backends: `/cpu/self`, and backends which rely upon the
513bcb2dfaeSJed Brown[OCCA](http://github.com/libocca/occa) package, such as `/cpu/occa`,
514bcb2dfaeSJed Brown`/gpu/occa`, and `/omp/occa`.
515bcb2dfaeSJed BrownIt also included several examples, in the `examples` folder:
516bcb2dfaeSJed BrownA standalone code that shows the usage of libCEED (with no external
517bcb2dfaeSJed Browndependencies) to apply the Laplace operator, `ex1`; an `mfem` example to perform BP1
518bcb2dfaeSJed Brown(with the application of the mass operator); and a `petsc` example to perform BP1
519bcb2dfaeSJed Brown(with the application of the mass operator).
520bcb2dfaeSJed Brown
521bcb2dfaeSJed BrownBackends available in this release:
522bcb2dfaeSJed Brown
52368e843eeSJed Brown| CEED resource (`-ceed`) | Backend                         |
52468e843eeSJed Brown|-------------------------|---------------------------------|
52568e843eeSJed Brown| `/cpu/self`             | Serial reference implementation |
52668e843eeSJed Brown| `/cpu/occa`             | Serial OCCA kernels             |
52768e843eeSJed Brown| `/gpu/occa`             | CUDA OCCA kernels               |
52868e843eeSJed Brown| `/omp/occa`             | OpenMP OCCA kernels             |
529bcb2dfaeSJed Brown
530bcb2dfaeSJed BrownExamples available in this release:
531bcb2dfaeSJed Brown
532bcb2dfaeSJed Brown| User code             | Example                           |
53368e843eeSJed Brown|-----------------------|-----------------------------------|
53468e843eeSJed Brown| `ceed`                | ex1 (scalar Laplace operator)     |
53568e843eeSJed Brown| `mfem`                | BP1 (scalar mass operator)        |
53668e843eeSJed Brown| `petsc`               | BP1 (scalar mass operator)        |
537bcb2dfaeSJed Brown```
538