Function main

tests/CIB/cib_double_shell.cpp:122–557

* For each run, the input filename and restart information (if needed) must * * be given on the command line. For non-restarted case, command line is: * * * * executable * *

 *                                                                             *
 *******************************************************************************/
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);

    // Several parts of the code (such as LDataManager) expect mesh files,
    // specified in the input file, to exist in the current working
    // directory. Since tests are run in temporary directories we need to regenerate these input
    // to work. We also create a petsc options file for CIB solvers.
    if (IBTK_MPI::getRank() == 0)
    {
        std::ifstream shell_in_vertex_stream(SOURCE_DIR "/shell_3d_in.vertex");
        std::ofstream shell_in_cwd("shell_3d_in.vertex");
        shell_in_cwd << shell_in_vertex_stream.rdbuf();

        std::ifstream shell_out_vertex_stream(SOURCE_DIR "/shell_3d_out.vertex");
        std::ofstream shell_out_cwd("shell_3d_out.vertex");
        shell_out_cwd << shell_out_vertex_stream.rdbuf();

        std::ifstream petsc_options_stream(SOURCE_DIR "/petsc_options.dat");
        std::ofstream petsc_options_cwd("petsc_options.dat");
        petsc_options_cwd << petsc_options_stream.rdbuf();
    }

    { // 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.
        Pointer<AppInitializer> app_initializer = new AppInitializer(argc, argv, "CIB.log");
        Pointer<Database> input_db = app_initializer->getInputDatabase();

        // some variants of this test use the SVD, some don't
        const auto inverse_type = IBAMR::string_to_enum<IBAMR::MobilityMatrixInverseType>(
            input_db->getStringWithDefault("mobility_inverse_type", "LAPACK_SVD"));

        // Read default Petsc options
        if (input_db->keyExists("petsc_options_file"))
        {
            std::string petsc_options_file = input_db->getString("petsc_options_file");
            PetscOptionsInsertFile(PETSC_COMM_WORLD, nullptr, petsc_options_file.c_str(), PETSC_TRUE);
        }

        // 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().isNull();

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

        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())
        {