def f(x, params): alpha = params return Numeric.log(alpha * x) / Numeric.sqrt(x)
def testLU_lndet(self): (lu, p, s) = LU_decomp(self.m1_4) lndet = LU_lndet(lu) result = fpcompare(lndet, numx.log(15.744), 3) self.assertTrue(result)
def testLU_lndet(self): (lu, p, s) = LU_decomp(self.m1_4) lndet = LU_lndet(lu) result = fpcompare(lndet, numx.log(15.744), 3) self.failUnless(result)
def f(x, params): alpha = params return Numeric.log(alpha *x) / Numeric.sqrt(x)