# Evgenii B. Rudnyi, http://MatrixProgramming.com
# the matrix and vector can be found among the files in the Matrix Programming group
# http://groups.google.com/group/matrixprogramming/files

import numpy
import scipy.linalg
import scipy.io
import time

mat = scipy.io.mmread('Test_matrix.mmf')
print 'matrix', mat.shape
vec = scipy.io.read_array('Test_vector.txt')
print 'vector', vec.shape

start = time.clock()
x = scipy.linalg.solve(mat, vec)
finish = time.clock()
print 'time for solution is ', finish - start,'s'
print 'residual', scipy.linalg.norm(numpy.dot(mat, x) - vec)