[Libmesh-users] NewmarkSystem with mesh refinement

[li p. <li...@ya...>]

Dear developers,
I tried to solve a dynamic problem with NewmarkSystem
and adaptive mesh refinement. In order to let the
system re-assemble vector and matrix, I commented
error() in NewmarkSystem::reinit(), and made
_finished_assemble as public in newmark_system.h . 
Then I put the follwoing code in ex8 after
"t_system.update_u_v_a()"

ErrorVector error;
KellyErrorEstimator error_estimator;
error_estimator.estimate_error(t_system, error);
mesh_refinement.flag_elements_by_error_fraction(error,
0.8,0.2,5);
mesh_refinement.refine_and_coarsen_elements();
t_system._finished_assemble = false;
equation_systems.reinit();

Could you please tell me, Is that everything what I
need to do?  

thanx

pan

__________________________________________________
Do You Yahoo!?
Tired of spam?  Yahoo! Mail has the best spam protection around 
http://mail.yahoo.com