On Thu, Feb 23, 2012 at 10:46 AM, Bojan Niceno <span dir="ltr"><<a href="mailto:bojan.niceno@psi.ch">bojan.niceno@psi.ch</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi all,<br>
<br>
I've never used a mailing list before, so I hope this message will reach PETSc users and experts and someone might be willing to help me. I am also novice in PETSc.<br>
<br>
I have developed an unstructured finite volume solver on top of PETSc libraries. In sequential, it works like a charm. For the parallel version, I do domain decomposition externally with Metis, and work out local and global numberings, as well as communication patterns between processor. (The latter don't seem to be needed for PETSc, though.) When I run my program in parallel, it also works, but I miss values in vectors' ghost points.<br>
<br>
I create vectors with command: VecCreate(PETSC_COMM_WORLD, &x);<br>
<br>
Is it possible to get the ghost values if a vector is created like this?<br></blockquote><div><br></div><div>I do not understand this question. By definition, "ghost values" are those not stored in the global vector.</div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
I have tried to use VecCreateGhost, but for some reason which is beyond my comprehension, PETSc goes berserk when it reaches the command: VecCreateGhost(PETSC_COMM_<u></u>WORLD, n, PETSC_DECIDE, nghost, ifrom, &x)<br>
</blockquote><div><br></div><div>I think you can understand that "berserk" tells me absolutely nothing. Error message? Stack trace? Did you try to run an</div><div>example which uses VecGhost?</div><div><br></div>
<div> Thanks,</div><div><br></div><div> Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Can anyone help me? Either how to reach ghost values for vector created by VecCreate, or how to use VecCreateGhost properly?<br>
<br>
<br>
Kind regards,<br>
<br>
Bojan<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener<br>