<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">Luc:</div><div class="gmail_quote">Run your code with option '-help |grep mumps', then you'll see what prefix should be used in your case with the mumps option </div><div class="gmail_quote"><span style="font-family:'Courier New',Courier,monospace">'-mat_mumps_icntl_14 30'.</span></div><div class="gmail_quote"><span style="font-family:'Courier New',Courier,monospace">You may try even larger </span><span style="font-family:'Courier New',Courier,monospace">icntl_14.</span></div><div class="gmail_quote"><span style="font-family:'Courier New',Courier,monospace"><br></span></div><div class="gmail_quote"><span style="font-family:'Courier New',Courier,monospace">Hong</span></div><div class="gmail_quote"><font face="Courier New, Courier, monospace"><br></font><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">
  

    
  
  <div bgcolor="#FFFFFF" text="#000000">
    Hi, I am using Petsc to solver a multiphysics problem and I have the
    following issue.<br>
    I partition my problem by declaring two fields:<br>
    <blockquote><font face="Courier New, Courier, monospace">-ksp_type
        gmres -pc_type fieldsplit -pc_fieldsplit_type schur
        -pc_fieldsplit_schur_factorization_type full
        -pc_fieldsplit_schur_precondition selfp -pc_fieldsplit_0_fields
        2,3 -pc_fieldsplit_1_fields 0,1<br>
      </font></blockquote>
    I want to solve the matrix representing field 0 with mumps so I pass
    the following following arguments:<br>
    <blockquote><font face="Courier New, Courier, monospace">-fieldsplit_0_ksp_type
        preonly -fieldsplit_0_pc_type lu -fieldsplit_0_pc_factor_mat_solver_package
        mumps<br>
      </font></blockquote>
    When I do this I get an error from mumps: INFO(1)=-9, INFO(2)=12532.
    This means that mumps main internal real workarray is too small and
    12532 are missing. To try to mitigate this I need to want to set
    mumps ICNTL(14)=30 (by default it is 20).<br>
    Reading Petsc documentation I find that I have to pass the following
    argument to my program:<br>
    <blockquote><font face="Courier New, Courier, monospace">-mat_mumps_icntl_14
        30<br>
      </font></blockquote>
    which does work fine when I work without fieldsplit but not when I
    use field split.<br>
    I also tried:<br>
    <blockquote><font face="Courier New, Courier, monospace">-fieldsplit_0_mat_mumps_icntl_14
        30<br>
      </font></blockquote>
    which does not work any better.<br>
    Any idea how I should pass the icntl_14 information to mumps in this
    case?<span class=""><font color="#888888"><br>
    <pre cols="72">-- 
Best,
Luc</pre>
  </font></span></div>

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