How do I use gaussian_pdf in Owl.Arr?

[LongPham7]

How do I use gaussian_pdf in Owl.Arr?

[ghennequin]

I don't get a segfault but I get nonsensical results (~1E-22). There must be a bug in Arr.gaussian_pdf. The scalar version Stats.gaussian_pdf works fine. I don't have time to chase this now but perhaps @jzstark can take a look?

[LongPham7]

Hi @ghennequin. It's interesting that you didn't get a segmentation fault. I think the low value in your result (~1E-22) makes sense because I accidentally set sigma to a very low value (0.1) in the above example.

Nonetheless, even if I use a larger (and more reasonable) value for sigma, I still get a segmentation fault. Specifically, the following simpler code produces a segmentation fault:

#require "owl";;
open Owl;;
let mu_arr = Arr.zeros [|2|];;
let sigma_arr = Arr.ones [|2|];;
let current_point = Arr.ones [|2|];;
let result_arr = Arr.gaussian_pdf ~mu:mu_arr ~sigma:sigma_arr current_point;;

My OS is Ubuntu.