projects
/
unres.git
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
corrections to changes form nostromo
[unres.git]
/
source
/
unres
/
src_MD-M
/
gradient_p.F
diff --git
a/source/unres/src_MD-M/gradient_p.F
b/source/unres/src_MD-M/gradient_p.F
index
b5ba82a
..
7295d2e
100644
(file)
--- a/
source/unres/src_MD-M/gradient_p.F
+++ b/
source/unres/src_MD-M/gradient_p.F
@@
-292,11
+292,9
@@
c If performing constraint dynamics, add the gradients of the constraint energy
do i=1,nres-3
gloc(i,icg)=gloc(i,icg)+dugamma(i)
enddo
do i=1,nres-3
gloc(i,icg)=gloc(i,icg)+dugamma(i)
enddo
- write(iout,*) "TU JESTEM?"
do i=1,nres-2
gloc(nphi+i,icg)=gloc(nphi+i,icg)+dutheta(i)
enddo
do i=1,nres-2
gloc(nphi+i,icg)=gloc(nphi+i,icg)+dutheta(i)
enddo
- write(iout,*) "TU JESTEM?"
endif
#ifdef TIMING
time01=MPI_Wtime()
endif
#ifdef TIMING
time01=MPI_Wtime()
@@
-310,7
+308,7
@@
cd write(iout,*) 'calling int_to_cart'
#ifdef DEBUG
write (iout,*) "gcart, gxcart, gloc before int_to_cart"
#endif
#ifdef DEBUG
write (iout,*) "gcart, gxcart, gloc before int_to_cart"
#endif
- do i=1,nct
+ do i=0,nct
do j=1,3
gcart(j,i)=gradc(j,i,icg)
gxcart(j,i)=gradx(j,i,icg)
do j=1,3
gcart(j,i)=gradc(j,i,icg)
gxcart(j,i)=gradx(j,i,icg)
@@
-351,7
+349,7
@@
C-------------------------------------------------------------------------
C
C Initialize Cartesian-coordinate gradient
C
C
C Initialize Cartesian-coordinate gradient
C
- do i=1,nres
+ do i=-1,nres
do j=1,3
gvdwx(j,i)=0.0D0
gradx_scp(j,i)=0.0D0
do j=1,3
gvdwx(j,i)=0.0D0
gradx_scp(j,i)=0.0D0
@@
-383,11
+381,25
@@
C
gradx(j,i,icg)=0.0d0
gscloc(j,i)=0.0d0
gsclocx(j,i)=0.0d0
gradx(j,i,icg)=0.0d0
gscloc(j,i)=0.0d0
gsclocx(j,i)=0.0d0
+ gliptranc(j,i)=0.0d0
+ gliptranx(j,i)=0.0d0
+ gradafm(j,i)=0.0d0
do intertyp=1,3
gloc_sc(intertyp,i,icg)=0.0d0
enddo
enddo
enddo
do intertyp=1,3
gloc_sc(intertyp,i,icg)=0.0d0
enddo
enddo
enddo
+c
+c Initialize the gradients of local restraints
+c
+ do i=1,nres
+ dutheta(i)=0.0d0
+ dugamma(i)=0.0d0
+ do j=1,3
+ duscdiff(j,i)=0.0d0
+ duscdiffx(j,i)=0.0d0
+ enddo
+ enddo
C
C Initialize the gradient of local energy terms.
C
C
C Initialize the gradient of local energy terms.
C