/* gdelay2.c - generalized circular delay with real-valued shift */

void gdelay2(D, c, q)
int D;
double c, *q;                             /* \(c\)=shift, \(q\)=offset index */
{
       *q -= c;                           /* decrement by \(c\) */

       if (*q < 0)  
              *q += D+1;

       if (*q > D)
              *q -= D+1;
}