I appreciate the clarification. I thought that global indexing was no longer present as the solution vector is distributed. I have one more doubt. I want to create a vector (H_vector) that stores some value for each Gauss point in the domain. For a serial problem, I was doing something like this:* H_vector = Vector<double> (8*(triangulation.n_active_cells()));*
*8 because the problem is in 3d, so I have 8 Gauss points per cell.* Now, for a MPI code, this H_vector would also be *LA::MPI::Vector*, and its size should be* 8*no. of cells owned by the processor.* How do I get the no.of cells owned by the processor? Thanks and regards Wasim On Thursday, April 6, 2023 at 9:24:14 PM UTC+5:30 Wolfgang Bangerth wrote: > On 4/6/23 06:02, Wasim Niyaz Munshi ce21d400 wrote: > > > > I don't have a solution_vector for a parallel code, but a > > locally_relevant_solution. I want to know that, given this > > locally_relevant_solution and the cell, how do I get the element_sol? > > The global_dof will not be helpful here, as the solution_vector is > > distributed across a number of processors. > > Daniel's question is correct, but to this specific point: A distributed > vector (and it's locally relevant incarnation) is still a global vector, > indexed by global degree of freedom numbers, and so the code remains > correct. > > Best > W. > > -- > ------------------------------------------------------------------------ > Wolfgang Bangerth email: bang...@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. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/afbbd383-46a9-40cc-9a02-d2b4347898b3n%40googlegroups.com.
