Wasim,
step-40 uses Q2 elements which means that you have 9 dofs/cell in 2D and
for 32 cells in each direction (33+32)^2=65^2=4225 dofs.
Best,
Daniel
On Fri, Apr 7, 2023 at 5:47 AM Wasim Niyaz Munshi ce21d400 <
ce21d...@smail.iitm.ac.in> wrote:
> Thank You, Prof. Bangerth.
> I have one more
Thank You, Prof. Bangerth.
I have one more doubt about the output of step-40.
For cycle 0, there are 1024 cells (32*32 mesh). So, we have 33 nodes each,
along x and y. So, shouldn't the number of DOFs be 1089 instead of 4225?
Regards
Wasim
On Thursday, April 6, 2023 at 11:27:32 PM UTC+5:30
On 4/6/23 10:18, Wasim Niyaz Munshi ce21d400 wrote:
How do I get the no.of cells owned by the processor?
Triangulation::n_locally_owned_active_cells().
Best
W.
--
Wolfgang Bangerth email:
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
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
Wasim,
The answer depends very much on what you actually want to do with that
solution vector.
Do you want a representation of the solution (assuming you are using Q1?
nodal elements) on a single process/all processes
or are you just interested in the partial solution on every process
separately?
Hello everyone.
I want to extract the element solution vector from the global solution once
the problem is solved in step-40. For a serial code, I would do something
like this:
*int i=0;*
*for (const auto vertex : cell->vertex_indices()) { int a =
(cell->vertex_dof_index(vertex, 0));