Hello,
I'm new to deal.ii so please excuse the ignorance I assume behind my question: (I'm playing with a most basic 1 cell elasticity model problem.) While iterating locally over the dofs inside a cell: Is there a quick / simple way to get the (real space) position of the degree of freedom by its local number? (like the way const std::vector< <http://www.dealii.org/6.2.1/doxygen/deal.II/classPoint.html> Point<spacedim> >& <http://www.dealii.org/6.2.1/doxygen/deal.II/classFEValuesBase.html> FEValuesBase< dim, spacedim >::get_quadrature_points ( ) const does for quadrature points.) I've already done some research on this problem but I just don't seem to find a simple method; all I could do would be to check the FE_Q object first as to get the degree and dimension of the shape function, and then conclude its position from the standard dof numbering; at least this works for dofs on vertices. Story behind it, just want to an inquiry on the vertices and assign a right hand side value to the proper dof, if true. I got it to work by just looping over the cell's vertices and concluding from their numbering the according dof, but I'm pretty sure there must be another way, looping over the local dofs (possibly without successively querying cell, dofhandler and triangulation.). After all, I would like to gain more insight on the library, as all the dof_handling and fe definitions are still a lot like a black box to me and I thought it might help to follow the flow of objects and data on a very basic model (likewise the way the data is presented in the gmv output). So I hope someone can help me or link me to a relevant post / documentation, as I certainly haven't found anything suitable yet. Greetings and thank you in advance Till
_______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
