* Check the gradient of the virtual-bond and SC vectors in the internal
* coordinates.
*
- aincr=1.0d-7
- aincr2=5.0d-8
+ aincr=1.5d-7
+ aincr2=2.5d-8
call cartder
write (iout,'(a)') '**************** dx/dalpha'
write (iout,'(a)')
integer uiparm(1)
double precision urparm(1)
external fdum
- r_cut=2.0d0
- rlambd=0.3d0
+c r_cut=2.0d0
+c rlambd=0.3d0
icg=1
nf=0
nfl=0
c call intcartderiv
c call checkintcartgrad
call zerograd
- aincr=8.0D-6
+ aincr=8.0D-7
write(iout,*) 'Calling CHECK_ECARTINT.'
nf=0
icall=0
c write(iout,'(2i5,2(a,f15.10))')i,j," etot",etot," etot1",etot1
dc(j,i)=ddc(j)-aincr
call chainbuild_cart
+C print *,c(j,i)
c call int_from_cart1(.false.)
if (.not.split_ene) then
call etotal(energia1(0))