scalarRectangularMatrix inside Foam::pinv(const tensor& t) uses memory allocation
Functionality to add/problem to solve
scalarRectangularMatrix is of variable size. Its use inside e.g. tensor makes a primitive suddenly allocate memory. This can give problems offloading.
Proposal
Have (equivalent of) RectangularMatrix templated on size?