Matrix vector product without computing the matrix

[tmigot]

Right now hprod, jprod and jtprod all need to compute the matrix (in sparce coo-format though) to compute the matrix-vector products.