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 > > --------- > > >
