diff --git a/src/prec/psb_dsplu.f90 b/src/prec/psb_dsplu.f90 index c6aca639..a1a1b487 100644 --- a/src/prec/psb_dsplu.f90 +++ b/src/prec/psb_dsplu.f90 @@ -170,10 +170,10 @@ contains end if do +! write(0,*)'++++',trw%infoa(psb_nnz_),size(trw%ia1),ktrw,k,i if (ktrw > trw%infoa(psb_nnz_)) exit if (trw%ia1(ktrw) > i) exit k = trw%ia2(ktrw) - ! write(0,*)'KKKKK',k if ((k < i).and.(k >= 1)) then l1 = l1 + 1 laspk(l1) = trw%aspk(ktrw) diff --git a/src/serial/psb_dspgtrow.f90 b/src/serial/psb_dspgtrow.f90 index 98c54508..2b2e627c 100644 --- a/src/serial/psb_dspgtrow.f90 +++ b/src/serial/psb_dspgtrow.f90 @@ -300,10 +300,11 @@ contains b%ia2(nzb+k) = (a%ia2(i)) endif enddo + nzb=nzb+k end if end if - b%infoa(psb_nnz_) = nzb + k + b%infoa(psb_nnz_) = nzb b%m = b%m+lrw-irw+1 b%k = max(b%k,a%k) end subroutine coo_dspgtrow diff --git a/test/pargen/RUNS/ppde.inp b/test/pargen/RUNS/ppde.inp index 65f32a65..8c11a1c2 100644 --- a/test/pargen/RUNS/ppde.inp +++ b/test/pargen/RUNS/ppde.inp @@ -1,8 +1,8 @@ 7 Number of entries below this -BICGSTAB Iterative method BICGSTAB CGS BICG BICGSTABL -2 Preconditioner ILU DIAGSC NONE +CGS Iterative method BICGSTAB CGS BICG BICGSTABL +4 Preconditioner ILU DIAGSC NONE 2 Number ov overlapping levels -COO A Storage format CSR COO JAD +CSR A Storage format CSR COO JAD 20 Domain size (acutal sistem is this**3) 1 Stopping criterion 80 MAXIT diff --git a/test/pargen/ppde90.f90 b/test/pargen/ppde90.f90 index 225e7d1d..7b6bedd2 100644 --- a/test/pargen/ppde90.f90 +++ b/test/pargen/ppde90.f90 @@ -126,7 +126,6 @@ program pde90 call dgamx2d(icontxt,'a',' ',ione, ione,t2,ione,t1,t1,-1,-1,-1) if (iam.eq.0) write(*,'("Overall matrix creation time : ",es10.4)')t2 if (iam.eq.0) write(*,'(" ")') - ! ! prepare the preconditioner. !