26 #include "grins_config.h"
36 #include "libmesh/exact_solution.h"
39 #ifdef GRINS_HAVE_GRVY
43 int main(
int argc,
char* argv[])
46 #ifdef GRINS_USE_GRVY_TIMERS
47 GRVY::GRVY_Timer_Class grvy_timer;
48 grvy_timer.Init(
"GRINS Timer");
55 std::cerr <<
"Error: Must specify libMesh input file." << std::endl;
60 std::string libMesh_input_filename = argv[1];
63 GetPot libMesh_inputfile( libMesh_input_filename );
65 #ifdef GRINS_USE_GRVY_TIMERS
66 grvy_timer.BeginTimer(
"Initialize Solver");
70 libMesh::LibMeshInit libmesh_init(argc, argv);
76 libmesh_init.comm() );
78 #ifdef GRINS_USE_GRVY_TIMERS
79 grvy_timer.EndTimer(
"Initialize Solver");
82 grins.attach_grvy_timer( &grvy_timer );
88 libMesh::Number qoi = grins.get_qoi_value( 0 );
91 const libMesh::Number exact_value = -0.5;
92 const libMesh::Number rel_error = std::fabs( (qoi - exact_value )/exact_value );
93 const libMesh::Number tol = 1.0e-11;
96 std::cerr <<
"Computed voriticity QoI mismatch greater than tolerance." << std::endl
97 <<
"Computed value = " << qoi << std::endl
98 <<
"Exact value = " << exact_value << std::endl
99 <<
"Relative error = " << rel_error << std::endl
100 <<
"Tolerance = " << tol << std::endl;
int main(int argc, char *argv[])