Hi, Bruno, Thanks for the reply. Of course, I am not expecting someone to go through my code and debug it for me. I post the codes just in case that someone is interested or someone can take it as a reference when they try to use FE_Nedelec.
What I am really expecting is that someone who has some experience of using FE_Nedelec to solve EM problems can give me some general advice. For example, can we use FE_Nedelec on a mesh imported from other mesh generators instead of the one generated from dealii itself? Or does dealii of the current version support this kind of mesh (obtained from Gmsh software and after transformed into hexahedrals) well? and so on and so forth. If someone can give me this kind of hint it will help a lot. Because the actual problem may be due to the not good support of non-rectangular meshes when using FE_Nedelec, and in that case my debugging is like a wast of time. Anyway, your suggestions make sense, I will try exact solutions of different forms and see when it works and when it does not work. Meanwhile, I will search previous works to see if anyone has met similar problems as me. Best Jianan 在 2017年12月20日星期三 UTC-5上午9:34:23，Bruno Turcksin写道： > > Jianan, > > On Wednesday, December 20, 2017 at 12:56:11 AM UTC-5, Jianan Zhang wrote: > >> >> I've spent like ten days finding the reason, and I am almost 100% sure >> that my exact solution is given correctly and the assembly is correct. The >> only thing I changed is the exact solution, and the difference between his >> and mine is that his is only z-component dependent while mine depend on >> both x and z. >> >> I post my codes and the vtk files for anyone who are interested. I am >> having a hard time on debugging this and any suggestions and help is much >> appreciated. >> > You cannot expect us to go through your code and debug it for you. We all > have our own codes to debug... Now you say that the only difference between > Ross' problem and yours is the dependence in both x and z. Now you have to > determine if the problem is in the x component or in the x-z coupling. What > happens if the solution depends only on x? only on y? What about y-z > instead of x-z? > > 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 [email protected]. For more options, visit https://groups.google.com/d/optout.
