LU factorization for 3d matrix

Hello,

I would need to solve P linear systems if possible without using loops. I was thinking about LU factorization which in my case is quite fast. I know how to approach the problem for a 2d matrix. But what should i do for the case of a 3d matrix M(N,N,:)?

thanks in advance!

Answer 1

[John D'Errico]

Convert it to a sparse block diagonal system. Then use backslash. No problem.

[rengell]

I will try that. Thanks!