lsqr_test Subroutine

public subroutine lsqr_test()

Unit test.

Arguments

None

Calls

proc~~lsqr_test~~CallsGraph proc~lsqr_test lsqrtest_module::lsqr_test proc~test lsqrtest_module::test_solver%test proc~lsqr_test->proc~test 4 4 proc~test->4 proc~acheck lsqr_module::lsqr_solver%acheck proc~test->proc~acheck proc~dcopy lsqpblas_module::dcopy proc~test->proc~dcopy proc~dnrm2 lsqpblas_module::dnrm2 proc~test->proc~dnrm2 proc~lsqr lsqr_module::lsqr_solver%LSQR proc~test->proc~lsqr proc~lstp lsqrtest_module::test_solver%lstp proc~test->proc~lstp proc~xcheck lsqr_module::lsqr_solver%xcheck proc~test->proc~xcheck proc~acheck->proc~dcopy proc~acheck->proc~dnrm2 aprod aprod proc~acheck->aprod proc~ddot lsqpblas_module::ddot proc~acheck->proc~ddot proc~dscal lsqpblas_module::dscal proc~acheck->proc~dscal proc~lsqr->proc~dcopy proc~lsqr->proc~dnrm2 proc~lsqr->aprod proc~d2norm lsqr_module::d2norm proc~lsqr->proc~d2norm proc~lsqr->proc~dscal proc~lstp->proc~dcopy proc~lstp->proc~dnrm2 proc~aprod1 lsqrtest_module::test_solver%aprod1 proc~lstp->proc~aprod1 proc~lstp->proc~dscal proc~hprod lsqrtest_module::hprod proc~lstp->proc~hprod proc~xcheck->proc~dcopy proc~xcheck->proc~dnrm2 proc~xcheck->aprod proc~xcheck->proc~dscal proc~aprod1->proc~hprod

Called by

proc~~lsqr_test~~CalledByGraph proc~lsqr_test lsqrtest_module::lsqr_test program~main main program~main->proc~lsqr_test

Source Code

    subroutine lsqr_test()

    integer :: iunit,nbar,nduplc,n,m,ndamp,npower
    real(wp) :: damp
    type(test_solver) :: solver

    open( newunit=iunit, file='LSQR.LIS', status='REPLACE' )

    solver%nout = iunit

    nbar   = 1000
    nduplc = 40

    m = 2*nbar
    n = nbar
    do ndamp = 2, 7
        npower = ndamp
        damp   = 10.0_wp**(-ndamp-6)
        call solver%test( m, n, nduplc, npower, damp )
    end do

    m = nbar
    n = nbar
    do ndamp = 2, 7
        npower = ndamp
        damp   = 10.0_wp**(-ndamp-6)
        call solver%test( m, n, nduplc, npower, damp )
    end do

    m = nbar
    n = 2*nbar
    do ndamp = 2, 7
        npower = ndamp
        damp   = 10.0_wp**(-ndamp-6)
        call solver%test( m, n, nduplc, npower, damp )
    end do

    close(iunit)

    end subroutine lsqr_test