sparskit.f90

 
module sparskit
    implicit none
 
    ! BLASSM.F
    interface
 
        subroutine amub(nrow, ncol, job, a, ja, ia, b, jb, ib, c, jc, ic, &
            nzmax, iw, ierr)
            use iso_fortran_env, only : int32, real64
            integer(int32), intent(in) :: nrow, ncol, job, nzmax
            integer(int32), intent(in) :: ja(*), ia(nrow+1), jb(*), ib(*)
            integer(int32), intent(out) :: jc(*), ic(*), iw(ncol), ierr
            real(real64), intent(in) :: a(*), b(*)
            real(real64), intent(out) :: c(*)
        end subroutine
 
        subroutine aplb(nrow, ncol, job, a, ja, ia, b, jb, ib, c, jc, ic, &
            nzmax, iw, ierr)
            use iso_fortran_env, only : int32, real64
            integer(int32), intent(in) :: nrow, ncol, job, nzmax
            integer(int32), intent(in) :: ja(*), ia(nrow+1), jb(*), ib(nrow+1)
            real(real64), intent(in) :: a(*), b(*)
            real(real64), intent(out) :: c(*)
            integer(int32), intent(out) :: jc(*), ic(nrow+1), iw(ncol), ierr
        end subroutine
 
        subroutine aplsb(nrow, ncol, a, ja, ia, s, b, jb, ib, c, jc, ic, &
            nzmax, iw, ierr)
            use iso_fortran_env, only : int32, real64
            integer(int32), intent(in) :: nrow, ncol, nzmax
            integer(int32), intent(in) :: ja(*), ia(nrow+1), jb(*), ib(nrow+1)
            integer(int32), intent(out) :: ierr
            integer(int32), intent(out) :: jc(*), ic(nrow+1), iw(ncol)
            real(real64), intent(in) :: s
            real(real64), intent(in) :: a(*), b(*)
            real(real64), intent(out) :: c(*)
        end subroutine
    end interface
 
    ! FORMATS.F
    interface
 
        subroutine csrcsc2(n, n2, job, ipos, a, ja, ia, a0, ja0, ia0)
            use iso_fortran_env, only : int32, real64
            integer(int32), intent(in) :: n, n2, job, ipos, ja(*), ia(n+1)
            integer(int32), intent(out) :: ja0(*), ia0(n2+1)
            real(real64), intent(in) :: a(*)
            real(real64), intent(out) :: a0(*)
        end subroutine
 
        subroutine bndcsr(n, abd, nabd, lowd, ml, mu, a, ja, ia, len, ierr)
            use iso_fortran_env, only : int32, real64
            integer(int32), intent(in) :: n, nabd, lowd, ml, mu, len
            real(real64), intent(in) :: abd(nabd,*)
            real(real64), intent(out) :: a(*)
            integer(int32), intent(out) :: ia(n+1), ja(*), ierr
        end subroutine
 
        subroutine csrmsr(n, a, ja, ia, ao, jao, wk, iwk)
            use iso_fortran_env, only : int32, real64
            integer(int32), intent(in) :: n
            integer(int32), intent(inout) :: ja(*), ia(n+1)
            integer(int32), intent(out) :: jao(*), iwk(n+1)
            real(real64), intent(inout) :: a(*)
            real(real64), intent(out) :: ao(*), wk(n)
        end subroutine
 
        subroutine msrcsr(n, a, ja, ao, jao, iao, wk, iwk)
            use iso_fortran_env, only : int32, real64
            integer(int32), intent(in) :: n, ja(*)
            integer(int32), intent(out) :: jao(*), iao(n+1), iwk(n+1)
            real(real64), intent(in) :: a(*)
            real(real64), intent(out) :: ao(*), wk(n)
        end subroutine
 
        subroutine coocsr(nrow, nnz, a, ir, jc, ao, jao, iao)
            use iso_fortran_env, only : int32, real64
            integer(int32), intent(in) :: nrow, nnz, jc(*)
            integer(int32), intent(inout) :: ir(*)
            real(real64), intent(in) :: a(*)
            integer(int32), intent(out) :: jao(*), iao(*)
            real(real64) :: ao(*)
        end subroutine
    end interface
 
    ! UNARY.F
    interface
 
        function getelm(i, j, a, ja, ia, iadd, sorted) result(rst)
            use iso_fortran_env, only : int32, real64
            integer(int32), intent(in) :: i, j, ia(*), ja(*)
            real(real64), intent(in) :: a(*)
            integer(int32), intent(out) :: iadd
            logical, intent(in) :: sorted
            real(real64) :: rst
        end function
 
        subroutine getdia(nrow, ncol, job, a, ja, ia, len, diag, idiag, ioff)
            use iso_fortran_env, only : int32, real64
            integer(int32), intent(in) :: nrow, ncol, job, ja(*), ia(*), ioff
            integer(int32), intent(out) :: len, idiag(*)
            real(real64), intent(in) :: a(*)
            real(real64), intent(out) :: diag(*)
        end subroutine
 
        subroutine csort(n, a, ja, ia, values)
            use iso_fortran_env, only : int32, real64
            integer(int32), intent(in) :: n
            real(real64), intent(inout) :: a(*)
            integer(int32), intent(inout) :: ja(*)
            integer(int32), intent(in) :: ia(*)
            logical, intent(in) :: values
        end subroutine
 
 
        subroutine clncsr(job, value2, nrow, a, ja, ia, indu, iwk)
            use iso_fortran_env, only : int32, real64
            integer(int32), intent(in) :: job, value2, nrow
            real(real64), intent(inout) :: a(*)
            integer(int32), intent(inout) :: ja(*), ia(*)
            integer(int32), intent(inout) :: indu(*), iwk(*)
        end subroutine
    end interface
 
    ! ILUT.F
    interface
 
        subroutine ilut(n, a, ja, ia, lfil, droptol, alu, jlu, ju, iwk, w, &
            jw, ierr)
            use iso_fortran_env, only : int32, real64
            integer(int32), intent(in) :: n, ja(*), ia(n+1), lfil, iwk
            integer(int32), intent(out) :: jlu(*), ju(n), jw(2*n), ierr
            real(real64), intent(in) :: a(*), droptol
            real(real64), intent(out) :: alu(*), w(n+1)
        end subroutine
 
        subroutine ilutp(n, a, ja, ia, lfil, droptol, permtol, mbloc, alu, &
            jlu, ju, iwk, w, jw, iperm, ierr)
            use iso_fortran_env, only : int32, real64
            integer(int32), intent(in) :: n, lfil, iwk, mbloc
            integer(int32), intent(inout) :: ja(*), ia(n+1)
            integer(int32), intent(out) :: jlu(*), ju(n), jw(2*n), &
                iperm(2*n), ierr
            real(real64), intent(in) :: droptol, permtol
            real(real64), intent(inout) :: a(*)
            real(real64), intent(out) :: alu(*), w(n+1)
        end subroutine
 
        subroutine ilud(n, a, ja, ia, alph, tol, alu, jlu, ju, iwk, w, jw, ierr)
            use iso_fortran_env, only : int32, real64
            integer(int32), intent(in) :: n, iwk, ja(*), ia(n+1)
            integer(int32), intent(out) :: jlu(*), ju(n), jw(2*n), ierr
            real(real64), intent(in) :: a(*), alph, tol
            real(real64), intent(out) :: alu(*), w(2*n)
        end subroutine
 
        subroutine iludp(n, a, ja, ia, alph, droptol, permtol, mbloc, alu, &
            jlu, ju, iwk, w, jw, iperm, ierr)
            use iso_fortran_env, only : int32, real64
            integer(int32), intent(in) :: n, iwk, mbloc
            integer(int32), intent(inout) :: ja(*), ia(n+1)
            integer(int32), intent(out) :: jlu(*), ju(n), jw(2*n), iperm(2*n), &
                ierr
            real(real64), intent(in) :: alph, droptol, permtol
            real(real64), intent(inout) :: a(*)
            real(real64), intent(out) :: alu(*), w(2*n)
        end subroutine
 
        subroutine pgmres(n, im, rhs, sol, vv, eps, maxits, iout, aa, ja, ia, &
            alu, jlu, ju, ierr)
            use iso_fortran_env, only : int32, real64
            integer(int32), intent(in) :: n, im, maxits, iout, ja(*), ia(n+1), &
                jlu(*), ju(n)
            integer(int32), intent(out) :: ierr
            real(real64), intent(in) :: aa(*), eps, alu(*)
            real(real64), intent(inout) :: rhs(n), sol(n)
            real(real64), intent(out) :: vv(n,*)
        end subroutine
 
        subroutine lusol(n, y, x, alu, jlu, ju)
            use iso_fortran_env, only : int32, real64
            integer(int32), intent(in) :: n, jlu(*), ju(*)
            real(real64), intent(in) :: y(n), alu(*)
            real(real64), intent(out) :: x(n)
        end subroutine
    end interface
end module