ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/group/trunk/OOPSE/libmdtools/do_Forces.F90
(Generate patch)

Comparing trunk/OOPSE/libmdtools/do_Forces.F90 (file contents):
Revision 648 by chuckv, Wed Jul 23 22:13:59 2003 UTC vs.
Revision 841 by mmeineke, Wed Oct 29 17:55:28 2003 UTC

# Line 4 | Line 4
4  
5   !! @author Charles F. Vardeman II
6   !! @author Matthew Meineke
7 < !! @version $Id: do_Forces.F90,v 1.24 2003-07-23 22:13:59 chuckv Exp $, $Date: 2003-07-23 22:13:59 $, $Name: not supported by cvs2svn $, $Revision: 1.24 $
7 > !! @version $Id: do_Forces.F90,v 1.34 2003-10-29 17:55:28 mmeineke Exp $, $Date: 2003-10-29 17:55:28 $, $Name: not supported by cvs2svn $, $Revision: 1.34 $
8  
9   module do_Forces
10    use force_globals
# Line 18 | Line 18 | module do_Forces
18    use reaction_field
19    use gb_pair
20    use vector_class
21 +  use eam
22 +  use status
23   #ifdef IS_MPI
24    use mpiSimulation
25   #endif
# Line 43 | Line 45 | contains
45    public :: do_force_loop
46    public :: setRlistDF
47  
48 + #ifdef PROFILE
49 +  real(kind = dp) :: forceTime
50 +  real(kind = dp) :: forceTimeInitial, forceTimeFinal
51 +  real(kind = dp) :: globalForceTime
52 +  real(kind = dp) :: maxForceTime
53 +  integer, save :: nloops = 0
54 + #endif
55 +
56   contains
57  
58    subroutine setRlistDF( this_rlist )
# Line 141 | Line 151 | contains
151            return
152         end if
153      endif
154 +
155 +
156 +    if (FF_uses_EAM) then
157 +         call init_EAM_FF(my_status)
158 +       if (my_status /= 0) then
159 +          write(*,*) "init_EAM_FF returned a bad status"
160 +          thisStat = -1
161 +          return
162 +       end if
163 +    endif
164 +
165 +
166      
167      if (FF_uses_GB) then
168         call check_gb_pair_FF(my_status)
# Line 161 | Line 183 | contains
183            return
184         endif
185      endif
186 +    
187  
188      havePolicies = .true.
189      if( haveRlist ) do_forces_initialized = .true.
190 <    
190 >
191    end subroutine init_FF
192    
193  
# Line 192 | Line 215 | contains
215      real( kind = DP ) :: pot_local
216      integer :: nrow
217      integer :: ncol
218 +    integer :: nprocs
219   #endif
220      integer :: nlocal
221      integer :: natoms    
# Line 221 | Line 245 | contains
245      nlocal = getNlocal()
246      natoms = nlocal
247   #endif
248 <  
248 >
249      call check_initialization(localError)
250      if ( localError .ne. 0 ) then
251 +       call handleError("do_force_loop","Not Initialized")
252         error = -1
253         return
254      end if
# Line 232 | Line 257 | contains
257      do_pot = do_pot_c
258      do_stress = do_stress_c
259  
260 +
261      ! Gather all information needed by all force loops:
262      
263   #ifdef IS_MPI    
# Line 248 | Line 274 | contains
274      endif
275      
276   #endif
277 <    
277 >
278 > !! Begin force loop timing:
279 > #ifdef PROFILE
280 >    call cpu_time(forceTimeInitial)
281 >    nloops = nloops + 1
282 > #endif
283 >  
284      if (FF_RequiresPrepairCalc() .and. SimRequiresPrepairCalc()) then
285         !! See if we need to update neighbor lists
286         call checkNeighborList(nlocal, q, listSkin, update_nlist)  
# Line 256 | Line 288 | contains
288         !! do_prepair_loop_if_needed
289         !! if_mpi_scatter_stuff_from_prepair
290         !! if_mpi_gather_stuff_from_prepair_to_main_loop
291 <
291 >    
292   !--------------------PREFORCE LOOP----------->>>>>>>>>>>>>>>>>>>>>>>>>>>
293   #ifdef IS_MPI
294      
# Line 332 | Line 364 | contains
364         neighborListSize = size(list)
365    
366         nlist = 0
367 <      
367 >
368         do i = 1, natoms-1
369            point(i) = nlist + 1
370            
# Line 344 | Line 376 | contains
376            
377  
378               if (rijsq < rlistsq) then
379 <                
379 >
380 >          
381                  nlist = nlist + 1
382                
383                  if (nlist > neighborListSize) then
# Line 369 | Line 402 | contains
402         point(natoms) = nlist + 1
403        
404      else !! (update)
405 <      
405 >  
406         ! use the list to find the neighbors
407         do i = 1, natoms-1
408            JBEG = POINT(i)
# Line 390 | Line 423 | contains
423      endif    
424   #endif
425      !! Do rest of preforce calculations
426 <   call do_preforce(nlocal,pot)
426 >    !! do necessary preforce calculations  
427 >    call do_preforce(nlocal,pot)
428 >   ! we have already updated the neighbor list set it to false...
429 >   update_nlist = .false.
430      else
431         !! See if we need to update neighbor lists for non pre-pair
432         call checkNeighborList(nlocal, q, listSkin, update_nlist)  
# Line 543 | Line 579 | contains
579   #endif
580      
581      ! phew, done with main loop.
582 <    
582 >
583 > !! Do timing
584 > #ifdef PROFILE
585 >    call cpu_time(forceTimeFinal)
586 >    forceTime = forceTime + forceTimeFinal - forceTimeInitial
587 > #endif
588 >
589 >
590   #ifdef IS_MPI
591      !!distribute forces
592    
# Line 657 | Line 700 | contains
700      endif
701  
702   #endif
703 <    
703 >
704 > #ifdef PROFILE
705 >    if (do_pot) then
706 >
707 > #ifdef IS_MPI
708 >
709 >      
710 >       call printCommTime()
711 >
712 >       call mpi_allreduce(forceTime,globalForceTime,1,MPI_DOUBLE_PRECISION, &
713 >            mpi_sum,mpi_comm_world,mpi_err)
714 >
715 >       call mpi_allreduce(forceTime,maxForceTime,1,MPI_DOUBLE_PRECISION, &
716 >            MPI_MAX,mpi_comm_world,mpi_err)
717 >      
718 >       call mpi_comm_size( MPI_COMM_WORLD, nprocs,mpi_err)
719 >      
720 >       if (getMyNode() == 0) then
721 >          write(*,*) "Total processor time spent in force calculations is: ", globalForceTime
722 >          write(*,*) "Total Time spent in force loop per processor is: ", globalforceTime/nprocs
723 >          write(*,*) "Maximum force time on any processor is: ", maxForceTime
724 >       end if
725 > #else
726 >       write(*,*) "Time spent in force loop is: ", forceTime
727 > #endif
728 >
729 >    
730 >    endif
731 >
732 > #endif
733 >
734 >
735 >    write(*,*) "In F, after: Tau = "
736 >    write(*,*) tau(1), tau(2), tau(3)
737 >    write(*,*) tau(4), tau(5), tau(6)
738 >    write(*,*) tau(7), tau(8), tau(9)
739 >    write(*,*) " "
740 >
741 >
742    end subroutine do_force_loop
743  
744    subroutine do_pair(i, j, rijsq, d, do_pot, do_stress, u_l, A, f, t, pot)
# Line 710 | Line 791 | contains
791         call getElementProperty(atypes, me_j, "is_DP", is_DP_j)
792        
793         if ( is_DP_i .and. is_DP_j ) then
713 –          
794            call do_dipole_pair(i, j, d, r, rijsq, pot, u_l, f, t, &
795                 do_pot, do_stress)
796            if (FF_uses_RF .and. SimUsesRF()) then
# Line 735 | Line 815 | contains
815  
816      if (FF_uses_GB .and. SimUsesGB()) then
817  
818 +
819         call getElementProperty(atypes, me_i, "is_GB", is_GB_i)
820         call getElementProperty(atypes, me_j, "is_GB", is_GB_j)
821        
# Line 780 | Line 861 | contains
861    
862     r = sqrt(rijsq)
863    
864 +
865   #ifdef IS_MPI
866     if (tagRow(i) .eq. tagColumn(j)) then
867        write(0,*) 'do_pair is doing', i , j, tagRow(i), tagColumn(j)
# Line 794 | Line 876 | contains
876     me_j = atid(j)
877    
878   #endif
879 <  
879 >    
880     if (FF_uses_EAM .and. SimUsesEAM()) then
881        call getElementProperty(atypes, me_i, "is_EAM", is_EAM_i)
882        call getElementProperty(atypes, me_j, "is_EAM", is_EAM_j)
# Line 802 | Line 884 | contains
884        if ( is_EAM_i .and. is_EAM_j ) &
885             call calc_EAM_prepair_rho(i, j, d, r, rijsq )
886     endif
805 –  end subroutine do_prepair
887  
888 + end subroutine do_prepair
889  
890  
891  
892 +
893    subroutine do_preforce(nlocal,pot)
894      integer :: nlocal
895      real( kind = dp ) :: pot
896  
897 <   if (FF_uses_EAM .and. SimUsesEAM()) then
898 <      call calc_EAM_preforce_Frho(nlocal,pot)
899 <   endif
897 >    if (FF_uses_EAM .and. SimUsesEAM()) then
898 >       call calc_EAM_preforce_Frho(nlocal,pot)
899 >    endif
900  
901  
902    end subroutine do_preforce
# Line 876 | Line 959 | contains
959      error = 0
960      ! Make sure we are properly initialized.
961      if (.not. do_forces_initialized) then
962 +       write(*,*) "Forces not initialized"
963         error = -1
964         return
965      endif
# Line 923 | Line 1007 | contains
1007  
1008   #endif
1009  
1010 +
1011 +    if (FF_uses_EAM .and. SimUsesEAM()) then
1012 +       call clean_EAM()
1013 +    endif
1014 +
1015 +
1016 +
1017 +
1018 +
1019      rf = 0.0_dp
1020      tau_Temp = 0.0_dp
1021      virial_Temp = 0.0_dp
# Line 1035 | Line 1128 | end module do_Forces
1128      doesit = FF_uses_RF
1129    end function FF_RequiresPostpairCalc
1130    
1131 + !! This cleans componets of force arrays belonging only to fortran
1132 +
1133   end module do_Forces

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines