--- rpl/lapack/lapack/dpotrf2.f 2016/08/27 15:34:35 1.2 +++ rpl/lapack/lapack/dpotrf2.f 2017/06/17 10:54:01 1.3 @@ -2,14 +2,14 @@ * * =========== DOCUMENTATION =========== * -* Online html documentation available at -* http://www.netlib.org/lapack/explore-html/ +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ * * Definition: * =========== * * RECURSIVE SUBROUTINE DPOTRF2( UPLO, N, A, LDA, INFO ) -* +* * .. Scalar Arguments .. * CHARACTER UPLO * INTEGER INFO, LDA, N @@ -17,7 +17,7 @@ * .. Array Arguments .. * REAL A( LDA, * ) * .. -* +* * *> \par Purpose: * ============= @@ -41,7 +41,7 @@ *> *> The subroutine calls itself to factor A11. Update and scale A21 *> or A12, update A22 then calls itself to factor A22. -*> +*> *> \endverbatim * * Arguments: @@ -94,22 +94,22 @@ * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Denver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. * -*> \date November 2015 +*> \date December 2016 * *> \ingroup doublePOcomputational * * ===================================================================== RECURSIVE SUBROUTINE DPOTRF2( UPLO, N, A, LDA, INFO ) * -* -- LAPACK computational routine (version 3.6.0) -- +* -- LAPACK computational routine (version 3.7.0) -- * -- LAPACK is a software package provided by Univ. of Tennessee, -- * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- -* November 2015 +* December 2016 * * .. Scalar Arguments .. CHARACTER UPLO @@ -126,7 +126,7 @@ PARAMETER ( ONE = 1.0D+0, ZERO = 0.0D+0 ) * .. * .. Local Scalars .. - LOGICAL UPPER + LOGICAL UPPER INTEGER N1, N2, IINFO * .. * .. External Functions .. @@ -189,7 +189,7 @@ IF ( IINFO.NE.0 ) THEN INFO = IINFO RETURN - END IF + END IF * * Compute the Cholesky factorization A = U**T*U * @@ -198,10 +198,10 @@ * Update and scale A12 * CALL DTRSM( 'L', 'U', 'T', 'N', N1, N2, ONE, - $ A( 1, 1 ), LDA, A( 1, N1+1 ), LDA ) + $ A( 1, 1 ), LDA, A( 1, N1+1 ), LDA ) * * Update and factor A22 -* +* CALL DSYRK( UPLO, 'T', N2, N1, -ONE, A( 1, N1+1 ), LDA, $ ONE, A( N1+1, N1+1 ), LDA ) CALL DPOTRF2( UPLO, N2, A( N1+1, N1+1 ), LDA, IINFO ) @@ -216,7 +216,7 @@ * * Update and scale A21 * - CALL DTRSM( 'R', 'L', 'T', 'N', N2, N1, ONE, + CALL DTRSM( 'R', 'L', 'T', 'N', N2, N1, ONE, $ A( 1, 1 ), LDA, A( N1+1, 1 ), LDA ) * * Update and factor A22