Menu

Commit [r23414]  Maximize  Restore  History

Critical fix for the math domain catching of model CR72.

This was discovered with the added 8 unit tests demonstrating edge case 'no Rex' failures.

This follows from the ideas in the post http://article.gmane.org/gmane.science.nmr.relax.devel/5858.
This is related to: task #7793: (https://gna.org/task/?7793) Speed-up of dispersion models.

This is to implement catching of math domain errors, before they occur.

When kex is large, ex: kex = 1e5, the values of:
etapos = eta_scale * sqrt(Psi + sqrt_psi2_zeta2) / cpmg_frqs

will exceed possible numerical representation.

The cathing of these occurences needed to be re-written.

tlinnet 2014-05-26

changed /branches/disp_speed/lib/dispersion/cr72.py
/branches/disp_speed/lib/dispersion/cr72.py Diff Switch to side-by-side view
Loading...
Want the latest updates on software, tech news, and AI?
Get latest updates about software, tech news, and AI from SourceForge directly in your inbox once a month.