hdiff output

r29560/multipot.f90 2015-12-03 17:30:12.019233755 +0000 r29559/multipot.f90 2015-12-03 17:30:12.411238903 +0000
296:     NDEGS = 3*TMP_NATOMS296:     NDEGS = 3*TMP_NATOMS
297: 297: 
298:     DO J1 = 1, NDEGS  ! Loop over all the atoms with this kind of potential298:     DO J1 = 1, NDEGS  ! Loop over all the atoms with this kind of potential
299:         TMP_X(J1) = X(DEGS_BY_POT(POTID,J1))299:         TMP_X(J1) = X(DEGS_BY_POT(POTID,J1))
300:     ENDDO300:     ENDDO
301: 301: 
302:     IF(DEBUG) THEN302:     IF(DEBUG) THEN
303:         WRITE(*,*) "Calling potential", POTTYPES(J1)303:         WRITE(*,*) "Calling potential", POTTYPES(J1)
304:         WRITE(*,*) "Degrees of freedom being used:"304:         WRITE(*,*) "Degrees of freedom being used:"
305:         WRITE(*,*) DEGS_BY_POT(POTID,:NDEGS)305:         WRITE(*,*) DEGS_BY_POT(POTID,:NDEGS)
306:     ENDIF 
307: 306: 
308:     ! The advantage of this slightly tortuous method is that we only need to call POT once. Compare with the pairwise307:     ! The advantage of this slightly tortuous method is that we only need to call POT once. Compare with the pairwise
309:     ! routine, where we need to make a function call for every pair of atoms in the system.308:     ! routine, where we need to make a function call for every pair of atoms in the system.
310:     CALL POT(TMP_NATOMS, TMP_X, POTPARAMS(POTID,:), TMP_E, TMP_G, TMP_HESS, GTEST, STEST)309:     CALL POT(TMP_NATOMS, TMP_X, POTPARAMS(POTID,:), TMP_E, TMP_G, TMP_HESS, GTEST, STEST)
311: 310: 
312:     ENERGY = ENERGY + TMP_E*POTSCALES(POTID)311:     ENERGY = ENERGY + TMP_E*POTSCALES(POTID)
313: 312: 
314:     ! Unpack the gradient, if it is being calculated313:     ! Unpack the gradient, if it is being calculated
315:     IF(GTEST) THEN314:     IF(GTEST) THEN
316:         DO J1 = 1,NDEGS315:         DO J1 = 1,NDEGS


r29560/potential.f 2015-12-03 17:30:12.219236381 +0000 r29559/potential.f 2015-12-03 17:30:12.603241424 +0000
 25:          USE MODCHARMM 25:          USE MODCHARMM
 26:          use PORFUNCS 26:          use PORFUNCS
 27:          use SDWATER, ONLY : SDPOTENTIAL, SDGRAD, SDHESS 27:          use SDWATER, ONLY : SDPOTENTIAL, SDGRAD, SDHESS
 28:          USE MCY, ONLY : MCYPOT=>POTENTIAL 28:          USE MCY, ONLY : MCYPOT=>POTENTIAL
 29:          USE BOWMANWATER, ONLY : BOWMANPOT 29:          USE BOWMANWATER, ONLY : BOWMANPOT
 30:          USE FINITE_DIFFERENCES 30:          USE FINITE_DIFFERENCES
 31:          USE MODAMBER9,only : ifswitch,goodstructure1,irespa,cisarray1,checkcistransalways,checkcistransalwaysdna, 31:          USE MODAMBER9,only : ifswitch,goodstructure1,irespa,cisarray1,checkcistransalways,checkcistransalwaysdna,
 32:      1   checkcistransalwaysrna 32:      1   checkcistransalwaysrna
 33:          ! hk286 33:          ! hk286
 34:          USE GENRIGID 34:          USE GENRIGID
 35:          USE MULTIPOT 
 36:          USE CHIRALITY, ONLY: CIS_TRANS_CHECK, CHIRALITY_CHECK 35:          USE CHIRALITY, ONLY: CIS_TRANS_CHECK, CHIRALITY_CHECK
 37:          USE AMBER12_INTERFACE_MOD, ONLY: AMBER12_ENERGY_AND_GRADIENT, POT_ENE_REC_C, AMBER12_NUM_HESS 36:          USE AMBER12_INTERFACE_MOD, ONLY: AMBER12_ENERGY_AND_GRADIENT, POT_ENE_REC_C, AMBER12_NUM_HESS
 38:          USE MODCUDALBFGS, ONLY: CUDA_ENEGRAD_WRAPPER 37:          USE MODCUDALBFGS, ONLY: CUDA_ENEGRAD_WRAPPER
 39:          ! use AMHGLOBALS 38:          ! use AMHGLOBALS
 40:  39: 
 41:          IMPLICIT NONE 40:          IMPLICIT NONE
 42:  41: 
 43:          DOUBLE PRECISION, DIMENSION(3*NATOMS) :: COORDS 42:          DOUBLE PRECISION, DIMENSION(3*NATOMS) :: COORDS
 44:          DOUBLE PRECISION ENERGY 43:          DOUBLE PRECISION ENERGY
 45:          DOUBLE PRECISION, DIMENSION(3*NATOMS) :: VNEW 44:          DOUBLE PRECISION, DIMENSION(3*NATOMS) :: VNEW
457:                   IF (DEBUG) PRINT *, ' potential> hessian created'456:                   IF (DEBUG) PRINT *, ' potential> hessian created'
458:                   ! PRINT '(A)','Hessian:'457:                   ! PRINT '(A)','Hessian:'
459:                   ! PRINT '(6G20.10)',HESS(1:NOPT,1:NOPT)458:                   ! PRINT '(6G20.10)',HESS(1:NOPT,1:NOPT)
460:                ENDIF459:                ENDIF
461:             ENDIF460:             ENDIF
462:             IF (PTEST) THEN461:             IF (PTEST) THEN
463:                WRITE(*,10) ' Energy for last cycle=',ENERGY462:                WRITE(*,10) ' Energy for last cycle=',ENERGY
464:                WRITE(ESTRING,10) ' Energy for last cycle=',ENERGY463:                WRITE(ESTRING,10) ' Energy for last cycle=',ENERGY
465:             ENDIF464:             ENDIF
466:             465:             
467:          ELSE IF(MULTIPOTT) THEN 
468:             CALL MULTIPOT_CALL(COORDS,VNEW,ENERGY,GTEST,SSTEST) 
469:             IF (PTEST) THEN 
470:                WRITE(*,10) ' potential> Energy for last cycle=',ENERGY,' ' 
471:                WRITE(ESTRING,10) 'Energy for last cycle=',ENERGY,' ' 
472:             ENDIF 
473: 466: 
474:          ELSE IF(EX1DT) THEN467:          ELSE IF(EX1DT) THEN
475:             468:             
476:             CALL example1Dpotential(COORDS,VNEW,ENERGY,GTEST,STEST)469:             CALL example1Dpotential(COORDS,VNEW,ENERGY,GTEST,STEST)
477:             470:             
478:             IF (PTEST) THEN471:             IF (PTEST) THEN
479:                WRITE(*,10) ' potential> Energy for last cycle=',ENERGY,' '472:                WRITE(*,10) ' potential> Energy for last cycle=',ENERGY,' '
480:                WRITE(ESTRING,10) 'Energy for last cycle=',ENERGY,' '473:                WRITE(ESTRING,10) 'Energy for last cycle=',ENERGY,' '
481:             ENDIF474:             ENDIF
482: 475: 


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0