Hi, Bruno, Hope you had a good holiday. I have done more checks on my code these days and found that the problem truly lies in the use of FE_Nedelec for non-rectangular meshes.
To illustrate, I tried to solve the problem with Ross's exact solution (but with a larger value of omega) on the mesh generated by Gmsh, Then the results obtained are obviously wrong. But when I try the same exact solution on the mesh created by dealii, the results are correct. This demonstrates that dealii of the current version truly has not so good support for non-rectangular meshes. So, can you help me contact the people working on the PR to try this out to see if my code works with the new implementation? Many thanks Jianan 在 2017年12月20日星期三 UTC-5下午5:35:11，Bruno Turcksin写道： > > Jianan, > > 2017-12-20 17:21 GMT-05:00 Jianan Zhang <adamzh...@gmail.com <javascript:> > >: > >> Thanks for the prompt reply. I will try the implementation as you >> suggested, hope I can make it work with it. Otherwise, I will just go with >> brick elements, although the EM structure I can model may be limited. >> > According to the comments in that PR it works on non-standard meshes. If > your code doesn't work, I would contact the people working on the PR to > make sure that your mesh works with the new implementation but don't > exclude too quickly that the bug is your code ;-) > > Best, > > Bruno > -- 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. For more options, visit https://groups.google.com/d/optout.
