Thank you for the reply and useful link, I shall give that one a try!
I have attached a much smaller system for you to take a look at. The outline of the file is as written in my initial post.
A bit of context for the problem, I am running a finite difference method solver for an 8 band k.p model in 3D. The exact problem we are trying to solve requires a very small step size - hence the large matrix.