Menu

Commit [r25396]  Maximize  Restore  History

Modified systemtest Relax_disp.verify_estimate_r2eff_err_compare_mc, to try the difference between the Jacobian.

The results are:

Printing the estimated R2eff error as function of estimation from Co-variance and number of Monte-Carlo simulations.
-2 0.069 0.081 0.085 0.092 0.085 0.074 0.083 0.069 0.066 0.074 0.025 0.035 0.018 0.016 sum= 0.874
-1 0.052 0.023 0.034 0.021 0.020 0.041 0.030 0.028 0.011 0.163 0.048 0.012 0.009 0.010 sum= 0.502
0 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 sum= 0.000
10 0.044 0.043 0.044 0.049 0.053 0.027 0.048 0.043 0.032 0.038 0.016 0.014 0.013 0.009 sum= 0.473
20 0.038 0.044 0.051 0.050 0.049 0.040 0.046 0.035 0.037 0.034 0.015 0.020 0.009 0.008 sum= 0.475
30 0.029 0.050 0.050 0.054 0.041 0.040 0.030 0.041 0.024 0.038 0.015 0.012 0.009 0.008 sum= 0.441
40 0.031 0.041 0.043 0.044 0.046 0.039 0.040 0.036 0.032 0.036 0.014 0.020 0.008 0.009 sum= 0.439

0.9 ++------+------+-------+-------+------+-------+-------+------+------++
+ A + + + R2eff error as function of MC number+**A*** +
0.8 ++ * ++
| * |
0.7 ++ * ++
| * |
0.6 ++ * ++
| * |
0.5 ++ A ++
| * A**************A*************** |
| * ** A**************A
0.4 ++ * * ++
| * ** |
0.3 ++ * * ++
| * ** |
0.2 ++ * * ++
| * ** |
0.1 ++ * * ++
+ *** + + + + + + + +
0 ++------A------+-------+-------+------+-------+-------+------+------++
-5 0 5 10 15 20 25 30 35 40

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

tlinnet 2014-08-28

changed /trunk/test_suite/system_tests/relax_disp.py
/trunk/test_suite/system_tests/relax_disp.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.