[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Getfem-users] two different "NonLin" terms in assembly - allowed?
From: |
Yves Renard |
Subject: |
Re: [Getfem-users] two different "NonLin" terms in assembly - allowed? |
Date: |
Mon, 29 Oct 2007 14:45:30 +0100 |
User-agent: |
KMail/1.9.5 |
On Friday 26 October 2007 14:14, Andriy Andreykiv wrote:
> Dear Yves,
>
> Should there be any problem to have two different "NonLin" terms in
> one "comp" call in generic assembly?
Hello Andriy,
No, of course, it should work. Could you give me a little compilable program
such that I can find the bug ?
Yves?
> I'm getting an error:
>
> An error has been detected !!! |
> ============================================
> Error in dal_static_stored_objects.cc, line 205 :
> This object is not stored : 0x8f3a28c typename: N6getfem13mat_elem_typeE
>
>
> The code is:
>
> getfem::NormalTerm norm(mf_lambda,deformed_mesh,ls,1);
> getfem::NormalTerm R_matrix(mf_lambda,deformed_mesh,ls,2);
> getfem::generic_assembly assem;
> assem.set(
> .
> .
>
> "KuuOm1=0.5*sym(comp(vGrad(#2).Grad(#3).Grad(#3).vBase(#2).NonLin$2
> (#1).NonLin$1(#1))(:,j,i,k,m,l,m,:,n,i,j,n).F(k).F(l));"
> );
> assem.push_mi(mim_line);
> assem.push_mf(mf_lambda); //scalar mf on line
> assem.push_mf(mf_U_line); // vectorial mf on line
> assem.push_mf(mf_interpolate); // interpolated scalar fem on line
> assem.push_nonlinear_term(&norm); // vector ..
> assem.push_nonlinear_term(&R_matrix); //matrix ..
> .
> .
> assem.assembly();
>
> I tested and when the above KuuOm1 expression has either only NonLin$1(#1),
> or NonLin2$(#1) it works fine. While combining them in one "comp" call
> gives the above error.
>
> Is there some remedy for this?
>
> Thanks in advance,
> Andriy
> __________________________________________
> Andriy Andreykiv (PhD, MSc)
> Delft University of Technology
> Faculty of Mechanical Engineering
> Fundamentals of Microsystems
> Mekelweg 2
> 2628 CD Delft
> The Netherlands
>
> E-mail : address@hidden
> Tel. : (31) 15 2786818
> Fax. : (31) 15 2789475
> www : http://www-tm.wbmt.tudelft.nl/~andrico
> private: (31) 6 47376804
>
> _______________________________________________
> Getfem-users mailing list
> address@hidden
> https://mail.gna.org/listinfo/getfem-users
--
Yves Renard (address@hidden) tel : (33) 04.72.43.87.08
Pole de Mathematiques, fax : (33) 04.72.43.85.29
INSA de Lyon, Universite de Lyon
20, rue Albert Einstein
69621 Villeurbanne Cedex, FRANCE
http://math.univ-lyon1.fr/~renard
---------