fractional power basis issue on restart

General discussion of the Cambridge quantum Monte Carlo code CASINO; how to install and setup; how to use it; what it does; applications.

fractional power basis issue on restart

Postby Vladimir_Konjkov » Sat Jun 02, 2018 4:38 pm

Hello CASINO developers

I've tried to use fractional power basis in the form r/(r+a) power in the Rank: [ 2, 0 ] term and in the Rank: [ 2, 1 ] for e-e distance only.
I successfully optimized JASTROW coefficients by EMIN, but when I try to restart calculation with optimized coefficients I've got unphysical E value.

example for B atom in attachment

Best Vladimir.

P.S.
same for 1/(r+a) power and r/(r^b+a) power in the Rank: [ 2, 1 ] for e-e distance. I tried to find a balance between the VMC optimization accuracy and the time of JASTROW evaluation in DMC calculation.
I wrote down a list of arbitrary combinations of JASTROW basic functions and cutoffs, then ran an automatic script to compute DMC timing for them all. Unfortunately, I have not even reached the middle of the list.
Attachments
b_8_b4.tgz
(115.96 KiB) Downloaded 73 times
In Soviet Russia Casino play you.
Vladimir_Konjkov
 
Posts: 123
Joined: Wed Apr 15, 2015 3:14 pm

Re: fractional power basis issue on restart

Postby Pablo_Lopez_Rios » Sun Jun 03, 2018 6:52 pm

Hi Vladimir,

I see that in your first run you get the warning:
Code: Select all
Warning: [REDO_GAUSSIAN_ELIMINATION] non-linear parameters produce 1 zero(es)
 in the coefficient matrix for the constraint equations.  We recommend that
 you stop this calculation, change the values of any basis parameters involved
 to non-integer values, all slightly different from each other, and re-run to
 avoid possible problems.  This run will continue in the hope that the problem
 will go away during optimization.  The affected equations and parameters are
 listed below.

Zero #1 encountered at equation #17, parameter #26.

which hints at a problem with solving the constraint equations at the initial parameter values.

I would suggest you set the initial values of a and L to different, non-integer values, e.g., 3.1, 3.2, 3.3 and 3.4 (or start them closer to their optimized values if you want). This ought to work around the issue.

Best,
Pablo
Hey there! I am using CASINO.
Pablo_Lopez_Rios
 
Posts: 40
Joined: Thu Jan 30, 2014 1:25 am

Re: fractional power basis issue on restart

Postby Vladimir_Konjkov » Mon Jun 04, 2018 1:49 am

Hello Pablo.

Oh magic! How did you know?

now I can compare the speed of the calculation, thanks.

Vladimir.
In Soviet Russia Casino play you.
Vladimir_Konjkov
 
Posts: 123
Joined: Wed Apr 15, 2015 3:14 pm


Return to The CASINO program

Who is online

Users browsing this forum: No registered users and 3 guests

cron