Two changes:

1. Avoid making a CPU banded matrix and avoid doing copies with
size N * (2 * K + 1)
2. Fix a bug in block banded LU.