Developments in preconditioned iterative methods with application to glacial isostatic adjustment models