Skip to content

Commit 95a761c

Browse files
langoulangou
authored
Merge pull request #195 from iyamazaki/aasen-2stage
Aasen 2stage
2 parents 9f572f6 + 6204cd9 commit 95a761c

14 files changed

Lines changed: 15 additions & 55 deletions

‎SRC/chesv_aa.f‎

Lines changed: 0 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -129,8 +129,6 @@
129129
*> The length of WORK. LWORK >= MAX(1,2*N,3*N-2), and for best
130130
*> performance LWORK >= MAX(1,N*NB), where NB is the optimal
131131
*> blocksize for CHETRF.
132-
*> for LWORK < N, TRS will be done with Level BLAS 2
133-
*> for LWORK >= N, TRS will be done with Level BLAS 3
134132
*>
135133
*> If LWORK = -1, then a workspace query is assumed; the routine
136134
*> only calculates the optimal size of the WORK array, returns

‎SRC/chesv_aa_2stage.f‎

Lines changed: 2 additions & 10 deletions
Original file line numberDiff line numberDiff line change
@@ -1,7 +1,5 @@
11
*> \brief <b> CHESV_AA_2STAGE computes the solution to system of linear equations A * X = B for HE matrices</b>
22
*
3-
* @precisions fortran c -> s d
4-
*
53
* =========== DOCUMENTATION ===========
64
*
75
* Online html documentation available at
@@ -203,16 +201,10 @@ SUBROUTINE CHESV_AA_2STAGE( UPLO, N, NRHS, A, LDA, TB, LTB,
203201
* ..
204202
*
205203
* =====================================================================
206-
* .. Parameters ..
207-
COMPLEX ZERO, ONE
208-
PARAMETER ( ZERO = ( 0.0E+0, 0.0E+0 ),
209-
$ ONE = ( 1.0E+0, 0.0E+0 ) )
210204
*
211205
* .. Local Scalars ..
212206
LOGICAL UPPER, TQUERY, WQUERY
213-
INTEGER I, J, K, I1, I2, TD
214-
INTEGER LDTB, LWKOPT, NB, KB, NT, IINFO
215-
COMPLEX PIV
207+
INTEGER LWKOPT
216208
* ..
217209
* .. External Functions ..
218210
LOGICAL LSAME
@@ -223,7 +215,7 @@ SUBROUTINE CHESV_AA_2STAGE( UPLO, N, NRHS, A, LDA, TB, LTB,
223215
EXTERNAL XERBLA
224216
* ..
225217
* .. Intrinsic Functions ..
226-
INTRINSIC CONJG, MIN, MAX
218+
INTRINSIC MAX
227219
* ..
228220
* .. Executable Statements ..
229221
*

‎SRC/chetrs_aa_2stage.f‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -204,7 +204,7 @@ SUBROUTINE CHETRS_AA_2STAGE( UPLO, N, NRHS, A, LDA, TB, LTB,
204204
*
205205
* Read NB and compute LDTB
206206
*
207-
NB = TB( 1 )
207+
NB = INT( TB( 1 ) )
208208
LDTB = LTB/N
209209
*
210210
IF( UPPER ) THEN

‎SRC/dsysv_aa_2stage.f‎

Lines changed: 2 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -203,16 +203,10 @@ SUBROUTINE DSYSV_AA_2STAGE( UPLO, N, NRHS, A, LDA, TB, LTB,
203203
* ..
204204
*
205205
* =====================================================================
206-
* .. Parameters ..
207-
DOUBLE PRECISION ZERO, ONE
208-
PARAMETER ( ZERO = ( 0.0D+0, 0.0D+0 ),
209-
$ ONE = ( 1.0D+0, 0.0D+0 ) )
210206
*
211207
* .. Local Scalars ..
212208
LOGICAL UPPER, TQUERY, WQUERY
213-
INTEGER I, J, K, I1, I2, TD
214-
INTEGER LDTB, LWKOPT, NB, KB, NT, IINFO
215-
DOUBLE PRECISION PIV
209+
INTEGER LWKOPT
216210
* ..
217211
* .. External Functions ..
218212
LOGICAL LSAME
@@ -223,7 +217,7 @@ SUBROUTINE DSYSV_AA_2STAGE( UPLO, N, NRHS, A, LDA, TB, LTB,
223217
EXTERNAL XERBLA
224218
* ..
225219
* .. Intrinsic Functions ..
226-
INTRINSIC CONJG, MIN, MAX
220+
INTRINSIC MAX
227221
* ..
228222
* .. Executable Statements ..
229223
*

‎SRC/dsytrs_aa_2stage.f‎

Lines changed: 1 addition & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -1,7 +1,5 @@
11
*> \brief \b DSYTRS_AA_2STAGE
22
*
3-
* @precisions fortran d -> s c
4-
*
53
* =========== DOCUMENTATION ===========
64
*
75
* Online html documentation available at
@@ -204,7 +202,7 @@ SUBROUTINE DSYTRS_AA_2STAGE( UPLO, N, NRHS, A, LDA, TB, LTB,
204202
*
205203
* Read NB and compute LDTB
206204
*
207-
NB = TB( 1 )
205+
NB = INT( TB( 1 ) )
208206
LDTB = LTB/N
209207
*
210208
IF( UPPER ) THEN

‎SRC/ssysv_aa_2stage.f‎

Lines changed: 5 additions & 13 deletions
Original file line numberDiff line numberDiff line change
@@ -1,7 +1,5 @@
11
*> \brief <b> SSYSV_AA_2STAGE computes the solution to system of linear equations A * X = B for SY matrices</b>
22
*
3-
* @generated from SRC/chesv_aasen_2stage.f, fortran c -> s, Tue Oct 31 11:22:31 2017
4-
*
53
* =========== DOCUMENTATION ===========
64
*
75
* Online html documentation available at
@@ -21,8 +19,8 @@
2119
* ===========
2220
*
2321
* SUBROUTINE SSYSV_AA_2STAGE( UPLO, N, NRHS, A, LDA, TB, LTB,
24-
* IPIV, IPIV2, B, LDB, WORK, LWORK,
25-
* INFO )
22+
* IPIV, IPIV2, B, LDB, WORK, LWORK,
23+
* INFO )
2624
*
2725
* .. Scalar Arguments ..
2826
* CHARACTER UPLO
@@ -203,16 +201,10 @@ SUBROUTINE SSYSV_AA_2STAGE( UPLO, N, NRHS, A, LDA, TB, LTB,
203201
* ..
204202
*
205203
* =====================================================================
206-
* .. Parameters ..
207-
REAL ZERO, ONE
208-
PARAMETER ( ZERO = ( 0.0E+0, 0.0E+0 ),
209-
$ ONE = ( 1.0E+0, 0.0E+0 ) )
210-
*
204+
* ..
211205
* .. Local Scalars ..
212206
LOGICAL UPPER, TQUERY, WQUERY
213-
INTEGER I, J, K, I1, I2, TD
214-
INTEGER LDTB, LWKOPT, NB, KB, NT, IINFO
215-
REAL PIV
207+
INTEGER LWKOPT
216208
* ..
217209
* .. External Functions ..
218210
LOGICAL LSAME
@@ -223,7 +215,7 @@ SUBROUTINE SSYSV_AA_2STAGE( UPLO, N, NRHS, A, LDA, TB, LTB,
223215
EXTERNAL XERBLA
224216
* ..
225217
* .. Intrinsic Functions ..
226-
INTRINSIC CONJG, MIN, MAX
218+
INTRINSIC MAX
227219
* ..
228220
* .. Executable Statements ..
229221
*

‎SRC/ssytrf_aa_2stage.f‎

Lines changed: 1 addition & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -1,7 +1,5 @@
11
*> \brief \b SSYTRF_AA_2STAGE
22
*
3-
* @generated from SRC/dsytrf_aa_2stage.f, fortran d -> s, Mon Oct 30 11:57:50 2017
4-
*
53
* =========== DOCUMENTATION ===========
64
*
75
* Online html documentation available at
@@ -29,7 +27,7 @@
2927
* ..
3028
* .. Array Arguments ..
3129
* INTEGER IPIV( * ), IPIV2( * )
32-
* REAL A( LDA, * ), TB( * ), WORK( * )
30+
* REAL A( LDA, * ), TB( * ), WORK( * )
3331
* ..
3432
*
3533
*> \par Purpose:

‎SRC/ssytrs_aa_2stage.f‎

Lines changed: 1 addition & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -1,7 +1,5 @@
11
*> \brief \b SSYTRS_AA_2STAGE
22
*
3-
* @generated from SRC/dsytrs_aa_2stage.f, fortran d -> s, Mon Oct 30 11:59:02 2017
4-
*
53
* =========== DOCUMENTATION ===========
64
*
75
* Online html documentation available at
@@ -204,7 +202,7 @@ SUBROUTINE SSYTRS_AA_2STAGE( UPLO, N, NRHS, A, LDA, TB, LTB,
204202
*
205203
* Read NB and compute LDTB
206204
*
207-
NB = TB( 1 )
205+
NB = INT( TB( 1 ) )
208206
LDTB = LTB/N
209207
*
210208
IF( UPPER ) THEN

‎SRC/zhesv_aa.f‎

Lines changed: 0 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -129,8 +129,6 @@
129129
*> The length of WORK. LWORK >= MAX(1,2*N,3*N-2), and for best
130130
*> performance LWORK >= max(1,N*NB), where NB is the optimal
131131
*> blocksize for ZHETRF.
132-
*> for LWORK < N, TRS will be done with Level BLAS 2
133-
*> for LWORK >= N, TRS will be done with Level BLAS 3
134132
*>
135133
*> If LWORK = -1, then a workspace query is assumed; the routine
136134
*> only calculates the optimal size of the WORK array, returns

‎TESTING/LIN/cdrvhe_aa_2stage.f‎

Lines changed: 0 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -1,7 +1,5 @@
11
*> \brief \b CDRVHE_AA_2STAGE
22
*
3-
* @precisions fortran c -> s d
4-
*
53
* =========== DOCUMENTATION ===========
64
*
75
* Online html documentation available at

0 commit comments

Comments
 (0)