Function main

tests/IBFE/explicit_ex8.cpp:313–717

 *******************************************************************************/

int
main(int argc, char** argv)
{
    // Initialize IBAMR and libraries. Deinitialization is handled by this object as well.
    IBTKInit ibtk_init(argc, argv, MPI_COMM_WORLD);
    const LibMeshInit& init = ibtk_init.getLibMeshInit();

    PetscOptionsSetValue(nullptr, "-ksp_rtol", "1e-16");
    PetscOptionsSetValue(nullptr, "-ksp_rtol", "1e-16");
    PetscOptionsSetValue(nullptr, "-stokes_ksp_atol", "1e-14");
    PetscOptionsSetValue(nullptr, "-stokes_ksp_rtol", "1e-14");

    { // cleanup dynamically allocated objects prior to shutdown
        // prevent a warning about timer initializations
        TimerManager::createManager(nullptr);

        // Parse command line options, set some standard options from the input
        // file, initialize the restart database (if this is a restarted run),
        // and enable file logging.
        Pointer<AppInitializer> app_initializer = new AppInitializer(argc, argv, "IB.log");
        Pointer<Database> input_db = app_initializer->getInputDatabase();

        // Get various standard options set in the input file.
        const bool dump_viz_data = app_initializer->dumpVizData();
        const int viz_dump_interval = app_initializer->getVizDumpInterval();
        const bool uses_visit = dump_viz_data && app_initializer->getVisItDataWriter();
#ifdef LIBMESH_HAVE_EXODUS_API
        const bool uses_exodus = dump_viz_data && !app_initializer->getExodusIIFilename().empty();
#else
        const bool uses_exodus = false;
        if (!app_initializer->getExodusIIFilename().empty())
        {
            plog << "WARNING: libMesh was compiled without Exodus support, so no "
                 << "Exodus output will be written in this program.\n";
        }
#endif
        const string block1_exodus_filename = app_initializer->getExodusIIFilename("block1");
        const string block2_exodus_filename = app_initializer->getExodusIIFilename("block2");
        const string beam_exodus_filename = app_initializer->getExodusIIFilename("beam");

        const bool dump_restart_data = app_initializer->dumpRestartData();
        const int restart_dump_interval = app_initializer->getRestartDumpInterval();
        const string restart_dump_dirname = app_initializer->getRestartDumpDirectory();

        // Create a simple FE mesh.
        const double dx = input_db->getDouble("DX");
        const double ds_block = input_db->getDouble("BLOCK_MFAC") * dx;
        const double ds_beam = input_db->getDouble("BEAM_MFAC") * dx;

        string block_elem_type = input_db->getStringWithDefault("BLOCK_ELEM_TYPE", "QUAD9");
        string beam_elem_type = input_db->getStringWithDefault("BEAM_ELEM_TYPE", "QUAD9");

        Mesh block1_mesh(init.comm(), NDIM);
        MeshTools::Generation::build_square(block1_mesh,
                                            ceil(0.5 / ds_block),
                                            ceil(0.5 / ds_block),
                                            0.0,
                                            0.5,