def test_orthogonal_project_to_b0(self) :
print "orthogonal projection to 1st vector basis"
print "========================================="
b1 = np.array([
[4,5,1],
[4,8,2],
[6,2,6]
])
b1_ = orthoproject(b1[0], b1[1:])
print b1_
b2 = np.array([
[-2, 7, 7, -5],
[3, -2, 6, -1],
[2, -8, -9, -7],
[8, -9, 6, -4],
])
b2_ = orthoproject(b2[0], b2[1:])
print b2_
def test_orthoproject_2(self) :
x = np.array([[4.0, 0.2, -0.6],
[0.9268292682926829, -1.853658536585366, 5.560975609756098]] )
x_ = orthoproject( x[0], x[1:] )
print x_
