Dear all,
I am using dealii for a CFD solver and the solution vector which is passed to the linear algebra solver is a BlockVector where block zero is the velocity and block 1 is the pressure. I would like to obtain the global DOF numbers of all the entries in the block vector; I would therefore like to know if the the indices of the solution vector (especially for block zero) corresponds to the global DOF numbers. The following lines (ignoring other details) illustrates the question I am asking: //solution vector BlockVector<double> newton_update; //the solution steps SolverFGMRES<BlockVector<double>> gmres(solver_control, vector_memory); gmres.solve(system_matrix, newton_update, system_rhs, *preconditioner); I want to know for instance if index "44" used in the snippet below will correspond to global dof number 44. newton_update.block(0)[44] In the case where the answer to the above question is no, could someone please suggest a way to obtain the global DOF number of each entry in the BlockVector? My ultimate goal is the obtain the DOF and real coordinates (x,y,z) of each entry in the block vector. Fortunately, I already obtained the DOFs and their corresponding coordinates in the real space using DoFTools::map_dofs_to_support_points(...) ; I want to get the DOFs for the BlockVector entries so I can match the the coordinates to the entries in the BlockVector. Any suggestion on a better way to do it will also be much appreciated. Thanks in advance for the help. Best, Isaac. -- 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. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/3c2f6af4-000c-44fa-bdea-d4edc6da3c71n%40googlegroups.com.
