Re: [Help-gsl] Question on Cholesky decomposition error - bspline penalty matrix

[Rhys Ulerich]

>> \Omega_{ij} = \int_0^1 { d^2B_i(x)/dx^2 * d^2 B_j(x)/dx^2 } dx

> *p += wl * gsl_pow_2(dB_i * dB_j);

Ha.  I apparently cheated.  Remove that gsl_pow_2 and I encounter the
same issues.

*shrugs*

- Rhys