hdiff output

r29867/MLPB3.f90 2016-03-16 18:32:39.070500492 +0000 r29866/MLPB3.f90 2016-03-16 18:32:39.274502590 +0000
182:       ENDDO182:       ENDDO
183: 183: 
184:       D2PCDBHID2=DPCDBHID(MLPOUTJ1)*(DYCDBHID(MLPOUTJ1)-DUMMY3) + PMLPOUTJ1*(D2YCDBHID2(MLPOUTJ1)-DUMMY1)184:       D2PCDBHID2=DPCDBHID(MLPOUTJ1)*(DYCDBHID(MLPOUTJ1)-DUMMY3) + PMLPOUTJ1*(D2YCDBHID2(MLPOUTJ1)-DUMMY1)
185: 185: 
186:       HESS(BOFFSET+1,BOFFSET+1)=HESS(BOFFSET+1,BOFFSET+1)+DPCDBHID(MLPOUTJ1)**2/PMLPOUTJ1**2-D2PCDBHID2/PMLPOUTJ1186:       HESS(BOFFSET+1,BOFFSET+1)=HESS(BOFFSET+1,BOFFSET+1)+DPCDBHID(MLPOUTJ1)**2/PMLPOUTJ1**2-D2PCDBHID2/PMLPOUTJ1
187: !187: !
188: !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!188: !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
189: !189: !
190: !     b hidden with w^2 beta gamma190: !     b hidden with w^2 beta gamma
191: !191: !
 192:       DUMMY3=0.0D0
 193:       DO J4=1,MLPOUT
 194:          DUMMY3=DUMMY3+PROB(J4)*DYCDBHID(J4)
 195:       ENDDO
 196:       DO J4=1,MLPOUT
 197:          DPCDBHID(J4)=PROB(J4)*(DYCDBHID(J4)-DUMMY2)
 198:       ENDDO
 199: 
192:       DUMMY1=0.0D0200:       DUMMY1=0.0D0
193:       DO J4=1,MLPOUT201:       DO J4=1,MLPOUT
194:          DUMMY1=DUMMY1+PROB(J4)*DYCDBHID(J4)202:          DUMMY1=DUMMY1+PROB(J4)*DYCDBHID(J4)
195:       ENDDO203:       ENDDO
196: 204: 
197:       DO J2=1,MLPHIDDEN ! J2 is beta205:       DO J2=1,MLPHIDDEN ! J2 is beta
198:          DO J3=1,MLPIN  ! J3 is gamma206:          DO J3=1,MLPIN  ! J3 is gamma
199:             DUMMY4=0.0D0207:             DUMMY4=0.0D0
200:             DO J4=1,MLPOUT208:             DO J4=1,MLPOUT
201:                DUMMY4=DUMMY4+PROB(J4)*(X(MLPOFFSET+(J4-1)*MLPHIDDEN+J2)+X(BOFFSET+2))209:                DUMMY4=DUMMY4+PROB(J4)*(X(MLPOFFSET+(J4-1)*MLPHIDDEN+J2)+X(BOFFSET+2))
203:             DO J4=1,MLPOUT211:             DO J4=1,MLPOUT
204:                DPCDW2BG(J4,J2,J3)=PROB(J4)*MLPDAT(J1,J3)*SECH2(J2)*(X(MLPOFFSET+(J4-1)*MLPHIDDEN+J2)+X(BOFFSET+2)-DUMMY4)212:                DPCDW2BG(J4,J2,J3)=PROB(J4)*MLPDAT(J1,J3)*SECH2(J2)*(X(MLPOFFSET+(J4-1)*MLPHIDDEN+J2)+X(BOFFSET+2)-DUMMY4)
205:             ENDDO213:             ENDDO
206:             DUMMY2=0.0D0214:             DUMMY2=0.0D0
207:             DUMMY3=0.0D0215:             DUMMY3=0.0D0
208:             DO J4=1,MLPOUT216:             DO J4=1,MLPOUT
209:                DUMMY2=DUMMY2+DYCDBHID(J4)*DPCDW2BG(J4,J2,J3)217:                DUMMY2=DUMMY2+DYCDBHID(J4)*DPCDW2BG(J4,J2,J3)
210:                DUMMY3=DUMMY3-2.0D0*PROB(J4)*MLPDAT(J1,J3)*SUMINPUTS*(X(MLPOFFSET+(J4-1)*MLPHIDDEN+J2)+X(BOFFSET+2))*SECH2(J2)*TANHSUM(J2)218:                DUMMY3=DUMMY3-2.0D0*PROB(J4)*MLPDAT(J1,J3)*SUMINPUTS*(X(MLPOFFSET+(J4-1)*MLPHIDDEN+J2)+X(BOFFSET+2))*SECH2(J2)*TANHSUM(J2)
211:             ENDDO219:             ENDDO
212: 220: 
213:             D2YCDBHIDDW2BG=-2.0D0*MLPDAT(J1,J3)*SUMINPUTS*(X(MLPOFFSET+(MLPOUTJ1-1)*MLPHIDDEN+J2)+X(BOFFSET+2))*SECH2(J2)*TANHSUM(J2)221:             D2YCDBHIDDW2BG=-2.0D0*PMLPOUTJ1*MLPDAT(J1,J3)*SUMINPUTS*(X(MLPOFFSET+(MLPOUTJ1-1)*MLPHIDDEN+J2)+X(BOFFSET+2))*SECH2(J2)*TANHSUM(J2) &
214: 222:   &         + DPCDW2BG(MLPOUTJ1,J2,J3)*(DYCDBHID(MLPOUTJ1)-DUMMY1)-PMLPOUTJ1*(DUMMY2+DUMMY3) 
215:             D2PCDBHIDDW2BG=DPCDW2BG(MLPOUTJ1,J2,J3)*(DYCDBHID(MLPOUTJ1)-DUMMY1)+PMLPOUTJ1*(D2YCDBHIDDW2BG-DUMMY2-DUMMY3)223:             HESS(BOFFSET+1,(J3-1)*MLPHIDDEN+J2)= &
216: 224:   &         HESS(BOFFSET+1,(J3-1)*MLPHIDDEN+J2) + DPCDBHID(MLPOUTJ1)*DPCDW2BG(MLPOUTJ1,J2,J3)/PMLPOUTJ1**2 - D2PCDBHIDDW2BG/PMLPOUTJ1
217:             HESS(BOFFSET+1,(J2-1)*MLPIN+J3)= & 
218:   &         HESS(BOFFSET+1,(J2-1)*MLPIN+J3) + DPCDBHID(MLPOUTJ1)*DPCDW2BG(MLPOUTJ1,J2,J3)/PMLPOUTJ1**2 - D2PCDBHIDDW2BG/PMLPOUTJ1 
219:          ENDDO225:          ENDDO
220:       ENDDO226:       ENDDO
221: !227: !
222: ! This block w^1 with w^1 is locally symmetric228: ! This block w^1 with w^1 is locally symmetric
223: !229: !
224:       DO J4=1,MLPOUT ! J4 is beta 230:       DO J4=1,MLPOUT ! J4 is beta 
225:          DO J2=1,MLPHIDDEN ! J2 is gamma231:          DO J2=1,MLPHIDDEN ! J2 is gamma
226:             DO K4=1,J4 ! K4 is alpha232:             DO K4=1,J4 ! K4 is alpha
227:                DO K2=1,MLPHIDDEN ! K2 is epsilon233:                DO K2=1,MLPHIDDEN ! K2 is epsilon
228:                   DUMMY1=0.0D0234:                   DUMMY1=0.0D0


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0