From aa376a22b73ddb5a59122112d842ebd2dae5b98c Mon Sep 17 00:00:00 2001 From: Salvatore Filippone Date: Tue, 2 May 2006 12:45:07 +0000 Subject: [PATCH] Fix bug in size estimate for CSR with diagonal matrices. --- src/serial/psb_cest.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/serial/psb_cest.f90 b/src/serial/psb_cest.f90 index 2434b251..7f84807a 100644 --- 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