Skip to content

Draft: Added support for simd assembly

Alexander Schell requested to merge alexander.schell/dune-pdelab:support_simd_assembly into master

This MR adds support for SIMD assembly of the DefaultLocalAssembler. It deals with the case where the DomainField and the RangeField of the gridoperator are SIMD and one wants multiple residuals to be assembled at once in a vectorized fashion.
The problem was that the local assembler's weight_ was not scalar anymore although it should be. The reason was that weight_ was defined to be of type RangeField which in that case was SIMD. To ensure that it remains scalar we changed its type to Real. And we (@christi) changed the type definition of Real to Simd::Scalar<RangeField> such that this type actually represents a real number and is not SIMD.

Merge request reports