main.C File Reference

Go to the source code of this file.

Functions
Number	initial_value (const Point &p, const Parameters &, const std::string &, const std::string &)
void	initialize (EquationSystems &es, const std::string &system_name)
int	main (int argc, char *argv[])

Function Documentation

initial_value()

Number initial_value	(	const Point &	p,
		const Parameters &	,
		const std::string &	,
		const std::string &
	)

Definition at line 44 of file main.C.

Referenced by initialize().

{
  return p(0)*p(0) + 1; // x^2 + 1
}

initialize()

void initialize	(	EquationSystems &	es,
		const std::string &	system_name
	)

Definition at line 52 of file main.C.

References libMesh::EquationSystems::get_system(), initial_value(), libMesh::EquationSystems::parameters, libMesh::System::project_solution(), libMesh::Real, libMesh::Parameters::set(), and libMesh::System::time.

Referenced by libMesh::ParmetisPartitioner::_do_repartition(), and main().

{
  ExplicitSystem & system = es.get_system<ExplicitSystem>(system_name);
  es.parameters.set<Real> ("time") = system.time = 0;
  system.project_solution(initial_value, nullptr, es.parameters);
}

main()

int main	(	int	argc,
		char *	argv[]
	)

Definition at line 60 of file main.C.

References libMesh::EquationSystems::add_system(), libMesh::System::add_variable(), libMesh::System::attach_init_function(), libMesh::MeshTools::Generation::build_cube(), libMesh::HEX8, libMesh::TriangleWrapper::init(), libMesh::EquationSystems::init(), initialize(), libMesh::TET4, libMesh::DTKSolutionTransfer::transfer(), libMesh::EquationSystems::update(), libMesh::System::variable(), and libMesh::MeshOutput< MT >::write_equation_systems().

{
  LibMeshInit init (argc, argv);

#if !defined(LIBMESH_TRILINOS_HAVE_DTK)
  // Skip this example (and use a different return code) if libMesh
  // was compiled without Trilinos+DTK support.
  libmesh_example_requires(false, "--enable-trilinos");

#else

  Mesh from_mesh(init.comm());
  MeshTools::Generation::build_cube(from_mesh, 4, 4, 4, 0, 1, 0, 1, 0, 1, HEX8);
  from_mesh.print_info();
  EquationSystems from_es(from_mesh);
  System & from_sys = from_es.add_system<ExplicitSystem>("From");
  unsigned int from_var = from_sys.add_variable("from");
  from_sys.attach_init_function(initialize);
  from_es.init();

  ExodusII_IO(from_mesh).write_equation_systems("from.e", from_es);

  Mesh to_mesh(init.comm());
  MeshTools::Generation::build_cube(to_mesh, 5, 5, 5, 0, 1, 0, 1, 0, 1, TET4);
  to_mesh.print_info();
  EquationSystems to_es(to_mesh);
  System & to_sys = to_es.add_system<ExplicitSystem>("To");
  unsigned int to_var = to_sys.add_variable("to");
  to_es.init();

  DTKSolutionTransfer dtk_transfer(init.comm());

  dtk_transfer.transfer(from_sys.variable(from_var), to_sys.variable(to_var));

  to_es.update();
  ExodusII_IO(to_mesh).write_equation_systems("to.e", to_es);

#endif

  return 0;
}