Invalid memory accesses in calc_coupling_elem_spin0and2_pure? about pspy HOT 2 CLOSED

I think a correct way to handle these cases is the following

        if (i2 < 1) then
            combin2 = 0d0
         else
            combin2 = fac_b * thrcofb(i2)
        end if

        if (i3 < 1) then
            combin3 = 0d0
         else
            combin3 = fac_c * thrcofc(i3)
        end if

        combin = thrcofa(i1) + combin2 + combin3

where combin2 and combin3 are real(8). I've added a test for pure BB mode with comparisons with NaMaster (see PR #82). The above code (not yet implemented in PR #82) is still "physically" correct but prevent access to undefined array values. Does it look fine for you ?

from pspy.

Looks great to me, thanks!

from pspy.