Re: [deal.II] FEValues and get_function_values when support points = quadrature points

2016-09-30 Thread Wolfgang Bangerth
On 09/30/2016 09:21 PM, Praveen C wrote: But,... if FEValues constructor could take an optional boolean argument where user can say whether support points == quadrature points, or not. For FESystem, one has to pass a vector of boolean. I think you already see all of the contortions one has to

Re: [deal.II] FEValues and get_function_values when support points = quadrature points

2016-09-30 Thread Praveen C
Hello Wolfgang, When I use a FESystem, I have to write a couple of lines of code to get function values. I thought if FEValues could figure this, it would be nice, and make my code compact. But making this automatic inside FEValues seems difficult, since all FE in an FESystem may not even have

Re: [deal.II] FEValues and get_function_values when support points = quadrature points

2016-09-30 Thread Wolfgang Bangerth
On 09/30/2016 08:59 PM, Praveen C wrote: In many cases, we use same nodes as support points for FE space and quadrature. In this case, FEValues.get_function_values can just read off the function values without performing a sum over all basis functions. Does FEValues actually recognise and make

[deal.II] FEValues and get_function_values when support points = quadrature points

2016-09-30 Thread Praveen C
Dear all In many cases, we use same nodes as support points for FE space and quadrature. In this case, FEValues.get_function_values can just read off the function values without performing a sum over all basis functions. Does FEValues actually recognise and make this optimization ? If not, is it

[deal.II] C1 elements for 4th order derivatives

2016-09-30 Thread Jinhyun Choo
Hello All, I'm interested in solving a PDE involving a 4th order derivative (specifically, the 4th order phase-field model proposed by Borden et al. 2014). I believe such equations require C1 continuous elements but this is

Re: [deal.II] Parallel hp::DoFHandler

2016-09-30 Thread Wolfgang Bangerth
On 09/30/2016 03:40 AM, Andreas Krämer wrote: can you already foresee when there will be a version of hp::DoFHandler that works with parallel::distributed::Triangulation? Andreas -- it's been on my list of things to do for a long time, and I think it has been on that of others as well. It's

Re: [deal.II] a few questions about matrix-free

2016-09-30 Thread Denis Davydov
Hi Martin, > On 30 Sep 2016, at 16:49, Martin Kronbichler > wrote: > > There are two possibilities for doing things. One is like you say, you loop > over the components in the input vector or create two FEEvaluation objects > and let them do similar things. Note

Re: [deal.II] a few questions about matrix-free

2016-09-30 Thread Martin Kronbichler
Hi Denis, I will put the description in the header of the file. PR underway. The best documentation gets written when explaining things to others. :-) > actually it’s not about components, it is about values at quadrature > points for a vector of solution. > Say I have two vectors representing

[deal.II] Re: Parallel hp::DoFHandler

2016-09-30 Thread Bruno Turcksin
Andreas, On Friday, September 30, 2016 at 5:40:43 AM UTC-4, Andreas Krämer wrote: > > can you already foresee when there will be a version of hp::DoFHandler > that works with parallel::distributed::Triangulation? > It is in our todo list but nobody has time to implement it right now. You will

Re: [deal.II] Re: how to get 3D structured mesh from rotated 2D mesh in dealii

2016-09-30 Thread Vinetou Incucuna
Hello Daniel, I have fixed the problem using Your proposed algorithm of mesh creation. The strange behaviour of velocity-pressure field was really caused by the bad quality of the mesh (possibly doubled vertices at the boundary between merged entities). I want to ask a question, which could be

Re: [deal.II] a few questions about matrix-free

2016-09-30 Thread Denis Davydov
Hi Martin, both of your posts helped a lot to clarify things. Thanks a lot! Would you mind adding the discussion on using MatrixFree with different operators to the documentation of MatrixFree class? Essentially copy-pasting of what you wrote in the two e-mails related to update flags and alike

[deal.II] Parallel hp::DoFHandler

2016-09-30 Thread Andreas Krämer
Dear developers, can you already foresee when there will be a version of hp::DoFHandler that works with parallel::distributed::Triangulation? Best, Andreas -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see