songtianlun/openfoam
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

tidying

Browse Source
This commit is contained in:
andy 2009-03-26 11:43:25 +00:00
parent 284f5d97f0
commit d504063ce1
1 changed files with 4 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
src/thermophysicalModels/radiation/radiationModel/fvDOM/fvDOM/fvDOM.C
View File

@ -132,8 +132,8 @@ Foam::radiation::fvDOM::fvDOM(const volScalarField& T)
{ {
if (mesh_.nSolutionD() == 3) //3D if (mesh_.nSolutionD() == 3) //3D
{ {
IRay_.setSize(4*nPhi_*nTheta_);
nRay_ = 4*nPhi_*nTheta_; nRay_ = 4*nPhi_*nTheta_;
IRay_.setSize(nRay_);
scalar deltaPhi = mathematicalConstant::pi/(2.0*nPhi_); scalar deltaPhi = mathematicalConstant::pi/(2.0*nPhi_);
scalar deltaTheta = mathematicalConstant::pi/nTheta_; scalar deltaTheta = mathematicalConstant::pi/nTheta_;
label i = 0; label i = 0;
@ -169,8 +169,8 @@ Foam::radiation::fvDOM::fvDOM(const volScalarField& T)
{ {
scalar thetai = mathematicalConstant::pi/2.0; scalar thetai = mathematicalConstant::pi/2.0;
scalar deltaTheta = mathematicalConstant::pi; scalar deltaTheta = mathematicalConstant::pi;
IRay_.setSize(4*nPhi_);
nRay_ = 4*nPhi_; nRay_ = 4*nPhi_;
IRay_.setSize(nRay_);
scalar deltaPhi = mathematicalConstant::pi/(2.0*nPhi_); scalar deltaPhi = mathematicalConstant::pi/(2.0*nPhi_);
label i = 0; label i = 0;
for (label m = 1; m <= 4*nPhi_; m++) for (label m = 1; m <= 4*nPhi_; m++)
@ -199,8 +199,8 @@ Foam::radiation::fvDOM::fvDOM(const volScalarField& T)
{ {
scalar thetai = mathematicalConstant::pi/2.0; scalar thetai = mathematicalConstant::pi/2.0;
scalar deltaTheta = mathematicalConstant::pi; scalar deltaTheta = mathematicalConstant::pi;
IRay_.setSize(2);
nRay_ = 2; nRay_ = 2;
IRay_.setSize(nRay_);
scalar deltaPhi = mathematicalConstant::pi; scalar deltaPhi = mathematicalConstant::pi;
label i = 0; label i = 0;
for (label m = 1; m <= 2; m++) for (label m = 1; m <= 2; m++)
@ -296,7 +296,7 @@ void Foam::radiation::fvDOM::calculate()
maxResidual = max(maxBandResidual, maxResidual); maxResidual = max(maxBandResidual, maxResidual);
} }
Info << "Radiation solver iter: " << radIter << endl; Info << "Radiation solver iter: " << radIter << endl;
} while(maxResidual > convergence_); } while(maxResidual > convergence_);