hdiff output

r31837/chiro.f90 2017-01-27 15:30:11.961009278 +0000 r31836/chiro.f90 2017-01-27 15:30:12.189012382 +0000
327:     double precision :: xi(3), xj(3)327:     double precision :: xi(3), xj(3)
328: 328: 
329:     double precision :: image_moljr(3)329:     double precision :: image_moljr(3)
330:     integer :: lshift, ushift, jshift330:     integer :: lshift, ushift, jshift
331:     logical :: newpair331:     logical :: newpair
332: 332: 
333: ! hk286333: ! hk286
334:     double precision :: epsilon_penalty, delta_penalty334:     double precision :: epsilon_penalty, delta_penalty
335:     integer :: power_penalty335:     integer :: power_penalty
336: 336: 
337:  
338:     epsilon_penalty = 1.0d20337:     epsilon_penalty = 1.0d20
339:     delta_penalty = 1.0d-3338:     delta_penalty = 1.0d-3
340:     power_penalty = 100339:     power_penalty = 100
341: 340: 
342:     energy = 0.0341:     energy = 0.0
343:     vt(:) = 0.0342:     vt(:) = 0.0
344:     if(gtest) grad(:) = 0.0343:     if(gtest) grad(:) = 0.0
345: 344: 
346:     ! update the values for all the molecules345:     ! update the values for all the molecules
347:     do moli_ind = 1, nmols346:     do moli_ind = 1, nmols
491:                 grad(moli%rind: moli%rind + 2) = grad(moli%rind: moli%rind + 2) + grad_contrib(1: 3)490:                 grad(moli%rind: moli%rind + 2) = grad(moli%rind: moli%rind + 2) + grad_contrib(1: 3)
492:                 grad(molj%rind: molj%rind + 2) = grad(molj%rind: molj%rind + 2) + grad_contrib(4: 6)491:                 grad(molj%rind: molj%rind + 2) = grad(molj%rind: molj%rind + 2) + grad_contrib(4: 6)
493:                 grad(moli%pind: moli%pind + 2) = grad(moli%pind: moli%pind + 2) + grad_contrib(7: 9)492:                 grad(moli%pind: moli%pind + 2) = grad(moli%pind: moli%pind + 2) + grad_contrib(7: 9)
494:                 grad(molj%pind: molj%pind + 2) = grad(molj%pind: molj%pind + 2) + grad_contrib(10: 12)493:                 grad(molj%pind: molj%pind + 2) = grad(molj%pind: molj%pind + 2) + grad_contrib(10: 12)
495:             end if494:             end if
496: 495: 
497:             ! increment pbc shift for inner loop496:             ! increment pbc shift for inner loop
498:             end do497:             end do
499:         end do498:         end do
500: 499: 
501:     end do 500:     end do
502:       
503:  ! 2D: set all z gradients to zero501:  ! 2D: set all z gradients to zero
504:     if(twod) then502:     if(twod) then
505:         do moli_ind = 1, nmols503:         do moli_ind = 1, nmols
506:             grad(molecules(moli_ind)%rind + 2) = 0.0504:             grad(molecules(moli_ind)%rind + 2) = 0.0
507:         end do505:         end do
508:     end if506:     end if
509: 507: 
510: end subroutine chiro508: end subroutine chiro
511: 509: 
512: subroutine chiro_output510: subroutine chiro_output


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0