I was writing a model and I used the Gaussian class. I had a weird bug that went away when I switched to boost::math::normal. I should write some tests to figure out why, and perhaps scrap the Gaussian class.
Weird, it seems to be working fine when I test Gaussian in isolation. Maybe I had a bug that I just couldn't see, but ended up correcting...closing for now.
I was writing a model and I used the Gaussian class. I had a weird bug that went away when I switched to boost::math::normal. I should write some tests to figure out why, and perhaps scrap the Gaussian class.