Type error from numpy in ufunc
1.22.0 version of numpy gives type error at these lines in mcsampler.py
Here is the full error message
Traceback (most recent call last):
File "/home/vinaya.valsan/.conda/envs/igwn-py38-clone_rapidpe_rift_devl/bin/integrate_likelihood_extrinsic_batchmode", line 7, in <module>
exec(compile(f.read(), __file__, 'exec'))
File "/home/vinaya.valsan/.conda/envs/igwn-py38-clone_rapidpe_rift_devl/src/rift/MonteCarloMarginalizeCode/Code/bin/integrate_likelihood_extrinsic_batchmode", line 1426, in <module>
res = analyze_event(P_list, indx, data_dict, psd_dict, fmax, opts)
File "/home/vinaya.valsan/.conda/envs/igwn-py38-clone_rapidpe_rift_devl/src/rift/MonteCarloMarginalizeCode/Code/bin/integrate_likelihood_extrinsic_batchmode", line 1236, in analyze_event
res, var, neff, dict_return = sampler.integrate(likelihood_function, *unpinned_params, **pinned_params)
File "/home/vinaya.valsan/.conda/envs/igwn-py38-clone_rapidpe_rift_devl/src/rift/MonteCarloMarginalizeCode/Code/RIFT/integrators/mcsampler.py", line 484, in integrate
p_s, p_prior, rv = self.draw(n, *self.params_ordered,**args_draw) # keep in order
File "/home/vinaya.valsan/.conda/envs/igwn-py38-clone_rapidpe_rift_devl/src/rift/MonteCarloMarginalizeCode/Code/RIFT/integrators/mcsampler.py", line 233, in draw
res.append((self.pdf[param](numpy.float64(cdf_rv)).astype(numpy.float64)/self._pdf_norm[param], self.prior_pdf[param](cdf_rv), cdf_rv))
File "/home/vinaya.valsan/.conda/envs/igwn-py38-clone_rapidpe_rift_devl/src/rift/MonteCarloMarginalizeCode/Code/RIFT/integrators/mcsampler.py", line 797, in <lambda>
uniform_samp_theta = lambda x,numpy=numpy: 0.5*numpy.sin(x)
The error is solved if I change x
to x.astype(float)
. But a more sensible solution is to make sure that the input to this function in the ILE script is in the correct float format.