getfem-users
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

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


From: Andriy Andreykiv
Subject: [Getfem-users] two different "NonLin" terms in assembly - allowed?
Date: Fri, 26 Oct 2007 14:14:15 +0200
User-agent: KMail/1.9.6

Dear Yves,

        Should there be any problem to have two different "NonLin" terms in 
one "comp" call in generic assembly?
        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



reply via email to

[Prev in Thread] Current Thread [Next in Thread]