Function main

examples/IIM/ex9/example.cpp:332–823

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

int
main(int argc, char* argv[])
{
    // Initialize libMesh, PETSc, MPI, and SAMRAI.
    LibMeshInit init(argc, argv);
    SAMRAI_MPI::setCommunicator(PETSC_COMM_WORLD);
    SAMRAI_MPI::setCallAbortInSerialInsteadOfExit();
    SAMRAIManager::startup();

    { // cleanup dynamically allocated objects prior to shutdown

        // 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.
        tbox::Pointer<AppInitializer> app_initializer = new AppInitializer(argc, argv, "IB.log");
        tbox::Pointer<tbox::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();

        const bool uses_exodus = dump_viz_data && !app_initializer->getExodusIIFilename().empty();

        const string viz_dump_dirname = app_initializer->getVizDumpDirectory();
        const string beam_filename = viz_dump_dirname + "/beam.ex2";
        const string bndry_beam_filename = viz_dump_dirname + "/bndry_beam.ex2";
        const string bndry_housing_filename = viz_dump_dirname + "/bndry_housing.ex2";

        const bool dump_postproc_data = app_initializer->dumpPostProcessingData();
        const int postproc_data_dump_interval = app_initializer->getPostProcessingDataDumpInterval();
        const string postproc_data_dump_dirname = app_initializer->getPostProcessingDataDumpDirectory();
        if (dump_postproc_data && (postproc_data_dump_interval > 0) && !postproc_data_dump_dirname.empty())
        {
            Utilities::recursiveMkdir(postproc_data_dump_dirname);
        }

        const bool dump_timer_data = app_initializer->dumpTimerData();
        const int timer_dump_interval = app_initializer->getTimerDumpInterval();

        // Create a simple FE mesh.

        string housing_elem_type = input_db->getStringWithDefault("HOUSING_ELEM_TYPE", "HEX8");
        string beam_elem_type = input_db->getString("BEAM_ELEM_TYPE");
        double mfac_b = input_db->getDouble("MFAC");
        const double n_cycles = input_db->getDouble("NCYCLE");
        if (beam_elem_type == "TET10" || beam_elem_type == "HEX27") mfac_b *= 2.0;
        dx = input_db->getDouble("DX");
        const double ds_b = mfac_b * dx;

        lambda_s = input_db->getDouble("LAMBDA_S");

        mu_s = input_db->getDouble("MU_S");
        bulk_mod = input_db->getDouble("BULK_MOD");
        shear_mod = input_db->getDouble("SHEAR_MOD");

        rho_s = input_db->getDouble("RHO_S");