XProblem.cpp

long double TImmersedBoundaryProblem::Coord(T3dNormalGrid *grid, const int index, const int type, const int data_axis)
{
    long double step;
    if (type == COORD_X) {
        step = grid->GetSeparator1().Separation[0];
    }
    else if (type == COORD_Y) {
        step = grid->GetSeparator2().Separation[0];
    }
    else if (type == COORD_Z) {
        step = grid->GetSeparator3().Separation[0];
    }
    else {
        throw "Incorrect type of axis";
    }

    return index * step + (type == data_axis ? 0.5 * step : 0);
}