source: altifloat/src/u_d.h

forward (YREAL ud, YREAL u_norm)
{
  //YS1=ud;             
  YS1=ud*u_norm*mask[Yi][Yj];
}

backward(YREAL ud, YREAL u_norm)
{
  /*if (Yi==9 && Yj==16)
    {
      printf("YG=%g\n",YG1_u_d(Yi,Yj));
      }*/
//YJ1I1=1;
 //YJ1I2=0;
  YJ1I1=u_norm*mask[Yi][Yj];
  YJ1I2=ud*mask[Yi][Yj];
}