minor - explicit casting of vector sizes

minor - update copyright headers

opt - reduce memory usage in CeedOperatorLinearAssembleQFunctionCore_Opt

rstr - use CeedSize for l_size

cpu - drop some empty TensorContract destructors

cpu - fall through now on INTERP/GRAD cases

cpu - fix element/block evec lengths

vec/qf - initial valid/borrowed/owned split for data (#853)

* vec/qf - initial valid/borrowed/owned split for data

* vec/qf - tidy logic for checking active/stale data

* minor - add missing NU

vec/qf - initial valid/borrowed/owned split for data (#853)

* vec/qf - initial valid/borrowed/owned split for data

* vec/qf - tidy logic for checking active/stale data

* minor - add missing NULL

* doc - explain VectorTakeArray update

* minor - update error messages

* test - update error message in junit/tap

* gpu - fix stray CeedScalar vs void for QFunctionContext

* vec/qf - clarify/simplify access logic

* vec - calloc host arrays when no value set to make empty

* style - minor

* style - minor

* minor - fix error messages

* vec/qf - move data validity checking to backend interface

* gpu - add missing sync error checking for qfcontext

* gpu - homogonize use of impl for backend data to reduce confusion

* vec - clarify access conditions

* python - update test for stricter vector access

* vec - minor fixes

* minor - fix ipython change

* vec - add missing declarations in ceed/backend.h

* ctx - mirror vector borrowed data check in ctx interface

* vec - add CeedVectorGetArrayWrite

* vec - consistent use of CeedVectorGetArray vs CeedVectorGetArrayWrite

* python - small vec fixes

* doc - describe vector data semantics

* magma - update restriction

* gpu - fix restr bug I added, need to sum into target

* magma - fix restriction bug

* cpu - fix restriction bug here too

* op - fix evec allocations

* julia - fix ElemRestriction for new vector access rules

* op - double check GetArray vs Read vs Write usage

* doc - small fix

* restr - clean up read/write logic for restr

* python - add vec.array_write

* magma - typo fix

show more ...

cpu - additional minor name updates

cpu - fix unclear variable names in operator impl

cpu - remove e_data array from backend operator data

op - turn magic number 16 into constant CEED_FIELD_MAX

opt - specalize tensor product for C=1 case

qf - cache more QF assembly components for CPU

op - add CeedOperatorLinearAssembleQFunctionBuildOrUpdate to reduce kernel calls in diagonal/full/fdm assembly

interface - refactor *GetFields to include number of fields

op - impl restriction only identity operators for CPU, error on code gen GPU backends

style - updating C variable style

Install install backend headers under include/ceed/

This makes it possible to distribute source plugins that provide
additional backends. It's also used in MFEM, perhaps temporarily.

Deprecate ceed

Install install backend headers under include/ceed/

This makes it possible to distribute source plugins that provide
additional backends. It's also used in MFEM, perhaps temporarily.

Deprecate ceed-backend.h, which was not previously installed, but some
users accessed it from an in-place build.

Also install CUDA and HIP headers that allow users to provide CUfunction
and hipFunction_t.

Co-authored-by: Jeremy L. Thompson <jeremy.thompson@colorado.edu>
Requested-by: Andrew T. Barker <barker29@llnl.gov>

show more ...

Full assembly of CeedOperator (#712)

* Add CeedOperatorLinearFullAssemble() to explicitly assemble a CeedOperator

The output format is a stream of (row, column, value) triples, which client
code

Full assembly of CeedOperator (#712)

* Add CeedOperatorLinearFullAssemble() to explicitly assemble a CeedOperator

The output format is a stream of (row, column, value) triples, which client
code can use to assemble a sparse matrix by adding the values in the
appropriate places.

* Add full assembly test for triangular mesh, t563-operator.c, which is analogous to t536-operator.c

* assembly: change interface to separate nonzero structure from values

* assembly: add potential backends for CeedOperatorLinearAssemble, CeedOperatorLinearAssemblySymbolic

* assembly: use CeedMatrixMultiply() instead of hand-rolled duplicated code

This also involves switching some arrays from column-major ordering to
row-major ordering.

* Remove CeedMatrixMultiply from ceed-impl.h because it is already in ceed-backend.h

* assembly: modify use of fallback resource

* operator - use new error codes in full assembly

* style - minor

* operator - switch default fallback to none

* operator - fix fallback for full assembly

* operator - fix fallback usage

* operator - fix fallback for digonal assembly

* avx - revert bug

* avx - revert bug

Co-authored-by: jeremylt <thompson.jeremy.luke@gmail.com>

show more ...

Error Handling Improvement [fix #696] (#691)

* Operator - add operator/qfunction field compatibility checks

* QFunction - do not allow adding fields to QFunction in use with an operator

* Exam

Error Handling Improvement [fix #696] (#691)

* Operator - add operator/qfunction field compatibility checks

* QFunction - do not allow adding fields to QFunction in use with an operator

* Examples - add some extra exclusion markers in ceed example

* error - add error enum

* error - update error enum names and numbering

* error - use CEED_ERROR_BACKEND in all backend errors

* error - begin classifying interface errors

* error - update backends to use CEED_ERROR_SUCCESS and CeedChkBackend

* error - use new errors in gallery

* error - add some unsaved modifications

* error - improve documentation

* error - define CEED_ERROR_SUCCESS in GPU JiT; we really should have a common header to pipe defines to the JiT code

* error - more error code editing

* error - fix error string

* operator - fix setting field qpts

* basis - add input/output dimension error checking

* python - move basis utility methods to ceed object, no basis required or used

* python - force exit with negative error code

* make style-py

* rust - initial work to add error handling logic

* rust - add ceed.resource method

* rust - add results for methods that may fail

* rust - also format doctests

* minor - drop unused CeedChk()

* error - rename terminal/nonterminal to major/minor

* rust - set ErrorStore as default errorhandler

* python - revert error handing change for python

* python - use success error code from C bindings

* error - only upgrade error code in backend if positive

show more ...

headers - clearify includes to not rely on transitive includes (#701)

* headers - clearify includes to not rely on transitive includes

* style - add header recommendations from 'include-what-you-

headers - clearify includes to not rely on transitive includes (#701)

* headers - clearify includes to not rely on transitive includes

* style - add header recommendations from 'include-what-you-use'

* style - apply 'include-what-you-use' changes to CUDA backends

* style - 'include-what-you-use' for hip backends

* style - drop ceed.h includes in gallery qf source

* docs - add dev notes for header files

* style - header style and alphabetize

show more ...

Add static library (libceed.a) [resolve #670]

We no longer use __attribute__((constructor)) to register backends and
gallery implementations because we can't ensure that the symbols are
linked into

Add static library (libceed.a) [resolve #670]

We no longer use __attribute__((constructor)) to register backends and
gallery implementations because we can't ensure that the symbols are
linked into applications that link the static library. We've switched to
creating CeedRegisterAll() and CeedQFunctionRegisterAll(), which are
called automatically by the library, and call weak symbols to register
all the backend/gallery implementations. This strategy was partly
motivated by not wanting to have preprocessor macros describing what is
available, and the associated need to recompile rather than just relink
when those macros change.

So we now have backends/ceed-backend-list.h that declares all the
backends wrapped in a macro. It is included by
backends/ceed-backend-weak.c to create weak definitions of all the
backends. In the makefile, we sort so this comes last when linking a
shared or static library, and thus these weak symbols will only be
picked up if they were not defined by the actual backend source files.
The same header is included (with different macro wrapping) in
interface/ceed-register.c, where CeedRegisterAll() is defined.

To add a new backend, one must do essentially the same registration
strategy as in the past, plus add one line to the common
ceed-backend-list.h.

show more ...

Jeremy/py debugging (#637)

* Ceed - only store error with parent

* Python - assert state is not in write status

* Vector - add reader check to destructor

* make style

* remove debugging

Jeremy/py debugging (#637)

* Ceed - only store error with parent

* Python - assert state is not in write status

* Vector - add reader check to destructor

* make style

* remove debugging assertion

* Python - minor documentation correctness

* Python - hold reference to data array for USE_POINTER to prevent gc

* Python - hold reference to data array for USE_POINTER to prevent gc for restriction and context

* Python - use None object for initial ref back to ceed objects

* Python - add USE_POINTER to imports as needed

* Junit - update valid error message for writable access lock

* Python - change exception raised test to help gc

show more ...

QFunction Context Data Object (#596)

* Ctx - create context object for QFunction context data

* Context - rename UserContext -> QFunctionContext

* Ctx - add lcov markers

* Ctx - fix leak in

QFunction Context Data Object (#596)

* Ctx - create context object for QFunction context data

* Context - rename UserContext -> QFunctionContext

* Ctx - add lcov markers

* Ctx - fix leak in identity QFunctions

* Hip/Cuda - rename sync functions for vector/context

* Tests - lcov marker update

* QFunction - drop unused function

* Python - fix copy-paste errors

* Ctx - update notes for Fortran usage

* Fortran - drop unneeded cast

Co-authored-by: Jed Brown <jed@jedbrown.org>

* Interface - use void* for SetData interfaces

* Make - use call quiet for NVCC

* Interface - use void* for GetData interfaces

* Make - add quiet call option for examples

* Makefile - create common makefile to reduce duplication/complexity in example makefiles

Co-authored-by: Jed Brown <jed@jedbrown.org>

show more ...