Running model in a HPC and would like to save intermediate outputs

As Martin says, you could probably benefit from sparse representations; I have a demo in the context of splines here