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?