numpy.linglg.singular_value_decomposition is defined as follows: def singular_value_decomposition(A, full_matrices=0): return svd(A, 0) Shouldn't that last line be return svd(A, full_matrices) Zach