relax SVN archive

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.

Merged revisions 21111 via svnmerge from
svn+ssh://bugman@.../svn/relax/trunk

........
r21111 | bugman | 2013-10-15 13:20:25 +0200 (Tue, 15 Oct 2013) | 5 lines

Fix for the new lib.linear_algebra.matrix_exponential.matrix_exponential() function.

This function now returns a numpy array type rather than matrix type.
........

Fix for the new lib.linear_algebra.matrix_exponential.matrix_exponential() function.

This function now returns a numpy array type rather than matrix type.

Merged revisions 21107-21109 via svnmerge from
svn+ssh://bugman@.../svn/relax/trunk

........
r21107 | bugman | 2013-10-15 11:31:40 +0200 (Tue, 15 Oct 2013) | 5 lines

Created the lib.check_types.is_complex() function.

This is used to determine if a number is a Python or numpy complex type.
........
r21108 | bugman | 2013-10-15 11:34:11 +0200 (Tue, 15 Oct 2013) | 5 lines

The lib.linear_algebra.matrix_exponential.matrix_exponential() function now uses lib.check_types.is_complex().

This fixes the function for complex matrices.
........
r21109 | bugman | 2013-10-15 11:45:29 +0200 (Tue, 15 Oct 2013) | 3 lines

Created a new unit test for lib.linear_algebra.matrix_exponential.matrix_exponential() for complex matrices.
........

Created a new unit test for lib.linear_algebra.matrix_exponential.matrix_exponential() for complex matrices.

The lib.linear_algebra.matrix_exponential.matrix_exponential() function now uses lib.check_types.is_complex().

This fixes the function for complex matrices.

Created the lib.check_types.is_complex() function.

This is used to determine if a number is a Python or numpy complex type.

Merged revisions 21103-21105 via svnmerge from
svn+ssh://bugman@.../svn/relax/trunk

........
r21103 | bugman | 2013-10-15 10:43:03 +0200 (Tue, 15 Oct 2013) | 6 lines

Created a unit test for the lib.linear_algebra.matrix_exponential module.

This module does not exist yet, but it will be used to replace the scipy.linalg.expm() function use
in the relaxation dispersion branch.
........
r21104 | bugman | 2013-10-15 11:07:19 +0200 (Tue, 15 Oct 2013) | 3 lines

Loosened the lib.linear_algebra.matrix_exponential.matrix_exponential() unit test checks.
........
r21105 | bugman | 2013-10-15 11:07:54 +0200 (Tue, 15 Oct 2013) | 5 lines

Implemented the lib.linear_algebra.matrix_exponential.matrix_exponential() function.

This handles square matrices in either complex or real form.
........

Implemented the lib.linear_algebra.matrix_exponential.matrix_exponential() function.

This handles square matrices in either complex or real form.

Loosened the lib.linear_algebra.matrix_exponential.matrix_exponential() unit test checks.