both Dirichlet boundary conditions and the constraints arising from
adaptively refined meshes can be imposed within the matrix-free
framework.
- \textcolor{red}{%
- \sout{The only restriction is that for two-dimensional meshes the
- finite element degree must be odd. There is no such restriction in
- three dimensions.}
- \emph{%
- WB: Can we explain this restriction? It seems odd (pun intended)...
- BT: No, we cannot. I have spent at least two full weeks trying to
- understand why it doesn't work, i.e., gives the wrong result, but I
- don't know. We decided at the time that even with the current
- limitation this was a useful addition.
- MM: I suggest that we then simply not mention the limitation and treat
- it as a known bug. As it is currently worded it can be easily
- interpreted as a fundamental limitation, i.\,e. ``FE with even degree
- are not admissible.''}}
+% \textcolor{red}{%
+% \sout{The only restriction is that for two-dimensional meshes the
+% finite element degree must be odd. There is no such restriction in
+% three dimensions.}
+% \emph{%
+% WB: Can we explain this restriction? It seems odd (pun intended)...
+% BT: No, we cannot. I have spent at least two full weeks trying to
+% understand why it doesn't work, i.e., gives the wrong result, but I
+% don't know. We decided at the time that even with the current
+% limitation this was a useful addition.
+% MM: I suggest that we then simply not mention the limitation and treat
+% it as a known bug. As it is currently worded it can be easily
+% interpreted as a fundamental limitation, i.\,e. ``FE with even degree
+% are not admissible.''}}
\item MPI matrix-free computations: using \texttt{LinearAlgebra::distributed::Vector}, the
matrix-free framework can scale to multiple GPUs by taking