hdiff output

r30054/finalio.f90 2016-03-16 18:33:35.519080969 +0000 r30053/finalio.f90 2016-03-16 18:33:35.727083106 +0000
 28:     USE AMHGLOBALS, ONLY:NMRES,IRES 28:     USE AMHGLOBALS, ONLY:NMRES,IRES
 29:     USE BGUPMOD 29:     USE BGUPMOD
 30:     USE PERMU 30:     USE PERMU
 31:     USE MODHESS, ONLY : HESS, MASSWT 31:     USE MODHESS, ONLY : HESS, MASSWT
 32:  32: 
 33:     IMPLICIT NONE 33:     IMPLICIT NONE
 34:  34: 
 35:     !   MCP 35:     !   MCP
 36:     INTEGER III, I3,  GLY_COUNT, ID, NUMCRD, NUMPRO, NCPHST, GETUNIT, AMHUNIT1, AMHUNIT2, LUNIT, NSYMOPS 36:     INTEGER III, I3,  GLY_COUNT, ID, NUMCRD, NUMPRO, NCPHST, GETUNIT, AMHUNIT1, AMHUNIT2, LUNIT, NSYMOPS
 37:     INTEGER J1, J2, J3, J4, J5, MYUNIT2, I1, NDUMMY, MYUNIT3, NC, NRBS1, NRBS2, MYUNIT4 37:     INTEGER J1, J2, J3, J4, J5, MYUNIT2, I1, NDUMMY, MYUNIT3, NC, NRBS1, NRBS2, MYUNIT4
 38:     DOUBLE PRECISION RBCOORDS(NRBSITES*3), DCOORDS(3*NATOMS), EDUMMY, ITDET 38:     DOUBLE PRECISION RBCOORDS(NRBSITES*3), DCOORDS(3*NATOMS), EDUMMY
 39:     DOUBLE PRECISION P3(3,3), P(3), DU(3), RMI(3,3), DRMI(3,3), PI, PHI, THT, CHI 39:     DOUBLE PRECISION P3(3,3), P(3), DU(3), RMI(3,3), DRMI(3,3), PI, PHI, THT, CHI
 40:     !DOUBLE PRECISION, ALLOCATABLE :: XCOORDS(:), YCOORDS(:) 40:     !DOUBLE PRECISION, ALLOCATABLE :: XCOORDS(:), YCOORDS(:)
 41:     CHARACTER(LEN=13) J1CHAR,J1CHAR2                  !for gay-berne output files 41:     CHARACTER(LEN=13) J1CHAR,J1CHAR2                  !for gay-berne output files
 42:     CHARACTER(LEN=6) ZSTR 42:     CHARACTER(LEN=6) ZSTR
 43:     CHARACTER(LEN=2)  DUMMYSTR 43:     CHARACTER(LEN=2)  DUMMYSTR
 44:     DOUBLE PRECISION EulerPhi,EulerPsi,EulerTheta,EulerThetadeg,EulerPhiDeg,EulerPsiDeg  ! Euler angles for ellipsoids of revolution 44:     DOUBLE PRECISION EulerPhi,EulerPsi,EulerTheta,EulerThetadeg,EulerPhiDeg,EulerPsiDeg  ! Euler angles for ellipsoids of revolution
 45:  45: 
 46:     DOUBLE PRECISION EPS2, RAD, HEIGHT,sumx,sumy,sumz,CM(3),SYMOPS(120,3,3) 46:     DOUBLE PRECISION EPS2, RAD, HEIGHT,sumx,sumy,sumz,CM(3),SYMOPS(120,3,3)
 47:     LOGICAL :: GTEST 47:     LOGICAL :: GTEST
 48:     COMMON /CAPS/ EPS2, RAD, HEIGHT 48:     COMMON /CAPS/ EPS2, RAD, HEIGHT
