FreeFOAM The Cross-Platform CFD Toolkit
GAMGSolverScalingFactor.C
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration |
5  \\ / A nd | Copyright (C) 1991-2010 OpenCFD Ltd.
6  \\/ M anipulation |
7 -------------------------------------------------------------------------------
8 License
9  This file is part of OpenFOAM.
10 
11  OpenFOAM is free software: you can redistribute it and/or modify it
12  under the terms of the GNU General Public License as published by
13  the Free Software Foundation, either version 3 of the License, or
14  (at your option) any later version.
15 
16  OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
17  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
18  FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
19  for more details.
20 
21  You should have received a copy of the GNU General Public License
22  along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
23 
24 \*---------------------------------------------------------------------------*/
25 
26 #include "GAMGSolver.H"
27 #include <OpenFOAM/vector2D.H>
28 
29 // * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * //
30 
31 Foam::scalar Foam::GAMGSolver::scalingFactor
32 (
33  scalarField& field,
34  const scalarField& source,
35  const scalarField& Acf,
36  const scalarField& D
37 ) const
38 {
39  scalar scalingFactorNum = 0.0;
40  scalar scalingFactorDenom = 0.0;
41 
42  forAll(field, i)
43  {
44  scalingFactorNum += source[i]*field[i];
45  scalingFactorDenom += Acf[i]*field[i];
46 
47  // While the matrix-multiply done for the scaling it is
48  // possible to perform a point-Jacobi smoothing operation cheaply
49  field[i] += (source[i] - Acf[i])/D[i];
50  }
51 
52  vector2D scalingVector(scalingFactorNum, scalingFactorDenom);
53  reduce(scalingVector, sumOp<vector2D>());
54  return scalingVector.x()/stabilise(scalingVector.y(), VSMALL);
55 }
56 
57 
58 Foam::scalar Foam::GAMGSolver::scalingFactor
59 (
60  scalarField& Acf,
61  const lduMatrix& A,
62  scalarField& field,
63  const FieldField<Field, scalar>& interfaceLevelBouCoeffs,
64  const lduInterfaceFieldPtrsList& interfaceLevel,
65  const scalarField& source,
66  const direction cmpt
67 ) const
68 {
69  A.Amul
70  (
71  Acf,
72  field,
73  interfaceLevelBouCoeffs,
74  interfaceLevel,
75  cmpt
76  );
77 
78  return scalingFactor
79  (
80  field,
81  source,
82  Acf,
83  A.diag()
84  );
85 }
86 
87 
88 // ************************ vim: set sw=4 sts=4 et: ************************ //