Open bruno-moraes opened 1 week ago
I will work on this now.
@fjaviersanchez , I could pinpoint more or less where the error is, but this is something I never fixed before. I assume it's numerical instabilities in the derivatives due to the small values of A_s compared to the other parameters. I won't be able to solve it now.
Thank you for checking and identifying the issue! I think this can be potentially solved when we add the stem derivative method. One potential solution is to normalize all parameters internally for the derivatives and then change back or change the step size dynamically, but will check first if the easy fix works.
TL;DR: A call to boltzmann.py within CCL is failing in the derivative calculation, presumably due to O(10^-9) values in A_s.
Explanation: When attempting a 3x2pt forecast, it runs and completes when the sigma8 parameter is used. If A_s is used instead, the Jupyter kernel dies after a few seconds when trying to calculate the Fisher matrix. Running the code as a python script from the command line, the error returned is:
The trace did not contain information about where this error occurred, so I ran
pdb
on the script (my first time, so not very skillfully) and found the spot where it is breaking: