Fixing #38: replaced interpolate( 0 ) with setToZero() again in inv diag

Closes #38

See merge request !21

Fixing #38: replaced interpolate( 0 ) with setToZero() again in inv diag kernel type. Also extended test such that the diagonal is assemble more than once (without this fix that failed).

Introduce NewTypes for Trial- and TestSpace

Closes #14

See merge request !20

introduce NewTypes for Trial- and TestSpace such that mypy can check they are not confused; fix trial/test mix-ups

Weak Nitsche freeslip integrands, improved and extended external parameter handling

See merge request !19

Not throwing exception because it is not a valid assumption that the free symbol must still exist after possible symbolic simplification.

* FE coeffs now passed with dicts instead of lists of tuples
* constant scalar and matrix coefficients can now easily be added with the `scalars` and `matrix` callables passed by the process_integrand() function
* added the coordinates as a symbol to the integrand signature

Refactoring forms

See merge request !16

Daniel Bauer authored 9 months ago and Nils Kohl committed 9 months ago

Co-authored-by: Daniel Bauer <daniel.j.bauer@fau.de>

Improving/fixing tabulation: by default now only a single table is constructed per integrand unless an optional name is supplied.

Refactored shear heating and SUPG forms. Extended integrand symbols by hessian stuff. Some refactoring of the Hessian FEM helper functions.

Fix static type check in frozen velocity form

See merge request !18

Custom deepcopy implementation for DoFSymbol since it seems to be required for dataclasses.asdict in integrand.py.

This commit moves the boilerplate code that has been part of every form to a dedicated place: the IntegrandSymbols class and the process_integrand() function (integrand.py). This change simplifies the construction of new operators/forms.

The refactoring is not complete (essentially, the files forms.py, forms_boundary.py, manifold_forms.py, etc. will become unnecessary). But for regression testing, the skeletons are kept for another while.

Concrete changes and additions:
* implemented the function process_integrand() that processes a symbolic integrand - see docstrings in integrand.py for details
* moved symbolic integrand expressions to hog/recipes/integrands - less code and better separation of concerns
* the TensorialVectorFunctionSpace can now be restricted to a single component - now the same integrands can be used for full vector forms and component-wise terms
* fixing test and trial mix-up in generate_all_operators.py (see #14)

Add general frozen velocity operator/form

See merge request !17

Boundary Integrals

See merge request !15

* cleaned up imports
* fixed comments
* printing space dimension in all ElementGeometry subclasses
* renamed vertex shuffle function
* moved substitution of last reference coordinate to add_boundary_integral()