Bug fixed

This commit is contained in:
Ricardo 2014-10-03 11:55:58 +01:00
parent 94b19242f0
commit 45f8751769

View file

@ -320,7 +320,8 @@ class HalfT(Prior):
return "hT(" + str(np.round(self.A)) + ', ' + str(np.round(self.nu)) + ')'
def lnpdf(self,theta):
return self.constant*(theta>0) -.5*(self.nu+1) * np.log( 1.+ (1./self.nu) * (theta/self.A)**2 )
stop
return (theta>0) * ( self.constant -.5*(self.nu+1) * np.log( 1.+ (1./self.nu) * (theta/self.A)**2 ) )
#theta = theta if isinstance(theta,np.ndarray) else np.array([theta])
#lnpdfs = np.zeros_like(theta)