hdiff output

r31173/finalio.f90 2016-09-20 11:30:09.861586529 +0100 r31172/finalio.f90 2016-09-20 11:30:10.617596638 +0100
 52:   DOUBLE PRECISION EulerPhi,EulerPsi,EulerTheta,EulerThetadeg,EulerPhiDeg,EulerPsiDeg  ! Euler angles for ellipsoids of revolution 52:   DOUBLE PRECISION EulerPhi,EulerPsi,EulerTheta,EulerThetadeg,EulerPhiDeg,EulerPsiDeg  ! Euler angles for ellipsoids of revolution
 53:    53:   
 54:   DOUBLE PRECISION EPS2, RAD, HEIGHT,sumx,sumy,sumz,CM(3),SYMOPS(120,3,3) 54:   DOUBLE PRECISION EPS2, RAD, HEIGHT,sumx,sumy,sumz,CM(3),SYMOPS(120,3,3)
 55:   LOGICAL :: GTEST 55:   LOGICAL :: GTEST
 56:   COMMON /CAPS/ EPS2, RAD, HEIGHT 56:   COMMON /CAPS/ EPS2, RAD, HEIGHT
 57:    57:   
 58:   CHARACTER(LEN=4) :: FPGRP 58:   CHARACTER(LEN=4) :: FPGRP
 59:   CHARACTER(LEN=6) :: CRMS 59:   CHARACTER(LEN=6) :: CRMS
 60:   CHARACTER(LEN=20) :: MYFILENAME2, ISTR, DBNUM, MYFILENAME3 60:   CHARACTER(LEN=20) :: MYFILENAME2, ISTR, DBNUM, MYFILENAME3
 61:   CHARACTER(LEN=15), ALLOCATABLE :: DBNAME(:) 61:   CHARACTER(LEN=15), ALLOCATABLE :: DBNAME(:)
 62:   CHARACTER(LEN=80) :: tempstring 
 63:    62:   
 64:   CHARACTER(LEN=6) :: J1_STRING, MYNODE_STRING 63:   CHARACTER(LEN=6) :: J1_STRING, MYNODE_STRING
 65:    64:   
 66:   !  AMH 65:   !  AMH
 67:   CHARACTER(LEN=3) :: RES_TYPE 66:   CHARACTER(LEN=3) :: RES_TYPE
 68:   CHARACTER(LEN=2) :: ATOM_TYPE 67:   CHARACTER(LEN=2) :: ATOM_TYPE
 69:   CHARACTER*1 COUNTTT 68:   CHARACTER*1 COUNTTT
 70:   INTEGER COUNTT 69:   INTEGER COUNTT
 71:   DOUBLE PRECISION  PPPCORD(NMRES*3*3,3,3,5) 70:   DOUBLE PRECISION  PPPCORD(NMRES*3*3,3,3,5)
 72:   EXTERNAL NUM_TO_CHAR 71:   EXTERNAL NUM_TO_CHAR
606:                 LEN('lowest.'//TRIM(ADJUSTL(MYNODE_STRING))), &605:                 LEN('lowest.'//TRIM(ADJUSTL(MYNODE_STRING))), &
607:                 LOGICAL(.FALSE.,KIND=1))606:                 LOGICAL(.FALSE.,KIND=1))
608:         ELSE 607:         ELSE 
609:            ! Pass the name, length of the string and whether or not to write a header.608:            ! Pass the name, length of the string and whether or not to write a header.
610:            CALL AMBER12_WRITE_XYZ(QMINP(J1,:), 'lowest', LEN('lowest'), LOGICAL(.FALSE.,KIND=1))609:            CALL AMBER12_WRITE_XYZ(QMINP(J1,:), 'lowest', LEN('lowest'), LOGICAL(.FALSE.,KIND=1))
611:         END IF610:         END IF
612:      ELSE IF (AMBERT) THEN611:      ELSE IF (AMBERT) THEN
613:         ! sf344> write out coordinates612:         ! sf344> write out coordinates
614:         COORDS1(1:3*NATOMS) = QMINP(J1,1:3*NATOMS)613:         COORDS1(1:3*NATOMS) = QMINP(J1,1:3*NATOMS)
615:         IF (DUMPSTRUCTURES) THEN614:         IF (DUMPSTRUCTURES) THEN
 615:            CALL AMBERFINALIO(j1,MYUNIT2,MYNODE+1,'0',0,COORDS1(1:3*NATOMS))
616:            WRITE(J1CHAR2,'(I3)') J1616:            WRITE(J1CHAR2,'(I3)') J1
617:            IF (MPIT) THEN617:            IF (MPIT) THEN
618:               CALL AMBERFINALIO(j1,20,MYNODE+1,tempstring,0,COORDS1(1:3*NATOMS)) 
619:               WRITE (ISTR, '(I10)') MYNODE+1618:               WRITE (ISTR, '(I10)') MYNODE+1
620:               MYUNIT3=GETUNIT()619:               MYUNIT3=GETUNIT()
621:               MYFILENAME2='coords.'//TRIM(ADJUSTL(J1CHAR2))//'.'//trim(adjustl(istr))620:               MYFILENAME2='coords.'//TRIM(ADJUSTL(J1CHAR2))//'.'//trim(adjustl(istr))
622:               OPEN(MYUNIT3,FILE=trim(adjustl(MYFILENAME2)), STATUS="unknown", form="formatted")621:               OPEN(MYUNIT3,FILE=trim(adjustl(MYFILENAME2)), STATUS="unknown", form="formatted")
623:            ELSE622:            ELSE
624:               CALL AMBERFINALIO(j1,MYUNIT2,MYNODE+1,tempstring,0,COORDS1(1:3*NATOMS)) 
625:               MYUNIT3=GETUNIT()623:               MYUNIT3=GETUNIT()
626:               MYFILENAME2='coords.'//TRIM(ADJUSTL(J1CHAR2))624:               MYFILENAME2='coords.'//TRIM(ADJUSTL(J1CHAR2))
627:               OPEN(MYUNIT3,FILE=trim(adjustl(MYFILENAME2)), STATUS="unknown", form="formatted")625:               OPEN(MYUNIT3,FILE=trim(adjustl(MYFILENAME2)), STATUS="unknown", form="formatted")
628:            END IF626:            END IF
629:            DO J2=1,NATOMS627:            DO J2=1,NATOMS
630:               WRITE(MYUNIT3,'(3F28.20)') QMINP(J1,3*(J2-1)+1),QMINP(J1,3*(J2-1)+2),QMINP(J1,3*(J2-1)+3)628:               WRITE(MYUNIT3,'(3F28.20)') QMINP(J1,3*(J2-1)+1),QMINP(J1,3*(J2-1)+2),QMINP(J1,3*(J2-1)+3)
631:            ENDDO629:            ENDDO
632:            CLOSE(MYUNIT3)630:            CLOSE(MYUNIT3)
633:         ELSE631:         ELSE
634:            DO I1=1,NATOMS632:            DO I1=1,NATOMS


r31173/modcudalbfgs.F90 2016-09-20 11:30:10.105589795 +0100 r31172/modcudalbfgs.F90 2016-09-20 11:30:10.873600081 +0100
201:         C_CUDAPOT = CUDAPOT201:         C_CUDAPOT = CUDAPOT
202:         C_DEBUG = DEBUG202:         C_DEBUG = DEBUG
203:         C_CUDATIMET = CUDATIMET203:         C_CUDATIMET = CUDATIMET
204:         C_MAXBFGS = MAXBFGS204:         C_MAXBFGS = MAXBFGS
205:         C_MAXERISE = MAXERISE205:         C_MAXERISE = MAXERISE
206:         C_ATOMRIGIDCOORDT = ATOMRIGIDCOORDT206:         C_ATOMRIGIDCOORDT = ATOMRIGIDCOORDT
207:         C_DEGFREEDOMS = DEGFREEDOMS207:         C_DEGFREEDOMS = DEGFREEDOMS
208:         C_NRIGIDBODY = NRIGIDBODY208:         C_NRIGIDBODY = NRIGIDBODY
209:         C_MAXSITE = MAXSITE209:         C_MAXSITE = MAXSITE
210:         C_COLDFUSIONLIMIT = COLDFUSIONLIMIT210:         C_COLDFUSIONLIMIT = COLDFUSIONLIMIT
211:         C_COLDFUSION = .FALSE. 
212:         C_DGUESS = DGUESS211:         C_DGUESS = DGUESS
213:         C_BQMAX = BQMAX212:         C_BQMAX = BQMAX
214:         C_FREEZE = FREEZE213:         C_FREEZE = FREEZE
215:         C_NFREEZE = NFREEZE214:         C_NFREEZE = NFREEZE
216:         PROJECT = .FALSE.215:         PROJECT = .FALSE.
217:         C_AACONVERGENCET = AACONVERGENCET216:         C_AACONVERGENCET = AACONVERGENCET
218: 217: 
219:         ! 'C_' prefix denotes those variables which have intent out or inout or are copies of those from common blocks/modules 218:         ! 'C_' prefix denotes those variables which have intent out or inout or are copies of those from common blocks/modules 
220:         CALL CUDA_LBFGS(N, C_XCOORDS, EPS, C_MFLAG, C_ENERGY, ITMAX, C_ITDONE, C_MAXBFGS, C_MAXERISE, C_RMS, C_CUDAPOT, & 219:         CALL CUDA_LBFGS(N, C_XCOORDS, EPS, C_MFLAG, C_ENERGY, ITMAX, C_ITDONE, C_MAXBFGS, C_MAXERISE, C_RMS, C_CUDAPOT, & 
221:                        C_DEBUG, C_CUDATIMET, ECALLS, C_COLDFUSION, C_COLDFUSIONLIMIT, C_DGUESS, MUPDATE, & 220:                        C_DEBUG, C_CUDATIMET, ECALLS, C_COLDFUSION, C_COLDFUSIONLIMIT, C_DGUESS, MUPDATE, & 


r31173/potential.f90 2016-09-20 11:30:10.353593118 +0100 r31172/potential.f90 2016-09-20 11:30:11.137603622 +0100
114: !      WRITE(88,'(G20.10)') (ACCRAT(J1), J1=1, NPAR)114: !      WRITE(88,'(G20.10)') (ACCRAT(J1), J1=1, NPAR)
115: !      WRITE(88,'(G20.10)') (STEP(J1), J1=1, NPAR)115: !      WRITE(88,'(G20.10)') (STEP(J1), J1=1, NPAR)
116: !      WRITE(88,'(G20.10)') (ASTEP(J1), J1=1, NPAR)116: !      WRITE(88,'(G20.10)') (ASTEP(J1), J1=1, NPAR)
117: !      WRITE(88,'(G20.10)') (OSTEP(J1), J1=1, NPAR)117: !      WRITE(88,'(G20.10)') (OSTEP(J1), J1=1, NPAR)
118: !      CALL SYSTEM('rm ssave')118: !      CALL SYSTEM('rm ssave')
119: !      STOP119: !      STOP
120: !   END IF120: !   END IF
121: 121: 
122: ! Count the number of potential calls.122: ! Count the number of potential calls.
123:    NPCALL=NPCALL+1123:    NPCALL=NPCALL+1
 124: 
124: 10 CONTINUE125: 10 CONTINUE
125: 126: 
126:    IF (MULTIPOTT) THEN127:    IF (MULTIPOTT) THEN
127:        IF (RIGIDINIT .AND. (.NOT. ATOMRIGIDCOORDT)) THEN128:        IF (RIGIDINIT .AND. (.NOT. ATOMRIGIDCOORDT)) THEN
128:            XRIGIDCOORDS(1:DEGFREEDOMS)=X(1:DEGFREEDOMS)129:            XRIGIDCOORDS(1:DEGFREEDOMS)=X(1:DEGFREEDOMS)
129:            CALL TRANSFORMRIGIDTOC(1, NRIGIDBODY, X, XRIGIDCOORDS)               130:            CALL TRANSFORMRIGIDTOC(1, NRIGIDBODY, X, XRIGIDCOORDS)               
130:        END IF131:        END IF
131: 132: 
132: !       CALL RAD(X, GRAD, EREAL, GRADT)133: !       CALL RAD(X, GRAD, EREAL, GRADT)
133:        CALL MULTIPOT_CALL(X, GRAD, EREAL, GRADT, SECT)134:        CALL MULTIPOT_CALL(X, GRAD, EREAL, GRADT, SECT)


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0