Method computeLagrangianFluidSource

src/IB/IBMethod.cpp:1053–1066

} // getInterpOperator

void
IBMethod::computeLagrangianFluidSource(const double data_time)
{
    if (!d_ib_source_fcn) return;
    const int coarsest_ln = 0;
    const int finest_ln = d_hierarchy->getFinestLevelNumber();
    for (int ln = coarsest_ln; ln <= finest_ln; ++ln)
    {
        if (d_n_src[ln] == 0) continue;
        std::fill(d_Q_src[ln].begin(), d_Q_src[ln].end(), 0.0);
        d_ib_source_fcn->computeSourceStrengths(d_Q_src[ln], d_hierarchy, ln, data_time, d_l_data_manager);
    }
    return;
} // computeLagrangianFluidSource

void
IBMethod::spreadFluidSource(const int q_data_idx,