>> I've written simple demo of curl-curl equation. It is based on >> eddy-currents phenomena in low conducting materials. So far, in the >> demo directory, there are no examples of this kind of PDE, where >> Nedelec elements are used. >> If you found it useful, please consider placing it in the repository. >> >> The patch which I'm sending, contains cpp and python code. Cpp works >> fine, but demo.py gives wrong results. I have no idea, whats wrong in >> the python code. Please, take a look, this is probably bug. >> > > (1) There seems to be a minus sign discrepancy between the two demos in > the EddyCurrents part.
You are right, there should be minus in both demos. > (2) In the python demo, replace > > bc = DirichletBC(P1, zero, DirichletBoundary()) > > by > bc = DirichletBC(PN, zero, DirichletBoundary()) Oh, good! How could I miss that. This is the bug which I was looking for. Now both examples works perfectly. Thank you very much. Do you want me to send improved patch? Or you can change this two characters? cheers, _______________________________________________ DOLFIN-dev mailing list [email protected] http://www.fenics.org/mailman/listinfo/dolfin-dev
