# test_simeq.py  simple test
from simeq import *

print "test simeq.py running"
a = [[0 for j in range(4)] for i in range(4)]
x = [0 for j in range(4)]
y = [0 for j in range(4)]
t = [0 for j in range(4)]
a[0][0] = 1.0
a[0][1] = 1.0
a[0][2] = 2.0
a[0][3] = 2.0
a[1][0] = 2.0
a[1][1] = 3.0
a[1][2] = 3.0
a[1][3] = 5.0
a[2][0] = 1.0
a[2][1] = 2.0
a[2][2] = 4.0
a[2][2] = 5.0
a[3][0] = 2.0
a[3][1] = 2.0
a[3][2] = 1.0
a[3][3] = 3.0
print "a="
print a
y[0] = 25.0
y[1] = 26.0
y[2] = 22.0
y[3] = 31.0
print "y="
print y
x = simeq(a,y)
print "x="
print x

for i in range(4):
    t[i] = 0.0
    for j in range(4):
        t[i]+=a[i][j]*x[j]

    print "i=",i,", err=",(y[i]-t[i])

print "test_simeq.py finished"