| 7 | { |
| 8 | |
| 9 | inline tmp<scalarField> Ek |
| 10 | ( |
| 11 | const scalar Ea, |
| 12 | const scalar k0, |
| 13 | const scalarField& k |
| 14 | ) |
| 15 | { |
| 16 | tmp<scalarField> tEk = Ea*pow(k/k0, 4.0)*exp(-2.0*sqr(k/k0)); |
| 17 | |
| 18 | /* |
| 19 | scalarField& Ekf = tEk(); |
| 20 | |
| 21 | label i; |
| 22 | forAll(Ekf, i) |
| 23 | { |
| 24 | if (k[i] < 2 || k[i] > 10) |
| 25 | { |
| 26 | Ekf[i] = 0.0; |
| 27 | } |
| 28 | } |
| 29 | */ |
| 30 | |
| 31 | return tEk; |
| 32 | } |
| 33 | |
| 34 | |
| 35 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // |