Home | History | Annotate | Download | only in lapack
      1 *> \brief \b ILACLR
      2 *
      3 *  =========== DOCUMENTATION ===========
      4 *
      5 * Online html documentation available at 
      6 *            http://www.netlib.org/lapack/explore-html/ 
      7 *
      8 *> \htmlonly
      9 *> Download ILACLR + dependencies 
     10 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/ilaclr.f"> 
     11 *> [TGZ]</a> 
     12 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/ilaclr.f"> 
     13 *> [ZIP]</a> 
     14 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/ilaclr.f"> 
     15 *> [TXT]</a>
     16 *> \endhtmlonly 
     17 *
     18 *  Definition:
     19 *  ===========
     20 *
     21 *       INTEGER FUNCTION ILACLR( M, N, A, LDA )
     22 * 
     23 *       .. Scalar Arguments ..
     24 *       INTEGER            M, N, LDA
     25 *       ..
     26 *       .. Array Arguments ..
     27 *       COMPLEX            A( LDA, * )
     28 *       ..
     29 *  
     30 *
     31 *> \par Purpose:
     32 *  =============
     33 *>
     34 *> \verbatim
     35 *>
     36 *> ILACLR scans A for its last non-zero row.
     37 *> \endverbatim
     38 *
     39 *  Arguments:
     40 *  ==========
     41 *
     42 *> \param[in] M
     43 *> \verbatim
     44 *>          M is INTEGER
     45 *>          The number of rows of the matrix A.
     46 *> \endverbatim
     47 *>
     48 *> \param[in] N
     49 *> \verbatim
     50 *>          N is INTEGER
     51 *>          The number of columns of the matrix A.
     52 *> \endverbatim
     53 *>
     54 *> \param[in] A
     55 *> \verbatim
     56 *>          A is array, dimension (LDA,N)
     57 *>          The m by n matrix A.
     58 *> \endverbatim
     59 *>
     60 *> \param[in] LDA
     61 *> \verbatim
     62 *>          LDA is INTEGER
     63 *>          The leading dimension of the array A. LDA >= max(1,M).
     64 *> \endverbatim
     65 *
     66 *  Authors:
     67 *  ========
     68 *
     69 *> \author Univ. of Tennessee 
     70 *> \author Univ. of California Berkeley 
     71 *> \author Univ. of Colorado Denver 
     72 *> \author NAG Ltd. 
     73 *
     74 *> \date April 2012
     75 *
     76 *> \ingroup complexOTHERauxiliary
     77 *
     78 *  =====================================================================
     79       INTEGER FUNCTION ILACLR( M, N, A, LDA )
     80 *
     81 *  -- LAPACK auxiliary routine (version 3.4.1) --
     82 *  -- LAPACK is a software package provided by Univ. of Tennessee,    --
     83 *  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
     84 *     April 2012
     85 *
     86 *     .. Scalar Arguments ..
     87       INTEGER            M, N, LDA
     88 *     ..
     89 *     .. Array Arguments ..
     90       COMPLEX            A( LDA, * )
     91 *     ..
     92 *
     93 *  =====================================================================
     94 *
     95 *     .. Parameters ..
     96       COMPLEX          ZERO
     97       PARAMETER ( ZERO = (0.0E+0, 0.0E+0) )
     98 *     ..
     99 *     .. Local Scalars ..
    100       INTEGER I, J
    101 *     ..
    102 *     .. Executable Statements ..
    103 *
    104 *     Quick test for the common case where one corner is non-zero.
    105       IF( M.EQ.0 ) THEN
    106          ILACLR = M
    107       ELSE IF( A(M, 1).NE.ZERO .OR. A(M, N).NE.ZERO ) THEN
    108          ILACLR = M
    109       ELSE
    110 *     Scan up each column tracking the last zero row seen.
    111          ILACLR = 0
    112          DO J = 1, N
    113             I=M
    114             DO WHILE((A(MAX(I,1),J).EQ.ZERO).AND.(I.GE.1))
    115                I=I-1
    116             ENDDO
    117             ILACLR = MAX( ILACLR, I )
    118          END DO
    119       END IF
    120       RETURN
    121       END
    122