Skip to content

Feature/loopy invert matrix

Marcel Koch requested to merge feature/loopy-invert-matrix into master

Moves the matrix inversion stuff from blockstructured/geometry.py to pdelab/tensor.py. This touches on the subject of #139 (closed).

You can get the name of the inverse of a matrix A with shape shape from

name = name_matrix_inverse(A, shape, visitor)

matrix_inverse() returns the inverse wrapped in a pymbolic variable, although I don't know how useful this is.

TODO:

  • handle inames correctly
  • what about tags? -> not really necessary
  • what about dependencies -> not really necessary
Edited by Marcel Koch

Merge request reports