[deal.II] Re: Extra layer around mesh

Dear Daniel, Here is the small example you asked for. If you run it as it is, it will work. If you can 1 to 2 in line 107, it crashes. Im run dealii-8.2.1, if that is needed. Best, Joel On Tuesday, November 1, 2016 at 11:44:48 AM UTC+1, Daniel Arndt wrote: > > Joel, > > I have code that runs

[deal.II] Re: Extra layer around mesh

se> >]* *The violated condition was: * *pairs1.size() == pairs2.size()* *The name and call sequence of the exception was:* *ExcMessage ("Unmatched faces on periodic boundaries")* *Additional Information: * *Unmatched faces on periodic boundaries* How can I fix it? Best, Joel On F

[deal.II] Re: Vector-valued gradient of solution vector

Dear Daniel, I did try it, but I cant just change the local_velocity_values from a std::vector> to std::vector. Then the fe_values[velocities].get_function_values function doesnt compile. So I commented that part out, the results you see is from the code as it is, I used the for

[deal.II] Re: Vector-valued gradient of solution vector

Dear Daniel, To clarify what I am searching for, I would like to get the vector-valued solution into a format that matches a scalar. Either as a std::vector> or a std::vector, the format doesn't matter as long as it matches a scalar dofs. So I can manipulate the vector

[deal.II] Re: Vector-valued gradient of solution vector

tem_to_component_index > <https://www.dealii.org/developer/doxygen/deal.II/classFiniteElement.html#a27220a135402b96c7e6eecbb04acda56>. > > See the tutorials for examples on how to use them. > > Regards, > J-P > > On Tuesday, September 13, 2016 at 2:05:56 PM UTC+2, Joel Davidss

[deal.II] Re: Vector-valued gradient of solution vector

gt; > and could then extract the solution based on your custom quadrature rule. > > Regards, > J-P > > On Tuesday, August 30, 2016 at 5:32:13 PM UTC+2, Joel Davidsson wrote: >> >> Dear Daniel, >> >> I was just wondering if their exist a similar function &

[deal.II] Re: nearest neighbor

Thanks Daniel, I had confused the dof_index with the vertex_number. Now It works as expected. Best, Joel On Wednesday, September 7, 2016 at 3:47:21 PM UTC+2, Daniel Arndt wrote: > > Joel, > > In GridTools::find_cells_adjacent_to_vertex(dof_handler,center_id) you > need to give the number of

[deal.II] nearest neighbor

Dear all, If I have the index to a dof, is there any easy way I can get indexes to the nearest neighbor? Thanks, Joel -- 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

[deal.II] Re: Vector-valued gradient of solution vector

Dear Daniel, Thank you for a fast answer. I think I got the matrix and right-hand side correct now, please see test.cc for the example code. What I want to from the start was to get the vector-valued gradient from a scalar field, but I could not get this to work. See section "Do not work" in

[deal.II] Re: Extra layer around mesh

s has been useful to you. Apologies in advance if I don't > respond to any further messages - from tomorrow I'll be away from my PC for > a while. > > Regards, > J-P > > On Thursday, August 4, 2016 at 3:33:29 PM UTC+2, Joel Davidsson wrote: >> >> D

[deal.II] Extra layer around mesh

Dear all, Is there some way that one can refine the mesh in the following way, see attached file: I looking for a way to increase the size of the grid but keeping the original nodes fixed. If you look in the attached file, you see that an extra "layer" has formed around the original mesh