(PC-12xx~14xx) qthsh Tanh-Sinh quadrature

[robve]

(04-12-2021 12:20 PM)Albert Chan Wrote:  To reduce noise, we should quit summing, if we actually hit the finite edge.
So this is the proposed patch:

Code:

      else {                // Exp-Sinh
        x = c + d/r;
        if (x == c) break;  // x hit finite edge
        y = f(x);
        if (isfinite(y)) q += y/w;
      }

And remove the guard x<b for Exp-Sinh and Sinh-Sinh:

Code:

        x = c + d*r;
        y = f(x);
        if (isfinite(y))       // if f(x) is finite, add to local sum
          q += y*w