[petsc-users] Regarding changing vector/matrix type at runtime

[Mohammad Mirzadeh]

Hi guys,

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:

1) Is it worth considering using pThread and/or CUSP versions of the
matrices for such problems?
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?
3) Are there equivalent versions of XXXWithArray functions for CUSP and
pThread types? I don't seem to find them in manual page.

Thanks
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120823/87b96b51/attachment.html>