Thomas,thank you, the vertices would have the same order as the velocity vector.
On Monday, December 5, 2016 at 3:19:20 AM UTC, thomas stephens wrote: > > Bernard, you probably want to take a look at Step 18 in the tutorial, > specifically the move_mesh() > <https://www.dealii.org/8.4.1/doxygen/deal.II/step_18.html#TopLevelmove_mesh>function. > > > > In short, you will want to loop through the cells of your triangulation > and call exactly cell->vertex(i) as you have already pointed out. If > you were to get the vertices all at once, what order would they be in > relative to your velocity vector? Looping over the cells at least gives you > a chance to call > cell->vertex_dof_index(v,d) so you can return the global dof index into > your velocity vector. > > > > > On Sunday, December 4, 2016 at 6:25:27 PM UTC-5, BENARD KIPCHUMBA wrote: >> >> Thanks Timo Heister, >> >> My interest is to move the vertices by some velocity and would therefore >> first ask for the vertex coordinates of each cell using cell->vertex(i) >> and then displace them by velocity. Is this the best approach? How can i >> get the vertices at once? Cheers >> >> On Sunday, December 4, 2016 at 5:57:23 PM UTC, Timo Heister wrote: >>> >>> Benard, >>> >>> the get_function_values is made to evaluate a finite element solution >>> at quadrature points within a cell. If you are just interested in >>> values at DoFs you can access the vector of coefficients directly: >>> Velocity_solution.block(0)[j] will return the value that corresponds >>> to the j-th DoF shape function. If your FE space is nodal, this will >>> be the function value at that support point. >>> >>> On Sun, Dec 4, 2016 at 12:37 PM, BENARD KIPCHUMBA <[email protected]> >>> wrote: >>> > >>> > I want to use fe_function_values to get values for the entire of >>> degrees and >>> > am using >>> > >>> > fe_values.get_function_values (Velocity_solution.block(0),x_value); >>> > >>> > where size of (Velocity_solution.block(0) is supposed to be >>> > dof_handler.n_dofs() but because x_value is defined only for a >>> particular >>> > cell, I get the error Dimension 0 not equal to 1. >>> > >>> > How do I get function values for the entire degrees of freedom? >>> > >>> > With thanks, >>> > >>> > Benard >>> > >>> > -- >>> > The deal.II project is located at >>> https://urldefense.proofpoint.com/v2/url?u=http-3A__www.dealii.org_&d=CwIBaQ&c=Ngd-ta5yRYsqeUsEDgxhcqsYYY1Xs5ogLxWPA_2Wlc4&r=4k7iKXbjGC8LfYxVJJXiaYVu6FRWmEjX38S7JmlS9Vw&m=cyAHZ9QFB-vto36D97sZUlYu6mcMYRJ7B7bnJKRrhKY&s=6RH5_9bnc-9bDAN-q9HgmchzIXlamTxSYoErARyQalA&e= >>> >>> > For mailing list/forum options, see >>> > >>> https://urldefense.proofpoint.com/v2/url?u=https-3A__groups.google.com_d_forum_dealii-3Fhl-3Den&d=CwIBaQ&c=Ngd-ta5yRYsqeUsEDgxhcqsYYY1Xs5ogLxWPA_2Wlc4&r=4k7iKXbjGC8LfYxVJJXiaYVu6FRWmEjX38S7JmlS9Vw&m=cyAHZ9QFB-vto36D97sZUlYu6mcMYRJ7B7bnJKRrhKY&s=s6d_WIfIWKuXFBgPJVVkbIM-tyuiNjEdJ9Ya5w-6LTs&e= >>> >>> > --- >>> > 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 [email protected]. >>> > For more options, visit >>> https://urldefense.proofpoint.com/v2/url?u=https-3A__groups.google.com_d_optout&d=CwIBaQ&c=Ngd-ta5yRYsqeUsEDgxhcqsYYY1Xs5ogLxWPA_2Wlc4&r=4k7iKXbjGC8LfYxVJJXiaYVu6FRWmEjX38S7JmlS9Vw&m=cyAHZ9QFB-vto36D97sZUlYu6mcMYRJ7B7bnJKRrhKY&s=PNMvimxd8ons5DaW7ZJtOfbZZlxH1A9te08pWeuM6lw&e= >>> >>> . >>> >>> >>> >>> -- >>> Timo Heister >>> http://www.math.clemson.edu/~heister/ >>> >> -- 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 [email protected]. For more options, visit https://groups.google.com/d/optout.
