Skip to content
Snippets Groups Projects

Fix: TrustRegionSolver: use blocked matrices for fufem scalar basis assembler

Merged Fix: TrustRegionSolver: use blocked matrices for fufem scalar basis assembler
fix/blocked-local-matrix into master
Merged Patrick Jaap requested to merge fix/blocked-local-matrix into master
Viewing commit f5f0bd78
Show latest version
Compare changes
  • Side-by-side
  • Inline
Files
1
  • f5f0bd78 
    valgrind complained here and it is right: in the internal loop of the
assembler we loop over localMatrix.N() and M(). Since dune-functions
uses flat matrices for the localMatrix the number of dof's is dim-times
too large for a powerBasis. In order to spare a double implementation of
the local assembler only for powerBasis, we can simply convert it in the
lambda. valgrind is now happy here.