On 11/29/2011 02:50 PM, Kent-Andre Mardal wrote:
> On 28 November 2011 23:26, Marie Rognes<email address hidden>
> wrote:
>
>> On 28. nov. 2011, at 23:15, Kent-Andre Mardal<email address hidden> wrote:
>>
>>> On 28 November 2011 22:07, Anders E. Johansen
>>> <email address hidden>wrote:
>>>
>>>> Kristian, I'm using the optimization options from the Hyperelasticity
>>>> demo. But with
>>>>
>>>> parameters["form_compiler"]["cpp_optimize"] = True
>>>>
>>>> FFC won't compile when computing 2PK (at least not in an hour). Using
>>>> the same bcs, parameters, etc. as in my first post, this code reads
>>>>
>>>> # Kinematics
>>>> I = Identity(3)
>>>> F = variable(I + grad(u))
>>>> C = variable(F.T*F)
>>>> E = 0.5*(C - I)
>>>> J = sqrt(det(C))
>>>>
>>>> # Material law
>>>> f=0; s=1; n=2
>>>> W = (bff*E[f,f]**2
>>>> + bxx*(E[n,n]**2 + E[s,s]**2 + E[n,s]**2)
>>>> + bfx*(E[f,n]**2 + E[n,f]**2 + E[f,s]**2 + E[s,f]**2))
>>>> psi = 0.5*K*(exp(W) - 1) + Cc*(J*ln(J) - J + 1)
>>>>
>>>> # 2PK
>>>> S = 2*diff(psi, C)
>>>>
>>>> # 1PK
>>>> P = F*S
>>>> #P = diff(psi, F)
>>>>
>>>> When computing 1PK directly, there is no problem at all.
>>>>
>>>> The reason I want to use FFC, is that SFC is not in active development,
>>>> and have lack of developers, users and features compared to FFC, I have
>>>> been told.
>>>>
>>>>
>>> :)
>>> We will make sure SFC works also in the future. But at present, SFC does
>>> what we need...
>>> If you have specific needs, let us know.
>>>
>> I tried filing a bug a while back but it did not seem to receive much
>> attention. Is there another channel for letting you know?
>>
>> --
>> Marie
>>
>>
> Sorry about that. Martin has fixed it now.
On 11/29/2011 02:50 PM, Kent-Andre Mardal wrote: "form_compiler" ]["cpp_ optimize" ] = True
> On 28 November 2011 23:26, Marie Rognes<email address hidden>
> wrote:
>
>> On 28. nov. 2011, at 23:15, Kent-Andre Mardal<email address hidden> wrote:
>>
>>> On 28 November 2011 22:07, Anders E. Johansen
>>> <email address hidden>wrote:
>>>
>>>> Kristian, I'm using the optimization options from the Hyperelasticity
>>>> demo. But with
>>>>
>>>> parameters[
>>>>
>>>> FFC won't compile when computing 2PK (at least not in an hour). Using
>>>> the same bcs, parameters, etc. as in my first post, this code reads
>>>>
>>>> # Kinematics
>>>> I = Identity(3)
>>>> F = variable(I + grad(u))
>>>> C = variable(F.T*F)
>>>> E = 0.5*(C - I)
>>>> J = sqrt(det(C))
>>>>
>>>> # Material law
>>>> f=0; s=1; n=2
>>>> W = (bff*E[f,f]**2
>>>> + bxx*(E[n,n]**2 + E[s,s]**2 + E[n,s]**2)
>>>> + bfx*(E[f,n]**2 + E[n,f]**2 + E[f,s]**2 + E[s,f]**2))
>>>> psi = 0.5*K*(exp(W) - 1) + Cc*(J*ln(J) - J + 1)
>>>>
>>>> # 2PK
>>>> S = 2*diff(psi, C)
>>>>
>>>> # 1PK
>>>> P = F*S
>>>> #P = diff(psi, F)
>>>>
>>>> When computing 1PK directly, there is no problem at all.
>>>>
>>>> The reason I want to use FFC, is that SFC is not in active development,
>>>> and have lack of developers, users and features compared to FFC, I have
>>>> been told.
>>>>
>>>>
>>> :)
>>> We will make sure SFC works also in the future. But at present, SFC does
>>> what we need...
>>> If you have specific needs, let us know.
>>>
>> I tried filing a bug a while back but it did not seem to receive much
>> attention. Is there another channel for letting you know?
>>
>> --
>> Marie
>>
>>
> Sorry about that. Martin has fixed it now.
Great, thanks!
--
Marie