40 for(i = 0; i < rows; i++) {
41 for(j = 0; j < bandwidth; j++) {
86 for(i = 0; i < rows; i++) {
87 for(j = 0; j < bandwidth; j++) {
96 for(i = 0; i < rows; i++) {
97 for(j = i; j < rows; j++) {
129 #pragma omp for schedule (static) private(i, j, tmp) 130 for (i = 0; i < rows; i++) {
132 for (j = 0; j < bandwidth; j++) {
134 tmp += A[i][j]*x[i + j];
141 for (i = 0; i < rows; i++) {
143 for (j = 1; j < bandwidth; j++) {
145 y[i + j] += A[i][j]*x[i];
double ** G_math_matrix_to_sband_matrix(double **A, int rows, int bandwidth)
Convert a symmetrix matrix into a symmetric band matrix.
double ** G_alloc_matrix(int, int)
Matrix memory allocation.
void G_math_Ax_sband(double **A, double *x, double *y, int rows, int bandwidth)
Compute the matrix - vector product of symmetric band matrix A and vector x.
double ** G_math_sband_matrix_to_matrix(double **A, int rows, int bandwidth)
Convert a symmetric band matrix into a symmetric matrix.