Few comments on current matrix-free implementation and missing features
authorkkormann <ekkormann@gmail.com>
Tue, 7 Oct 2014 19:26:31 +0000 (12:26 -0700)
committerkkormann <ekkormann@gmail.com>
Tue, 7 Oct 2014 19:26:31 +0000 (12:26 -0700)

index 9e47f709fd18dbeb806269c16835a25696b7c304..d89c1b43942414e79807ddfe611caef087f9eab1 100644 (file)
@@ -31,6 +31,14 @@ most algorithms become memory-bound, we must consider retiring
   product structure
   - XFEM
   - Overlapping nonmatching grids
+- Note: *hp* adaptivity is supported by the `FEEvaluation` class.
+## Design and limitations of the `FEEvaluation` class
+The matrix-free implementation uses the class `FEEvaluation` instead of `FEValues`. Different from the `FEValues` class, information about the degrees of freedom and the mapping from physical to unit cell is precomputed and stored. The degrees of freedom are numbered MPI-local and the constraints are incorporated into the DoF list in order to provide direct array access into vectors.
+This enables e.g. fast evaluation of the shape functions at the quadrature points. On the other hand, only values, gradients, etc. that are precomputed can be accessed while the `FEValues` object enables on-the-fly computations of more general evaluations.
+The main task would therefore be to merge the different functionalities of the two classes.
 # Tensor product quadrature
@@ -123,3 +131,4 @@ well. For these we need truncated tensor products.
   numbering of the shape functions, which is the inversion of a
   similar process when we generate the functions. Should we reconsider
   the ordering in `DoFHandler`?
+- `FEEvaluation` currently knows the elements `FE_Q`, `FE_DGQ`, `FE_DGP` and `FE_Q_DG0`.
\ No newline at end of file

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams

Typeset in Trocchi and Trocchi Bold Sans Serif.