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

Bug-fix:

Browse source 
- Account for solution singularity in the CG solver
This commit is contained in:
Sandeep Menon 2011-05-19 12:31:03 -04:00
parent 2001db95f6
commit af58d024c3
1 changed files with 11 additions and 2 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

13
src/dynamicMesh/meshMotion/mesquiteMotionSolver/mesquiteMotionSolver.C
View file

@ -2886,7 +2886,7 @@ label mesquiteMotionSolver::CG
)
{
// Local variables
scalar alpha, beta, rho, rhoOld, residual;
scalar alpha, beta, rho, rhoOld, residual, wApA;
label maxIter = x.size(), iter = 0;
reduce(maxIter, sumOp<label>());
@ -2915,7 +2915,16 @@ label mesquiteMotionSolver::CG
{
A(p,w);
alpha = rho / dot(p,w);
wApA = dot(p,w);
if ((mag(wApA)/norm) < VSMALL)
{
Info<< " Solution singularity.";
return iter;
}
alpha = rho / wApA;
forAll (x, i)
{