<div dir="ltr"><div>subramanya:<br></div><div class="gmail_extra"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">





<div>
<div>Hi,<br>
I figured out that block size was missing. However I am not quite sure what the block size means.
<br></div></div></blockquote><div> </div><div>number of fields.</div><div>See</div><div> ./ex19 -snes_monitor_short -ksp_monitor_short -pc_type fieldsplit -pc_fieldsplit_block_size 2 -pc_fieldsplit_0_fields 0,1 -pc_fieldsplit_1_fields 0,1 -pc_fieldsplit_type multiplicative -snes_view   -fieldsplit_pc_type lu -da_refine 1 -ksp_type fgmres -h |grep fieldsplit_block_size</div>
<div>  -pc_fieldsplit_block_size <-1>: Blocksize that defines number of fields (PCFieldSplitSetBlockSize)</div><div><br></div><div>Hong</div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div><div>
<br>
--- Original Message ---<br>
<br>
From: "Hong Zhang" <<a href="mailto:hzhang@mcs.anl.gov" target="_blank">hzhang@mcs.anl.gov</a>><br>
Sent: December 28, 2013 12:06 AM<br>
To: "subramanya sadasiva" <<a href="mailto:potaman@outlook.com" target="_blank">potaman@outlook.com</a>><br>
Cc: <a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>, <a href="mailto:libmesh-users@lists.sourceforge.net" target="_blank">libmesh-users@lists.sourceforge.net</a><br>
Subject: Re: [petsc-users] fieldsplit error<br>
<br>
</div><div><div class="h5">
<div>
<div dir="ltr">Add
<div>-pc_fieldsplit_block_size 2</div>
<div><br>
</div>
<div>See example </div>
<div>petsc/src/snes/examples/tutorials/ex19.c:</div>
<div>./ex19 -snes_monitor_short -ksp_monitor_short -pc_type fieldsplit -pc_fieldsplit_block_size 2 -pc_fieldsplit_0_fields 0,1 -pc_fieldsplit_1_fields 0,1 -pc_fieldsplit_type multiplicative -snes_view   -fieldsplit_pc_type lu -da_refine 1 -ksp_type fgmres<br>

</div>
<div><br>
</div>
<div>(see runex19_8 in petsc/src/snes/examples/tutorials/makefile)</div>
<div><br>
</div>
<div>Hong</div>
<div><br>
</div>
<div>
<div><br>
<div>On Fri, Dec 27, 2013 at 9:25 PM, subramanya sadasiva <span dir="ltr">
<<a href="mailto:potaman@outlook.com" target="_blank">potaman@outlook.com</a>></span> wrote:<br>
<blockquote style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>
<div dir="ltr">
<div>Hi, </div>
<div>I have a nonlinear system with 6 variables grouped in two sets of 3. I am trying to get a fieldsplit preconditioner working . I have the following options set up </div>
<div>
<div><br>
</div>
<div>-pc_type fieldsplit</div>
<div>-pc_fieldsplit_type multiplicative</div>
<div>-fieldsplit_0_fields 0,2,4</div>
<div>-fieldsplit_1_fields 1,3,5</div>
<div>-fieldsplit_0_ksp cg</div>
<div>-fieldsplit_1_ksp cg</div>
<div>-fieldsplit_0_pc hypre</div>
<div>-fieldsplit_1_pc hypre</div>
</div>
<div><br>
</div>
<div>However, I get the following error. My installation of libmesh has block matrices disabled, which I've seen to cause problems with fieldsplit. Any ideas what might cause this error?</div>
<div>Thanks, </div>
<div>Subramanya </div>
<div><br>
</div>
<div><br>
</div>
<div><br>
</div>
<div>[0]PETSC ERROR: --------------------- Error Message ------------------------------------</div>
<div>[0]PETSC ERROR: Petsc has generated inconsistent data!</div>
<div>[0]PETSC ERROR: Unhandled case, must have at least two fields, not 1!</div>
<div>[0]PETSC ERROR: ------------------------------------------------------------------------</div>
<div>[0]PETSC ERROR: Petsc Release Version 3.4.3, unknown </div>
<div>[0]PETSC ERROR: See docs/changes/index.html for recent updates.</div>
<div>[0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.</div>
<div>[0]PETSC ERROR: See docs/index.html for manual pages.</div>
<div><br>
</div>
</div>
</div>
</blockquote>
</div>
<br>
</div>
</div>
</div>
</div>
</div></div></div>

</blockquote></div><br></div></div>