| 389 | else |
| 390 | { |
| 391 | scalarField dd((d & binDir_) - binMin_); |
| 392 | |
| 393 | forAll(dd, i) |
| 394 | { |
| 395 | label bini = min(max(floor(dd[i]/binDx_), 0), force_[0].size() - 1); |
| 396 | |
| 397 | force_[0][bini] += fN[i]; |
| 398 | force_[1][bini] += fT[i]; |
| 399 | force_[2][bini] += fP[i]; |
| 400 | moment_[0][bini] += Md[i]^fN[i]; |
| 401 | moment_[1][bini] += Md[i]^fT[i]; |
| 402 | moment_[2][bini] += Md[i]^fP[i]; |
| 403 | } |
| 404 | } |
| 405 | } |
| 406 |
no test coverage detected