hdiff output

r30265/commons.f90 2016-04-14 11:30:11.825630136 +0100 r30264/commons.f90 2016-04-14 11:30:12.737642657 +0100
 75:      &                 PEPSILON1(3),PSCALEFAC1(2),PSCALEFAC2(2),PYA11(3),PYA21(3),PYA12(3),PYA22(3), & 75:      &                 PEPSILON1(3),PSCALEFAC1(2),PSCALEFAC2(2),PYA11(3),PYA21(3),PYA12(3),PYA22(3), &
 76:      &                 PEPSILONATTR(2),PSIGMAATTR(2), PYOVERLAPTHRESH, PYCFTHRESH, LJSITECOORDS(3), & 76:      &                 PEPSILONATTR(2),PSIGMAATTR(2), PYOVERLAPTHRESH, PYCFTHRESH, LJSITECOORDS(3), &
 77:      &                 MSTART,MFINISH,MBSTART1,MBFINISH1,MBSTART2,MBFINISH2,MBHEIGHT1,MBHEIGHT2,ME1,ME2,ME3, & 77:      &                 MSTART,MFINISH,MBSTART1,MBFINISH1,MBSTART2,MBFINISH2,MBHEIGHT1,MBHEIGHT2,ME1,ME2,ME3, &
 78:      &                 BSPTQMAX, BSPTQMIN, PFORCE, CSMNORM, CSMGUIDENORM, CSMEPS, PERCCUT, & 78:      &                 BSPTQMAX, BSPTQMIN, PFORCE, CSMNORM, CSMGUIDENORM, CSMEPS, PERCCUT, &
 79:      &                 LOWESTE, PERTSTEP, GCPLUS, & 79:      &                 LOWESTE, PERTSTEP, GCPLUS, &
 80:      &                 KINT, INTFREEZETOL, IMSEPMIN, IMSEPMAX, CONCUTABS, CONCUTFRAC, & 80:      &                 KINT, INTFREEZETOL, IMSEPMIN, IMSEPMAX, CONCUTABS, CONCUTFRAC, &
 81:      &                 LPDGEOMDIFFTOL, INTCONFRAC, MAXCONE, INTRMSTOL, BFGSTSTOL, ORBITTOL, & 81:      &                 LPDGEOMDIFFTOL, INTCONFRAC, MAXCONE, INTRMSTOL, BFGSTSTOL, ORBITTOL, &
 82:      &                 INTCONSTRAINTTOL, INTCONSTRAINTDEL, RBCUTOFF, INTCONSTRAINTREP, INTCONSTRAINREPCUT, & 82:      &                 INTCONSTRAINTTOL, INTCONSTRAINTDEL, RBCUTOFF, INTCONSTRAINTREP, INTCONSTRAINREPCUT, &
 83:      &                 INTLJTOL, INTLJDEL, INTLJEPS, REPCON, INTDGUESS, CHECKREPCUTOFF, INTMINFAC, FREEZETOL, & 83:      &                 INTLJTOL, INTLJDEL, INTLJEPS, REPCON, INTDGUESS, CHECKREPCUTOFF, INTMINFAC, FREEZETOL, &
 84:      &                 LOCALPERMCUT, LOCALPERMCUT2, INTCONCUT, QCIRADSHIFT, MLPLAMBDA, & 84:      &                 LOCALPERMCUT, LOCALPERMCUT2, INTCONCUT, QCIRADSHIFT, MLPLAMBDA, &
 85:      &                 CAPSIDRHO,CAPSIDEPS,SIGMAPENT,RADPENT,SIGMAHEX,RADHEX,SIGMAPH, KLIM, SCA 85:      &                 CAPSIDRHO,CAPSIDEPS,SIGMAPENT,RADPENT,SIGMAHEX,RADHEX,SIGMAPH
 86:  86: 
 87:       LOGICAL DEBUG, TARGET, MORSET, CUTT, SEEDT, CENT, TSALLIST, FREEZECORE, NEWJUMP, RENORM, CAPSID, FREEZE, & 87:       LOGICAL DEBUG, TARGET, MORSET, CUTT, SEEDT, CENT, TSALLIST, FREEZECORE, NEWJUMP, RENORM, CAPSID, FREEZE, &
 88:      &        OTPT, LJMFT, STRANDT, PAHT, SWT, MSTRANST, STOCKT, STICKYT, BLNT, MYSDT, FREEZERES, CENTXY, & 88:      &        OTPT, LJMFT, STRANDT, PAHT, SWT, MSTRANST, STOCKT, STICKYT, BLNT, MYSDT, FREEZERES, CENTXY, &
 89:      &        MSORIGT, SQUEEZET, PERIODIC, SCT, MSCT, MGUPTAT, RESIZET, TIP, RIGID, CALCQT, MPIT, JMT, LJCOULT, SETCENT, & 89:      &        MSORIGT, SQUEEZET, PERIODIC, SCT, MSCT, MGUPTAT, RESIZET, TIP, RIGID, CALCQT, MPIT, JMT, LJCOULT, SETCENT, &
 90:      &        SORTT, HIT, SAVEQ, PARALLELT, FIXD, RKMIN, BSMIN, PERMDIST, PERMOPT, BSWL, BSPT, BSPTRESTART, & 90:      &        SORTT, HIT, SAVEQ, PARALLELT, FIXD, RKMIN, BSMIN, PERMDIST, PERMOPT, BSWL, BSPT, BSPTRESTART, &
 91:      &        SYMMETRIZE, SYMMETRIZECSM, PRINT_PTGRP, PRINT_MINDATA, DUMPT, NEON, ARGON, P46, NORESET, TABOOT, EVSTEPT, PACHECO, DL_POLY, QUCENTRE, & 91:      &        SYMMETRIZE, SYMMETRIZECSM, PRINT_PTGRP, PRINT_MINDATA, DUMPT, NEON, ARGON, P46, NORESET, TABOOT, EVSTEPT, PACHECO, DL_POLY, QUCENTRE, &
 92:      &        STAR, PLUS, TWOPLUS, GROUND, DIPOLE, DFTBT, DFTBCT, SW, SUPERSTEP, EAMLJT, PBGLUET, TRACKDATAT, & 92:      &        STAR, PLUS, TWOPLUS, GROUND, DIPOLE, DFTBT, DFTBCT, SW, SUPERSTEP, EAMLJT, PBGLUET, TRACKDATAT, &
 93:      &        EAMALT, ALGLUET, MGGLUET, GUPTAT, LJATT, FST, DECAY, COOP, FIXBIN, GAUSST, QUENCHDOS, FIXDIHEFLAG, & 93:      &        EAMALT, ALGLUET, MGGLUET, GUPTAT, LJATT, FST, DECAY, COOP, FIXBIN, GAUSST, QUENCHDOS, FIXDIHEFLAG, &
 94:      &        FRAUSIT, ANGST, SELFT, STEPOUT, WENZEL, THRESHOLDT, THOMSONT, MULLERBROWNT, CHARMMENERGIES, & 94:      &        FRAUSIT, ANGST, SELFT, STEPOUT, WENZEL, THRESHOLDT, THOMSONT, MULLERBROWNT, CHARMMENERGIES, &
 95:      &        PROJ, RGCL2, TOSI, WELCH, AXTELL, AMBER, FIXIMAGE, BINARY, SHIFTCUT, ARNO, TUNNELT, TWOD, &  95:      &        PROJ, RGCL2, TOSI, WELCH, AXTELL, AMBER, FIXIMAGE, BINARY, SHIFTCUT, ARNO, TUNNELT, TWOD, & 
