Menu

SVN archive Commit Log


Commit Date  
[r25330] by tlinnet

Tried to implement the Jacobian from C-code.

This though also report errors which are to small.

Maybe some scaling is wrong.

task #7822(https://gna.org/task/index.php?7822): Implement user function to estimate R2eff and associated errors for exponential curve fitting.

2014-08-27 09:29:24 Tree
[r25329] by tlinnet

Implemented the first try to compute the Variance of R2eff and i0, by the co-variance.

This uses the Jacobian matrix.
The errors calculated, are though way to small compared 2000 Monte-Carlo simulations.

task #7822(https://gna.org/task/index.php?7822): Implement user function to estimate R2eff and associated errors for exponential curve fitting.

2014-08-27 09:12:50 Tree
[r25328] by bugman

Merged revisions 25317-25319,25321-25327 via svnmerge from
svn+ssh://bugman@.../svn/relax/trunk

........
r25317 | bugman | 2014-08-26 19:38:01 +0200 (Tue, 26 Aug 2014) | 8 lines

Added a script and log file for calculating the numerical gradient for an exponential curve.

This uses the data at http://thread.gmane.org/gmane.science.nmr.relax.devel/6807/focus=6840 and
calculates the Hessian using the numdifftools.Hessian object construct and obtain the matrix, both
at the minimum and at a point away from the minimum. The values will be used to construct a unit
test to check the C module implementation.
........
r25318 | bugman | 2014-08-26 19:49:36 +0200 (Tue, 26 Aug 2014) | 3 lines

Fixes for the Hessian.py script for numerical integrating the Hessian for an exponential curve.
........
r25319 | bugman | 2014-08-26 19:51:12 +0200 (Tue, 26 Aug 2014) | 7 lines

Implemented two unit tests to check the Hessian of the target_functions.relax_fit.d2func() function.

This compares the calculated Hessian to the numerically integrated values from the
test_suite/shared_data/curve_fitting/numeric_gradient/Hessian.py script, showing that the d2func()
function is implemented correctly.
........
r25321 | tlinnet | 2014-08-27 08:50:56 +0200 (Wed, 27 Aug 2014) | 3 lines

Modified profiling script, but it seems that the dfunc from target_functions.relax_fit does not work.

task #7822(https://gna.org/task/index.php?7822): Implement user function to estimate R2eff and associated errors for exponential curve fitting.
........
r25322 | tlinnet | 2014-08-27 09:19:45 +0200 (Wed, 27 Aug 2014) | 4 lines

Modified estimate R2eff module, to use C-code.

But systemtest Relax_disp.test_estimate_r2eff_error shows that the Jacobian is not correctly implemented
to be called in minfx.
........
r25323 | bugman | 2014-08-27 09:47:13 +0200 (Wed, 27 Aug 2014) | 7 lines

Created an initial test suite data directory for a mixed R1rho + CPMG dispersion analysis.

The generate.py script will be extended in the future to generate both synthetic R1rho and CPMG data
for a common exchange process. Such a data combination should show some minor flaws in the current
design of the dispersion analysis and will help to solve these.
........
r25324 | bugman | 2014-08-27 10:01:54 +0200 (Wed, 27 Aug 2014) | 7 lines

Improvements to the pipe_control.minimise.reset_min_stats() function.

The minimise statistics resetting is now more elegantly implemented. And the sim_index keyword
argument is accepted by the function and individual Monte Carlo simulation elements can now be
reset.
........
r25325 | tlinnet | 2014-08-27 10:05:03 +0200 (Wed, 27 Aug 2014) | 3 lines

Modified wrapper function for curve_fit, to only change to list type, if the type is a ndarray.

task #7822(https://gna.org/task/index.php?7822): Implement user function to estimate R2eff and associated errors for exponential curve fitting.
........
r25326 | bugman | 2014-08-27 10:10:57 +0200 (Wed, 27 Aug 2014) | 6 lines

The model-free reset_min_stats() function has been replaced with the pipe_control.minimise version.

The specific_analyses.model_free.optimisation.reset_min_stats() function has been deleted and
instead the pipe_control.minimise version is being used.
........
r25327 | bugman | 2014-08-27 10:45:28 +0200 (Wed, 27 Aug 2014) | 11 lines

The minimisation statistics are now being reset for all analysis types.

The minimise.calcululate, minimise.grid_search, and minimise.execute user functions now all reset
the minimisation statistics for either the model or the Monte Carlo simulations prior to performing
any optimisation. This is required for both parallelised grid searches and repetitive optimisation
schemes to allow the result to overwrite an old result in all situations, as sometimes the original
chi-squared value is lower and the new result hence is rejected.

The equivalent functionality from the model-free minimise() API method has been removed.
........

2014-08-27 08:54:27 Tree
[r25327] by bugman

The minimisation statistics are now being reset for all analysis types.

The minimise.calcululate, minimise.grid_search, and minimise.execute user functions now all reset
the minimisation statistics for either the model or the Monte Carlo simulations prior to performing
any optimisation. This is required for both parallelised grid searches and repetitive optimisation
schemes to allow the result to overwrite an old result in all situations, as sometimes the original
chi-squared value is lower and the new result hence is rejected.

The equivalent functionality from the model-free minimise() API method has been removed.

2014-08-27 08:45:28 Tree
[r25326] by bugman

The model-free reset_min_stats() function has been replaced with the pipe_control.minimise version.

The specific_analyses.model_free.optimisation.reset_min_stats() function has been deleted and
instead the pipe_control.minimise version is being used.

2014-08-27 08:10:57 Tree
[r25325] by tlinnet

Modified wrapper function for curve_fit, to only change to list type, if the type is a ndarray.

task #7822(https://gna.org/task/index.php?7822): Implement user function to estimate R2eff and associated errors for exponential curve fitting.

2014-08-27 08:05:03 Tree
[r25324] by bugman

Improvements to the pipe_control.minimise.reset_min_stats() function.

The minimise statistics resetting is now more elegantly implemented. And the sim_index keyword
argument is accepted by the function and individual Monte Carlo simulation elements can now be
reset.

2014-08-27 08:01:54 Tree
[r25323] by bugman

Created an initial test suite data directory for a mixed R1rho + CPMG dispersion analysis.

The generate.py script will be extended in the future to generate both synthetic R1rho and CPMG data
for a common exchange process. Such a data combination should show some minor flaws in the current
design of the dispersion analysis and will help to solve these.

2014-08-27 07:47:13 Tree
[r25322] by tlinnet

Modified estimate R2eff module, to use C-code.

But systemtest Relax_disp.test_estimate_r2eff_error shows that the Jacobian is not correctly implemented
to be called in minfx.

2014-08-27 07:19:45 Tree
[r25321] by tlinnet

Modified profiling script, but it seems that the dfunc from target_functions.relax_fit does not work.

task #7822(https://gna.org/task/index.php?7822): Implement user function to estimate R2eff and associated errors for exponential curve fitting.

2014-08-27 06:50:56 Tree
Older >
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.