ffc team mailing list archive

[Patrick Riesen <priesen@xxxxxxxxxxxxxxx>]

Nope, won't do differnt with Zero().
Here is a simple complete program:
-----------------------------------------------------------------
cell = interval
FE = FiniteElement("CG", cell, 2)

v = TestFunction(FE)
U = TrialFunction(FE)

u = Coefficient(FE)

K = Constant(cell)
c = Constant(cell)

e = as_matrix([[0.5*u.dx(0)**2, 0.5*u.dx(0)],
               [0.5*u.dx(0)   ,     Zero()]])

e_dev = e - 0.5*tr(e)*Identity(2)

Ie = 0.5*tr(e_dev*e_dev)
beta = K*exp(-c*Ie)

L = v.dx(0)*beta*0.5*u.dx(0)*dx    # ok
a = derivative(L, u, U)            # it fails here
-----------------------------------------------------------------


regards,
patrick

On 11/17/2010 02:43 PM, Anders Logg wrote:
> I don't know what that means, but I assume the 0.0 is the problem.
>
> Does it work if you write Zero()?
>
> --
> Anders
>
>
>
> On Wed, Nov 17, 2010 at 02:36:54PM +0100, Patrick Riesen wrote:
> > hello,
> > i'm trying to update my code for the most recent fenics development
> > versions.
> >
> > compiling+demos for ffc/ufl/dolfin etc. is all ok.
> >
> > now my code has a problem, i have the following in the ufl-file:
> >
> > ****
> > e = as_matrix([[0.5*u.dx(0)**2, 0.5*u.dx(0)],
> >                 [0.5*u.dx(0)   ,        0.0]])
> >
> > e_dev = e - 0.5*tr(e)*Identity(2)
> >
> > Ie = 0.5*tr(e_dev*e_dev)
> > beta = K*exp(-c*Ie)
> >
> > ****
> >
> > note:
> > the variable u is a coefficient from a mixed element as
> > uw = Coefficient(ME)
> > u, w = split(uw)
> >
> >
> > compiling with ffc now terminates after the compiler stage 1 and gives:
> >
> > All components of a list tensor most provide same arguments
> >
> > *** FFC: All components of a list tensor most provide same arguments
> > *** FFC: To get more information about this error, rerun FFC with --verbose.
> >
> > what does this mean? how do i need to fix my syntax above?
> >
> > regards,
> > patrick
> >
> > _______________________________________________
> > Mailing list: https://launchpad.net/~ffc
> > Post to     : ffc@xxxxxxxxxxxxxxxxxxx
> > Unsubscribe : https://launchpad.net/~ffc
> > More help   : https://help.launchpad.net/ListHelp