Thankyou! that works :) On Saturday, 18 March 2017 12:03:55 UTC+5:30, Sumedh Yadav wrote: > > Hello, > I need to do some post-processing on computed solution data (say > Vector<double> Solution(dof_handler.n_dofs()) is the solution vector). For > post-processing field (say Vector<double> Curvature(triangulation.n_dofs()) > is the scalar field to be computed) I decided to have one entry per cell. > Once I calculate this entry in 'cell iterator loop' I have to transfer it > to global vector Curvature based on some index of the cell. I am just > drawing the analogy from what we do during the assembling of system_matrix > and/or system_rhs. There we use cell iterator functionality > get_dof_indices() to get the global indices of DoFs that are of on the > current cell. But I am not able to find analogous method to get an > cell-wise index (we can't use cell iterator return value since it is not > unsigned int) which can be used to transfer cell data to the global vector > Curvature(triangulation.n_dofs()). Please help me to figure it out. > Secondly, I want to use this global vector > Curvature(triangulation.n_dofs()) in a DataOut object essentially to add > this vector to the object via add_data_vector(Curvature, "curvature") > method. I believe (read in another question thread in the group) the > DataOut object will not have any problems accepting a vector which has size > triangulation.n_dofs() but not dof_handler.n_dofs() since I initialize the > DataOut object with dof_handler? I meant that it would judiciously > distinguish this vector from other vectors which have size > dof_handler.n_dofs() and output it correctly? >
-- 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.
