Skip to content

WIP: Implement `forEach` matrix loops for rows and entries

Patrick Jaap requested to merge patrick.jaap/dune-istl:feature/matrix-iterator into master

The routine multiIndexedForEach can be used to loop over any blocked container on a scalar level. For a meaningful application the matrix blocks have to be aligned.

TODOs:

  • Improve naming of the routines?
  • Documentation!
  • There is one further restriction on the blocks: If there is an empty sparse row (e.g. in BCRSMatrix), we assume all descending blocks to be of static size. Otherwise the number of scalar rows cannot be extracted efficiently.
Edited by Patrick Jaap

Merge request reports