Switched to btgs solver to increase speed
Handle cases when input or pmf is zero
Unit test plnorm and plnorm_lowtrunc