Did you try running with -on_error_attach_debugger and and using your debugger to try to figure out where your code is segfaulting?<br><br>Max<div><br><div class="gmail_quote">On Thu, Feb 23, 2012 at 10:07 AM, Bojan Niceno <span dir="ltr"><<a href="mailto:bojan.niceno@psi.ch">bojan.niceno@psi.ch</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div bgcolor="#FFFFFF" text="#000000">
Thanks Ju. I studied this case carefully, and it seems clear to
me. When I apply the same techniques in my code, I get error
messages I sent to my reply to Matthew.<br>
<br>
<br>
Cheers,<br>
<br>
<br>
Bojan<div><div class="h5"><br>
<br>
<br>
On 2/23/2012 6:46 PM, Ju LIU wrote:
<blockquote type="cite"><br>
<br>
<div class="gmail_quote">2012/2/23 Bojan Niceno <span dir="ltr"><<a href="mailto:bojan.niceno@psi.ch" target="_blank">bojan.niceno@psi.ch</a>></span><br>
<blockquote class="gmail_quote" style="margin:0pt 0pt 0pt 0.8ex;border-left:1px solid rgb(204,204,204);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>
<br>
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_WORLD, n, PETSC_DECIDE,
nghost, ifrom, &x)<br>
<br>
Can anyone help me? Either how to reach ghost values for
vector created by VecCreate, or how to use VecCreateGhost
properly?<br>
<br>
</blockquote>
<div><a href="http://www.mcs.anl.gov/petsc/petsc-current/src/vec/vec/examples/tutorials/ex9.c.html" target="_blank">http://www.mcs.anl.gov/petsc/petsc-current/src/vec/vec/examples/tutorials/ex9.c.html</a>
could be helpful.<br>
<br>
</div>
<blockquote class="gmail_quote" style="margin:0pt 0pt 0pt 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<br>
Bojan<br>
</blockquote>
</div>
<br>
</blockquote>
<br>
<br>
</div></div><span class="HOEnZb"><font color="#888888"><div>-- <br>
<img src="cid:part1.08030804.07060605@psi.ch" border="0"></div>
</font></span></div>
</blockquote></div><br></div>