<div>Hello all,</div><div><br></div><div>In section 4.2 of the manual, it is stated that the only option for solving linear systems having the same preconditioner is to call KSPSolve() multiple times.  However, when dealing with such a system, one can take advantage of DGEMM when updating the vectors as a &#39;multivector&#39;.  Has anyone attempted to create a multivector object in PETSc?  If so, what is the interface and where can I find the code?  If not, can anyone speak to the difficulties?  </div>
<div><br></div><div>Thank you,</div><div>Henc</div><div><br></div>