109:      &        HARMONICDONTMOVE, DUMPUNIQUE, FREEZESAVE, TBP, RBSYMT, PTMCDUMPSTRUCT, PTMCDUMPENERT, PYCOLDFUSION, MONITORT,&109:      &        HARMONICDONTMOVE, DUMPUNIQUE, FREEZESAVE, TBP, RBSYMT, PTMCDUMPSTRUCT, PTMCDUMPENERT, PYCOLDFUSION, MONITORT,&
110:      &        CHARMMDFTBT, PERMINVOPT, BLOCKMOVET, MAXERISE_SET, PYT, BINARY_EXAB, CHIROT, SANDBOXT, &110:      &        CHARMMDFTBT, PERMINVOPT, BLOCKMOVET, MAXERISE_SET, PYT, BINARY_EXAB, CHIROT, SANDBOXT, &
111:      &        RESERVOIRT, DISTOPT, ONEDAPBCT, ONEDPBCT, TWODAPBCT, TWODPBCT, THREEDAPBCT, THREEDPBCT, RATIOT, &111:      &        RESERVOIRT, DISTOPT, ONEDAPBCT, ONEDPBCT, TWODAPBCT, TWODPBCT, THREEDAPBCT, THREEDPBCT, RATIOT, &
112:      &        PTRANDOM, PTINTERVAL, PTSINGLE, PTSETS, CHEMSHIFT, CHEMSHIFT2, CSH, DEBUGss2029, UNIFORMMOVE, RANSEEDT, &112:      &        PTRANDOM, PTINTERVAL, PTSINGLE, PTSETS, CHEMSHIFT, CHEMSHIFT2, CSH, DEBUGss2029, UNIFORMMOVE, RANSEEDT, &
113:      &        TTM3T, NOINVERSION, RIGIDCONTOURT, UPDATERIGIDREFT, HYBRIDMINT, COMPRESSRIGIDT, MWFILMT, &113:      &        TTM3T, NOINVERSION, RIGIDCONTOURT, UPDATERIGIDREFT, HYBRIDMINT, COMPRESSRIGIDT, MWFILMT, &
114:      &        SUPPRESST, MFETT, POLIRT, QUIPT, SWPOTT, MWPOTT, REPMATCHT, GLJT, MLJT, READMASST, SPECMASST, NEWTSALLIST, &114:      &        SUPPRESST, MFETT, POLIRT, QUIPT, SWPOTT, MWPOTT, REPMATCHT, GLJT, MLJT, READMASST, SPECMASST, NEWTSALLIST, &
115:      &        PHI4MODELT, CUDAT, CUDATIMET, AMBER12T, ENERGY_DECOMPT, NEWMOVEST, DUMPMINT, MBPOLT, MOLECULART, GCBHT, SEMIGRAND_MUT, USEROT, & 115:      &        PHI4MODELT, CUDAT, CUDATIMET, AMBER12T, ENERGY_DECOMPT, NEWMOVEST, DUMPMINT, MBPOLT, MOLECULART, GCBHT, SEMIGRAND_MUT, USEROT, & 
116:      &        SAVEMULTIMINONLY, GRADPROBLEMT, INTLJT, CONDATT, QCIPERMCHECK, &116:      &        SAVEMULTIMINONLY, GRADPROBLEMT, INTLJT, CONDATT, QCIPERMCHECK, &
117:      &        INTCONSTRAINTT, INTFREEZET, CHECKCONINT, CONCUTABST, CONCUTFRACT, INTERPCOSTFUNCTION, &117:      &        INTCONSTRAINTT, INTFREEZET, CHECKCONINT, CONCUTABST, CONCUTFRACT, INTERPCOSTFUNCTION, &
118:      &        RBAAT, FREEZENODEST, DUMPINTEOS, DUMPINTXYZ, QCIPOTT, QCIPOT2T, INTSPRINGACTIVET, LPERMDIST, LOCALPERMDIST, QCIRADSHIFTT, &118:      &        RBAAT, FREEZENODEST, DUMPINTEOS, DUMPINTXYZ, QCIPOTT, QCIPOT2T, INTSPRINGACTIVET, LPERMDIST, LOCALPERMDIST, QCIRADSHIFTT, &
119:      &        MLP3T, MKTRAPT, MLPB3T, MULTIPOTT, QCIAMBERT, MLPNEWREG, DJWRBT, STEALTHYT119:      &        MLP3T, MKTRAPT, MLPB3T, MULTIPOTT, QCIAMBERT, MLPNEWREG, DJWRBT
120: !120: !
121:       DOUBLE PRECISION, ALLOCATABLE :: SEMIGRAND_MU(:) 121:       DOUBLE PRECISION, ALLOCATABLE :: SEMIGRAND_MU(:) 
122:       DOUBLE PRECISION, ALLOCATABLE:: ATMASS(:)122:       DOUBLE PRECISION, ALLOCATABLE:: ATMASS(:)
123:       DOUBLE PRECISION, ALLOCATABLE:: SPECMASS(:) 123:       DOUBLE PRECISION, ALLOCATABLE:: SPECMASS(:) 
124: 124: 
125: ! csw34> FREEZEGROUP variables125: ! csw34> FREEZEGROUP variables
126: !126: !
127:       INTEGER :: GROUPCENTRE127:       INTEGER :: GROUPCENTRE
128:       DOUBLE PRECISION :: GROUPRADIUS128:       DOUBLE PRECISION :: GROUPRADIUS
129:       CHARACTER (LEN=2) :: FREEZEGROUPTYPE129:       CHARACTER (LEN=2) :: FREEZEGROUPTYPE


