CFD/hypre
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

bug fix; nonsquare rap (#581)

Browse Source 
This PR fixes a corner case of the RAP routine for RAP matrix that is globally square but not locally.
This commit is contained in:
Ruipeng Li 2022-03-03 21:26:17 -08:00 committed by GitHub
parent 04af9a4cd9
commit ebd6eb88c3
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
1 changed files with 2 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
src/parcsr_ls/par_rap.c
View File

@ -234,7 +234,8 @@ hypre_BoomerAMGBuildCoarseOperatorKT( hypre_ParCSRMatrix *RT,
num_nz_cols_A = num_cols_diag_A + num_cols_offd_A;
n_coarse_RT = hypre_ParCSRMatrixGlobalNumCols(RT);
if (n_coarse != n_coarse_RT)
if (n_coarse != n_coarse_RT || num_cols_diag_RT != num_cols_diag_P)
{
square = 0;
}