CArl_FETI_solution.cpp

Go to the documentation of this file.

#include "CArl_FETI_solution.h"

int main(int argc, char** argv) {

    // --- Initialize libMesh
    libMesh::LibMeshInit init(argc, argv);

    // Do performance log?
    libMesh::PerfLog perf_log("Main program");

    // libMesh's C++ / MPI communicator wrapper
    libMesh::Parallel::Communicator& WorldComm = init.comm();

    // Number of processors and processor rank.
    int rank = WorldComm.rank();
    int nodes = WorldComm.size();

    // --- Set up inputs

    // Command line parser
    GetPot command_line(argc, argv);

    // File parser
    GetPot field_parser;

    // If there is an input file, parse it to get the parameters. Else, parse the command line
    std::string input_filename;
    if (command_line.search(2, "--inputfile", "-i")) {
        input_filename = command_line.next(input_filename);
        field_parser.parse_input_file(input_filename, "#", "\n", " \t\n");
    } else {
        field_parser = command_line;
    }

    carl::feti_set_sol_params input_params;
    get_input_params(field_parser, input_params);

    // Object containing the FETI operations
    carl::FETI_Operations feti_op(WorldComm,input_params.scratch_folder_path);

    // --- Define if the rb modes will be used or not
    feti_op.using_rb_modes(input_params.bUseRigidBodyModes);

    // Read the decoupled solutions, K_i * u_0,i  = F_i 
    feti_op.read_decoupled_solutions();

    // Read the solutions of K_i * x_i(FINAL)  = C_i^t * phi(FINAL)
    feti_op.read_ext_solver_output();

    // Read the rb modes correction, 'RB_corr(FINAL)'
    if(input_params.bUseRigidBodyModes)
    {
        feti_op.read_rb_corr();
    }

    // Calculate the solution
    feti_op.calculate_coupled_solution();

    // Export it (finallly!)
    feti_op.export_coupled_solution(input_params.output_folder);

    return 0;
}