<html>
<head>
<meta content="text/html; charset=utf-8" http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<meta http-equiv="content-type" content="text/html; charset=utf-8">
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</font><br>
</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</font><br>
</blockquote>
<font face="Courier New, Courier, monospace"> </font> 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 <a
href="http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MATSOLVERMUMPS.html">documentation</a>
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</font><br>
</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</font><br>
</blockquote>
which does not work any better.<br>
Any idea how I should pass the icntl_14 information to mumps in this
case?<br>
<pre class="moz-signature" cols="72">--
Best,
Luc</pre>
<blockquote cite="mid:546229AF.1010808@columbi.edu" type="cite"> </blockquote>
<br>
</body>
</html>