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

added check for rhs = 0 , if rhs = 0 the routine sets x = 0 and returns

Browse Source
This commit is contained in:
ulrikey 1999-05-07 18:25:34 +00:00
parent 4bf083aba1
commit da359d1c62
2 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
parcsr_linear_solvers/gmres.c
View File

@ -224,7 +224,7 @@ hypre_GMRESSolve(void *gmres_vdata,
/* initialize first term of hessenberg system */
rs[0] = r_norm;
if (r_norm == 0.0)
if (r_norm == 0.0 || b_norm == 0)
{
hypre_PCGCopyVector(b,x);
ierr = 0;

2
parcsr_ls/gmres.c
View File

@ -224,7 +224,7 @@ hypre_GMRESSolve(void *gmres_vdata,
/* initialize first term of hessenberg system */
rs[0] = r_norm;
if (r_norm == 0.0)
if (r_norm == 0.0 || b_norm == 0)
{
hypre_PCGCopyVector(b,x);
ierr = 0;