Non-square operator full assembly (#1316)

* Consistency formatting for operator tests

* Add (failing) test for non-square operator full assembly

* Add support for CPU-based full assembly of no

Non-square operator full assembly (#1316)

* Consistency formatting for operator tests

* Add (failing) test for non-square operator full assembly

* Add support for CPU-based full assembly of non-square CeedOperators

* Fix full assembly of identity quadrature functions and operators with CEED_BASIS_NONE

Also adds a unit test for full assembly which fails prior to the changes in this commit.

* Minor refactor to improve coverage

show more ...

Merge pull request #1270 from CEED/jeremy/constantinople

Rename CEED_BASIS_COLLOCATED to CEED_BASIS_NONE

interface - rename CEED_BASIS_COLLOCATED to CEED_BASIS_NONE

Merge pull request #1197 from sebastiangrimberg/sjg/style-whitespace-fix

Minor style consistency updates

Formatting fixes for tests including consistent include order

IWYU fixes (#1182)

* iwyu - include fixes

* iwyu - silence some iwyu output

* minor - clearer macro names

* iwyu - fix suggestion of "ceed/ceed.h" externally

* iwyu - lighter petsc heade

IWYU fixes (#1182)

* iwyu - include fixes

* iwyu - silence some iwyu output

* minor - clearer macro names

* iwyu - fix suggestion of "ceed/ceed.h" externally

* iwyu - lighter petsc headers

* iwyu - ceed/ceed.h -> ceed.h

* iwyu - cuda/hip include fixes

show more ...

Merge pull request #1167 from CEED/jeremy/drop-i

Consistency in c tests

style - consistency in c tests

Switch to clang-format (#1051)

* style - switch to clang-format

* ci - use newer libxsmm

* action - update format action

* format - consistent use of {} for multi-line if/for

* make - re

Switch to clang-format (#1051)

* style - switch to clang-format

* ci - use newer libxsmm

* action - update format action

* format - consistent use of {} for multi-line if/for

* make - remove stray newline

* make - simpler 'make format' target

* ci - use newer libxsmm

* doc - minor release note claification

* minor - minor fix

* minor - minor fix

* minor - minor fix

* minor - minor fix

* make format

* format - less aggressive alignment rules

* tidy - check for argument name mismatches

* fix newline

* format - mirror Ratel update to .clang-format

* fix merge error

* fix merge conflict

* fix merge error

* drop style in .phony list

* Update .clang-format

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

* apply updated format

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

show more ...

Merge pull request #1002 from CEED/jeremy/format-ceed-types

Add CeedInt_FMT

test - update tests to use CeedInt_FMT

Merge pull request #913 from CEED/jeremy/coo-ptrdiff

Create CeedSize as ptrdiff_t

ceed - create type CeedSize

pc - use PetscCount/ptrdiff_t for num COO entries

Merge pull request #877 from CEED/jeremy/field-names-for-people

Make qf/op field names human readable

style - make qf/op field names human readable

Allow CeedScalar to be single precision (#788)

One can modify `ceed.h` to include `ceed-f32.h` and then use single precision. This is tested for C in CI and has been tested by developers with Rust,

Allow CeedScalar to be single precision (#788)

One can modify `ceed.h` to include `ceed-f32.h` and then use single precision. This is tested for C in CI and has been tested by developers with Rust, Julia, and Python. This interface is evolving and should be considered experimental at this time (thus lack of automated build support).

* Introduce CeedScalarType enum

* WIP changes to allow different definitions of CeedScalar

* Introduce new header files for float and double

* Only use avx tensor contract and MAGMA non-tensor basis if CeedScalar is double

* WIP changes to allow CeedScalar to be float

* WIP start trying to adjust test tolerances for float or double

* fix typos in comments

* install ceed-f32/64 headers

* Fix missing casts for hipMAGMA element restrictions

* make CeedQFunctionContextGetContextSize available for Python bindings

* Changes to Python bindings to allow CeedScalar to be float

* WIP adjust Python tests for float or double

* make style

* remove QFunctionContextGetContextSize from backend header

* Use quotes instead of <> in include statement

* Remove unncessary includes

* Update tolerances for tests

* [Julia] allow CeedScalar to be Float32

* [Julia] Use Preferences instead of custom build configuration

# Conflicts:
# julia/LibCEED.jl/src/C.jl

* [Makefile] Change definition of CC_VENDOR so it works with cross-compilation

* [Julia] Use Preferences in CI

# Conflicts:
# .github/workflows/julia-test-with-style.yml

* [Julia] Update docs about preferences

* [Julia] Add test/Project.toml workaround for Preferences

* Add CeedGetScalarType to get the type of CeedScalar at runtime

* [Julia] Move functions from Ceed.jl to LibCEED.jl

* [Julia] Add support for getting library path and scalar type at runtime

* [Julia] Minor change to checking if CUDA is loaded

* [Julia] Check correct CeedScalar types in basis functions

* [Julia] Fix tests comparing with output file

* [Julia] Change devtests to use CeedScalar instead of Float64

* Update test 402 so context will be same size in double or float

* Update tolerances for ceed examples

* [Julia] CUDA fixes

* remove unused variable in t208

* SchurDecomposition: do not compute tau on final iteration

* Update tolerances for some basis tests (for single precision)

* Make style

* Python style fixes for basis test

* Add single precision output for t300 and t320 and adjust checks; skip t541 in single

* Add LCOV exclusions after moving to new line

* fix spacing

* Python: make CEED_EPSILON available as libceed.EPSILON

* Python: optional parameter to specify different output file for test comparison

* Python: update tests' use of EPSILON and change test_300 output file for single precision

* Python: add convenience function for getting dtype corresponding to CeedScalar

* rust - add single precision support

* [Julia] Fall back on Float64 if CeedGetScalarType is not available

* [Julia] style

* Adjust tolerance for t301

* xsmm - add single precision support

* avx - add single precision support

* Add initial single precision support for MAGMA non-tensor basis

* Skip t300 and t320 in single precision; revert Python t300 changes

* Revert output changes for t300 and t320 in junit

* [Julia] Changes to autogenerated bindings for mixed precision

* [Julia] style

* [Julia] Check scalar type when changing libceed library path

The check is also performed when the package is loaded. This prevents having to
restart the Julia session twice

* [Julia] Require JLLWrappers version 1.3

This is needed to use Preferences to change the library path

* Add documentation page for precision development

Co-authored-by: Will Pazner <will.e.p@gmail.com>

* Cleanup from merge: remove old README

* Return CEED_ALIGN to backend.h

* Make Fortran compiler (FC) optional; empty skips Fortran tests

Use in Python and Rust builds, which may not have a Fortran compiler
installed and thus would produce confusing output.

* Add single precision CI test for Noether

Co-authored-by: Jeremy L Thompson <jeremy@jeremylt.org>

Co-authored-by: Will Pazner <will.e.p@gmail.com>
Co-authored-by: Jeremy L Thompson <jeremy@jeremylt.org>
Co-authored-by: Jed Brown <jed@jedbrown.org>

show more ...

Merge pull request #735 from CEED/jeremy/style-guide

Style Unification

style - updating C variable style

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 ...