Mentions légales du service

Skip to content
Snippets Groups Projects
Commit 9224a468 authored by Mathieu Faverge's avatar Mathieu Faverge
Browse files

Add compilation of the kernels

parent 3b69d027
No related branches found
No related tags found
1 merge request!3Add ztpqrt, and ztpgqrt functions
...@@ -73,6 +73,8 @@ set(ZSRC ...@@ -73,6 +73,8 @@ set(ZSRC
core_zsyrk.c core_zsyrk.c
core_zsyssq.c core_zsyssq.c
core_zsytf2_nopiv.c core_zsytf2_nopiv.c
core_ztpqrt.c
core_ztpmqrt.c
core_ztradd.c core_ztradd.c
core_ztrasm.c core_ztrasm.c
core_ztrmm.c core_ztrmm.c
......
...@@ -148,8 +148,7 @@ int CORE_ztpmqrt( MORSE_enum side, MORSE_enum trans, ...@@ -148,8 +148,7 @@ int CORE_ztpmqrt( MORSE_enum side, MORSE_enum trans,
MORSE_Complex64_t *B, int LDB, MORSE_Complex64_t *B, int LDB,
MORSE_Complex64_t *WORK ) MORSE_Complex64_t *WORK )
{ {
int m1; int m1, n1, ldwork;
int n1;
/* Check input arguments */ /* Check input arguments */
if ((side != MorseLeft) && (side != MorseRight)) { if ((side != MorseLeft) && (side != MorseRight)) {
...@@ -160,7 +159,7 @@ int CORE_ztpmqrt( MORSE_enum side, MORSE_enum trans, ...@@ -160,7 +159,7 @@ int CORE_ztpmqrt( MORSE_enum side, MORSE_enum trans,
if ( side == MorseLeft ) { if ( side == MorseLeft ) {
m1 = K; m1 = K;
n1 = N; n1 = N;
ldwork = ib; ldwork = IB;
} }
else { else {
m1 = M; m1 = M;
......
...@@ -102,12 +102,6 @@ int CORE_ztpqrt( int M, int N, int L, int IB, ...@@ -102,12 +102,6 @@ int CORE_ztpqrt( int M, int N, int L, int IB,
MORSE_Complex64_t *T, int LDT, MORSE_Complex64_t *T, int LDT,
MORSE_Complex64_t *WORK ) MORSE_Complex64_t *WORK )
{ {
static MORSE_Complex64_t zone = 1.0;
static MORSE_Complex64_t zzero = 0.0;
MORSE_Complex64_t alpha;
int i, ii, sb;
#if !defined(NDEBUG) #if !defined(NDEBUG)
/* Check input arguments */ /* Check input arguments */
if (M < 0) { if (M < 0) {
...@@ -144,7 +138,7 @@ int CORE_ztpqrt( int M, int N, int L, int IB, ...@@ -144,7 +138,7 @@ int CORE_ztpqrt( int M, int N, int L, int IB,
if ((M == 0) || (N == 0) || (IB == 0)) if ((M == 0) || (N == 0) || (IB == 0))
return MORSE_SUCCESS; return MORSE_SUCCESS;
if ( L == O ) { if ( L == 0 ) {
CORE_ztsqrt( M, N, IB, A, LDA, B, LDB, T, LDT, WORK, WORK+N ); CORE_ztsqrt( M, N, IB, A, LDA, B, LDB, T, LDT, WORK, WORK+N );
} }
else if (L == M) { else if (L == M) {
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment