<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
Thanks for the reply Matthew<br>
Your suggestion to use MatGetValues helped me to figure out my
problem.<br>
<br>
The legacy part of my software decided I didn't have enough memory
to solve the problem in one shot, so it split it into two pieces,
unbeknownst to me. So my A matrix was being partially loaded (hence
the unfilled diagonals)<br>
<br>
<br>
<div class="moz-cite-prefix">On 08/16/2013 11:34 AM, Matthew Knepley
wrote:<br>
</div>
<blockquote
cite="mid:CAMYG4GmOc7yE7bRNG_ELfPjfj6PR0fXWeMTzuawUuk6Qf2=m6w@mail.gmail.com"
type="cite">
<div dir="ltr">On Fri, Aug 16, 2013 at 10:09 AM, Bill McGrory <span
dir="ltr"><<a moz-do-not-send="true"
href="mailto:mcgrory@aerosoftinc.com" target="_blank">mcgrory@aerosoftinc.com</a>></span>
wrote:<br>
<div class="gmail_extra">
<div class="gmail_quote">
<blockquote class="gmail_quote" style="margin:0px 0px 0px
0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">More
information pertaining to the referenced post.<br>
<br>
I am using default parameters for my KSP solver, so my
routines are pretty simple I don't pass any options in
through the command line, so this is all default stuff
here.<br>
<br>
I create, and fill A and b, and then call<br>
MatAssemblyBegin(A,MAT_FINAL_ASSEMBLY);<br>
<br>
MatAssemblyEnd(A,MAT_FINAL_ASSEMBLY); <br>
<br>
VecAssemblyBegin(b);<br>
VecAssemblyEnd(b);<br>
<br>
VecSet(x,0.);<br>
<br>
// Solve A x = b<br>
<br>
KSPCreate(PETSC_COMM_WORLD,&ksp);<br>
<br>
KSPSetOperators(ksp,A,A,DIFFERENT_NONZERO_PATTERN);<br>
<br>
KSPSetFromOptions(ksp);<br>
<br>
KSPSolve(ksp,b,x);<br>
<br>
The missing diagonals, I see are in the preconditioner,
not the original Matrix, so I thought I would check my
original.<br>
</blockquote>
<div><br>
</div>
<div>From the code above, this makes no sense. You do not
have a separate preconditioner matrix.</div>
<div><br>
</div>
<blockquote class="gmail_quote" style="margin:0px 0px 0px
0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">When
I make a call to MatMissingDiagonal, after assembling my
matrix, I get the error telling me that MatMissingDiagonal
is not supported for a mpibaij matrix.<br>
<br>
Do I have any alternative for querying my assembled
matrix?<br>
</blockquote>
<div><br>
</div>
<div>You can always use MatGetValues().</div>
<div><br class="">
Matt</div>
<div> </div>
<blockquote class="gmail_quote" style="margin:0px 0px 0px
0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">Thanks
again<span class=""><font color="#888888"><br>
Bill<br>
<br>
<br>
</font></span></blockquote>
</div>
<br>
<br clear="all">
<div><br>
</div>
</div>
</div>
</blockquote>
Very true statement in my case below. Ah, but to know what you are
taking for granted.<br>
<blockquote
cite="mid:CAMYG4GmOc7yE7bRNG_ELfPjfj6PR0fXWeMTzuawUuk6Qf2=m6w@mail.gmail.com"
type="cite">
<div dir="ltr">
<div class="gmail_extra">-- <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</div>
</div>
</blockquote>
<br>
</body>
</html>