ffc team mailing list archive

[Marie Rognes <meg@xxxxxxxxx>]

On 17. nov. 2010 17:50, Patrick Riesen wrote:
> Ok, i wanted to file a bug report just now, but i got two
>
> 'time out errors' from launchpad services. i'll try again later.
>
> btw. launchpad told me about a related (fixed) bug:
>     
> #601965 compute_form_arities crashes in ListTensor
>
>

That bug was related to the same functionality, but a different problem
(which is fixed).
So, feel free to open a new one.

--
Marie


> anyway, many thanks and best regards,
> patrick
>
> On 11/17/2010 05:28 PM, Marie Rognes wrote:
>> On 17. nov. 2010 16:15, Patrick Riesen wrote:
>>> 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
>>> -----------------------------------------------------------------
>>>
>>>
>>
>>
>> Thanks. I see what is going on, but do not quite know how to fix it
>> immediately.
>> Please add a bug report (for UFL, which is where the problem is), and
>> I'll think about how to fix it.
>>
>> Here is a simpler example btw:
>>
>> from ufl.algorithms import compute_form_arities
>>
>> cell = interval
>> FE = FiniteElement("CG", cell, 2)
>> u = Coefficient(FE)
>>
>> e = as_matrix([[u, 0.0],
>>                 [0.0, 0.0]])
>> L = tr(e)*dx
>>
>> U = TrialFunction(FE)
>> a = derivative(L, u, U)
>> compute_form_arities(a)
>>
>> -- 
>> Marie
>>
>>
>> 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
>>>
>>>
>>> _______________________________________________
>>> Mailing list: https://launchpad.net/~ffc
>>> Post to     : ffc@xxxxxxxxxxxxxxxxxxx
>>> Unsubscribe : https://launchpad.net/~ffc
>>> More help   : https://help.launchpad.net/ListHelp
>>
>