This is the mail archive of the gsl-discuss@sources.redhat.com mailing list for the GSL project.
Index Nav: | [Date Index] [Subject Index] [Author Index] [Thread Index] | |
---|---|---|
Message Nav: | [Date Prev] [Date Next] | [Thread Prev] [Thread Next] |
Other format: | [Raw text] |
Hi, thank you for the analysis. Apparently, my analytical calculation is able to guess the minimum reasonably good, so I need only 20000 iterations to find the minimum.
However, the strange thing is if I take the final minimized vector and run it through the conjugate gradient routine again, the energy (my function to minimize) keeps going down for a few more 1000s iterations (albeit it change only the 4th significant digit). Do you know what's wrong?
I follow the example in the GSL document and set my success criteria as
status = gsl_multimin_test_gradient (s->gradient, 1e-20);
The initial setup of the vector is
T = gsl_multimin_fdfminimizer_conjugate_pr; s = gsl_multimin_fdfminimizer_alloc (T, 3*nparts); gsl_multimin_fdfminimizer_set (s, &energy, x, 0.001, 1e-12);
For clarity, I used dimensionless units so ALL the nearest neighbor distances and the coupling constants are of the order of unity. nparts is the number of particles which is about 100000.
Could you give some advises? Thanks a lot in advance.
Toan
Index Nav: | [Date Index] [Subject Index] [Author Index] [Thread Index] | |
---|---|---|
Message Nav: | [Date Prev] [Date Next] | [Thread Prev] [Thread Next] |