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

[emece67]

OK, I see now. I've changed it to:

Code:

      try:
        fpl = f(bpa2 + bma2*r)
      except ArithmeticError:
        fpl = 0;
      p = fpl*w if isnormal(fpl) else 0
      try:
        fmi = f(bpa2 + bma2/r) if expsinh else f(bpa2 - bma2*r)
      except ArithmeticError:
        fmi = 0
      tnfe += 2
      p += (fmi/w if expsinh else fmi*w) if isnormal(fmi) else 0
      p *= ch
      wsl += p

which seems more readable to me. I've committed it to GitHub.

Regards.