Re: [deal.II] determine global dof index of a point
There is also DoFTools::map_support_points_to_dofs, see https://www.dealii.org/current/doxygen/deal.II/namespaceDoFTools.html#a3540ceb577e65414bde1b6b14808da2c Best, Daniel On Thu, May 18, 2023 at 11:15 AM Wolfgang Bangerth wrote: > On 5/17/23 12:00, Mathieu wrote: > > > > I apologize for my simple question: > > Given a point, which is also a vertex (support point of FE_Q element), > in real > > coordinates, is there a function that allows me to query the global dof > index > > associated with that support point? > > > > I know all the cells at which these support point lives. > > One could then loop over the vertices of anyone of these cells, and call > > cell->vertex_dof_index(vertex,...) > > if vertex is the one that matches the support point. > > > > But I was wondering if there is a direct way to get the dof index of a > vertex > > based on the coordinates in real space. > > You can use DoFTools::map_dofs_to_support_points(). This gives you the map > "in > the wrong direction", but would probably still be useful for your case. > > Best > W. > > -- > > Wolfgang Bangerth email: bange...@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/78768f0a-b5a4-b689-adeb-e6a95e6da1a8%40colostate.edu > . > -- 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/CAOYDWbLt1EKRRu-H6T3Sxk2fYm5UZYquuBc_Mb_QuEJ_SvYw%2Bg%40mail.gmail.com.
[deal.II] deal.II Newsletter #254
Hello everyone! This is deal.II newsletter #254. It automatically reports recently merged features and discussions about the deal.II finite element library. ## Below you find a list of recently proposed or merged features: #15239: Fix the number of pyramid lines and one pyramid table. (proposed by drwells) https://github.com/dealii/dealii/pull/15239 #15238: Initialize and check an assertion in the NOX solver. (proposed by bangerth) https://github.com/dealii/dealii/pull/15238 #15237: Properly treat exceptions in a user callback in NOX. (proposed by bangerth; merged) https://github.com/dealii/dealii/pull/15237 #15236: Minor updates to the NOX solver documentation. (proposed by bangerth; merged) https://github.com/dealii/dealii/pull/15236 #15235: Fix Tpetra for newer Trilinos (proposed by masterleinad) https://github.com/dealii/dealii/pull/15235 #15233: fix bug in find_active_cell_around_point (proposed by richardschu) https://github.com/dealii/dealii/pull/15233 #15231: improve documentation for PETSc SNES (proposed by tjhei; merged) https://github.com/dealii/dealii/pull/15231 #15230: Patterns: Fix parsing lists of lists (proposed by masterleinad; merged) https://github.com/dealii/dealii/pull/15230 #15229: MatrixFree: Avoid setting face vectorization group boundary in serial (proposed by kronbichler) https://github.com/dealii/dealii/pull/15229 #15228: Avoid compiler warning (proposed by kronbichler; merged) https://github.com/dealii/dealii/pull/15228 #15226: Tpetra is not functional without sync() function (proposed by sebproell) https://github.com/dealii/dealii/pull/15226 #15223: Convert ARKode interfaces to conform to our callback policy. (proposed by bangerth) https://github.com/dealii/dealii/pull/15223 #15222: Minor doc improvement in step-6. (proposed by bangerth; merged) https://github.com/dealii/dealii/pull/15222 #15221: Update output of step-77. (proposed by bangerth; merged) https://github.com/dealii/dealii/pull/15221 #15220: fix some compiler warnings (proposed by richardschu; merged) https://github.com/dealii/dealii/pull/15220 #15219: fixed some compiler warnings in tria.cc (proposed by richardschu) https://github.com/dealii/dealii/pull/15219 #15218: Add a function for computing distances to unit cells. (proposed by drwells) https://github.com/dealii/dealii/pull/15218 #15217: Make MatrixFree evaluation kernels more flexible w.r.t. number type (proposed by kronbichler) https://github.com/dealii/dealii/pull/15217 #15216: SUNDIALS: remove outdated copy functionality (proposed by sebproell; merged) https://github.com/dealii/dealii/pull/15216 #15215: Face path for FEPointEvaluation (proposed by bergbauer; merged) https://github.com/dealii/dealii/pull/15215 #15214: Fix step-77 NOX test (proposed by sebproell; merged) https://github.com/dealii/dealii/pull/15214 #15212: Fix Kokkos compilation with GCC 13. (proposed by drwells; merged) https://github.com/dealii/dealii/pull/15212 #15211: FEPointEvaluation: Use template argument consistently (proposed by bergbauer; merged) https://github.com/dealii/dealii/pull/15211 #15210: Fix step-77: zero out residual before assembly (proposed by sebproell; merged) https://github.com/dealii/dealii/pull/15210 #15209: Integrate DG into non-nested MG (proposed by jh66637) https://github.com/dealii/dealii/pull/15209 #15208: Fix typo in documentation (proposed by kronbichler; merged) https://github.com/dealii/dealii/pull/15208 #15207: Tensor product kernels: Use the right type in some templated functions (proposed by kronbichler; merged) https://github.com/dealii/dealii/pull/15207 #15206: Fix warnings with gcc-13.1 (proposed by masterleinad; merged) https://github.com/dealii/dealii/pull/15206 #15205: CMake: create a compile_commands.json symlink (proposed by tamiko; merged) https://github.com/dealii/dealii/pull/15205 #15204: Adds create_with_same_size function (proposed by MarcelKoch) https://github.com/dealii/dealii/pull/15204 #15203: Updates Ginkgo solver wrappers (proposed by MarcelKoch) https://github.com/dealii/dealii/pull/15203 #15202: Adds wrappers for Ginkgo matrices (proposed by MarcelKoch) https://github.com/dealii/dealii/pull/15202 #15200: Use Kokkos to launch kernels in CUDAWrappers::MatrixFree (proposed by Rombur) https://github.com/dealii/dealii/pull/15200 #15199: Allow finding user-installed HDF5 (proposed by masterleinad; merged) https://github.com/dealii/dealii/pull/15199 #15197: Add an abstract base class ReadVector purely for vector access. (proposed by drwells) https://github.com/dealii/dealii/pull/15197 #15196: Avoid variables named 'import', part 2. (proposed by bangerth; merged) https://github.com/dealii/dealii/pull/15196 #15195: [WIP] Avoid the name 'import' as a function name. (proposed by bangerth) https://github.com/dealii/dealii/pull/15195 #15194: KINSOL: do not copy vectors internally (proposed by sebproell; merged) https://github.com/dealii/dealii/pull/15194 #15193:
Re: [deal.II] determine global dof index of a point
On 5/17/23 12:00, Mathieu wrote: I apologize for my simple question: Given a point, which is also a vertex (support point of FE_Q element), in real coordinates, is there a function that allows me to query the global dof index associated with that support point? I know all the cells at which these support point lives. One could then loop over the vertices of anyone of these cells, and call cell->vertex_dof_index(vertex,...) if vertex is the one that matches the support point. But I was wondering if there is a direct way to get the dof index of a vertex based on the coordinates in real space. You can use DoFTools::map_dofs_to_support_points(). This gives you the map "in the wrong direction", but would probably still be useful for your case. Best W. -- Wolfgang Bangerth email: bange...@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/78768f0a-b5a4-b689-adeb-e6a95e6da1a8%40colostate.edu.
