working dsymv in RHS function
working sparse matrix multiply