Creating a new vector based on a already distributed set of Vecs
Matthew Knepley
knepley at gmail.com
Sun Mar 23 15:12:33 CDT 2008
2008/3/23 Vijay S. Mahadevan <vijay.m at gmail.com>:
> Hi all,
>
> I am sure the subject line wasn't too explanatory and so here goes.
>
> I have a system which already has a Vec object created and distributed on
> many processors based on domain decomposition. I could have several such
> systems, each with a solution vector and a residual vector pertaining to its
> physics.
>
> Now, can I create a global vector which would just be pointers to the
> already existing vectors but from PETSc's point of view, appears to be a
> globally valid, (and possibly weirdly) distributed vector. Such an option
> will save the memory needed for the global vector and eliminates errors as
> to synchronization of the solution, residuals for the systems.
I do not understand what you want. If you want a vector with the same layout,
use VecDuplicate(). If you want the same data, just use the Vec itself.
Thanks,
Matt
> I was reading the documentation and found the PetscMap data structure. But I
> am not entirely sure if this is what I am looking for.
>
> I hope that makes sense. I would appreciate any help you can provide to
> point me in the right direction.
>
> Cheers,
> Vijay
>
> No virus found in this outgoing message.
> Checked by AVG.
> Version: 7.5.519 / Virus Database: 269.21.8/1339 - Release Date: 3/22/2008
> 4:43 PM
>
>
>
--
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which
their experiments lead.
-- Norbert Wiener
More information about the petsc-users
mailing list