<div dir="ltr"><div><div>Thank you, Jed. "-snes_compare_explicit" worked with a nested matrix. I'll try to use MatGetLocalSubMatrix in my programm.<br><br></div>Best,<br></div>nori<br></div><div class="gmail_extra">
<br><br><div class="gmail_quote">On Thu, Apr 17, 2014 at 4:56 PM, Jed Brown <span dir="ltr"><<a href="mailto:jed@jedbrown.org" target="_blank">jed@jedbrown.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div class="">Norihiro Watanabe <<a href="mailto:norihiro.w@gmail.com">norihiro.w@gmail.com</a>> writes:<br>
<br>
> Hi, I wonder whether "-snes_type test" option can work also with a nested<br>
> matrix. I got the following error in my program and am not sure if there is<br>
> a problem in my codes or just PETSc doesn't support it.<br>
<br>
</div>I believe -snes_compare_explicit should convert if necessary (it<br>
provides similar functionality to -snes_type test).<br>
<div class=""><br>
> Testing hand-coded Jacobian, if the ratio is<br>
> O(1.e-8), the hand-coded Jacobian is probably correct.<br>
> Run with -snes_test_display to show difference<br>
> of hand-coded and finite difference Jacobian.<br>
> [0]PETSC ERROR: --------------------- Error Message<br>
> ------------------------------------<br>
> [0]PETSC ERROR: No support for this operation for this object type!<br>
> [0]PETSC ERROR: Mat type nest!<br>
<br>
</div>We strongly recommend writing your code in a way that is agnostic to<br>
matrix format (i.e., assemble using MatGetLocalSubMatrix).  Then just<br>
use a monolithic format for comparisons like this and switch to MatNest<br>
when trying to get the best performance out of a fieldsplit<br>
preconditioner.<br>
</blockquote></div><br><br clear="all"><br>-- <br>Norihiro Watanabe
</div>