Fix bug in size estimate for CSR with diagonal matrices.

19 years ago · aa376a22b7
parent 208952a9a5
commit aa376a22b7
1 changed files with 1 additions and 1 deletions
--- a/src/serial/psb_cest.f90
+++ b/src/serial/psb_cest.f90
@ -81,7 +81,7 @@ subroutine psb_cest(afmt, m,n,nnz, lia1, lia2, lar, iup, info)
      lar = nnz
    else if(afmt.eq.'CSR') then
      lia1 = nnz
-      lia2 = nnz
+      lia2 = max(nnz,m+1)
      lar = nnz
    else
      info = 136