50 void matrix_init(double *&A, int &n, const char *fname);
57 void matrix_write(double *A, int n, const char *fname, bool is_triangular = false);
70 static double ***create_tile_array(double *A, int n, int b) {  in create_tile_array()  argument
71     const int p = n / b;  in create_tile_array()
84                     temp_block[T_j * b + T_i] = A[A_j * n + A_i];  in create_tile_array()
94 static void collapse_tile_array(double ***tile, double *A, int n, int b) {  in collapse_tile_array()  argument
95     const int p = n / b;  in collapse_tile_array()
103                     A[A_j * n + A_i] = temp_block[T_j * b + T_i];  in collapse_tile_array()
124     bool check_if_valid(double *A0, double *C, double *A, int n) {  in check_if_valid()  argument
129         for (int i = 0; i < n; ++i) {  in check_if_valid()
130             for (int j = i + 1; j < n; ++j) {  in check_if_valid()
131                 A0[j * n + i] = 0.;  in check_if_valid()
135         dgemm(&transa, &transb, &n, &n, &n, &alpha, A0, &n, A0, &n, &beta, C, &n);  in check_if_valid()
137         for (int j = 0; j < n; ++j) {  in check_if_valid()
138             for (int i = 0; i < n; ++i) {  in check_if_valid()
139                 const double epsilon = std::abs(A[j * n + i] * 0.1);  in check_if_valid()
141                 if (std::abs(C[j * n + i] - A[j * n + i]) > epsilon) {  in check_if_valid()
146                            C[j * n + i],  in check_if_valid()
149                            A[j * n + i]);  in check_if_valid()
152                            A[j * n + i] - epsilon,  in check_if_valid()
153                            C[j * n + i],  in check_if_valid()
154                            A[j * n + i] + epsilon);  in check_if_valid()
165     double operator()(double *A, int n, int b, int trials) {  in operator ()()  argument
168         double *A0 = (double *)calloc(sizeof(double), n * n);  in operator ()()
169         double *C = (double *)calloc(sizeof(double), n * n);  in operator ()()
173                 double ***tile = create_tile_array(A, n, b);  in operator ()()
175                 func(tile, n, b);  in operator ()()
178                 collapse_tile_array(tile, A0, n, b);  in operator ()()
181                 memcpy(A0, A, sizeof(double) * n * n);  in operator ()()
183                 func(A0, n, b);  in operator ()()
190             if (!g_benchmark_run && !check_if_valid(A0, C, A, n)) {  in operator ()()
212                n,  in operator ()()
223     virtual void func(void *ptr, int n, int b) = 0;
269     virtual void func(void *ptr, int n, int b) {  in func()  argument
271         const int p = n / b;  in func()
292     virtual void func(void *ptr, int n, int /* b */) {  in func()  argument
294         int lda = n;  in func()
297         dpotrf(&uplo, &n, A, &lda, &info);  in func()
469     virtual void func(void *ptr, int n, int b) {  in func()  argument
475         const int p = n / b;  in func()
545     virtual void func(void *ptr, int n, int b) {  in func()  argument
548         const int p = n / b;  in func()