1 // CLASS_GETRF_NOPIVOT
2 //
call_routines(struct RoutineArgs * args)3 void call_routines(struct RoutineArgs *args){
4
5 // unpack args
6
7 // routine call
8 //
9 BLASFEO(ROUTINE)(
10 args->m, args->n,
11 args->sA_po, args->ai, args->aj,
12 args->sD, args->di, args->dj
13 );
14
15 REF(BLASFEO(ROUTINE))(
16 args->m, args->n,
17 args->rA_po, args->ai, args->aj,
18 args->rD, args->di, args->dj,
19 );
20
21 }
22
print_routine(struct RoutineArgs * args)23 void print_routine(struct RoutineArgs *args){
24 // unpack args
25
26 printf("%s\n", string(ROUTINE));
27 printf(
28 "Solving A[%d:%d,%d:%d] = P * LU[%d:%d,%d:%d]\n",
29 args->ai, args->m, args->aj, args->n,
30 args->di, args->m, args->dj, args->n
31 );
32
33 }
34
print_routine_matrices(struct RoutineArgs * args)35 void print_routine_matrices(struct RoutineArgs *args)
36 {
37 printf("\nPrint A:\n");
38 blasfeo_print_xmat_debug(args->m, args->n, args->sA_po, args->ai, args->aj, 0, 0, 0);
39 print_xmat_debug(args->m, args->n, args->rA_po, args->ai, args->aj, 0, 0, 0);
40
41 printf("\nPrint LU:\n");
42 blasfeo_print_xmat_debug(args->m, args->n, args->sD, args->ai, args->aj, 0, 0, 0);
43 print_xmat_debug(args->m, args->n, args->rD, args->ai, args->aj, 0, 0, 0);
44 }
45
46
set_test_args(struct TestArgs * targs)47 void set_test_args(struct TestArgs *targs)
48 {
49 targs->nis = 9;
50 targs->njs = 9;
51 targs->nks = 9;
52
53 targs->alphas = 1;
54 }
55