The noise-scaling maps were not being extracted correctly, because the self.chi_squared_map was being altered when we extract them from a ci_fit:
@property
def noise_scaling_map_of_ci_regions(self):
return self.ci_data_fit.chinj.ci_regions_from_array(array=self.chi_squared_map)
I fixed this by copying the array.
@property
def noise_scaling_map_of_ci_regions(self):
return self.ci_data_fit.chinj.ci_regions_from_array(array=self.chi_squared_map.copy())
We need to make sure this sort of thing isn't happening anywhere else we extract arrays from ci_frame.