Add missing include / use std::invoke_result

See merge request !372

Clang master (futur 15) warns about it

Caused compilation error in DuMuX

[CI] Run the 'codespell' tool

See merge request !371

This will warn about misspellings in the source code.

derivative of `DiscreteGlobalBasisFunction`

See merge request !366

Carsten Gräser authored Jun 30, 2022 and Ansgar Burchardt committed Jul 01, 2022

The new namespace is dedicated to implemenation details that should
not be generically ignore by doxygen, e.g. because the contained
classes are mixins whole members should show up in the documentation
of the derived class.

Carsten Gräser authored Jun 30, 2022 and Ansgar Burchardt committed Jul 01, 2022

This mixin base class does not do anything with this member.
Hence it's easier to not pass it as template parameter, but
directly store it in the derived class.

Carsten Gräser authored Jun 30, 2022 and Ansgar Burchardt committed Jun 30, 2022

* Unify algorithm for mapping shape function range and coefficient
  to the function range. Effectively the new implementation produuces
  the same results but is much easier to understand and describe.
* Move `DiscreteGlobalBasisFunctionBase` to implementation name space.
* Factor out access to nodeToRangeMap` and assignment of different
  container types.
* Uniformly initialize all vectors/coefficients with `istlVectorBackend(c)=0;`
* Don't use `operator+=` when we actually assign a zero-initialized value.
* Swap loop order. Having the statically sized loop inside it is more
  likely to be unrolled.

Ansgar Burchardt authored Jun 15, 2022 and Ansgar Burchardt committed Jun 30, 2022

The `DiscreteGlobalBasisFunctionDerivative` evaluates the derivative
of the `DiscreteGlobalBasisFunction` with respect to global
coordinates.

This implementation is based on !188.

The `LocalFunction` now holds a shared pointer for the struct
containing the required data.

Ansgar Burchardt authored Jun 15, 2022 and Ansgar Burchardt committed Jun 27, 2022

This allows the LocalFunction to exist even when the associated
DiscreteGlobalBasisFunction has been removed.

It furthermore is expected to ease addition of a derivative: the
derivative's global and local functions have to share the same data
and it is easier to use if there are no shared lifetime requirements
between all of these.

[bugfix] Define hash_value for StaticMultiIndex

See merge request !370

We already defined the macro, but the actual hash function was missing.

[cleanup] Simplify index computation in LagrangeDGPreBasis

See merge request !369

With this patch the element offset is computed once per element
and not once per basis function. Using `if constexpr` instead
of `switch` also makes this slightly more readable (no need to
escape `switch` with `continue`).

Add missing include of hash.hh

See merge request !368

Improve multi index type in global basis

See merge request !329

* We now separate the stored and exported multi index type
  `MultiIndexStorage` and `MultiIndex` in the `LocalView`.
* `MultiIndexStorage` will always have dynamic size and may need
  a larger capacity compared to `MultiIndex`s size.
* If it is clear that the multi indices generated by the `PreBasis`
  have a static length, the`MultiIndex` has static size, while
  the `MultiIndexStorage` is an `OverflowArray<MultiIndex, ...>`.
  Otherwise the `MultiIndex` is a `ReservedVector`.
* If the size of multi indices is guaranteed to be one, then the
  `MultiIndex` can be cast to the plain index to avoid using `[0]`
  in simple cases.

This represents a statically sized multiindex.
It adds a static `size()` methods and a cast to the first
entry for `size()==1`.
`FlatMultiIndex<T>` is now an alias for `StaticMultiIndex<T,1>`.

It mostly behaves like `Dune::ReservedVector` but derives from a statically
sized array base class. This allows to have temporary dynamic size but cast
to the result to the statically sized base class.

Specialize FlatVectorBackend for std::array

See merge request !367

This adds a specialization for `std::array` and some improvements
to hand out size statically. Notice that this only works if
the entries of the `array` have a statically known size.
Ideally one would implement a recursive `FlatVectorBackend`
but I can't get this to compile (with a 'use of ... before decution
of auto' error).

[example][cleanup] Use geometry.jacobian(...) directly

See merge request !365

Simplify examples by replacing the chain rule hack
```cpp
gradient = referenceGradient * transpose(jacobianTransposed);
```
by plain use of
```cpp
gradient = referenceGradient * jacobian;
```

add `bound()` member function to LocalView and LocalFunction

Closes #72

See merge request !360

Carsten Gräser authored May 24, 2022 and Ansgar Burchardt committed Jun 15, 2022

Surprisingly the test worked before, although it does a concept check
for the interface. The reason is, that the interface of the type erasure
class contains the `bound()` method, so the substitution in the concept
check works. But this just checks if the expression is valid when
instantiating the template. It does not check if the instantiated
template compiles. Hence we really need to call the method in the check.

Ansgar Burchardt authored May 20, 2022 and Ansgar Burchardt committed Jun 15, 2022

It might be confusing to have the function still bound after calling
`unbind()`. So reset the geometry to have `bound()` return false
afterwards.

Ansgar Burchardt authored May 20, 2022 and Ansgar Burchardt committed Jun 15, 2022

Also implement the method for other local functions besides
DiscreteGlobalBasisFunction.

Ansgar Burchardt authored May 20, 2022 and Ansgar Burchardt committed Jun 15, 2022

Also implement it in SubspaceLocalView in addition to DefaultLocalView.