Is there a good way to explicitly compute the nodal data and save them? I use LinearImplicitSystem to solve the Stokes equation for velocities and pressure. I want to use the computed nodal velocities to explicitly compute the displacements at each node using Euler's method or higher order method. Also, I want to use the displacements and velocities to explicitly compute other physical parameters, such as deformation rate, for post-processing purpose.
Right now, I explicitly compute and save them in NumericVector whose length is equal to the number of nodes and put them into the system using Syetem::add_vector. If I need the value inside a local element, I just interpolate them using shape functions. However, I think there is a more efficient way using ExplicitSystem or other classes. If I am doing something unnecessary, could you tell me a better way to do it? Is there any examples? Thanks, Yusuke ------------------------------------------------------------------------------ Live Security Virtual Conference Exclusive live event will cover all the ways today's security and threat landscape has changed and how IT managers can respond. Discussions will include endpoint security, mobile security and the latest in malware threats. http://www.accelacomm.com/jaw/sfrnl04242012/114/50122263/ _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
