libmesh_apply_solution_homogeneous.cpp File Reference

#include "libmesh_apply_solution_homogeneous.h"

Go to the source code of this file.

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

Function Documentation

int main	(	int	argc,
		char **	argv
	)

Definition at line 3 of file libmesh_apply_solution_homogeneous.cpp.

                                {

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

    // --- 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::libmesh_apply_solution_input_params input_params;
    carl::get_input_params(field_parser, input_params);

    // Check libMesh installation dimension
    const unsigned int dim = 3;

    libmesh_example_requires(dim == LIBMESH_DIM, "3D support");

    // - Parallelized meshes A
    libMesh::Mesh system_mesh(WorldComm, dim);
    system_mesh.read(input_params.input_mesh);
    system_mesh.prepare_for_use();

    // Set the equation systems object
    libMesh::EquationSystems equation_systems(system_mesh);

    // Add linear elasticity and stress
    libMesh::LinearImplicitSystem& elasticity_system
                                        = add_elasticity(equation_systems);
    add_stress(equation_systems);

    // Initialize the equation systems
    equation_systems.init();

    // Homogeneous physical properties
    set_homogeneous_physical_properties(equation_systems, input_params.physical_params_file);

    // Read the solution vector
    libMesh::PetscVector<libMesh::Real> * sol_vec_ptr = libMesh::cast_ptr<libMesh::PetscVector<libMesh::Real> * >(elasticity_system.solution.get());
    carl::read_PETSC_vector(sol_vec_ptr->vec(),input_params.input_vector, WorldComm.get());

    // Close it and update
    elasticity_system.solution->close();
    elasticity_system.update();

    // Calculate the stress
    compute_stresses(equation_systems);

    // Export solution
#ifdef LIBMESH_HAVE_EXODUS_API
        libMesh::ExodusII_IO exo_io_interface(system_mesh, /*single_precision=*/true);

        std::set<std::string> system_names;
        system_names.insert("Elasticity");
        exo_io_interface.write_equation_systems(input_params.output_mesh,equation_systems,&system_names);
        exo_io_interface.write_element_data(equation_systems);
#endif

    return 0;
}