Antony Lewis
Antony Lewis
Thanks, PRs welcome. Would need a bit of care with multiple massive neutrino models etc.
Looks far too large to be physical, effect of (inaccurate linear) reion temp change is very small. (though this is crazy small scale, not sure I've looked much beyond 500)...
Possibly forcing WantCls is causing minuscule time steps that are then numerically unstable?
I think the difference comes from WantCl=F turning off reionization, so it really is just an effect of reionization on the growth via the sound speed. Obviously CAMB can't calculate...
As shown in PR (https://github.com/cmbant/CAMB/pull/180), looks like the sound speed was going negative, which is why the sign of the effect is different to what we expected.
May be a compatibility option that would allow this. However, I doubt it will work with that compiler in general -the Python interface has only been tested with ifort and...
Do the forutils unit tests pass? If not, probably compiler bug/difference. Or it's interpretting the dummy object pointer argument to dverk in subroutines.f90 differently.
Which version gcc and what error? Have not tested with valgrind recently; but camb unit tests do have a test for memory leaks.
Not super suprising, what the code is doing here is very hacky.
This is true, though the optimizer will normally apply short-circuit evaluation. Easy enough to change e.g. ``` if (allocated(ajl)) then if (any(ubound(ajl) < [num_xx, max_ix])) deallocate(ajl, ajlpr, ddajlpr) end if...