LU for large band matrices

[henrikandresen]

Hi All

I need a solver for a bandmatrix with a band of approx. 1000. There is only 1 solution required, so LU or direct solution are both almost identical in performance.

The tricky part is that the matrix is too large to fit in the memory of a single GPU, so I need some sort of data streaming between the system memory and GPU memory.

Is this something that can be done using CULA?

Best regards

Henrik Andresen

[john]

This might not be an easy fit for CULA. What is the long dimension of the matrix?

[henrikandresen]

The long dimension will be something in the 10^5 to 10^6.

Cheers

Henrik Andresen

[john]

That's a tricky one since it's just over the size of GPU memory. We tend to call this "out of core" operation, which is unfortunately not something that's slated for the banded solvers any time soon.