projects
/
unres.git
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
ENDMDL cluster_wham/src
[unres.git]
/
source
/
unres
/
src_MD
/
checkder_p.F
diff --git
a/source/unres/src_MD/checkder_p.F
b/source/unres/src_MD/checkder_p.F
index
4d0379e
..
05ea860
100644
(file)
--- a/
source/unres/src_MD/checkder_p.F
+++ b/
source/unres/src_MD/checkder_p.F
@@
-284,8
+284,9
@@
C Check the gradient of the energy in Cartesian coordinates.
c call intcartderiv
c call checkintcartgrad
call zerograd
c call intcartderiv
c call checkintcartgrad
call zerograd
- aincr=1.0D-5
- write(iout,*) 'Calling CHECK_ECARTINT.'
+c aincr=1.0D-4
+ aincr=checkgrad_inc
+ write(iout,*) 'Calling CHECK_ECARTINT.',aincr
nf=0
icall=0
call geom_to_var(nvar,x)
nf=0
icall=0
call geom_to_var(nvar,x)
@@
-395,6
+396,7
@@
c call int_from_cart1(.false.)
if (.not.split_ene) then
call etotal(energia1(0))
etot1=energia1(0)
if (.not.split_ene) then
call etotal(energia1(0))
etot1=energia1(0)
+c write (iout,*) "i",i," etot",etot," etot1",etot1
else
!- split gradient
call etotal_long(energia1(0))
else
!- split gradient
call etotal_long(energia1(0))
@@
-441,6
+443,7
@@
c write (iout,*)
if (.not.split_ene) then
call etotal(energia1(0))
etot1=energia1(0)
if (.not.split_ene) then
call etotal(energia1(0))
etot1=energia1(0)
+c write (iout,*) "i",i," etot",etot," etot1",etot1
else
!- split gradient
call etotal_long(energia1(0))
else
!- split gradient
call etotal_long(energia1(0))
@@
-686,7
+689,7
@@
cd write (iout,'(i3,1pe14.4)') (i,gana(i),i=1,nvar)
call etotal(energia2(0))
etot2=energia2(0)
gg(i)=(etot2-etot1)/aincr
call etotal(energia2(0))
etot2=energia2(0)
gg(i)=(etot2-etot1)/aincr
- write (iout,*) i,etot1,etot2
+c write (iout,*) i,etot1,etot2
x(i)=xi
enddo
write (iout,'(/2a)')' Variable Numerical Analytical',
x(i)=xi
enddo
write (iout,'(/2a)')' Variable Numerical Analytical',