<html>
<head>
<meta content="text/html; charset=windows-1252"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<p>Further to this, I debugged into PetscObjectTypeCompare().</p>
<p>It's returning PETSC_FALSE because the type name of the variable
mat is "seqbaij" (and I guess would normally be something like
"mpibaij" if running in parallel outside gdb), which doesn't match
with MATMFFD (= "mffd").</p>
<p>Does that mean I should be setting the matrix type to MATMFFD to
get this to work?<br>
</p>
<p>- Adrian<br>
</p>
<br>
<div class="moz-cite-prefix">On 24/06/16 16:35, Adrian Croucher
wrote:<br>
</div>
<blockquote
cite="mid:ad72d07a-8660-8b63-e6f0-d1cc4c1bc9f9@auckland.ac.nz"
type="cite">
<meta content="text/html; charset=windows-1252"
http-equiv="Content-Type">
<div class="moz-cite-prefix">On 24/06/16 15:58, Adrian Croucher
wrote:<br>
</div>
<blockquote
cite="mid:8deba8e3-d969-fe53-038b-b7447168a05d@auckland.ac.nz"
type="cite">
<meta content="text/html; charset=windows-1252"
http-equiv="Content-Type">
<br>
I will try debugging into MatMFFDSetType() and seeing what it's
doing when I pass in "ds".<br>
<br>
</blockquote>
<br>
I did this. In MatMFFDSetType() (mffd.c: 116), it executes <br>
<br>
ierr =
PetscObjectTypeCompare((PetscObject)mat,MATMFFD,&match);CHKERRQ(ierr);<br>
<br>
which returns match = PETSC_FALSE. The next line<br>
<br>
if (!match) PetscFunctionReturn(0);<br>
<br>
means it just returns out of there without doing anything.<br>
<br>
- Adrian<br>
</blockquote>
<br>
<pre class="moz-signature" cols="72">--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: <a class="moz-txt-link-abbreviated" href="mailto:a.croucher@auckland.ac.nz">a.croucher@auckland.ac.nz</a>
tel: +64 (0)9 923 84611
</pre>
</body>
</html>