[SciPy-user] cholesky decomposition fails

[Robert Kern]

Chris Fonnesbeck wrote:

> I did find a bug in basic_lite, in the cholesky_decomposition
> function, however.
> 
> --> 120     return Numeric.transpose(MLab.triu(a,k=0)).copy()
>         global Numeric.transpose = <function transpose at 0x444830>
>         global MLab.triu = undefined
>         a = array([[ 1.,  0.],
>        [ 0.,  1.]])
>         global k.copy = undefined
>     121
>     122
> 
> NameError: global name 'MLab' is not defined
> 
> Apparently, MLab is not imported by this module.

>From the top of basic_lite.py:

import scipy.base as Numeric
import scipy.lib.lapack_lite as lapack_lite
MLab = Numeric
from scipy.base import asarray, multiply
import math

Are you sure you have the latest checkout from SVN?

-- 
Robert Kern
rkern at ucsd.edu

"In the fields of hell where the grass grows high
 Are the graves of dreams allowed to die."
  -- Richard Harter