r30265/keywords.f 2016-04-14 11:30:12.017632776 +0100 r30264/keywords.f 2016-04-14 11:30:12.933645345 +0100
5743: !           CALL READF(XX)5743: !           CALL READF(XX)
5744: !           SQUEEZER=XX5744: !           SQUEEZER=XX
5745: !        ENDIF5745: !        ENDIF
5746: !        IF (NITEMS.GT.3) THEN5746: !        IF (NITEMS.GT.3) THEN
5747: !           CALL READF(XX)5747: !           CALL READF(XX)
5748: !           SQUEEZED=XX5748: !           SQUEEZED=XX
5749: !        ENDIF5749: !        ENDIF
5750: 5750: 
5751:       ELSE IF (WORD.EQ.'STAR') THEN5751:       ELSE IF (WORD.EQ.'STAR') THEN
5752:          STAR=.TRUE.5752:          STAR=.TRUE.
5753:  
5754:       ELSE IF (WORD.EQ.'STEALTHY') THEN 
5755:          STEALTHYT=.TRUE. 
5756:          CALL READF(KLIM) 
5757:          IF (NITEMS.GT.2) THEN 
5758:             CALL READF(SCA) 
5759:          ELSE 
5760:             SCA=1 
5761:          END IF 
5762: !5753: !
5763: ! Read in the maximum initial step size, factor for determining angular5754: ! Read in the maximum initial step size, factor for determining angular
5764: ! moves, and for rigid bodies the angular step size and the size of the5755: ! moves, and for rigid bodies the angular step size and the size of the
5765: ! blocks for Cartesian and angular moves.5756: ! blocks for Cartesian and angular moves.
5766: !5757: !
5767: ! For parallel runs different values can be used for different runs by5758: ! For parallel runs different values can be used for different runs by
5768: ! adding additional "STEP" lines to the data file. Otherwise the5759: ! adding additional "STEP" lines to the data file. Otherwise the
5769: ! parameters for subsequent parallel runs are set to the values for the5760: ! parameters for subsequent parallel runs are set to the values for the
5770: ! first one.5761: ! first one.
5771: !5762: !


