exaFOAM/foam-extend4.1-coherent-io
Archived
1
1
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

Preconditioning bugfix by Alexander Monakov

Browse source
This commit is contained in:
Dominik Christ 2013-01-07 17:55:06 +00:00
parent 8616d140e0
commit d25102187a
1 changed files with 3 additions and 2 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

5
src/lduSolvers/lduSolver/bicgStabSolver/bicgStabSolver.C
View file

@ -172,8 +172,9 @@ Foam::lduSolverPerformance Foam::bicgStabSolver::solve
s[i] = r[i] - alpha*v[i];
}
// Execute preconditioning transpose
preconPtr_->preconditionT(sh, s, cmpt);
// Execute preconditioning
// Bug fix, Alexander Monakov, 11/Jul/2012
preconPtr_->precondition(sh, s, cmpt);
matrix_.Amul(t, sh, coupleBouCoeffs_, interfaces_, cmpt);
omega = gSumProd(t, s)/gSumProd(t, t);