When I do pm.logp(Y_obs, test_val).eval() I don’t get any NaN values.
Is there anything the smc sampler could struggle with that I am not aware of, like HalfNormal distributions, matrix output of the function I am using, pytensor.dimshuffle() or pytensor.switch()?