Re: [Getfem-users] two different "NonLin" terms in assembly - allowed?

[Yves Renard]

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

---------