Hi to all,
I am trying to learn libmesh but I still have some problems.
Right now I want to solve the diffusion equation but after every time
step I would like to look at the results to evaluate them.
Unfortunately I couldn't figure out how to do so.
I have just tried to use something which I found in
exact_error_estimator.C Here my slightly modified code:
AutoPtr<MeshFunction> fine_values;
AutoPtr<NumericVector<Number> >
fine_soln=NumericVector<Number>::build();
const System& fine_system = es.get_system(system.name());
std::vector<Number> global_soln;
// system fine_system.update_global_solution(global_soln);
fine_soln->init (global_soln.size(), true, SERIAL);
(*fine_soln) = global_soln;
fine_values = AutoPtr<MeshFunction> (new MeshFunction(es, *fine_soln,
fine_system.get_dof_map(), fine_system.variable_number("u")));
fine_values->init();
const std::vector<Point>& fe_points = fe->get_xyz();
DenseVector<Number> evaluation;
fine_values.operator(fe_points,time,&evaluation);
When I try to compile this I get the following error message: error:
expected type-specifier before ‘(’ token which I don't understand
because in my opinien I specified the types. However, I am not ever sure
if I can use this to get where I want.
thank you very much, Robert
------------------------------------------------------------------------------
Increase Visibility of Your 3D Game App & Earn a Chance To Win $500!
Tap into the largest installed PC base & get more eyes on your game by
optimizing for Intel(R) Graphics Technology. Get started today with the
Intel(R) Software Partner Program. Five $500 cash prizes are up for grabs.
http://p.sf.net/sfu/intelisp-dev2dev
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
