From: Rob S. <rs...@MI...> - 2010-03-26 16:33:09
|
> The matvec behavior is the intended behavior. The whole idea of the > PysparseMatrix class is that it is "easy" and "natural" to use. > Amongst other things, that means you can say y=A*x instead of > A.matvec(x,y). The savvy user will notice (as you did) that ll_mat > objects are more efficient and will use those. The typical user will > use the PysparseMatrix class. That creates a rather broken API. Having an inplace version of matvec is fine, but please at least call it "matvec_inplace" or something, in both ll_mat and jdsym, to make it clear that it's not supposed to be matvec. The current situation is very confusing. -- Rob |