Method reconstructVariables

src/IB/IBFECentroidPostProcessor.cpp:148–327

/////////////////////////////// PROTECTED ////////////////////////////////////

void
IBFECentroidPostProcessor::reconstructVariables(double data_time)
{
    EquationSystems* equation_systems = d_fe_data_manager->getEquationSystems();
    MeshBase& mesh = equation_systems->get_mesh();
    const int dim = mesh.mesh_dimension();
    std::unique_ptr<QBase> qrule = QBase::build(QGAUSS, NDIM, CONSTANT);

    // Set up all system data required to evaluate the mesh functions.
    FEDataInterpolation fe(dim, d_fe_data_manager->getFEData());
    fe.attachQuadratureRule(qrule.get());
    fe.evalQuadraturePoints();

    auto& X_system = equation_systems->get_system<System>(d_fe_data_manager->getCurrentCoordinatesSystemName());
    NumericVector<double>& X_data = *(X_system.current_local_solution);
    copy_and_synch(*X_system.solution, X_data);
    std::vector<int> vars(NDIM);
    for (unsigned int d = 0; d < NDIM; ++d) vars[d] = d;
    const size_t X_sys_idx = fe.registerInterpolatedSystem(X_system, vars, vars, &X_data);

    const size_t num_scalar_vars = d_scalar_var_systems.size();
    std::vector<const DofMap*> scalar_var_dof_maps(num_scalar_vars);
    std::vector<std::vector<unsigned int>> scalar_var_dof_indices(num_scalar_vars);
    std::vector<std::vector<size_t>> scalar_var_fcn_system_idxs(num_scalar_vars);
    for (unsigned int k = 0; k < num_scalar_vars; ++k)
    {
        scalar_var_dof_maps[k] = &d_scalar_var_systems[k]->get_dof_map();
        fe.setupInterpolatedSystemDataIndexes(
            scalar_var_fcn_system_idxs[k], d_scalar_var_system_data[k], equation_systems);
    }

    const size_t num_vector_vars = d_vector_var_systems.size();
    std::vector<const DofMap*> vector_var_dof_maps(num_vector_vars);
    std::vector<std::vector<std::vector<unsigned int>>> vector_var_dof_indices(num_vector_vars);
    std::vector<std::vector<size_t>> vector_var_fcn_system_idxs(num_vector_vars);
    for (unsigned int k = 0; k < num_vector_vars; ++k)
    {
        vector_var_dof_maps[k] = &d_vector_var_systems[k]->get_dof_map();
        vector_var_dof_indices[k].resize(d_vector_var_dims[k]);
        fe.setupInterpolatedSystemDataIndexes(
            vector_var_fcn_system_idxs[k], d_vector_var_system_data[k], equation_systems);
    }

    const size_t num_tensor_vars = d_tensor_var_systems.size();
    std::vector<const DofMap*> tensor_var_dof_maps(num_tensor_vars);
    std::vector<boost::multi_array<std::vector<unsigned int>, 2>> tensor_var_dof_indices(num_tensor_vars);
    std::vector<std::vector<size_t>> tensor_var_fcn_system_idxs(num_tensor_vars);
    for (unsigned int k = 0; k < num_tensor_vars; ++k)
    {
        tensor_var_dof_maps[k] = &d_tensor_var_systems[k]->get_dof_map();
        using array_type = boost::multi_array<std::vector<unsigned int>, 2>;
        array_type::extent_gen extents;
        tensor_var_dof_indices[k].resize(extents[d_tensor_var_dims[k]][d_tensor_var_dims[k]]);
        fe.setupInterpolatedSystemDataIndexes(
            tensor_var_fcn_system_idxs[k], d_tensor_var_system_data[k], equation_systems);
    }

    fe.init();