<div dir="ltr">Hi guys,<div><br></div><div>I've added a small function to only solve the linear system in parallel for problems that are not big enough to justify parallelizing the whole thing but would still take some time to solve in serial -- so far It's been handy. To prevent extra copy from my own format to petsc's, I'm using VecCreateMPIWithArray and MatCreateMPIWithSplitArray functions and they work with MPI on very few processes (like 4 or so). I have a couple of questions:</div>

<div><br></div><div>1) Is it worth considering using pThread and/or CUSP versions of the matrices for such problems?</div><div>2) If so, what would be a better approach. Using hardwired vector types in the code or using generic VecCreate functions to be able to change the type at run-time?</div>

<div>3) Are there equivalent versions of XXXWithArray functions for CUSP and pThread types? I don't seem to find them in manual page.</div><div><br></div><div>Thanks</div></div>