24 static void posdef_gen(double *A, int n) {  in posdef_gen()  argument
32     memset(A, 0, sizeof(double) * n * n);  in posdef_gen()
51     cblas_dgemm(CblasColMajor, CblasNoTrans, CblasNoTrans, n, n, n, 1, L, n, LT, n, 0, A, n);  in posdef_gen()
58 void matrix_init(double *&A, int &n, const char *fname) {  in matrix_init()  argument
73         A = (double *)calloc(sizeof(double), n * n);  in matrix_init()
76                 if (fscanf(fp, "%lf ", &A[i * n + j]) <= 0) {  in matrix_init()
81                     A[j * n + i] = A[i * n + j];  in matrix_init()
88         A = (double *)calloc(sizeof(double), n * n);  in matrix_init()
89         posdef_gen(A, n);  in matrix_init()
94 void matrix_write(double *A, int n, const char *fname, bool is_triangular = false) {  in matrix_write()  argument
108                 fprintf(fp, "%lf ", A[j * n + i]);  in matrix_write()
112                     fprintf(fp, "%lf ", A[i * n + j]);  in matrix_write()
129                     fprintf(fp, "%lf ", A[i * n + j]);  in matrix_write()