Comment on initialisation of xCorr in ILU smoother

This commit is contained in:
Hrvoje Jasak 2017-06-16 12:47:28 +01:00
parent 7c4511894b
commit 058c2b425d

View file

@ -121,6 +121,12 @@ public:
residual_[i] = b[i] - residual_[i];
}
// Note.
// For the ILU smoother it is not necessary to initialise
// xCorr to zero, because it gets over-written by the product
// but for more complex smoothers it is
// HJ, 15/Jun/2017
// xCorr_ = pTraits<Type>::zero;
precon_.precondition(xCorr_, residual_);
// Add correction to x