Function build_meter_web

src/IB/IBInstrumentPanel.cpp:230–303

! * \brief Build a local mesh database entry corresponding to a meter web. */

 * \brief Build a local mesh database entry corresponding to a meter web.
 */
void
build_meter_web(DBfile* dbfile,
                std::string& dirname,
                const boost::multi_array<Point, 2>& X_web,
                const boost::multi_array<Vector, 2>& dA_web,
                const int timestep,
                const double simulation_time)
{
    const int npoints = static_cast<int>(X_web.num_elements());

    std::vector<float> block_X(NDIM * npoints);
    std::vector<float> block_dA(NDIM * npoints);

    for (unsigned int m = 0, i = 0; m < X_web.shape()[0]; ++m)
    {
        for (unsigned int n = 0; n < X_web.shape()[1]; ++n, ++i)
        {
            // Get the coordinate and normal vector data.
            for (unsigned int d = 0; d < NDIM; ++d)
            {
                block_X[d * npoints + i] = static_cast<float>(X_web[m][n][d]);
                block_dA[d * npoints + i] = static_cast<float>(dA_web[m][n][d]);
            }
        }
    }

    // Set the working directory in the Silo database.
    if (DBSetDir(dbfile, dirname.c_str()) == -1)
    {
        TBOX_ERROR("IBInstrumentPanel::build_meter_web():\n"
                   << "  Could not set directory " << dirname << std::endl);
    }

    // Write out the variables.
    int cycle = timestep;
    auto time = static_cast<float>(simulation_time);
    double dtime = simulation_time;

    static const int MAX_OPTS = 3;
    DBoptlist* optlist = DBMakeOptlist(MAX_OPTS);
    DBAddOption(optlist, DBOPT_CYCLE, &cycle);
    DBAddOption(optlist, DBOPT_TIME, &time);
    DBAddOption(optlist, DBOPT_DTIME, &dtime);

    const char* meshname = "mesh";
    std::vector<float*> coords(NDIM);
    for (unsigned int d = 0; d < NDIM; ++d)
    {
        coords[d] = &block_X[d * npoints];
    }

    int ndims = NDIM;

    DBPutPointmesh(dbfile, meshname, ndims, &coords[0], npoints, DB_FLOAT, optlist);

    const char* varname = "scaled_normal";
    std::vector<float*> vars(NDIM);
    for (unsigned int d = 0; d < NDIM; ++d)