Merge pull request #5 from jalvesz/dev

merge from dev
jalvesz · May 11, 2024 · cea1950 · cea1950
2 parents bf71378 + d05b1dc
commit cea1950
Show file tree

Hide file tree

Showing 3 changed files with 184 additions and 209 deletions.
diff --git a/README.md b/README.md
@@ -49,9 +49,8 @@ call dense2coo( dense , coo )
 ### Sparse Matrix-Vector product
 (available) Matrix vector products are interfaced by the procedure
 ```fortran
-call matvec( Mat , vec_x, vec_y ) ! vec_y = Mat * vec_x
+call matvec( Mat , vec_x, vec_y ) ! vec_y = vec_y + Mat * vec_x
 ```
-**NOTE**: The matvec in this lib do not initialize `vec_y=0`. Why? enable easy recycling of data to do `vec_y=vec_y+Mat*vec_x`. If a clean vector is needed make sure to put a `vec_y=0` just before calling matvec.
 
 | Matrix | full | symmetric |
 |--------|-------|------------------|