*** empty log message ***

psblas3-type-indexed
Alfredo Buttari 19 years ago
parent 0e235de552
commit b878b6dc8f

@ -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)

@ -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

@ -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

@ -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.
!

Loading…
Cancel
Save