Function compute_flow_correction

src/IB/IBInstrumentPanel.cpp:188–224

} // init_meter_elements

double
compute_flow_correction(const boost::multi_array<Vector, 1>& U_perimeter,
                        const Vector& U_centroid,
                        const boost::multi_array<Point, 1>& X_perimeter,
                        const Point& X_centroid)
{
    double U_dot_dA = 0.0;
#if (NDIM == 2)
    TBOX_ERROR("no support for 2D flow meters at this time!\n");
    NULL_USE(U_perimeter);
    NULL_USE(U_centroid);
    NULL_USE(X_perimeter);
    NULL_USE(X_centroid);
#endif
#if (NDIM == 3)
    const int num_perimeter_nodes = static_cast<int>(X_perimeter.shape()[0]);
    for (int m = 0; m < num_perimeter_nodes; ++m)
    {
        const Vector& U_perimeter0(U_perimeter[m]);
        const Point& X_perimeter0(X_perimeter[m]);

        const Vector& U_perimeter1(U_perimeter[(m + 1) % num_perimeter_nodes]);
        const Point& X_perimeter1(X_perimeter[(m + 1) % num_perimeter_nodes]);

        // Compute the linear interpolation of the velocity at the center of the
        // triangle.
        const Vector U = (U_perimeter0 + U_perimeter1 + U_centroid) / 3.0;

        // Compute the area weighted normal to the triangle.
        const Vector dA = 0.5 * (X_centroid - X_perimeter0).cross(X_centroid - X_perimeter1);

        // Compute the contribution to U_dot_dA.
        U_dot_dA += U.dot(dA);
    }
#endif
    return U_dot_dA;
} // compute_flow_correction

#if defined(IBAMR_HAVE_SILO)
/*!