Open keflavich opened 4 years ago
I have been having this same issue running Ubuntu 16.04. Is there a workaround available?
No, we couldn't figure out what was causing it! Any insights you can offer would be very welcome. Maybe try different compilers?
I did manage to avoid the singular matrix error by editing line 952 of core.py:
frac_level_diff = level_diff/self.level_population
Here there was a divide by zero problem as self.level_population
had a number of elements equal to zero (which may be part of the original problem). I removed the elements equal to zero and got some output but it did not match the example output provided on this site for the following input:
R = pyradex.Radex(collider_densities={'oH2':900,'pH2':100}, column=1e16, species='co', temperature=20)
Tlvg = R(escapeProbGeom='lvg')
So in short I am still having major problems with it.
cc @low-sky
We have an issue where RADEX is producing nonsense results and crashing on ubuntu, but not on mac, for identical inputs.
results in