Hi all, I have orthonormal set of vectors B = [b_0, b_1,..., b_k-1], b_i from R^n (k may be less than n), and vector a from R^n What is most efficient way in numpy to get r from R^n and c_0, ..., c_k-1 from R: a = c_0*b_0+...+c_k-1*b_k-1 + r (r is rest) Thank you in advance, D.