cur_resid[2 * i + 1, :] = cur_zerob[2 * i + 1, :] + cur_resp[i, :]
cur_ref[i, :] = (cur_refall[2 * i, :] - cur_refall[2 * i + 1, :]) / 2.0
+ print 'cur_offs'
+ print cur_offs
print cur_zerob
print cur_resid
#print cur_resp
U = cur_resp.T
C = np.dot(Y, linalg.inv(U))
+ print 'C value'
print C
+ print 'Y value'
print Y
print np.dot(C, U)
for i in range(0, rows):