#include <stdlib.h>#include <stdio.h>#include <string.h>#include <math.h>#include <stdbool.h>#include <assert.h>#include <omp.h>#include "../lah_intf.h"#include "./lah_crs.h"#include "lin_alg_intf.h"
Functions | |
| int | lar_allocate_SM_and_LV_crs (int Nrdof_glob, int Max_SM_size, int *Nroffbl, int **L_offbl) |
| int | lar_create_solver_structures_accel_crs (int Matrix_id) |
| int | lar_get_crs_data_crs (int Matrix_id, int *Nrdof_glob_p, int *Nnz_p, int **Crs_col_ind_p, int **Crs_row_ptr_p, double **Crs_val_p, double **Rhs_p) |
| int | lar_initialize_SM_and_LV_crs (int Matrix_id, int Scope) |
| double | lar_get_storage_crs (int Matrix_id) |
| int | lar_fill_assembly_table_int_ent_crs (int Matrix_id, int Nr_dof_bl, int *L_bl_id, int *L_bl_nrdof, int *Assembly_table_int_ent) |
| int | lar_assemble_SM_and_LV_with_table_crs (int Matrix_id, int Scope, int Nr_dof_bl, int *Assembly_table_int_ent, int *Local_to_global_int_ent, double *Stiff_mat, double *Rhs_vect, char *Rewr_dofs) |
| int | lar_assemble_SM_and_LV_crs (int Matrix_id, int Scope, int Nr_dof_bl, int *L_bl_id, int *L_bl_nrdof, double *Stiff_mat, double *Rhs_vect, char *Rewr_dofs) |
| int | lar_allocate_preconditioner_crs (int Matrix_id, int Precon, int ILU_k) |
| int | lar_fill_preconditioner_crs (int Matrix_id) |
| int | lar_free_preconditioner_crs (int Matrix_id) |
| int | lar_free_SM_and_LV_crs (int Matrix_id) |
| void | lar_compute_residual_crs (int Matrix_id, int Use_rhs, int Ini_zero, int Ndof, double *X, double *B, double *V) |
| void | lar_perform_BJ_or_GS_iterations_crs (int Matrix_id, int Use_rhs, int Ini_zero, int Nr_prec, int Ndof, double *V, double *B) |
| void | lar_perform_rhsub_crs (int Matrix_id, int Ndof, double *V, double *B) |
| int | lar_get_SM_and_LV_crs_from_crs (int Matrix_id, int offset, int **crs_row, int **crs_col, double **crs_val, double **rhs) |
| int | lar_get_SM_and_LV_coo_from_crs (int Matrix_id, int *nnz, int offset, int **coo_row, int **coo_col, double **coo_val, double **rhs) |
Variables | |
| int | itv_nr_crs_matrices = 0 |
| itt_crs_matrices | itv_crs_matrices [LAC_MAX_MATRICES] |
| int lar_allocate_preconditioner_crs | ( | int | Matrix_id, | |
| int | Precon, | |||
| int | ILU_k | |||
| ) |
-------------------------------------------------------- lar_allocate_preconditioner_crs - to allocate space for preconditioner ---------------------------------------------------------

| int lar_allocate_SM_and_LV_crs | ( | int | Nrdof_glob, | |
| int | Max_SM_size, | |||
| int * | Nroffbl, | |||
| int ** | L_offbl | |||
| ) |


| int lar_assemble_SM_and_LV_crs | ( | int | Matrix_id, | |
| int | Scope, | |||
| int | Nr_dof_bl, | |||
| int * | L_bl_id, | |||
| int * | L_bl_nrdof, | |||
| double * | Stiff_mat, | |||
| double * | Rhs_vect, | |||
| char * | Rewr_dofs | |||
| ) |

| int lar_assemble_SM_and_LV_with_table_crs | ( | int | Matrix_id, | |
| int | Scope, | |||
| int | Nr_dof_bl, | |||
| int * | Assembly_table_int_ent, | |||
| int * | Local_to_global_int_ent, | |||
| double * | Stiff_mat, | |||
| double * | Rhs_vect, | |||
| char * | Rewr_dofs | |||
| ) |

