PLASMA  2.4.5
PLASMA - Parallel Linear Algebra for Scalable Multi-core Architectures
 All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Macros Groups
core_clansy.c
Go to the documentation of this file.
1 
17 #include <lapacke.h>
18 #include "common.h"
19 
20 /***************************************************************************/
25 #if defined(PLASMA_HAVE_WEAK)
26 #pragma weak CORE_clansy = PCORE_clansy
27 #define CORE_clansy PCORE_clansy
28 #endif
29 void CORE_clansy(int norm, int uplo, int N,
30  PLASMA_Complex32_t *A, int LDA,
31  float *work, float *normA)
32 {
33  *normA = LAPACKE_clansy_work(
34  LAPACK_COL_MAJOR,
35  lapack_const(norm), lapack_const(uplo),
36  N, A, LDA, work);
37 }
38 
39 /***************************************************************************/
42 void QUARK_CORE_clansy(Quark *quark, Quark_Task_Flags *task_flags,
43  int norm, int uplo, int N,
44  PLASMA_Complex32_t *A, int LDA, int szeA,
45  int szeW, float *result)
46 {
47  szeW = max(1, szeW);
48  DAG_CORE_LANSY;
49  QUARK_Insert_Task(quark, CORE_clansy_quark, task_flags,
50  sizeof(PLASMA_enum), &norm, VALUE,
51  sizeof(PLASMA_enum), &uplo, VALUE,
52  sizeof(int), &N, VALUE,
53  sizeof(PLASMA_Complex32_t)*szeA, A, INPUT,
54  sizeof(int), &LDA, VALUE,
55  sizeof(float)*szeW, NULL, SCRATCH,
56  sizeof(float), result, OUTPUT,
57  0);
58 }
59 
60 /***************************************************************************/
63 #if defined(PLASMA_HAVE_WEAK)
64 #pragma weak CORE_clansy_quark = PCORE_clansy_quark
65 #define CORE_clansy_quark PCORE_clansy_quark
66 #endif
67 void CORE_clansy_quark(Quark *quark)
68 {
69  float *normA;
70  int norm;
71  int uplo;
72  int N;
73  PLASMA_Complex32_t *A;
74  int LDA;
75  float *work;
76 
77  quark_unpack_args_7(quark, normA, norm, uplo, N, A, LDA, work);
78  *normA = LAPACKE_clansy_work(
79  LAPACK_COL_MAJOR,
80  lapack_const(norm), lapack_const(uplo),
81  N, A, LDA, work);
82 }
83 
84 /***************************************************************************/
87 void QUARK_CORE_clansy_f1(Quark *quark, Quark_Task_Flags *task_flags,
88  PLASMA_enum norm, PLASMA_enum uplo, int N,
89  PLASMA_Complex32_t *A, int LDA, int szeA,
90  int szeW, float *result,
91  float *fake, int szeF)
92 {
93  szeW = max(1, szeW);
94  DAG_CORE_LANSY;
95  QUARK_Insert_Task(quark, CORE_clansy_f1_quark, task_flags,
96  sizeof(PLASMA_enum), &norm, VALUE,
97  sizeof(PLASMA_enum), &uplo, VALUE,
98  sizeof(int), &N, VALUE,
99  sizeof(PLASMA_Complex32_t)*szeA, A, INPUT,
100  sizeof(int), &LDA, VALUE,
101  sizeof(float)*szeW, NULL, SCRATCH,
102  sizeof(float), result, OUTPUT,
103  sizeof(float)*szeF, fake, OUTPUT | GATHERV,
104  0);
105 }
106 
107 /***************************************************************************/
110 #if defined(PLASMA_HAVE_WEAK)
111 #pragma weak CORE_clansy_f1_quark = PCORE_clansy_f1_quark
112 #define CORE_clansy_f1_quark PCORE_clansy_f1_quark
113 #endif
114 void CORE_clansy_f1_quark(Quark *quark)
115 {
116  float *normA;
117  int norm;
118  int uplo;
119  int N;
120  PLASMA_Complex32_t *A;
121  int LDA;
122  float *work;
123  float *fake;
124 
125  quark_unpack_args_8(quark, norm, uplo, N, A, LDA, work, normA, fake);
126  *normA = LAPACKE_clansy_work(
127  LAPACK_COL_MAJOR,
128  lapack_const(norm),
129  lapack_const(uplo),
130  N, A, LDA, work);
131 }