r30265/potential.f90 2016-04-14 11:30:12.353637383 +0100 r30264/potential.f90 2016-04-14 11:30:13.121647930 +0100
1123: 1123: 
1124:    ELSE IF (JMT) THEN1124:    ELSE IF (JMT) THEN
1125:       CALL JMEC(NATOMS, X, GRAD, EREAL, GRADT, SECT)1125:       CALL JMEC(NATOMS, X, GRAD, EREAL, GRADT, SECT)
1126: 1126: 
1127:    ELSE IF (DF1T) THEN1127:    ELSE IF (DF1T) THEN
1128:       CALL DF1GRAD(X, NATOMS, GRAD, EREAL, GRADT, BOXLX, BOXLY)1128:       CALL DF1GRAD(X, NATOMS, GRAD, EREAL, GRADT, BOXLX, BOXLY)
1129: 1129: 
1130:    ELSE IF (ONEDAPBCT) THEN1130:    ELSE IF (ONEDAPBCT) THEN
1131:       CALL ENERGY_1D_APBC(X, GRAD, EREAL, GRADT, SECT)1131:       CALL ENERGY_1D_APBC(X, GRAD, EREAL, GRADT, SECT)
1132: 1132: 
1133:    ELSE IF (STEALTHYT) THEN 
1134:       CALL STEALTHY(X, GRAD, EREAL, GRADT, SECT) 
1135:  
1136:    ELSE IF (ONEDPBCT) THEN1133:    ELSE IF (ONEDPBCT) THEN
1137:       WRITE(*, *) 'This potential option is commented out.'1134:       WRITE(*, *) 'This potential option is commented out.'
1138:       STOP1135:       STOP
1139: !      CALL ENERGY_1D_PBC(X, GRAD, EREAL, GRADT, SECT)1136: !      CALL ENERGY_1D_PBC(X, GRAD, EREAL, GRADT, SECT)
1140:    ELSE IF (TWODAPBCT) THEN1137:    ELSE IF (TWODAPBCT) THEN
1141:       WRITE(*, *) 'This potential option is commented out.'1138:       WRITE(*, *) 'This potential option is commented out.'
1142:       STOP1139:       STOP
1143: !      CALL ENERGY_2D_APBC(X, GRAD, EREAL, GRADT, SECT)            1140: !      CALL ENERGY_2D_APBC(X, GRAD, EREAL, GRADT, SECT)            
1144:    ELSE IF (TWODPBCT) THEN1141:    ELSE IF (TWODPBCT) THEN
1145:       WRITE(*, *) 'This potential option is commented out.'1142:       WRITE(*, *) 'This potential option is commented out.'


