| 122 | } |
| 123 | |
| 124 | void |
| 125 | advect() |
| 126 | { |
| 127 | #ifdef _MSC_VER |
| 128 | concurrency::parallel_for_each |
| 129 | #else |
| 130 | __gnu_parallel::for_each |
| 131 | #endif |
| 132 | (positions.begin(), positions.end(), [&](std::array<Real, 3>& x) |
| 133 | { |
| 134 | std::array<Real, 3> v = enright_velocity_field(x); |
| 135 | x[0] += 0.005 * v[0]; |
| 136 | x[1] += 0.005 * v[1]; |
| 137 | x[2] += 0.005 * v[1]; |
| 138 | } |
| 139 | ); |
| 140 | } |
| 141 | |
| 142 | int main(int argc, char* argv[]) |
| 143 | { |
no test coverage detected