cuIBM/helpers_8cu_source.html

 #include "helpers.h"


 real dhRoma(real x, real h)
 {
     real r = fabs(x)/h;

     if(r>1.5)
         return 0.0;
     else if(r>0.5 && r<=1.5)
         return 1.0/(6*h)*( 5.0 - 3.0*r - sqrt(-3.0*(1-r)*(1-r) + 1.0) );
     else
         return 1.0/(3*h)*( 1.0 + sqrt(-3.0*r*r + 1.0) );
 } // dhRoma


 real delta(real x, real y, real h)
 {
     return dhRoma(x, h) * dhRoma(y, h);
 } // delta
real
double real
Is a float or a double depending on the machine precision.
Definition: types.h:116

delta
real delta(real x, real y, real h)
Two-dimensional discrete delta function.
Definition: helpers.cu:44

dhRoma
real dhRoma(real x, real h)
Discrete delta function from Roma et al. (1999).
Definition: helpers.cu:22

helpers.h
Declaration of the discrete delta function.