63 RootDeltaT(
sqrt(DeltaT)),
74 scalar sqrKupper =
sqr(Kupper);
75 scalar sqrKlower =
sqr(Klower) + SMALL;
80 if ((sqrK =
magSqr(K[i])) < sqrKupper && sqrK > sqrKlower)
82 UOfield[i] = Scale*Sigma*WeinerProcess();
104 scalar sqrKupper =
sqr(Kupper);
105 scalar sqrKlower =
sqr(Klower) + SMALL;
110 if ((sqrK =
magSqr(K[i])) < sqrKupper && sqrK > sqrKlower)
114 (1.0 - Alpha*DeltaT)*UOfield[i]
115 + Scale*Sigma*WeinerProcess();
119 Info<<
" Number of forced K = " << count <<
nl;