Re: [Getfem-users] generic_as​​​sembly of complex_ty​​​pe

[Yves Renard]

Sorry, and of course the same for the result vector :

assem3.push_vec(gmm::real_part(B));
assem3.push_vec(gmm::imag_part(B));

Yves.


On mardi 3 janvier 2012, Yves Renard wrote:
> Dear wan rong Joung,
> 
> Unfortunately, the generic assembly of Getfem++ do not accept directly a
> complex type vector or matrix. You have to push both the real and imaginary
> part as follows :
> 
> assem3.push_data(gmm::real_part(CCGik));
> assem3.push_data(gmm::real_part(Ui));
> assem3.push_data(gmm::imag_part(CCGik));
> assem3.push_data(gmm::imag_part(Ui));
> 
> and take it into account in your assembly formula of course.
> 
> Yves.
> 
> On mardi 3 janvier 2012, wan rong Joung wrote:
> > Dear all,
> > 
> > I have used “GetFem++” for several month and I appreciate it very much.
> > 
> > I have some problem about <complex_type>. I set some parameters for
> > complex_type(without brick). However I get some error from -------
> > 
> > ------------------à getfem::generic_assembly
> > 
> > 
> > 
> > The code following:
> > 
> > std::vector<scalar_type> H,CCG,CCGK2;
> > 
> > std::vector<complex_type> CCGik,CCGikUi,Ui;
> > 
> >  getfem::generic_assembly assem3;
> >  
> >  assem3.push_mi(mim);
> >  
> >  assem3.push_mf(mf_u);
> >  
> >  assem3.push_mf(mf_data);
> >  
> >  assem3.push_data(CCGik); // bug
> >  
> >  assem3.push_data(Ui); // bug
> >  
> >  assem3.push_vec(B);
> >  
> >  assem3.set("a=data$1(#2);"
> >  
> >             "b=data$2(#2);"
> >             
> >             "V(#1)+=2*comp(Base(#1).Base(#2))(:,k).a(k).b(k);");
> >  
> >  assem3.assembly(INCIDENT_BOUNDARY);
> > 
> > This caused the following error:
> > 
> > gmm::strongest_value_type<gmm::cs_vector_ref<const double*, const
> > unsigned int*, 0>, std::vector<std::complex<double> > >::value_type' to
> > 'bgeot::scalar_type' in assignment|
> > 
> > const std::complex<double>' to 'bgeot::scalar_type' in assignment|
> > 
> > 
> > 
> > I don’t understand this error information. Any comment for you will be
> > highly appreciated.


-- 

  Yves Renard (address@hidden)       tel : (33) 04.72.43.87.08
  Pole de Mathematiques, INSA-Lyon             fax : (33) 04.72.43.85.29
  20, rue Albert Einstein
  69621 Villeurbanne Cedex, FRANCE
  http://math.univ-lyon1.fr/~renard

---------