| void lar_compute_residual_crs | ( | int | Matrix_id, | |
| int | Use_rhs, | |||
| int | Ini_zero, | |||
| int | Ndof, | |||
| double * | X, | |||
| double * | B, | |||
| double * | V | |||
| ) |

| int lar_create_solver_structures_accel_crs | ( | int | Matrix_id | ) |
--------------------------------------------------------- lar_create_solver_structures_accel_crs - utility to create data structures on GPU ---------------------------------------------------------


| int lar_fill_assembly_table_int_ent_crs | ( | int | Matrix_id, | |
| int | Nr_dof_bl, | |||
| int * | L_bl_id, | |||
| int * | L_bl_nrdof, | |||
| int * | Assembly_table_int_ent | |||
| ) |
----------------------------------------------------------- lar_fill_assembly_table_int_ent_crs - to fill a part of the global assembly table related to one integration entity, for which lists of DOF blocks (their global positions) are provided ------------------------------------------------------------

| int lar_fill_preconditioner_crs | ( | int | Matrix_id | ) |
-------------------------------------------------------- lar_fill_preconditioner_crs - to fill preconditioner ---------------------------------------------------------

| int lar_free_preconditioner_crs | ( | int | Matrix_id | ) |
-------------------------------------------------------- lar_free_preconditioner_crs - to free space for a block structure ---------------------------------------------------------

| int lar_free_SM_and_LV_crs | ( | int | Matrix_id | ) |

| int lar_get_crs_data_crs | ( | int | Matrix_id, | |
| int * | Nrdof_glob_p, | |||
| int * | Nnz_p, | |||
| int ** | Crs_col_ind_p, | |||
| int ** | Crs_row_ptr_p, | |||
| double ** | Crs_val_p, | |||
| double ** | Rhs_p | |||
| ) |
--------------------------------------------------------- lar_get_crs_data_crs - utility to get CRS parameters and pointers ---------------------------------------------------------

| int lar_get_SM_and_LV_coo_from_crs | ( | int | Matrix_id, | |
| int * | nnz, | |||
| int | offset, | |||
| int ** | coo_row, | |||
| int ** | coo_col, | |||
| double ** | coo_val, | |||
| double ** | rhs | |||
| ) |

| int lar_get_SM_and_LV_crs_from_crs | ( | int | Matrix_id, | |
| int | offset, | |||
| int ** | crs_row, | |||
| int ** | crs_col, | |||
| double ** | crs_val, | |||
| double ** | rhs | |||
| ) |

| double lar_get_storage_crs | ( | int | Matrix_id | ) |
-------------------------------------------------------- lar_get_storage_crs - to compute storage of SM, LV and preconditioner ---------------------------------------------------------

| int lar_initialize_SM_and_LV_crs | ( | int | Matrix_id, | |
| int | Scope | |||
| ) |
-------------------------------------------------------- lar_initialize_SM_and_LV_crs - to initialize stiffness matrix and/or load vector ---------------------------------------------------------

| void lar_perform_BJ_or_GS_iterations_crs | ( | int | Matrix_id, | |
| int | Use_rhs, | |||
| int | Ini_zero, | |||
| int | Nr_prec, | |||
| int | Ndof, | |||
| double * | V, | |||
| double * | B | |||
| ) |
-------------------------------------------------------- lar_perform_BJ_or_GS_iterations_crs - to perform one iteration of block Gauss-Seidel or block Jacobi algorithm: v_out = v_in + M^-1 * ( b - A * v_in ) where M^-1 results from stored preconditioner matrix and the algorithm ---------------------------------------------------------

| void lar_perform_rhsub_crs | ( | int | Matrix_id, | |
| int | Ndof, | |||
| double * | V, | |||
| double * | B | |||
| ) |
-------------------------------------------------------- lar_perform_rhsub_crs - to perform forward reduction and back-substitution for ILU preconditioning ---------------------------------------------------------

| itt_crs_matrices itv_crs_matrices[LAC_MAX_MATRICES] |
| int itv_nr_crs_matrices = 0 |
1.6.1