Bartoz Would it be possible to send a tar.gz of the demos to me so that I can also have a look without messing around in the source tree.
Also, just keep reminding the guys about your patch - I know they are sometimes a little busy :) Cheers Evan On Wed, Apr 1, 2009 at 11:39 PM, Bartosz Sawicki <[email protected]>wrote: > Even if core team doesn't like my demo, could some point me thats wrong > with the following code. It's so simple and similar to the cpp version, that > I completely don't understand why the results are different. > > =============== > from dolfin import * > > mesh = UnitSphere(10) > > PN = FunctionSpace(mesh, "Nedelec", 0) > P1 = VectorFunctionSpace(mesh, "CG", 1) > > v0 = TestFunction(PN) > u0 = TrialFunction(PN) > v1 = TestFunction(P1) > u1 = TrialFunction(P1) > > dBdt = Function(P1, ("0.0", "0.0", "1.0")) > zero = Function(P1, ("0.0", "0.0", "0.0")) > T = Function(PN) > J = Function(P1) > > class DirichletBoundary(SubDomain): > def inside(self, x, on_boundary): > return on_boundary > > bc = DirichletBC(P1, zero, DirichletBoundary()) > > # Eddy currents equation (using potential T) > Teqn = (dot(rot(v0), rot(u0))*dx, -dot(v0, dBdt)*dx) > Tproblem = VariationalProblem( Teqn[0], Teqn[1], bc) > T = Tproblem.solve() > > # Current density equation > Jeqn = (dot(v1, u1)*dx, dot(v1, curl(T))*dx) > Jproblem = VariationalProblem( Jeqn[0], Jeqn[1]) > J = Jproblem.solve() > > plot(J) > > interactive() > ============== > > cheers, > BArtek > > > > > > Evan Lezar wrote: > >> That is true. And more examples are never a bad thing :) >> >> On 3/31/09, Bartosz Sawicki <[email protected]> wrote: >> >>> Yes, you are right, but your waveguide is eigenvalue problem and need >>> SLEPc to be installed. I think that it would be nice to have some basic >>> example which doesn't create extra requirements. >>> Just for dolfin beginners and for testing purposes. >>> >>> BArtosz >>> >>> >>> Evan Lezar wrote: >>> >>>> Hi >>>> >>>> The waveguide demo (pde/waveguide) uses the curl-curl formulation and >>>> Nedelec elements for electromagnetic eigenvalue problems. I am working >>>> on >>>> a >>>> couple of other demos for full wave electromagnetic problems. >>>> >>>> Evan >>>> >>>> >>>> 2009/3/30 Bartosz Sawicki <[email protected]> >>>> >>>> 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. >>>>> >>>>> cheers, >>>>> BArtosz >>>>> >>>>> _______________________________________________ >>>>> DOLFIN-dev mailing list >>>>> [email protected] >>>>> http://www.fenics.org/mailman/listinfo/dolfin-dev >>>>> >>>>> >>>>> >>> >> >
_______________________________________________ DOLFIN-dev mailing list [email protected] http://www.fenics.org/mailman/listinfo/dolfin-dev
