Re: [ quantlib-Patches-3047358 ] fixed gaussian orthogonal polynomial

Posted by Bojan Nikolic on
URL: http://quantlib.414.s1.nabble.com/quantlib-Patches-3047358-fixed-gaussian-orthogonal-polynomial-tp9136p9137.html


This probably doesn't work because the logic was reversed here:

,----
|         if (std::fabs(denom)<QL_EPSILON) {
|             if (std::fabs(num)<QL_EPSILON) {
|                 QL_FAIL("can't compute b_k for jacobi integration\n");
`----

I think it should be something like

,----
|         if (std::fabs(denom)<QL_EPSILON) {
|             if (std::fabs(num)>QL_EPSILON) {
|                 QL_FAIL("can't compute b_k for jacobi integration\n");
`----

This test exists because l'Hopital rule doesn't work if the nominator is
finite, so an error is signalled if denom is close to zero while num is
_not_ close to zero.

Best,
Bojan

--
Bojan Nikolic          ||          http://www.bnikolic.co.uk/ql/

------------------------------------------------------------------------------
Start uncovering the many advantages of virtual appliances
and start using them to simplify application deployment and
accelerate your shift to cloud computing.
http://p.sf.net/sfu/novell-sfdev2dev
_______________________________________________
QuantLib-dev mailing list
[hidden email]
https://lists.sourceforge.net/lists/listinfo/quantlib-dev