<p dir="ltr">MATNEST is merely an optimization, and not usually the most important one. Get everything working for PC fieldsplit first and then we can see if NEST makes sense.</p>
<p dir="ltr">  Thanks</p>
<p dir="ltr">      Matt</p>
<div class="gmail_quote">On Dec 6, 2014 7:56 AM, "linjing bo" <<a href="mailto:francium87@hotmail.com">francium87@hotmail.com</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">


<div><div dir="ltr">The Problem in short : How to manually set IS for MATNEST instead of calling MatNestGetISs in fortran?<div><br></div><div>Problems in detail: I'm trying to solve a 3D vector field in Fortran. Base on manual page, a efficient way is to use MATNEST with 3x3 block configuration. So I follow the ex70.c and use MatNestGetISs to set the Vector and other things. But no symbol is found in library. After searching the mailing list I notice the Fortran support is fixed recently in 3.5. But for capability issue ( most of the Supercomputer I can access only support optimized version 3.4.4 and older ), using the latest 3.5 is not a option. So I'm wondering is there a way to manually set IS to overcome this issue, or good suggestion other than MATNEST of solving a multiphysics problem?</div><div><br></div><div>Some Trial: I tried to use ISCreateBlock to create a block IS with 3 blocks on vector and use it instead of MatNestGetISs, But PETSc report that the local Matrix and Vector size are not match:</div><div><br></div><div><div>[6]PETSC ERROR: MatMult() line 2172 in /tmp/petsc-3.4.4/src/mat/interface/matrix.c Mat mat,Vec y: local dim 7320 7319</div></div><div><br></div>                                         </div></div>
</blockquote></div>