meshplot.C File Reference

Go to the source code of this file.

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

Function Documentation

main()

int main	(	int	argc,
		char **	argv
	)

Definition at line 28 of file meshplot.C.

References libMesh::TriangleWrapper::init(), mesh, libMesh::out, libMesh::MeshBase::print_info(), libMesh::MeshBase::read(), libMesh::EquationSystems::READ_ADDITIONAL_DATA, libMesh::EquationSystems::READ_BASIC_ONLY, libMesh::EquationSystems::READ_DATA, libMesh::EquationSystems::READ_HEADER, libMesh::MeshBase::write(), and libMesh::NameBasedIO::write_equation_systems().

{
  LibMeshInit init(argc, argv);

  Mesh mesh(init.comm());
  EquationSystems es(mesh);

  if (argc < 3 || argc > 4)
    libmesh_error_msg
      ("Usage: " << argv[0] <<
       " inputmesh [inputsolution] outputplot");

  LOG_CALL("mesh.read()", "main", mesh.read(argv[1]));

  libMesh::out << "Loaded mesh " << argv[1] << std::endl;
  mesh.print_info();

  if (argc > 3)
    {
      std::string solnname = argv[2];

      LOG_CALL("es.read()", "main",
               es.read(solnname,
                       EquationSystems::READ_HEADER |
                       EquationSystems::READ_DATA |
                       EquationSystems::READ_ADDITIONAL_DATA |
                       EquationSystems::READ_BASIC_ONLY));

      libMesh::out << "Loaded solution " << solnname << std::endl;
      es.print_info();
    }

  std::string outputname(argv[argc-1]);
  {
    LOG_SCOPE("write_equation_systems()", "main");
    if ((outputname.find(".xda") != std::string::npos) ||
        (outputname.find(".xdr") != std::string::npos))
      {
        mesh.write("mesh-"+outputname);
        es.write("soln-"+outputname);
      }
    else
      NameBasedIO(mesh).write_equation_systems (outputname, es);
  }

  libMesh::out << "Wrote output " << outputname << std::endl;

  return 0;
}