Profile the scaling behaviour for larger models.

[FredDeCeuster]

Profile the scaling behaviour for larger models (e.g. 256 x 256 x 256).
Memory consumption seems rather large (50 GB), and computations seem unreasonably slow (~10 minutes).
Profile seems to suggest that a nonzero function is the culprit.