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, it does not solve the system of linear equations. It substitutes all zeros into the solution vector and then makes solution(0) = 1. 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), looks just as what I would expect. The 2-nd order, however, looks just like the 0-th order. I opened the files in spreadsheet and indeed – the shape functions of the 0-th and 2-nd order are the same. I ruled out the possibility of mistake (the program prints fe.degree) and began worrying. I think VectorTools::point_value is to blame probably… Is there a replacement for it? Best, John On Monday, April 19, 2021 at 4:59:25 AM UTC+2 Wolfgang Bangerth wrote: > 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 > > < > https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Fwww.dealii.org%2Fcurrent%2Fdoxygen%2Fdeal.II%2FclassDoFHandler.html&data=04%7C01%7CWolfgang.Bangerth%40colostate.edu%7C39054f37290545641e5608d902d486d0%7Cafb58802ff7a4bb1ab21367ff2ecfc8b%7C0%7C0%7C637543934089426552%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&sdata=cfNPiGUOk%2BsIlssgAev6wBAmMiD%2BsZHYrl32blPab0E%3D&reserved=0> > > > > and nodal vector fe_function at the given point point, and return the > (vector) > > value of this function through the last argument. > > John -- have you tried? Oftentimes, just giving it a try is going to > answer > your question :-) > Best > W. > > > -- > ------------------------------------------------------------------------ > Wolfgang Bangerth email: [email protected] > 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 [email protected]. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/3cfc8137-d48e-4b98-8870-3dab4f64365en%40googlegroups.com.
