[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Getfem-users] [Getfem-commits] What is the best way to add paramete
From: |
Tetsuo Koyama |
Subject: |
Re: [Getfem-users] [Getfem-commits] What is the best way to add parameter for each direction in term? |
Date: |
Tue, 12 Feb 2019 11:47:19 +0900 |
Dear Yves and Konstantinos
Than you for your comment. I'll try it and comment later.
(Sorry, I mistook the sending mail address.)
Best Regards,
Tetsuo
2019年2月12日(火) 2:56 Konstantinos Poulios <address@hidden>:
>
> 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
>>
>> ---------
>>
>>