My main interest is sparse matrix-vector and matrix-transpose-vector multiplication of the form y=y+AA'x. Is there any library that performs y=y+AA'x efficiently? I have investigated SPARSKIT, MKL, ARPACK,... but I was not able to find the related routine. Isn't the y=y+AA'x operation widely used?
Library that performs sparse matrix-vector and matrix-transpose-vector multiplication
sparse;matrix