# Evgenii B. Rudnyi, http://Evgenii.Rudnyi.Ru

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)