Re: [deal.II] Re: Why does the L2-norm increase when the number of nodes becomes too large?

> On 21-Dec-2017, at 8:04 AM, Praveen C wrote: > > > >> On 20-Dec-2017, at 11:52 PM, jie liu > > wrote: >> >> Hi Luca, >> >> On Wednesday, December 20, 2017 at 5:20:43 PM UTC+1, Luca Heltai wrote: >> >> What’s the tolerance of your solver? Is it a direct or an

Re: [deal.II] Re: Why does the L2-norm increase when the number of nodes becomes too large?

> On 20-Dec-2017, at 11:52 PM, jie liu wrote: > > Hi Luca, > > On Wednesday, December 20, 2017 at 5:20:43 PM UTC+1, Luca Heltai wrote: > > What’s the tolerance of your solver? Is it a direct or an iterative solver? > If it is an iterative solver, is the tolerance right around 1e-12? > > Th

Re: [deal.II] Re: Problem with the use of FE_Nedelec with a mesh file generated from Gmsh

Bruno, Okay, I will do more experiments and checks on my code and Ross' code to double confirm that the bug is not my code. Once I confirm this, it will be great if you can contact the people working on the PR to see whether my mesh can work with the new implementation. I'll contact you again

Re: [deal.II] Re: Problem with the use of FE_Nedelec with a mesh file generated from Gmsh

Jianan, 2017-12-20 17:21 GMT-05:00 Jianan Zhang : > 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

Re: [deal.II] Re: Problem with the use of FE_Nedelec with a mesh file generated from Gmsh

Hi, Bruno, 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. Best Jianan 在 2017年12月20日星期三 UTC-5下午3:26:21，Bruno Turcksin写道： > > Jian

[deal.II] Re: Problem with Dirichlet BCs for FE_Nedelec for 3D?

Dear Ross, I am using dealii to solve the curl-curl equation: curl(mu^(-1)curl(E)) + (-omega^2*epsilon+j*omega*sigma)*E=0, with boundary conditions: n x E = n X F , (Dirichlet) or n X (curl(E)) = n X (curl(F)) (Neumann) where

Re: [deal.II] Re: Problem with the use of FE_Nedelec with a mesh file generated from Gmsh

Jianan, 2017-12-20 14:54 GMT-05:00 Jianan Zhang : > > 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

[deal.II] Re: Problem with the use of FE_Nedelec with a mesh file generated from Gmsh

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 s

Re: [deal.II] Re: Why does the L2-norm increase when the number of nodes becomes too large?

Hi Luca, On Wednesday, December 20, 2017 at 5:20:43 PM UTC+1, Luca Heltai wrote: > > What’s the tolerance of your solver? Is it a direct or an iterative > solver? If it is an iterative solver, is the tolerance right around 1e-12? > > Thank you for your reply. I adopted both iterative and direct

Re: [deal.II] Re: Why does the L2-norm increase when the number of nodes becomes too large?

M… I’m actually not sure it’s a bug… The point where the error starts growing is at 10e-12. This is basically machine precision, considering the fact that you are computing the L2 error. After you have reached machine precision, any addition you make to that, is just roundoff error addin

Re: [deal.II] DOF Ordering in DofHandler

Hello, First of all thank you so much for all the help. I printed coordinates and global dof no. for each node from DofHandler. I observed that global dof numbers are not the same as the global numbers in my input mesh file. So I compared the two and developed an association between old and new

[deal.II] Re: Problem with the use of FE_Nedelec with a mesh file generated from Gmsh

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 differe

[deal.II] Re: Why does the L2-norm increase when the number of nodes becomes too large?

Dear Jie, I'm with Bruno, you definately have a bug in your code. Maybe the application of the Dirichlet boundary is incorrect. When you use the 1d grid generation the left and right boundary nodes are coloured by different numbers. This is a different behaviour compared to 2d or 3d grid genera