Menu

Commit [r21113]  Maximize  Restore  History

Replaced all usage of scipy.linalg.expm() with lib.linear_algebra.matrix_exponential.matrix_exponential().

This is for the functions of the lib.dispersion package used for the relaxation dispersion numeric
solution models. The change eliminates a bug in the scipy function which uses the Pade
approximation which fails horribly for the complex part of the matrix. The real part looks good,
but the complex part looks to have nasty truncation artefacts which is propagated and amplified
through the Bloch-McConnell equations.

bugman 2013-10-15

changed /branches/relax_disp/lib/dispersion/mq_ns_cpmg_2site.py
changed /branches/relax_disp/lib/dispersion/ns_cpmg_2site_3d.py
changed /branches/relax_disp/lib/dispersion/ns_cpmg_2site_expanded.py
changed /branches/relax_disp/lib/dispersion/ns_cpmg_2site_star.py
changed /branches/relax_disp/lib/dispersion/ns_r1rho_2site.py
/branches/relax_disp/lib/dispersion/mq_ns_cpmg_2site.py Diff Switch to side-by-side view
Loading...
/branches/relax_disp/lib/dispersion/ns_cpmg_2site_3d.py Diff Switch to side-by-side view
Loading...
/branches/relax_disp/lib/dispersion/ns_cpmg_2site_expanded.py Diff Switch to side-by-side view
Loading...
/branches/relax_disp/lib/dispersion/ns_cpmg_2site_star.py Diff Switch to side-by-side view
Loading...
/branches/relax_disp/lib/dispersion/ns_r1rho_2site.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.