D:/demona/2.0/manual/doc/2-31set/processor/ilu/basicilu0CRS.f90

Go to the documentation of this file.

subroutine basicilu0CRS(val_A,column_indicator_A,row_pointer_A,n,valSize,diag)

implicit none
integer n


integer i,j,k,valSize

integer e,a,a1,b,c,c1,d,g,f

real(8) val_A(valSize)
integer column_indicator_A(valSize),row_pointer_A(n+1),diag(n)



do i=2,N


    do k=1,row_pointer_A(i+1)-row_pointer_A(i)

                !print*,k
        
        e=0
        a=row_pointer_A(i)
        a1=row_pointer_A(i+1)
        b=a+(k-1)

                !print*,b,column_indicator_A(b)

                !stop




        if (column_indicator_A(b).lt.i) then
            c=row_pointer_A(column_indicator_A(b))
            c1=row_pointer_A(column_indicator_A(b)+1)
            do d=c,c1-1            
                if (column_indicator_A(d).eq.column_indicator_A(b)) then
                    val_A(b)=val_A(b)/val_A(d)
                    e=e+1
                    exit
                end if
            end do
        end if 
        if (e.eq.0) then
            exit
        else
            g=b
            do b=a+1+(k-1),a1-1                 
                do f=c,c1-1                     
                    if (column_indicator_A(f).eq.column_indicator_A(b)) then
                        val_A(b)=val_A(b)-val_A(g)*val_A(f)
                        exit
                    end if                                     
                end do
            end do
        end if 
    end do
end do


































!do k=1,n-1


!end do









!do k=1,n-1
!    do i=k+1,n
!        if (a(i,k) .ne. 0.d0) then
!            a(i,k) = a(i,k) / a(k,k); 
!            do j=k+1,n
!                if (a(i,j) .ne. 0.d0 ) then
!                    a(i,j) = a(i,j) -a(i,k)*a(k,j)
!                end if
!            end do
!        end if
!    end do
!end do





end

---