dolfin team mailing list archive

[Bartosz Sawicki <sawickib@xxxxxxxxxxxxx>]

tgz has just been sent to your mailbox.

I have already submitted several patches. So I know that we have to be 
patient. When next release is close, I don't expect quick answers. :)

cheers,
BArtek


Evan Lezar wrote:
> 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 <sawickib@xxxxxxxxxxxxx>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 <sawickib@xxxxxxxxxxxxx> 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 <sawickib@xxxxxxxxxxxxx>
> > > > >
> > > > >  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
> > > > > > DOLFIN-dev@xxxxxxxxxx
> > > > > > http://www.fenics.org/mailman/listinfo/dolfin-dev