meshavg.C File Reference

Go to the source code of this file.

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

Variables
unsigned char	dim = 2

Function Documentation

main()

int main	(	int	argc,
		char **	argv
	)

Definition at line 32 of file meshavg.C.

References dim, libMesh::index_range(), libMesh::TriangleWrapper::init(), libMesh::out, libMesh::EquationSystems::READ_BASIC_ONLY, libMesh::EquationSystems::READ_DATA, libMesh::EquationSystems::READ_HEADER, and libMesh::NameBasedIO::write_equation_systems().

{
  LibMeshInit init(argc, argv);

  Mesh mesh1(init.comm(), dim);
  EquationSystems es1(mesh1);

  libMesh::out << "Usage: " << argv[0]
               << " outputsolution mesh firstsolution [secondsolution ...]" << std::endl;

  mesh1.read(argv[2]);
  libMesh::out << "Loaded mesh " << argv[2] << std::endl;

  mesh1.print_info();

  es1.read(argv[3],
           EquationSystems::READ_HEADER |
           EquationSystems::READ_DATA |
           EquationSystems::READ_BASIC_ONLY);
  libMesh::out << "Loaded first solution " << argv[3] << std::endl;

  es1.print_info();

  std::vector<std::string> sysnames;
  std::vector<std::unique_ptr<NumericVector<libMesh::Number>>> summed_solutions;

  for (unsigned int s = 0; s != es1.n_systems(); ++s)
    {
      sysnames.push_back(es1.get_system(s).name());
      summed_solutions.emplace_back(es1.get_system(s).solution->clone());
    }

  for (int i=4; i < argc; ++i)
    {
      Mesh mesh2(init.comm(), dim);
      EquationSystems es2(mesh2);
      mesh2.read(argv[2]);
      es2.read(argv[i],
               EquationSystems::READ_HEADER |
               EquationSystems::READ_DATA |
               EquationSystems::READ_BASIC_ONLY);
      libMesh::out << "Loaded next solution " << argv[i] << std::endl;

      for (auto s : index_range(sysnames))
        {
          libmesh_error_msg_if(!es2.has_system(sysnames[s]),
                               "EquationSystems object does not have " << sysnames[s]);

          (*summed_solutions[s]) += *es2.get_system(s).solution;
        }
    }

  int n_solutions = argc - 3;

  for (auto s : index_range(sysnames))
    {
      (*summed_solutions[s]) /= n_solutions;
      es1.get_system(s).solution->swap(*summed_solutions[s]);
      es1.get_system(s).solution->close();
    }

  NameBasedIO(mesh1).write_equation_systems (argv[1], es1);

  return 0;
}

Variable Documentation

dim

unsigned char dim = 2

Definition at line 30 of file meshavg.C.

Referenced by main().