Just in case anyone is looking at this thread in future, I moved to using the Python implementation of SLEPc (Slepc4py) to solve my large sparse matrices, as I couldn’t get the Julia wrapper to work 
. On the upside, performance is quite good! I can now, quite quickly, solve very large systems now without running into major memory requirements.
Thanks all for your help!