>
>
>> You're sure you need your own from-scratch GUI here? You definitely
>> can't make do with Paraview, VisIt, or some such?
>>
>> Well ... it is possible, of course. I had this simple GUI previously
>>
>
> Okay; if it's not "from-scratch" then that does change the calculation
> a bit.
>
> and being able to interact with the mesh is quite useful for
>> debugging.
>>
>
> I can't argue with that. High on my "wishlist that I will never ever
> get to myself" is support for live Paraview+simulation integration -
> I think their API for that is called Catalyst or Insitu?
>
Precisely. I need only a simple interface right now.
I see value in being able to output internals, move simulation forth and
back, experiment remeshing on-the-gly with variable parameters ...
With a more stable code - and needing more sophisticated plots - I will
move to paraview or similar.
But I have no energy to learn a new API right now (learning CGAL and
Libmesh was already quite consuming).
Potentially, I can do that later on. Thanks for the suggestion.
Yeah, I forgot that the sys.point_value(0,pt,elem) overload is really
> only intended for use on local elements. If all you care about is
> nodal values for now, then I'd use sys.solution->localize_to_one to
> make a serial copy of your solution vector on proc 0 too, then inside
> your loop query that copy (using node->dof_number(0,0,0) to get the
> dof index for the first system's first variable) for the value.
>
I did it - but still could not get it working.
I added this call on the master process (id=0):
vector<double> _sys_solution; /// in .hxx
sys.solution->localize_to_one( _sys_solution );
The system got stuck, with all four processes with 100% processor usage.
Below is my main routine, I am listing here because I am not sure the
non-master processes are handled correctly (simply returning zero).
69 sys.solve();
72 MeshSerializer serialize(lm_mesh, true, true);
74 if ( lm_mesh.processor_id() == 0) {
77 FEMainWindow *win = new FEMainWindow(lm_mesh, sys);
78 win->show();
79 return app.exec();
80 } else
81 return 0;
82 }
------------------------------------------------------------------------------
Check out the vibrant tech community on one of the world's most
engaging tech sites, Slashdot.org! http://sdm.link/slashdot
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
