Watch this Author's files
How can I adjust the solver such that I can solve the regular Poisson equation?
Works flawlessly! Good job!
[SOLUTION]: To be able to read the par file correctly, add the following codes in par_read_header mfile:
RescaleIntercept = vald(12);
RescaleSlope = vald(13);
ScaleSlope = vald(14);
right after the line:
Then also add the following codes:
info.RescaleIntercept = RescaleIntercept;
info.RescaleSlope = RescaleSlope;
info.ScaleSlope = ScaleSlope;
after the line:
Another bug is that, par_read_header function does not really read 'RescaleSlope' and 'RescaleIntercept' and simply assumes 1 and 0, respectively. So it does not really matter whether we set real to true or not.