Debugging Functions

Functions
void	Force_debug (double ave_dSdpi[3], double *dSdpi)
	Gets average for each generator of the Lie algebra.

Detailed Description

Function Documentation

Force_debug()

void Force_debug	(	double	ave_dSdpi[3],
		double *	dSdpi )

Gets average for each generator of the Lie algebra.

Parameters

[out]	ave_dSdpi	The average
[in]	dSdpi	What's being averaged. Can be the momentum either

Postcondition

ave_dSdpi overwritten with the averages.

Definition at line 19 of file integrate.c.

                                                   {
#ifdef USE_GPU
   ave_dSdpi[0]=cureduce_sum_d(dSdpi,kvol*ndim,0);
   ave_dSdpi[1]=cureduce_sum_d(dSdpi+kvol*ndim,kvol*ndim,1);
   ave_dSdpi[2]=cureduce_sum_d(dSdpi+2*kvol*ndim,kvol*ndim,2);
   cudaDeviceSynchronise();
#else
   ave_dSdpi[0]=0; ave_dSdpi[1]=0; ave_dSdpi[2]=0;
   for(unsigned int i=0;i<kvol*ndim;i++){
      ave_dSdpi[0]+=dSdpi[i];
      ave_dSdpi[1]+=dSdpi[i+kvol*ndim];
      ave_dSdpi[2]+=dSdpi[i+2*kvol*ndim];
   }
#endif
   ave_dSdpi[0]/=(ndim*kvol); ave_dSdpi[1]/=(ndim*kvol); ave_dSdpi[2]/=(ndim*kvol);
   return;
}

References cudaDeviceSynchronise, cureduce_sum_d(), kvol, and ndim.

Here is the call graph for this function:

Here is the caller graph for this function: