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 <[email protected]>: > > 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 <[email protected]> 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 <[email protected]>: >> >> 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 ([email protected]) tel : (33) 04.72.43.87.08 >> INSA-Lyon >> 20, rue Albert Einstein >> 69621 Villeurbanne Cedex, FRANCE >> http://math.univ-lyon1.fr/~renard >> >> --------- >> >>
