C++ (Cpp) FlowField::getRHSの例

コード例 #1

ファイル: RHSStencil.cpp プロジェクト: cqql/turbulence

void RHSStencil::apply(FlowField& flowField, int i, int j) {
  flowField.getRHS().getScalar(i, j) =
      1.0 / _parameters.timestep.dt *
      ((flowField.getFGH().getVector(i, j)[0] -
        flowField.getFGH().getVector(i - 1, j)[0]) /
           _parameters.meshsize->getDx(i, j) +
       (flowField.getFGH().getVector(i, j)[1] -
        flowField.getFGH().getVector(i, j - 1)[1]) /
           _parameters.meshsize->getDy(i, j));
}