zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

COMP: avoid ambiguous construct from tmp - randomProcesses

Browse Source
This commit is contained in:
Mark Olesen
2010-12-17 17:12:38 +01:00
parent 8de19df4a7
commit f4a5bea2b9
2 changed files with 9 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files

9
src/randomProcesses/noise/noiseFFT.C
View File

@ -34,6 +34,7 @@ License
Foam::scalar Foam::noiseFFT::p0 = 2e-5;
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
// Construct from pressure field
@ -226,9 +227,9 @@ Foam::graph Foam::noiseFFT::meanPf
scalarField MeanPf(N/2, 0.0);
scalarField Hwf = Hanning(N);
scalarField Hwf(Hanning(N));
for (label wi=0; wi<nw; wi++)
for (label wi=0; wi<nw; ++wi)
{
MeanPf += Pf(Hwf*window(N, wi));
}
@ -272,9 +273,9 @@ Foam::graph Foam::noiseFFT::RMSmeanPf
scalarField RMSMeanPf(N/2, 0.0);
scalarField Hwf = Hanning(N);
scalarField Hwf(Hanning(N));
for (label wi=0; wi<nw; wi++)
for (label wi=0; wi<nw; ++wi)
{
RMSMeanPf += sqr(Pf(Hwf*window(N, wi)));
}

4
src/randomProcesses/turbulence/turbGen.C
View File

@ -60,12 +60,14 @@ Foam::vectorField Foam::turbGen::U()
s = Ek(Ea, k0, mag(K))*s;
complexVectorField up =
complexVectorField up
(
fft::reverseTransform
(
ComplexField(cos(constant::mathematical::twoPi*rndPhases)*s,
sin(constant::mathematical::twoPi*rndPhases)*s),
K.nn()
)
);
return ReImSum(up);