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

Reversed order of update does not work for scalar matrix. Reverting

Browse source
This commit is contained in:
Hrvoje Jasak 2016-03-18 10:56:42 +00:00
parent 501d28c568
commit e5e0ac01c3
1 changed files with 8 additions and 14 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

22
src/foam/matrices/lduMatrix/lduMatrix/lduMatrixATmul.C
View file

@ -44,13 +44,6 @@ void Foam::lduMatrix::Amul
// HJ, 5/Nov/2007
Ax = 0;
// AmulCore must be additive to account for initialisation step
// in ldu interfaces. HJ, 6/Nov/2007
AmulCore(Ax, x);
// Note: changed order of interface update: init after core
// HJ, 14/Mar/2016
// Initialise the update of coupled interfaces
initMatrixInterfaces
(
@ -61,6 +54,10 @@ void Foam::lduMatrix::Amul
cmpt
);
// AmulCore must be additive to account for initialisation step
// in ldu interfaces. HJ, 6/Nov/2007
AmulCore(Ax, x);
// Update coupled interfaces
updateMatrixInterfaces
(
@ -132,13 +129,6 @@ void Foam::lduMatrix::Tmul
// HJ, 5/Nov/2007
Tx = 0;
// TmulCore must be additive to account for initialisation step
// in ldu interfaces. HJ, 6/Nov/2007
TmulCore(Tx, x);
// Note: changed order of interface update: init after core
// HJ, 14/Mar/2016
// Initialise the update of coupled interfaces
initMatrixInterfaces
(
@ -149,6 +139,10 @@ void Foam::lduMatrix::Tmul
cmpt
);
// TmulCore must be additive to account for initialisation step
// in ldu interfaces. HJ, 6/Nov/2007
TmulCore(Tx, x);
// Update coupled interfaces
updateMatrixInterfaces
(