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.

Reply via email to