Re: [deal.II] point_value FE_Nedelec

2021-04-22 Thread Wolfgang Bangerth
On 4/22/21 7:11 PM, John Smith wrote: I have made a silly mistake. Thank you for pointing out this mistake! The shape functions indeed look fine! I apologize for wasting time your time like that … No worries -- we've all made these sorts of mistakes! Best W. --

Re: [deal.II] point_value FE_Nedelec

2021-04-22 Thread John Smith
Dear Wolfgang, I have made a silly mistake. Thank you for pointing out this mistake! The shape functions indeed look fine! I apologize for wasting time your time like that … Best wishes, John On Friday, April 23, 2021 at 1:21:42 AM UTC+2 Wolfgang Bangerth wrote: > > > For the second

Re: [deal.II] point_value FE_Nedelec

2021-04-22 Thread Wolfgang Bangerth
For the second question: You are absolutely right. I do expect two edge elements of different degree to be different. Right, but that does not mean that every shape function of FE_Nedelec(1) must be different from all of the FE_Nedelec(0) shape functions. Some of them may in fact be the

Re: [deal.II] point_value FE_Nedelec

2021-04-19 Thread John Smith
Dear Wolfgang, I do not think it is a bug. I think it may very well be that VectorTools::point_value it is not meant to be used with FE_Nedelec. The first question: in this instance I do not use DataOut at all. I use VectorTools::point_value to get data point-by-point and save it into a

Re: [deal.II] point_value FE_Nedelec

2021-04-19 Thread Wolfgang Bangerth
On 4/18/21 10:50 PM, John Smith wrote: Then I have fed this solution together with the dof handler to the VectorTools::point_value on a grid of regularly spaced points. The results were written into a *.csv file. Then I have made the vector plot out of it. The 0-th order, FE_Nedelec<2> fe(0),

Re: [deal.II] point_value FE_Nedelec

2021-04-18 Thread John Smith
Yes, I have. Thus the question. I would like to have a simple silly test to see that the FE_Nedelec elements are what I expect them to be. So, I made a simple program similar to step-3 that sets up a test. The program uses FE_Nedelec<2>. It does all the steps similarly to step-3. However,

Re: [deal.II] point_value FE_Nedelec

2021-04-18 Thread Wolfgang Bangerth
On 4/18/21 7:43 PM, John Smith wrote: Is VectorTools::point_value supposed to work with FE_Nedelec? The documentation looks very promising: Evaluate a possibly vector-valued finite element function defined by the given DoFHandler

[deal.II] point_value FE_Nedelec

2021-04-18 Thread John Smith
Hello, Is VectorTools::point_value supposed to work with FE_Nedelec? The documentation looks very promising: Evaluate a possibly vector-valued finite element function defined by the given DoFHandler and nodal vector