The CASINO forum

Vladimir_Konjkov wrote: Yes, I also wanted to improve significantly your inelegant code, and compiled profiling version that gave me 'empty' gmon.out.

FFLAGS_prof = -pg -no-pie fix generation of the gmon.out
Also mpif90 is the deprecated OpenMPI wrapper in favor of mpifort.

Hello Mike. In the gm2_bf SUBROUTINE evaluates the elements of the prim array in CASE of g-orbitals, it seems to me that in this function such calculations are not required. I apologize, I just automatically made the code review. Some guy tried to do that once to the Gaussian evaluators (he was talk...

Hello all It is well known that for effective DMC calculation one need, among other things, to choose the optimal time-step. Fort all-electron atomic systems, the linear regime to occur for time steps less than of the order τ= 1/(3Z^2), where Z is the largest atomic number occurring in the system. h...

Hello Mike

Additionally try to change the next line from dg(4,0)/dx, dg(4,0)/dy, dg(4,0)/dz code:

t14=6.d0*xx*(alphay2+fouralphaz2-1.d0) to t14=6.d0*xx*(alphay2-fouralphaz2-1.d0)

This helps the 'g-ane' example to pass the the numerical check of the analytic derivatives test.

Best, Vladimir.

Dear Vladimir, Thanks for the report. Do you have an example that fails? (The MDET examples in CASINO/examples/TEST seem to work OK, at least with the gfortran and NAG compilers.) Best wishes, Neil. Hello Neil. My example is in the attachment. I'm still using the old version, it works completely. V...

Hello CASINO developers. I've found after upgrade to v2.13.673 from v2.13.639 that newest version fails when running with ANY multi determinant WFN in parallel mod with message: --Job's stderr-- [vladimir-Kubuntu-16:17932] *** An error occurred in MPI_Bcast [vladimir-Kubuntu-16:17932] *** reported b...

I think that unit tests for gaussian WFN gradients, second derivatives and laplacians are strongly required and it is desirable to follow the DRY principle.

Hello Mike. It seems to me that all remaining analytic derivatives except first three are right. All first three are wrong similar way. Because g(4,0) = 3/8 * x^4 + 3/4 * x^2y^2 - 3 * x^2z^2 + 3/8 * y^4 - 3 * y^2z^2 + z^4 is symmetric upon x,y permutation, so dg(4,0)/dx + dg(4,0)/dy also should be b...

Hi Vladimir, OK - I've now had a look at this. You're correct that there is a problem in the analytic first derivatives of g functions. One can see that - for example - the 'g-ane' example in examples/generic/gauss_dfg/gaussian/g-ane fails the numerical check of the analytic derivatives enabled by ...

Hi Vladimir, I also found that the pgra(27) has a wrong sign, should be pgra(27)=840.d0*alpha*xyz*(yy-xx) Can you triple-check that before I have a look? M. OK let's look at the following piece if code: t16=xx*(3.d0+twoalphay2) t17=2.d0*alphax4+yy pgra(25)=420.d0*y*(t16-t17) t14=2.d0*alphay4-3.d0*y...