64 RootDeltaT(
sqrt(DeltaT)),
67 Alpha(UOdict.get<scalar>(
"UOalpha")),
68 Sigma(UOdict.get<scalar>(
"UOsigma")),
69 Kupper(UOdict.get<scalar>(
"UOKupper")),
70 Klower(UOdict.get<scalar>(
"UOKlower")),
71 Scale((Kupper - Klower)*
pow(scalar(Mesh.size()), 1.0/
vector::dim))
75 scalar sqrKupper =
sqr(Kupper);
76 scalar sqrKlower =
sqr(Klower) + SMALL;
81 if ((sqrK =
magSqr(
K[i])) < sqrKupper && sqrK > sqrKlower)
83 UOfield[i] = Scale*Sigma*WeinerProcess();
100 scalar sqrKupper =
sqr(Kupper);
101 scalar sqrKlower =
sqr(Klower) + SMALL;
106 if ((sqrK =
magSqr(
K[i])) < sqrKupper && sqrK > sqrKlower)
110 (1.0 - Alpha*DeltaT)*UOfield[i]
111 + Scale*Sigma*WeinerProcess();
115 Info<<
" Number of forced K = " << count <<
nl;
CGAL::Exact_predicates_exact_constructions_kernel K
Calculate the wavenumber vector field corresponding to the space vector field of a finite volume mesh...
const complexVectorField & newField()
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
dimensionedSymmTensor sqr(const dimensionedVector &dv)
messageStream Info
Information stream (stdout output on master, null elsewhere)
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
dimensionedScalar sqrt(const dimensionedScalar &ds)
Vector< complex > complexVector
A Vector of complex values with 'scalar' precision.
dimensioned< typename typeOfMag< Type >::type > magSqr(const dimensioned< Type > &dt)
constexpr char nl
The newline '\n' character (0x0a)
#define forAll(list, i)
Loop across all elements in list.
A non-counting (dummy) refCount.