237:               IF (ALLOCATED(HESS)) DEALLOCATE(HESS)237:               IF (ALLOCATED(HESS)) DEALLOCATE(HESS)
238:               ALLOCATE(HESS(3*NATOMS,3*NATOMS))238:               ALLOCATE(HESS(3*NATOMS,3*NATOMS))
239:               IF(NSPECIES(0)>1) & ! ds656> this also sets ATMASS239:               IF(NSPECIES(0)>1) & ! ds656> this also sets ATMASS
240:                    CALL SET_ATOMLISTS(QMINT(J1,1:NATOMS),1)240:                    CALL SET_ATOMLISTS(QMINT(J1,1:NATOMS),1)
241:               CALL POTENTIAL(QMINP(J1,1:3*NATOMS), EVALS, EDUMMY, &241:               CALL POTENTIAL(QMINP(J1,1:3*NATOMS), EVALS, EDUMMY, &
242:                    .TRUE., .TRUE.) ! EVALS is dummy for gradient here242:                    .TRUE., .TRUE.) ! EVALS is dummy for gradient here
243:               IF(ABS(EDUMMY-QMIN(J1)) > ECONV) THEN243:               IF(ABS(EDUMMY-QMIN(J1)) > ECONV) THEN
244:                  WRITE(MYUNIT,'(A)') &244:                  WRITE(MYUNIT,'(A)') &
245:                       'finalio> WARNING: Unexpected energy value!'245:                       'finalio> WARNING: Unexpected energy value!'
246:               ENDIF246:               ENDIF
247:               CALL MASSWT() ! Weigh the Hessian by mass. 247:               CALL MASSWT() 
248:               !248:               !
249:               WRITE(MYUNIT, '(A,I6)') &249:               WRITE(MYUNIT, '(A,I6)') &
250:                    'finalio> Using DSYEV to calculate Hessian eigenvalues for minimum ', J1250:                    'finalio> Using DSYEV to calculate Hessian eigenvalues for minimum ', J1
251:               CALL DSYEV('V','U',3*NATOMS,HESS,3*NATOMS,EVALS,&251:               CALL DSYEV('V','U',3*NATOMS,HESS,3*NATOMS,EVALS,&
252:                    WORK,9*NATOMS,INFO)252:                    WORK,9*NATOMS,INFO)
253:               WRITE(MYUNIT,'(A)',ADVANCE='NO') 'finalio> Eigenvalues:' 
254:               DO J2=1,3*NATOMS 
255:                  WRITE(MYUNIT,'(1X, E12.5E2)',ADVANCE='NO') EVALS(J2) 
256:               ENDDO 
257:               WRITE(MYUNIT,*) 
258:               ! Log prod. of frequencies is log prod. of sqrt of evals 253:               ! Log prod. of frequencies is log prod. of sqrt of evals 
259:               LOG_PROD = &254:               LOG_PROD = &
260:                    SUM(DLOG(EVALS(NUM_ZERO_EVS+1:3*NATOMS)))255:                    SUM(DLOG(EVALS(NUM_ZERO_EVS+1:3*NATOMS)))
261:               WRITE(MYUNIT,'(A,G20.10)') &256:               WRITE(MYUNIT,'(A,G20.10)') &
262:                    'finalio> Log product of non-zero eigenvalues: ', &257:                    'finalio> Log product of frequencies= ', &
263:                    LOG_PROD258:                    LOG_PROD
264:               CALL MYINERTIA(QMINP(J1,1:3*NATOMS),ITDET)259:               WRITE(MYUNIT2, '(2F20.10,I6)') &
265:               WRITE(MYUNIT,'(A,G20.10)') &260:                    EDUMMY, LOG_PROD, NSYMOPS
266:                    'finalio> Determinant of inertia tensor: ', ITDET 
267:               WRITE(MYUNIT2, '(2F20.10,I6,F20.10)') & 
268:                    EDUMMY, LOG_PROD, NSYMOPS, ITDET 
269:               !261:               !
270:            ENDIF262:            ENDIF
271:            !263:            !
272:         ELSE ! <ds656264:         ELSE ! <ds656
273:            WRITE(MYUNIT2,10) J1, QMIN(J1), FF(J1), NPCALL_QMIN(J1)265:            WRITE(MYUNIT2,10) J1, QMIN(J1), FF(J1), NPCALL_QMIN(J1)
274: 10         FORMAT('Energy of minimum ',I6,'=',G20.10, &266: 10         FORMAT('Energy of minimum ',I6,'=',G20.10, &
275:                 ' first found at step ',I8,' after ',I20,' function calls')267:                 ' first found at step ',I8,' after ',I20,' function calls')
276:         ENDIF268:         ENDIF
277:         !269:         !
278:         IF (MSORIGT.OR.FRAUSIT) THEN270:         IF (MSORIGT.OR.FRAUSIT) THEN


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0