from numpy import *
import analysis.cMontecarlo as ACM, time

N = 2000
t = 0
for c in range (1):
    W = 0.1 * abs (random.rand (N, N)) + identity (N)
    s = abs (random.rand (N))
    t1 = time.time ()
    try:
        i, W_sh, m, no_iter = ACM.rec_solve (W, s)
        t2 = time.time ()
        t += t2-t1
        print t, no_iter
    except RuntimeError, msg:
        print msg
print t