Re: [Getfem-users] [Getfem-commits] What is the best way to add parameter for each direction in term?

[Konstantinos Poulios]

I guess componentwise multiplication

getfem::add_linear_term(model, mim, "(rho.*u).Test_u", TOP);

with rho as a vector should work as well.

BR

Kostas

On Mon, Feb 11, 2019 at 5:27 PM Yves Renard <address@hidden> wrote:

Dear Tetsuo,

If I understand well, you need probably a matrix, not a vector

model.add_fixed_size_data("rho", bgeot::multi_index(2, 2));
model.set_real_variable("rho")[0] = rho_for_x;
model.set_real_variable("rho")[2] = rho_for_y;

getfem::add_linear_term(model, mim, "(rho*u).Test_u", TOP);


Best Regards,

Yves




Le 11/02/2019 à 09:15, Tetsuo Koyama a écrit :
> Dear GetFEM++ Project
>
> I'm sorry, code was wrong.
>
> model.add_fixed_size_data("rho", 2);
> model.set_real_variable("rho")[0] = rho_for_x;
> model.set_real_variable("rho")[1] = rho_for_y;
>
> getfem::add_linear_term(model, mim, "rho*u.Test_u", TOP);
>
> Could you anser it? Thank you.
>
> Best Regard Tetsuo
>
> 2019年2月11日(月) 17:03 Tetsuo Koyama <address@hidden>:
>> Dear GetFEM++
>>
>> What is the best way to add parameter for each space direction in term?
>> I want to use fixed_size_data as a coefficient of the weak form
>> language, because parameter of each direction is different.
>> Is it work if the dimension of rho is equal to the dimension of mesh
>> in following code?
>>
>> Code of C++:
>> model.add_fixed_size_data("rho", rho);
>> getfem::add_linear_term(model, mim, "rho*u.Test_u", TOP);
>>
>> Best Regard Tetsuo


--

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

---------