r30265/stealthy.f90 2016-04-14 11:30:12.545640020 +0100 r30264/stealthy.f90 2016-04-14 11:30:13.313650568 +0100
  1: SUBROUTINE STEALTHY(XS,TIDU,ENERGY,GTEST,STEST)  1: svn: E195012: Unable to find repository location for 'svn+ssh://svn.ch.private.cam.ac.uk/groups/wales/trunk/GMIN/source/stealthy.f90' in revision 30264
  2:     USE COMMONS 
  3:     USE MODHESS 
  4:     IMPLICIT NONE 
  5:     LOGICAL GTEST,STEST 
  6:     INTEGER M, M1, M2, M3, J1, J2, I1, I2, L1, L2, KI!, KNUM, GNUM 
  7:     DOUBLE PRECISION XS(3*NATOMS), KX(3), KY(3), KZ(3), KN(3), TIDU(3*NATOMS), & 
  8:                      CKR(NATOMS), SKR(NATOMS) 
  9:     DOUBLE PRECISION PI, KR ,C, S, ENERGY, VF, IM, KM 
 10:  
 11: !    IF (ALLOCATED(HESS)) DEALLOCATE(HESS) 
 12: !    ALLOCATE(HESS(3*NATOMS, 3*NATOMS)) 
 13:  
 14: !   initialize 
 15:     PI=3.141592653589793 
 16:     KM=KLIM*KLIM 
 17:     KI=DINT(KLIM) 
 18:     ENERGY=0 
 19:  
 20:     IF (GTEST) THEN 
 21:        TIDU=0 
 22:     END IF 
 23:  
 24:     IF (STEST) THEN 
 25:        HESS=0 
 26:     END IF 
 27:  
 28: !    KNUM=0 
 29: !    GNUM=0 
 30:  
 31: !   get the initial box 
 32:  
 33: !   calculate the wave vector 
 34:     KX(1)=2*PI/BOXLX 
 35:     KX(2)=0 
 36:     KX(3)=0 
 37:     KY(1)=0 
 38:     KY(2)=2*PI/BOXLY 
 39:     KY(3)=0 
 40:     KZ(1)=0 
 41:     KZ(2)=0 
 42:     KZ(3)=2*PI/BOXLZ 
 43:     VF=BOXLX*BOXLY*BOXLZ 
 44: !   WRITE(*,*) KI(1), KI(2), KI(3) 
 45:  
 46: !   periodic boundary condition 
 47:     IF (PERIODIC) THEN 
 48:        DO J1=1, NATOMS 
 49:           J2=3*(J1-1) 
 50:           XS(J2+1)=XS(J2+1) - BOXLX*DNINT(XS(J2+1)/BOXLX) 
 51:           XS(J2+2)=XS(J2+2) - BOXLY*DNINT(XS(J2+2)/BOXLY) 
 52:           XS(J2+3)=XS(J2+3) - BOXLZ*DNINT(XS(J2+3)/BOXLZ) 
 53:        END DO 
 54:     END IF 
 55:  
 56: !   sum of k 
 57:     DO M1=0, KI 
 58:        DO M2=-KI, KI 
 59:           DO M3=-KI, KI 
 60:  
 61:              IF (M1==0.AND.M2>=0.AND.M3<0) CYCLE 
 62:              IF (M1==0.AND.M2<=0.AND.M3<=0) CYCLE 
 63:  
 64:              M=M1*M1+M2*M2+M3*M3 
 65:              IF (M>KM.OR.M==0) CYCLE 
 66:  
 67: !             KNUM=KNUM+1 
 68:              C=0 
 69:              S=0 
 70:              KN(1)=KX(1)*M1 + KY(1)*M2 + KZ(1)*M3 
 71:              KN(2)=KX(2)*M1 + KY(2)*M2 + KZ(2)*M3 
 72:              KN(3)=KX(3)*M1 + KY(3)*M2 + KZ(3)*M3 
 73:  
 74: !            calculation of Re[n(k)] and Im[n(k)] 
 75:              DO I1=1, NATOMS 
 76:                 I2=3*(I1-1) 
 77:                 KR=KN(1)*XS(I2+1) + KN(2)*XS(I2+2) + KN(3)*XS(I2+3) 
 78:                 CKR(I1)=COS(KR) 
 79:                 SKR(I1)=SIN(KR) 
 80:                 C=C+CKR(I1) 
 81:                 S=S+SKR(I1) 
 82:              END DO 
 83:  
 84: !            E=n(k)^2/Vf 
 85:              ENERGY=ENERGY+ SCA*(C*C+S*S)/VF 
 86:  
 87: !            gradient 
 88:              IF (GTEST) THEN 
 89: !                GNUM=GNUM+1 
 90:                 DO L1=1, NATOMS 
 91:                    L2=3*(L1-1) 
 92:                    IM=( C*SKR(L1) - S*CKR(L1) )*SCA 
 93:                    TIDU(L2+1)=TIDU(L2+1) - 2*KN(1)*IM/VF 
 94:                    TIDU(L2+2)=TIDU(L2+2) - 2*KN(2)*IM/VF 
 95:                    TIDU(L2+3)=TIDU(L2+3) - 2*KN(3)*IM/VF 
 96:                 END DO 
 97:              END IF 
 98:  
 99: !            hessian matrix 
