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
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
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
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()
149                            A[j * n + i]);  in check_if_valid()
152                            A[j * n + i] - epsilon,  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
173                 double ***tile = create_tile_array(A, n, b);  in operator ()()
181                 memcpy(A0, A, sizeof(double) * n * n);  in operator ()()
190             if (!g_benchmark_run && !check_if_valid(A0, C, A, n)) {  in operator ()()
293         double *A = (double *)ptr;  in func()  local
297         dpotrf(&uplo, &n, A, &lda, &info);  in func()
690     double *A = nullptr;  in main()  local
693     matrix_init(A, g_n, g_input_file_name);  in main()
699         matrix_write(A, g_n, s.c_str());  in main()
705             (*alg)(A, g_n, g_b, g_num_trials);  in main()
713             (*alg)(A, g_n, g_b, g_num_trials);  in main()
721     free(A);  in main()