<div dir="ltr">Hi,<div><br></div><div>I am trying to do something apparently really simple but with no success.</div><div><br></div><div>I need to perform a matrix-vector multiplication x = B f , where the length of x is bigger than the length of f (or viceversa). Thus, B cannot be created using <font face="monospace, monospace">DMCreateMatrix</font>.</div><div><br></div><div>Both x and f are obtained from different DMs, the smaller covering only a subdomain of the larger. The application is to apply a control f to a system, e.g. \dot{x} = A x + B f.</div><div><br></div><div>The problem is, when running on more than one core, the vector x is not organized as I would expect (everything works fine on a single core). </div><div><br></div><div>I attach a short example where B is intended to map f to the interior of x.</div><div><br></div><div><font face="monospace, monospace">mpirun -n 1 ./test -draw_pause -1</font>     works fine while</div><div><font face="monospace, monospace">mpirun -n 2 ./test -draw_pause -1</font>     shows the problem</div><div><br></div><div>I have not found any example with non square matrices in the src folder, any help is very welcome.</div><div><br></div><div>Thanks for your time,</div><div><br></div><div>Gianluca</div></div>