<div dir="ltr"><div>Hello, <br></div><div><br></div><div>I am looking for a way to compute two matrix vector products simultaneously, where the matrix is the same for both products (i.e. Ax and Ay). I know some numerical libraries have the concept of "block 
vectors" or "basis vectors" and was wondering if there is anything similar in PETSc. Right now I am just using two MatMults sequentially, but I'm hoping there are probably better ways to do this. My matrices are MPIAIJ if that makes any difference.<br></div><div><br></div><div>Best,<br></div><div>Tyler<br></div></div>