Saturation with no discontinuity handling

Saturation The first piece of code does not handle discontinuity. It simply affects the correct value to the output y, depending on the current value of x.

Copy

void saturation_(int *n, double *y, double *x, double rp[2])
{
   double xmax, xmin;

   xmax = rp[0];
   xmin = rp[1];

   if (*x > xmax) /* Limited to higher value. */
   {
      *y = xmax;
   }
   else if (*x < xmin) /* Limited to lower value. */
   {
      *y = xmin;
   }
   else /* Unlimited. */
   {
      *y = *x;
   }
}