Either that or create an inverse mapping that maps back from the current to the reference configuration. This method could potentially introduce to some numerical error though.
On Thursday, March 30, 2017 at 6:57:12 PM UTC+2, Wolfgang Bangerth wrote: > > > Rajat, > > > I am updating my mesh after every time-increment. So, cell->vertex(v) > > contains the current coordinates of the mesh. I know the coordinates of > > the cell (and whole mesh) at t = 0, I wanted > > to know how I can initialize the fe_values object wrt to the coordinates > > of the cell at t = 0. > > > > If possible, this will enable me to use > > fe_values->get_function_gradients() and fe_values->shape_grad() to get > > the gradients wrt to reference mesh. > > You can't. The FEValues classes only consider the mesh as it is right > now, not how it was a while ago. > > But you could consider leaving the mesh as it is (i.e., in reference > coordinates) and using MappingQEulerian to describe the displacement. > Then, if you use an FEValues without a mapping, you get everything with > regard to the original (reference) configuration. If you use an FEValues > with the mapping, you get the current configuration. > > Best > W. > > -- > ------------------------------------------------------------------------ > Wolfgang Bangerth email: bange...@colostate.edu > www: http://www.math.colostate.edu/~bangerth/ > -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.
