InsightSoftwareConsortium/ITK — pull request #4 by daniel-perry (2014-07-24 00:30:37+00:00)
Title: vnl_sparse_matrix::mult - size of q buffer should actually be (this->rows())*pcols.

I ran into a subtle bug in vnl_sparse_matrix::mult(), where the code assumes the p and q buffers are the same size (prows*pcols). However the size of q should really be (this->rows())*pcols.

ie, in matrix multiplication an (a,b) size matrix multiplied with a (b,c) size matrix results in a (a,c) size matrix, not (b,c).