dolfin team mailing list archive

[Hatef Monajemi <hmonajem@xxxxxxxxxxxxxxxxxxx>]

Kristian Oelgaard wrote:
> Quoting Hatef Monajemi :
> 
> > Hi again.
> >
> > I ran into something strange in ffc which I have no answer for that.
> > referring to the following the following piece of code, I am 
defining a
> > Function named "Tau0" , however when in the main program I define
> > "L.Tau0" it gives me a strange error saying that "  ‘class
> > SWELinearForm’ has no member named ‘Tau0’ "  . Is it related 
to the use
> > of MixedElement???
> >
> > _________________________________________________________
> >
> > P2 = VectorElement("Lagrange", "triangle", 1)
> > P1 = FiniteElement("Lagrange", "triangle", 1)
> > P0 = FiniteElement("Discontinuous Lagrange", "triangle", 0)
> > TH = MixedElement([P2 ,P1, P0])
> >
> >
> > (v,w,q) = TestFunctions(TH)
> > (U,P,Tau) = TrialFunctions(TH)
> >
> >
> > U0=Function(P2)
> > P0=Function(P1)
> > Tau0=Function(P0)
> 
> 1) I think it is less strange when you realise that you have used the 
function
> P0 to initialise Tau0.

Yes, You are right and thanks for noticing that. However strange is 
that FFC in not reporting it as a bug!!!

> 2) Your example form does not compile because k and H are not defined.

I deleted some part of the code in order to make the email shorter. Of 
course I defined them in the program. 
> 
> 3) This post should probably have been sent to the FFC mailing list.

ffc did not give me any error. I got the error while doing "scons". so 
I thought that would be a better idea to send it to dolfin-dev.

Thanks anyway for noticing that mistake.
Hatef:)

> Kristian
> 
> > #Defining Normal components
> > n=FacetNormal("triangle")
> >
> >
> > a1 = (dot(U,v)+ dot(grad(P),v))*dx
> > a2 = (P*w - dot(U,grad(w)))*dx
> > a3 = Tau*q*dx
> >
> > L1 = (dot(U0,v)- dot(grad(P0),v)-.5*k*Tau0*dot(U0,v))*dx
> > L2 = (P0*w+dot(U0,grad(w)))*dx +w*H*k*dot(n,U0)*ds
> > L3 = Tau0*q*dx
> >
> > a=a1+a2+a3
> > L=L1+L2+L3
> > _____________________________________________
> >
> >
> >
> >
> >
> >
> >
> >
> >
> > Thanks,
> > Hatef
> >
> >
> > Hatef Monajemi wrote:
> > > Thanks Andy. That worked!
> > >
> > > Andy Ray Terrel wrote:
> > > > On Mon, Mar 16, 2009 at 7:32 PM, Andy Ray Terrel  wrote:
> > > > > 2009/3/16 Hatef Monajemi :
> > > > >> Hi,
> > > > >>
> > > > >> I would assume more than 2 elemnts should be O.K for ffc.
> > However,
> > > I am
> > > > >> getting some error while trying the following code:
> > > > >>
> > > > >> -----------------------------------
> > > > >>
> > > > >> P2 = VectorElement("Lagrange", "triangle", 1)
> > > > >> P1 = FiniteElement("Lagrange", "triangle", 1)
> > > > >> P3 = FiniteElement("Discontinuous Lagrange", "triangle", 0)
> > > > >> TH = (P2 + P1)+ P3
> > > > > ^^^^^^^^^^^^^^^^^^^^^
> > > > >
> > > > > To fix the trace back replace the above line with
> > > > > TH = MixedElement(P2,P1,P3)
> > > >
> > > > opps copied the wrong line =D, use:
> > > >
> > > > TH = MixedElement([P2,P1,P3])
> > > >
> > > >
> > > > >
> > > > > The above code creates a Mixed element like so ((v,w), q)
> > > > >
> > > > > -- Andy
> > > > >
> > > > >>
> > > > >> (v,w,q) = TestFunctions(TH)
> > > > >> (U,P,Tau) = TrialFunctions(TH)
> > > > >>
> > > > >>
> > > > >>
> > > > >> a1 =a1(U,P;v)
> > > > >> a2 = a2(U,P;w)
> > > > >> a3 = a3(Tau;q)
> > > > >>
> > > > >> a=a1+a2+a3
> > > > >> ----------------------------------------------
> > > > >>
> > > > >> does anyone know what is wrong here?
> > > > >>
> > > > >> This is the error I am geting:
> > > > >>
> > > > >>  ****need more than 2 values to unpack
> > > > >>
> > > > >> I want to basically add a constant number on every element
> > besides
> > > velocity
> > > > >> and pressure on nodal points , That's why I used the P3 type 
of
> > > > >> FiniteElement.
> > > > >>
> > > > >>
> > > > >>
> > > > >>
> > > > >> Thanks ,
> > > > >>
> > > > >> H! at! ! ef
> > > > >>
> > > > >> _______________________________________________
> > > > >> DOLFIN-dev mailing list
> > > > >> DOLFIN-dev@xxxxxxxxxx
> > > > >> http://www.fenics.org/mailman/listinfo/dolfin-dev
> > > > >>
> > > > >>
> > > > >
> > >
> > > --------------------------------------
> > > Hatef Monajemi
> > > Graduate Student
> > > Dept of Civil and Environmental Eng
> > > 2032 MC Building, Carleton University
> > > 1125 Colonel By Drive, Ottawa
> > > Ontario, Canada, K1S 5B6
> > > _______________________________________________
> > > DOLFIN-dev mailing list
> > > DOLFIN-dev@xxxxxxxxxx
> > > http://www.fenics.org/mailman/listinfo/dolfin-dev
> >
> > --------------------------------------
> >
> > Hatef Monajemi
> > Graduate Student
> > Dept of Civil and Environmental Eng
> > 2032 MC Building, Carleton University
> > 1125 Colonel By Drive, Ottawa
> > Ontario, Canada, K1S 5B6
> > _______________________________________________
> > DOLFIN-dev mailing list
> > DOLFIN-dev@xxxxxxxxxx
> > http://www.fenics.org/mailman/listinfo/dolfin-dev
> >
> 
>

--------------------------------------

Hatef Monajemi
Graduate Student
Dept of Civil and Environmental Eng
2032 MC Building, Carleton University 
1125 Colonel By Drive, Ottawa 
Ontario, Canada, K1S 5B6