100:              IF (STEST) THEN 
101:                 CALL STHESS(XS, KN, VF, CKR, C, SKR, S) 
102:              END IF 
103:  
104:           END DO 
105:        END DO 
106:     END DO 
107:  
108: !    WRITE(*,*) KNUM 
109: !    WRITE(*,*) GNUM 
110: !   change the box back to cubic 
111:  
112: END SUBROUTINE 
113:  
114: SUBROUTINE STHESS(X, K, V, CK, CSUM, SK, SSUM) 
115:    USE COMMONS 
116:    USE MODHESS 
117:    IMPLICIT NONE 
118:    INTEGER I1, I2, J1, J2 
119:    DOUBLE PRECISION X(3*NATOMS), K(3), RIJ(3), CK(NATOMS), SK(NATOMS) 
120:    DOUBLE PRECISION KRIJ, V, CSUM, SSUM, CI, SI, CF, CSKR 
121:  
122:    DO I1=1, NATOMS 
123:       I2=3*(I1-1) 
124:  
125:       CI=CSUM-CK(I1) 
126:       SI=SSUM-SK(I1) 
127:       CF=( CK(I1)*CI + SK(I1)*SI )*SCA 
128:  
129:       HESS(I2+1, I2+1)=HESS(I2+1, I2+1) + 2*K(1)*K(1)/V * CF 
130:       HESS(I2+1, I2+2)=HESS(I2+1, I2+2) + 2*K(1)*K(2)/V * CF 
131:       HESS(I2+1, I2+3)=HESS(I2+1, I2+3) + 2*K(1)*K(3)/V * CF 
132:       HESS(I2+2, I2+2)=HESS(I2+2, I2+2) + 2*K(2)*K(2)/V * CF 
133:       HESS(I2+2, I2+3)=HESS(I2+2, I2+3) + 2*K(2)*K(3)/V * CF 
134:       HESS(I2+3, I2+3)=HESS(I2+3, I2+3) + 2*K(3)*K(3)/V * CF 
135:  
136:       DO J1=I1+1, NATOMS 
137:          J2=3*(J1-1) 
138:  
139:          RIJ(1)=X(I2+1) - X(J2+1) 
140:          RIJ(2)=X(I2+2) - X(J2+2) 
141:          RIJ(3)=X(I2+3) - X(J2+3) 
142:  
143:          KRIJ=K(1)*RIJ(1) + K(2)*RIJ(2) + K(3)*RIJ(3) 
144:          CSKR= COS(KRIJ)*SCA 
145:  
146:          HESS(I2+1, J2+1)=HESS(I2+1, J2+1) + 2*K(1)*K(1)/V * CSKR 
147:          HESS(I2+1, J2+2)=HESS(I2+1, J2+2) + 2*K(1)*K(2)/V * CSKR 
148:          HESS(I2+1, J2+3)=HESS(I2+1, J2+3) + 2*K(1)*K(3)/V * CSKR 
149:  
150:          HESS(I2+2, J2+1)=HESS(I2+2, J2+1) + 2*K(2)*K(1)/V * CSKR 
151:          HESS(I2+2, J2+2)=HESS(I2+2, J2+2) + 2*K(2)*K(2)/V * CSKR 
152:          HESS(I2+2, J2+3)=HESS(I2+2, J2+3) + 2*K(2)*K(3)/V * CSKR 
153:  
154:          HESS(I2+3, J2+1)=HESS(I2+3, J2+1) + 2*K(3)*K(1)/V * CSKR 
155:          HESS(I2+3, J2+2)=HESS(I2+3, J2+2) + 2*K(3)*K(2)/V * CSKR 
156:          HESS(I2+3, J2+3)=HESS(I2+3, J2+3) + 2*K(3)*K(3)/V * CSKR 
157:  
158:       END DO 
159:    END DO 
160:  
161:    DO I1=1, NATOMS 
162:       I2=3*(I1-1) 
163:  
164:       HESS(I2+2, I2+1)=HESS(I2+1, I2+2) 
165:       HESS(I2+3, I2+1)=HESS(I2+1, I2+3) 
166:       HESS(I2+3, I2+2)=HESS(I2+2, I2+3) 
167:  
168:       DO J1=I1+1, NATOMS 
169:          J2=3*(J1-1) 
170:  
171:          HESS(J2+1, I2+1)=HESS(I2+1, J2+1) 
172:          HESS(J2+2, I2+1)=HESS(I2+1, J2+2) 
173:          HESS(J2+3, I2+1)=HESS(I2+1, J2+3) 
174:  
175:          HESS(J2+1, I2+2)=HESS(I2+2, J2+1) 
176:          HESS(J2+2, I2+2)=HESS(I2+2, J2+2) 
177:          HESS(J2+3, I2+2)=HESS(I2+2, J2+3) 
178:  
179:          HESS(J2+1, I2+3)=HESS(I2+3, J2+1) 
180:          HESS(J2+2, I2+3)=HESS(I2+3, J2+2) 
181:          HESS(J2+3, I2+3)=HESS(I2+3, J2+3) 
182:  
183:       END DO 
184:    END DO 
185:  
186: END SUBROUTINE 


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0