Commit 540d53b0 authored by BOUCHERIE Raphael's avatar BOUCHERIE Raphael

test for zgels_param pass

parent a58864d1
......@@ -157,9 +157,9 @@ void morse_pzgeqrf_param( const libhqr_tree_t *qrtree, MORSE_desc_t *A, MORSE_de
m = tiles[j];
p = qrtree->currpiv(qrtree, k, m);
tempmm == A->mt-1 ? A->m-m*A->mb : A->mb;
ldam = BLKLDD(A, m);
tempmm = m == A->mt-1 ? A->m-m*A->mb : A->mb;
ldap = BLKLDD(A, p);
ldam = BLKLDD(A, m);
/* Tiles killed is a TS */
if(qrtree->gettype(qrtree, k, m) == 0){
......
......@@ -58,6 +58,7 @@
#include "coreblas/include/coreblas_s.h"
#include "coreblas/include/coreblas_zc.h"
#include "coreblas/include/coreblas_ds.h"
#include <assert.h>
#ifdef __cplusplus
extern "C" {
......@@ -66,8 +67,10 @@ extern "C" {
/** ****************************************************************************
* Coreblas Error
**/
#define coreblas_error(k, str) fprintf(stderr, "%s: Parameter %d / %s\n", __func__, k, str)
#define coreblas_error(k, str) do { \
fprintf(stderr, "%s: Parameter %d / %s\n", __func__, k, str) ; \
assert(0); \
} while(0)
/** ****************************************************************************
* CBlas enum
**/
......
......@@ -31,6 +31,8 @@
#include "runtime/quark/include/quark_blas.h"
#include "runtime/quark/include/core_blas_dag.h"
#define QUARK_Insert_Task QUARK_Execute_Task
#include "control/common.h"
typedef struct quark_option_s {
......
......@@ -64,7 +64,7 @@ set(ZSRC
# LAPACK
##################
testing_zgels.c
testing_zgeqrf_param.c
testing_zgels_param.c
#testing_zgesv.c
testing_zgesv_incpiv.c
#testing_zgetri.c
......@@ -141,6 +141,7 @@ if(NOT CHAMELEON_SIMULATION)
${CBLAS_LIBRARIES}
${BLAS_LIBRARIES}
${HWLOC_LIBRARIES}
${LIBHQR_LIBRARIES}
${EXTRA_LIBRARIES}
)
......
......@@ -248,6 +248,9 @@ int main (int argc, char **argv)
else if ( strcmp(func, "GESV_INCPIV") == 0 ) {
info += testing_zgesv_incpiv( argc, argv );
}
else if ( strcmp(func, "GELS_PARAM") == 0 ) {
info += testing_zgels_param( argc, argv );
}
/* else if ( strcmp(func, "GESV") == 0 ) { */
/* info += testing_zgesv( argc, argv ); */
/* } */
......
......@@ -92,6 +92,7 @@ int testing_zgeadd(int argc, char **argv);
int testing_zposv(int argc, char **argv);
int testing_zgels(int argc, char **argv);
int testing_zgels_param(int argc, char **argv);
int testing_zgesv(int argc, char **argv);
int testing_zgesv_incpiv(int argc, char **argv);
......
......@@ -280,16 +280,12 @@ int testing_zgels(int argc, char **argv)
/* Initialize A1 and A2 */
LAPACKE_zlarnv_work(IONE, ISEED, LDAxN, A1);
for (i = 0; i < M; i++)
for (j = 0; j < N; j++)
A2[LDA*j+i] = A1[LDA*j+i] ;
LAPACKE_zlacpy_work(LAPACK_COL_MAJOR, 'A', M, N, A1, LDA, A2, LDA );
/* Initialize B1 and B2 */
memset(B2, 0, LDB*NRHS*sizeof(MORSE_Complex64_t));
LAPACKE_zlarnv_work(IONE, ISEED, LDBxNRHS, B1);
for (i = 0; i < M; i++)
for (j = 0; j < NRHS; j++)
B2[LDB*j+i] = B1[LDB*j+i] ;
LAPACKE_zlacpy_work(LAPACK_COL_MAJOR, 'A', M, NRHS, B1, LDB, B2, LDB );
/* MORSE ZGELS */
MORSE_zgels(MorseNoTrans, M, N, NRHS, A2, LDA, T, B2, LDB);
......
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment