| 58 |
|
#define __FORTRAN90 |
| 59 |
|
#include "brains/fSimulation.h" |
| 60 |
|
#include "UseTheForce/fSwitchingFunction.h" |
| 61 |
+ |
#include "UseTheForce/DarkSide/fElectrostaticSummationMethod.h" |
| 62 |
|
|
| 63 |
|
type (simtype), public, save :: thisSim |
| 64 |
|
|
| 566 |
|
doesit = thisSim%SIM_uses_RF |
| 567 |
|
end function SimRequiresPostpairCalc |
| 568 |
|
|
| 569 |
< |
! Function returns true if the simulation has this atype |
| 569 |
> |
! Function returns true if the simulation has this atype |
| 570 |
|
function SimHasAtype(thisAtype) result(doesit) |
| 571 |
|
logical :: doesit |
| 572 |
|
integer :: thisAtype |
| 595 |
|
end if |
| 596 |
|
SimHasAtypeMap = .false. |
| 597 |
|
end if |
| 598 |
< |
! Loop through the local atoms and grab the atypes present |
| 598 |
> |
! Loop through the local atoms and grab the atypes present |
| 599 |
|
do me_i = 1,nLocal |
| 600 |
|
SimHasAtypeMap(atid(me_i)) = .true. |
| 601 |
|
end do |
| 602 |
< |
! For MPI, we need to know all possible atypes present in simulation on all |
| 603 |
< |
! processors. Use LOR operation to set map. |
| 602 |
> |
! For MPI, we need to know all possible atypes present in |
| 603 |
> |
! simulation on all processors. Use LOR operation to set map. |
| 604 |
|
#ifdef IS_MPI |
| 605 |
< |
call mpi_allreduce(SimHasAtypeMap, SimHasAtypeMap, nAtypes, mpi_logical, & |
| 606 |
< |
MPI_LOR, mpi_comm_world, mpiErrors) |
| 607 |
< |
#endif |
| 607 |
< |
|
| 605 |
> |
call mpi_allreduce(SimHasAtypeMap, SimHasAtypeMap, nAtypes, & |
| 606 |
> |
mpi_logical, MPI_LOR, mpi_comm_world, mpiErrors) |
| 607 |
> |
#endif |
| 608 |
|
end subroutine createSimHasAtype |
| 609 |
< |
|
| 609 |
> |
|
| 610 |
|
subroutine InitializeSimGlobals(thisStat) |
| 611 |
|
integer, intent(out) :: thisStat |
| 612 |
|
integer :: alloc_stat |
