code1.c

   1 //kernel to test the convergence of roots
   2 __global__
   3 void TestConvergence_kernel(int size,cuDoubleComplex *Z,
   4 cuDoubleComplex *ZPrec,double *Error)
   5 { int i= blockIdx.x*blockDim.x+ threadIdx.x;
   6
   7 if(i<size) {
   8    d_Error[i]=cuCmodule(cuCdiv(cuCsub(d_Z[i],d_ZPrec[i]), d_Z[i]));
   9
  10         }
  11 }