calc.py

import os, sys
from math import *
import numpy
from numpy import *
import scipy.linalg
#from scipy.linalg import *
from numpy.linalg import *
from numpy.random import *

#Matrix = numpy.matrix

#import scipy
#from scipy import *
#from scipy.linalg import *
from time import time

#import numarray
#import Numeric
#if os.getenv('NUMERIX') == 'numarray':
#    from numarray.matrix import Matrix
#else:
#    from scipy.linalg import *
#    from Matrix import Matrix

#def adj(X):
#    return X.H # conj(transpose(X))

#def inv(X):
#    return X.I # Matrix(numpy.linalg.inv(M))

#def eigvals(X):
#    print "eigvals:",
#    X = asarray_chkfinite(X)
#    geev, = scipy.linalg.lapack.get_lapack_funcs(('geev',),(X,))
#    print geev.module_name
#    return scipy.linalg.eigvals(X)
#    return numpy.linalg.eigvals(X)

#def integral(x,y):
#    return .5*sum((x[1:]-x[:-1])*(y[1:]+y[:-1]))

#def norm(X):
#    v = asarray(X).flatten()
#    return abs(dot(v,conj(v))**.5)

#def vdot(a,b):
#    sum = 0.0
#    for i in range(len(a)):
#        sum += a[i]*b[i]
#    return sum

#def gcd(a,b):
#    if a>b:
#        a,b = b,a
#    while a != 0:
#        a,b = b%a,a
#    return b

#def lcm(a,b):
#    return a*b/gcd(a,b)

class timer:
    def __init__(self,msg):
        print "%s..."%msg,
        sys.stdout.flush()
        self.starttime = time()

    def stop(self):
        print " %g s"%(time()-self.starttime)