enddo
x=0.0d0
+#ifdef MPI
time00=MPI_Wtime()
+#else
+ time00=tcpu()
+#endif
c Compute the stochastic forces acting on bodies. Store in force.
do i=nnt,nct-1
sig=stdforcp(i)
force(j,i+nres)=anorm_distr(x,sig2,lowb2,highb2)
enddo
enddo
+#ifdef MPI
time_fsample=time_fsample+MPI_Wtime()-time00
+#else
+ time_fsample=time_fsample+tcpu()-time00
+#endif
c Compute the stochastic forces acting on virtual-bond vectors.
do j=1,3
ff(j)=0.0d0
c------------------------------------------------------------------
subroutine setup_fricmat
implicit real*8 (a-h,o-z)
+#ifdef MPI
include 'mpif.h'
+#endif
include 'DIMENSIONS'
include 'COMMON.VAR'
include 'COMMON.CHAIN'
if (nfgtasks.gt.1) then
if (fg_rank.eq.0) then
c The matching BROADCAST for fg processors is called in ERGASTULUM
+#ifdef MPI
time00=MPI_Wtime()
+#else
+ time00=tcpu()
+#endif
call MPI_Bcast(10,1,MPI_INTEGER,king,FG_COMM,IERROR)
time_Bcast=time_Bcast+MPI_Wtime()-time00
c print *,"Processor",myrank,
& myginv_ng_count,MPI_DOUBLE_PRECISION,king,FG_COMM,IERROR)
time_scatter=time_scatter+MPI_Wtime()-time00
#ifdef TIMING
+#ifdef MPI
time_scatter_fmat=time_scatter_fmat+MPI_Wtime()-time00
+#else
+ time_scatter_fmat=time_scatter_fmat+tcpu()-time00
+#endif
#endif
do i=1,dimen
do j=1,2*my_ng_count