TP-openmp/BE_OpenMP_2019/lu_tasks/lu_seq.c
2023-06-22 20:19:48 +02:00

32 lines
558 B
C

#include "trace.h"
#include "common.h"
/* This is a sequential routine for the LU factorization of a square
matrix in block-columns */
void lu_seq(Matrix A, info_type info){
int i, j;
trace_init();
for(i=0; i<info.NB; i++){
/* Do the panel */
panel(A[i], i, info);
for(j=i+1; j<info.NB; j++){
/* Do all the correspondint updates */
update(A[i], A[j], i, j, info);
}
}
/* Do row permutations resulting from the numerical pivoting */
backperm(A, info);
trace_dump("trace_seq.svg");
return;
}