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

Fix chebyshev smoother for singular problems (#657)

Browse Source 
See PR's description for additional info
This commit is contained in:
Victor A. Paludetto Magri 2022-06-22 17:47:09 -07:00 committed by GitHub
parent b58585e0f0
commit 850fd47d07
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
1 changed files with 7 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
src/parcsr_ls/par_cheby.c
View File

@ -68,7 +68,7 @@ hypre_ParCSRRelax_Cheby_Setup(hypre_ParCSRMatrix *A, /* matrix to relax
hypre_CSRMatrix *A_diag = hypre_ParCSRMatrixDiag(A); hypre_CSRMatrix *A_diag = hypre_ParCSRMatrixDiag(A);
HYPRE_Real theta, delta; HYPRE_Real theta, delta;
HYPRE_Real den; HYPRE_Real den;
HYPRE_Real upper_bound = 0.0, lower_bound = 0.0; HYPRE_Real upper_bound, lower_bound;
HYPRE_Int num_rows = hypre_CSRMatrixNumRows(A_diag); HYPRE_Int num_rows = hypre_CSRMatrixNumRows(A_diag);
HYPRE_Real *coefs = NULL; HYPRE_Real *coefs = NULL;
HYPRE_Int cheby_order; HYPRE_Int cheby_order;
@ -89,18 +89,18 @@ hypre_ParCSRRelax_Cheby_Setup(hypre_ParCSRMatrix *A, /* matrix to relax
/* we are using the order of p(A) */ /* we are using the order of p(A) */
cheby_order = order - 1; cheby_order = order - 1;
if (min_eig >= 0.0) if (max_eig <= 0.0)
{
upper_bound = min_eig * 1.1;
lower_bound = max_eig - (max_eig - upper_bound) * fraction;
}
else
{ {
/* make sure we are large enough - Adams et al. 2003 */ /* make sure we are large enough - Adams et al. 2003 */
upper_bound = max_eig * 1.1; upper_bound = max_eig * 1.1;
/* lower_bound = max_eig/fraction; */ /* lower_bound = max_eig/fraction; */
lower_bound = (upper_bound - min_eig) * fraction + min_eig; lower_bound = (upper_bound - min_eig) * fraction + min_eig;
} }
else if (max_eig <= 0.0)
{
upper_bound = min_eig * 1.1;
lower_bound = max_eig - (max_eig - upper_bound) * fraction;
}
/* theta and delta */ /* theta and delta */
theta = (upper_bound + lower_bound) / 2; theta = (upper_bound + lower_bound) / 2;