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

[wan rong Joung]

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.

--

Wan Rong Joung

Department of Marine Environmental Engineering
National Kaohsiung Marine University
Kaohsiung, Taiwan