level.py 文件源码

def sum(self,x):
        """  sum(x) """
        t0 = time() 
#        MPI.COMM_WORLD.Barrier()
        nbduplicates = (self.np0*self.mp0)/(self.np*self.mp)
        local_sum = computesum(self.msk,x,self.nh)
        t1 = time()
        self.time['sum']+=t1-t0
        self.ncalls['sum']+=1
        global_sum=array(MPI.COMM_WORLD.allgather(local_sum))
        t0 = time()
        self.time['reduce']+=t0-t1
        self.ncalls['reduce']+=1

        return global_sum.sum() / nbduplicates

#----------------------------------------