mld2p4-2:

Docs updates for DIAG/Gauss-Seidel

docs/html/node12.html

@@ -726,7 +726,7 @@ $w = y_1$;
 \begin{tabbing}
 \quad \=\quad...
 ...= y_l+r_l$\\
-\textbf{endfor} \ [1mm]
+\textbf{endfor}  [1mm]
 $w = y_1$;
 \end{tabbing}}
 \end{minipage}}">

docs/html/node18.html

@@ -191,7 +191,7 @@ according to their needs.
 
 <P>
 <BR><P></P>
-<DIV ALIGN="CENTER"><A NAME="1267"></A>
+<DIV ALIGN="CENTER"><A NAME="1271"></A>
 <TABLE>
 <CAPTION><STRONG>Table 2:</STRONG>
 Parameters defining the type of multi-level preconditioner.
@@ -214,11 +214,11 @@ Parameters defining the type of multi-level preconditioner.
 </TR>
 <TR><TD ALIGN="LEFT" VALIGN="TOP" WIDTH=142><code>mld_smoother_type_</code> <BR><code>SMOOTHER_TYPE</code></TD>
 <TD ALIGN="LEFT"><code>character(len=*)</code></TD>
-<TD ALIGN="LEFT" VALIGN="TOP" WIDTH=57><TT>'DIAG'</TT>    <TT>'BJAC'</TT>    <TT>'AS'</TT></TD>
+<TD ALIGN="LEFT" VALIGN="TOP" WIDTH=57><TT>'JACOBI'</TT>    <TT>'BJAC'</TT>    <TT>'AS'</TT></TD>
 <TD ALIGN="LEFT"><TT>'AS'</TT></TD>
 <TD ALIGN="LEFT" VALIGN="TOP" WIDTH=198>Basic predefined one-level preconditioner
-                         (i.e. smoother): diagonal, 
-                           block Jacobi, AS.</TD>
+                         (i.e. smoother): Jacobi (with diagonal local
+                           solver), block Jacobi, AS.</TD>
 </TR>
 <TR><TD ALIGN="LEFT" VALIGN="TOP" WIDTH=142><code>mld_smoother_pos_</code> <BR><code>SMOOTHER_POS</code></TD>
 <TD ALIGN="LEFT"><code>character(len=*)</code></TD>
@@ -236,7 +236,7 @@ Parameters defining the type of multi-level preconditioner.
 
 <P>
 <BR><P></P>
-<DIV ALIGN="CENTER"><A NAME="1269"></A>
+<DIV ALIGN="CENTER"><A NAME="1273"></A>
 <TABLE>
 <CAPTION><STRONG>Table 3:</STRONG>
 Parameters defining the one-level preconditioner used as smoother.
@@ -277,21 +277,25 @@ Parameters defining the one-level preconditioner used as smoother.
 </TR>
 <TR><TD ALIGN="LEFT" VALIGN="TOP" WIDTH=99><code>mld_sub_solve_</code> <BR><code>SUB_SOLVE</code></TD>
 <TD ALIGN="LEFT"><code>character(len=*)</code></TD>
-<TD ALIGN="LEFT" VALIGN="TOP" WIDTH=91><TT>'ILU'</TT>   <TT>'MILU'</TT>   <TT>'ILUT'</TT> 
+<TD ALIGN="LEFT" VALIGN="TOP" WIDTH=91><TT>'DIAG'</TT>  
+                           <TT>'GS'</TT>  
+                           <TT>'ILU'</TT>   
+                           <TT>'MILU'</TT>   <TT>'ILUT'</TT> 
                              <TT>'UMF'</TT>   <TT>'SLU'</TT></TD>
 <TD ALIGN="LEFT"><TT>'ILU'</TT></TD>
-<TD ALIGN="LEFT" VALIGN="TOP" WIDTH=142>Predefined local solver: ILU(<IMG
+<TD ALIGN="LEFT" VALIGN="TOP" WIDTH=142>Predefined local solver: diagonal scaling,
+                           Gauss-Seidel, ILU(<IMG
  WIDTH="13" HEIGHT="31" ALIGN="MIDDLE" BORDER="0"
  SRC="img35.png"
- ALT="$p$">),
-                         MILU(<IMG
+ ALT="$p$">),  MILU(<IMG
  WIDTH="13" HEIGHT="31" ALIGN="MIDDLE" BORDER="0"
  SRC="img35.png"
- ALT="$p$">), ILU(<IMG
+ ALT="$p$">),
+                           ILU(<IMG
  WIDTH="27" HEIGHT="31" ALIGN="MIDDLE" BORDER="0"
  SRC="img36.png"
- ALT="$p,t$">), LU from UMFPACK, LU  from SuperLU 
-                           (plus triangular solve).</TD>
+ ALT="$p,t$">), LU from UMFPACK, LU  from
+                           SuperLU (plus triangular solve).</TD>
 </TR>
 <TR><TD ALIGN="LEFT" VALIGN="TOP" WIDTH=99><code>mld_sub_fillin_</code> <BR><code>SUB_FILLIN</code></TD>
 <TD ALIGN="LEFT"><code>integer</code></TD>
@@ -328,6 +332,27 @@ Parameters defining the one-level preconditioner used as smoother.
                            reordering according to the global numbering of the rows and columns of
                            the whole matrix.</TD>
 </TR>
+<TR><TD ALIGN="LEFT" VALIGN="TOP" WIDTH=99><code>mld_solver_sweeps_</code> <BR><code>SOLVER_SWEEPS</code></TD>
+<TD ALIGN="LEFT"><code>integer</code></TD>
+<TD ALIGN="LEFT" VALIGN="TOP" WIDTH=91>Any&nbsp;int.&nbsp;num.&nbsp;<IMG
+ WIDTH="31" HEIGHT="31" ALIGN="MIDDLE" BORDER="0"
+ SRC="img92.png"
+ ALT="$\ge 1$"></TD>
+<TD ALIGN="LEFT">1</TD>
+<TD ALIGN="LEFT" VALIGN="TOP" WIDTH=142>Number of sweeps for iterative local solver
+                           (currently only Gauss-Seidel).</TD>
+</TR>
+<TR><TD ALIGN="LEFT" VALIGN="TOP" WIDTH=99><code>mld_solver_eps_</code> <BR><code>SOLVER_EPS</code></TD>
+<TD ALIGN="LEFT"><code>real</code></TD>
+<TD ALIGN="LEFT" VALIGN="TOP" WIDTH=91>Any&nbsp;real&nbsp;num.</TD>
+<TD ALIGN="LEFT">0</TD>
+<TD ALIGN="LEFT" VALIGN="TOP" WIDTH=142>Stopping tolerance for iterative local solver
+                           (currently only Gauss-Seidel); if <IMG
+ WIDTH="31" HEIGHT="31" ALIGN="MIDDLE" BORDER="0"
+ SRC="img93.png"
+ ALT="$\le0$">, then
+                           perform prespecified number of iterations.</TD>
+</TR>
 </TABLE>
 </DIV>
 </TD></TR>
@@ -337,7 +362,7 @@ Parameters defining the one-level preconditioner used as smoother.
 
 <P>
 <BR><P></P>
-<DIV ALIGN="CENTER"><A NAME="1271"></A>
+<DIV ALIGN="CENTER"><A NAME="1275"></A>
 <TABLE>
 <CAPTION><STRONG>Table 4:</STRONG>
 Parameters defining the aggregation algorithm.
@@ -381,28 +406,28 @@ Parameters defining the aggregation algorithm.
 <TD ALIGN="LEFT"><code>real(</code><I>kind_parameter</I><code>)</code></TD>
 <TD ALIGN="LEFT" VALIGN="TOP" WIDTH=68>Any&nbsp;real&nbsp;num. <IMG
  WIDTH="56" HEIGHT="36" ALIGN="MIDDLE" BORDER="0"
- SRC="img92.png"
+ SRC="img94.png"
  ALT="$\in [0, 1]$"></TD>
 <TD ALIGN="LEFT" VALIGN="TOP" WIDTH=68>0.05</TD>
 <TD ALIGN="LEFT" VALIGN="TOP" WIDTH=170>Threshold <IMG
  WIDTH="13" HEIGHT="16" ALIGN="BOTTOM" BORDER="0"
- SRC="img93.png"
+ SRC="img95.png"
  ALT="$\theta$"> in the aggregation algorithm.</TD>
 </TR>
 <TR><TD ALIGN="LEFT" VALIGN="TOP" WIDTH=142><code>mld_aggr_scale_</code> <BR><code>AGGR_SCALE</code></TD>
 <TD ALIGN="LEFT"><code>real(</code><I>kind_parameter</I><code>)</code></TD>
 <TD ALIGN="LEFT" VALIGN="TOP" WIDTH=68>Any&nbsp;real&nbsp;num. <IMG
  WIDTH="56" HEIGHT="36" ALIGN="MIDDLE" BORDER="0"
- SRC="img92.png"
+ SRC="img94.png"
  ALT="$\in [0, 1]$"></TD>
 <TD ALIGN="LEFT" VALIGN="TOP" WIDTH=68>1.0</TD>
 <TD ALIGN="LEFT" VALIGN="TOP" WIDTH=170>Scale factor applied to the threshold going
  from level <IMG
  WIDTH="33" HEIGHT="16" ALIGN="BOTTOM" BORDER="0"
- SRC="img94.png"
+ SRC="img96.png"
  ALT="$ilev$"> to level <IMG
  WIDTH="63" HEIGHT="34" ALIGN="MIDDLE" BORDER="0"
- SRC="img95.png"
+ SRC="img97.png"
  ALT="$ilev+1$">.</TD>
 </TR>
 <TR><TD ALIGN="LEFT" VALIGN="TOP" WIDTH=142><code>mld_aggr_omega_alg_</code>  <BR><code>AGGR_OMEGA_ALG</code></TD>
@@ -440,7 +465,7 @@ Parameters defining the aggregation algorithm.
  -->
 <IMG
  WIDTH="113" HEIGHT="39" ALIGN="MIDDLE" BORDER="0"
- SRC="img96.png"
+ SRC="img98.png"
  ALT="$4/(3\rho(D^{-1}A))$"></TD>
 <TD ALIGN="LEFT" VALIGN="TOP" WIDTH=170>Damping parameter <IMG
  WIDTH="16" HEIGHT="18" ALIGN="BOTTOM" BORDER="0"
@@ -452,7 +477,7 @@ Parameters defining the aggregation algorithm.
                            otherwise it is computed by the library, using the
                            selected estimate of the spectral radius <IMG
  WIDTH="73" HEIGHT="39" ALIGN="MIDDLE" BORDER="0"
- SRC="img97.png"
+ SRC="img99.png"
  ALT="$\rho(D^{-1}A)$"> of
                            <IMG
  WIDTH="50" HEIGHT="21" ALIGN="BOTTOM" BORDER="0"
@@ -468,7 +493,7 @@ Parameters defining the aggregation algorithm.
 
 <P>
 <BR><P></P>
-<DIV ALIGN="CENTER"><A NAME="1274"></A>
+<DIV ALIGN="CENTER"><A NAME="1278"></A>
 <TABLE>
 <CAPTION><STRONG>Table 5:</STRONG>
 Parameters defining the coarse-space correction at the coarsest
@@ -526,7 +551,7 @@ level.</CAPTION>
 <TD ALIGN="LEFT"><code>integer</code></TD>
 <TD ALIGN="LEFT" VALIGN="TOP" WIDTH=91>Any&nbsp;int.&nbsp;num.&nbsp;<IMG
  WIDTH="32" HEIGHT="31" ALIGN="MIDDLE" BORDER="0"
- SRC="img98.png"
+ SRC="img100.png"
  ALT="$&gt; 0$"></TD>
 <TD ALIGN="LEFT">4</TD>
 <TD ALIGN="LEFT" VALIGN="TOP" WIDTH=142>Number of Block-Jacobi sweeps when 'BJAC' is used as

docs/html/node20.html

@@ -69,14 +69,14 @@ This routine computes <!-- MATH
  -->
 <IMG
  WIDTH="117" HEIGHT="39" ALIGN="MIDDLE" BORDER="0"
- SRC="img99.png"
+ SRC="img101.png"
  ALT="$y = op(M^{-1})  x$">, where <IMG
  WIDTH="23" HEIGHT="15" ALIGN="BOTTOM" BORDER="0"
  SRC="img60.png"
  ALT="$M$"> is a previously built
 preconditioner, stored into <code>p</code>, and <IMG
  WIDTH="22" HEIGHT="31" ALIGN="MIDDLE" BORDER="0"
- SRC="img100.png"
+ SRC="img102.png"
  ALT="$op$">
 denotes the preconditioner itself or its transpose, according to
 the value of <code>trans</code>.
@@ -106,7 +106,7 @@ and hence it is completely transparent to the user.
 <TR><TD ALIGN="LEFT" VALIGN="TOP" WIDTH=34>&nbsp;</TD>
 <TD ALIGN="LEFT" VALIGN="TOP" WIDTH=340>The local part of the vector <IMG
  WIDTH="14" HEIGHT="18" ALIGN="BOTTOM" BORDER="0"
- SRC="img101.png"
+ SRC="img103.png"
  ALT="$x$">. Note that <I>type</I> and   
                 <I>kind_parameter</I> must be chosen according
                 to the real/complex, single/double precision version of MLD2P4 under use.</TD>
@@ -117,7 +117,7 @@ and hence it is completely transparent to the user.
 <TR><TD ALIGN="LEFT" VALIGN="TOP" WIDTH=34>&nbsp;</TD>
 <TD ALIGN="LEFT" VALIGN="TOP" WIDTH=340>The local part of the vector <IMG
  WIDTH="14" HEIGHT="31" ALIGN="MIDDLE" BORDER="0"
- SRC="img102.png"
+ SRC="img104.png"
  ALT="$y$">. Note that <I>type</I> and
                 <I>kind_parameter</I> must be chosen according
                 to the real/complex, single/double precision version of MLD2P4 under use.</TD>
@@ -144,28 +144,28 @@ and hence it is completely transparent to the user.
  -->
 <IMG
  WIDTH="132" HEIGHT="39" ALIGN="MIDDLE" BORDER="0"
- SRC="img103.png"
+ SRC="img105.png"
  ALT="$op(M^{-1}) = M^{-1}$">;
                 if <code>trans</code> = <code>'T','t'</code> then <!-- MATH
  $op(M^{-1}) = M^{-T}$
  -->
 <IMG
  WIDTH="135" HEIGHT="40" ALIGN="MIDDLE" BORDER="0"
- SRC="img104.png"
+ SRC="img106.png"
  ALT="$op(M^{-1}) = M^{-T}$">
                 (transpose of <IMG
  WIDTH="48" HEIGHT="39" ALIGN="MIDDLE" BORDER="0"
- SRC="img105.png"
+ SRC="img107.png"
  ALT="$M^{-1})$">;  if <code>trans</code> = <code>'C','c'</code> then <!-- MATH
  $op(M^{-1}) = M^{-C}$
  -->
 <IMG
  WIDTH="136" HEIGHT="40" ALIGN="MIDDLE" BORDER="0"
- SRC="img106.png"
+ SRC="img108.png"
  ALT="$op(M^{-1}) = M^{-C}$">
                 (conjugate transpose of <IMG
  WIDTH="48" HEIGHT="39" ALIGN="MIDDLE" BORDER="0"
- SRC="img105.png"
+ SRC="img107.png"
  ALT="$M^{-1})$">.</TD>
 </TR>
 <TR><TD ALIGN="LEFT" VALIGN="TOP" WIDTH=34><code>work</code></TD>

docs/mld2p4-2.1-guide.pdf

@@ -3361,7 +3361,7 @@ endstream
 endobj
 435 0 obj
 <<
-/Length 4494      
+/Length 4559      
 >>
 stream
 0 g 0 G
@@ -3489,7 +3489,7 @@ q
 []0 d 0 J 0.398 w 0 0 m 0 40.648 l S
 Q
 BT
-/F44 10.9091 Tf 563.081 143.921 Td [('DIAG')]TJ 0 -13.549 Td [('BJAC')]TJ 0 -13.549 Td [('AS')]TJ
+/F44 10.9091 Tf 563.081 143.921 Td [('JACOBI')]TJ 0 -13.549 Td [('BJAC')]TJ 0 -13.549 Td [('AS')]TJ
 ET
 q
 1 0 0 1 625.752 112.758 cm
@@ -3503,7 +3503,7 @@ q
 []0 d 0 J 0.398 w 0 0 m 0 40.648 l S
 Q
 BT
-/F15 10.9091 Tf 695.229 143.921 Td [(Basic)-434(p)1(rede\014ned)-434(one-lev)28(el)-433(preconditioner)]TJ 0 -13.549 Td [(\050i.e.)-534(sm)-1(o)-27(other\051:)-847(diagonal,)-584(blo)-28(c)28(k)-535(Jacobi,)]TJ 0 -13.549 Td [(AS.)]TJ
+/F15 10.9091 Tf 695.229 143.921 Td [(Basic)-434(p)1(rede\014ned)-434(one-lev)28(el)-433(preconditioner)]TJ 0 -13.549 Td [(\050i.e.)-357(smo)-27(other\051:)-491(Jacobi)-357(\050with)-356(diagonal)-356(lo-)]TJ 0 -13.549 Td [(cal)-333(solv)27(er\051,)-333(blo)-28(c)28(k)-333(Jacobi,)-334(AS)1(.)]TJ
 ET
 q
 1 0 0 1 899.632 112.758 cm
@@ -3572,7 +3572,7 @@ endstream
 endobj
 440 0 obj
 <<
-/Length 10037     
+/Length 12420     
 >>
 stream
 0 g 0 G
@@ -3588,385 +3588,471 @@ q
 0 g 0 G
 0 g 0 G
 q
-1 0 0 1 14.002 196.767 cm
+1 0 0 1 14.002 265.708 cm
 []0 d 0 J 0.398 w 0 0 m 565.57 0 l S
 Q
 q
-1 0 0 1 14.002 183.018 cm
+1 0 0 1 14.002 251.96 cm
 []0 d 0 J 0.398 w 0 0 m 0 13.549 l S
 Q
 1 0 0 1 -292.625 -715.095 cm
 BT
-/F44 10.9091 Tf 312.605 902.178 Td [(what)]TJ
+/F44 10.9091 Tf 312.605 971.119 Td [(what)]TJ
 ET
 q
-1 0 0 1 417.795 898.113 cm
+1 0 0 1 417.795 967.055 cm
 []0 d 0 J 0.398 w 0 0 m 0 13.549 l S
 Q
 BT
-/F41 10.9091 Tf 423.773 902.178 Td [(d)22(a)67(t)67(a)-378(type)]TJ
+/F41 10.9091 Tf 423.773 971.119 Td [(d)22(a)67(t)67(a)-378(type)]TJ
 ET
 q
-1 0 0 1 536.768 898.113 cm
+1 0 0 1 536.768 967.055 cm
 []0 d 0 J 0.398 w 0 0 m 0 13.549 l S
 Q
 BT
-/F44 10.9091 Tf 542.746 902.178 Td [(val)]TJ
+/F44 10.9091 Tf 542.746 971.119 Td [(val)]TJ
 ET
 q
-1 0 0 1 639.432 898.113 cm
+1 0 0 1 639.432 967.055 cm
 []0 d 0 J 0.398 w 0 0 m 0 13.549 l S
 Q
 BT
-/F41 10.9091 Tf 645.41 902.178 Td [(def)89(a)22(ul)67(t)]TJ
+/F41 10.9091 Tf 645.41 971.119 Td [(def)89(a)22(ul)67(t)]TJ
 ET
 q
-1 0 0 1 718.51 898.113 cm
+1 0 0 1 718.51 967.055 cm
 []0 d 0 J 0.398 w 0 0 m 0 13.549 l S
 Q
 BT
-/F41 10.9091 Tf 724.488 902.178 Td [(comments)]TJ
+/F41 10.9091 Tf 724.488 971.119 Td [(comments)]TJ
 ET
 q
-1 0 0 1 872.198 898.113 cm
+1 0 0 1 872.198 967.055 cm
 []0 d 0 J 0.398 w 0 0 m 0 13.549 l S
 Q
 q
-1 0 0 1 306.627 897.914 cm
+1 0 0 1 306.627 966.855 cm
 []0 d 0 J 0.398 w 0 0 m 565.57 0 l S
 Q
 q
-1 0 0 1 306.627 870.616 cm
+1 0 0 1 306.627 939.558 cm
 []0 d 0 J 0.398 w 0 0 m 0 27.098 l S
 Q
 BT
-/F44 10.9091 Tf 312.605 888.23 Td [(mld_sub_ovr_)]TJ 0 -13.549 Td [(SUB_OVR)]TJ
+/F44 10.9091 Tf 312.605 957.172 Td [(mld_sub_ovr_)]TJ 0 -13.55 Td [(SUB_OVR)]TJ
 ET
 q
-1 0 0 1 417.795 870.616 cm
+1 0 0 1 417.795 939.558 cm
 []0 d 0 J 0.398 w 0 0 m 0 27.098 l S
 Q
 BT
-/F44 10.9091 Tf 423.773 888.23 Td [(integer)]TJ
+/F44 10.9091 Tf 423.773 957.172 Td [(integer)]TJ
 ET
 q
-1 0 0 1 536.768 870.616 cm
+1 0 0 1 536.768 939.558 cm
 []0 d 0 J 0.398 w 0 0 m 0 27.098 l S
 Q
 BT
-/F15 10.9091 Tf 542.746 888.23 Td [(an)28(y)-333(in)27(t.)-333(n)28(um.)]TJ/F25 10.9091 Tf 67.879 0 Td [(\025)]TJ/F15 10.9091 Tf 11.515 0 Td [(0)]TJ
+/F15 10.9091 Tf 542.746 957.172 Td [(an)28(y)-333(in)27(t.)-333(n)28(um.)]TJ/F25 10.9091 Tf 67.879 0 Td [(\025)]TJ/F15 10.9091 Tf 11.515 0 Td [(0)]TJ
 ET
 q
-1 0 0 1 639.432 870.616 cm
+1 0 0 1 639.432 939.558 cm
 []0 d 0 J 0.398 w 0 0 m 0 27.098 l S
 Q
 BT
-/F15 10.9091 Tf 645.41 888.23 Td [(1)]TJ
+/F15 10.9091 Tf 645.41 957.172 Td [(1)]TJ
 ET
 q
-1 0 0 1 718.51 870.616 cm
+1 0 0 1 718.51 939.558 cm
 []0 d 0 J 0.398 w 0 0 m 0 27.098 l S
 Q
 BT
-/F15 10.9091 Tf 724.488 888.23 Td [(Num)28(b)-28(er)-333(of)-334(o)28(v)28(erlap)-333(la)28(y)27(ers.)]TJ
+/F15 10.9091 Tf 724.488 957.172 Td [(Num)28(b)-28(er)-333(of)-334(o)28(v)28(erlap)-333(la)28(y)27(ers.)]TJ
 ET
 q
-1 0 0 1 872.198 870.616 cm
+1 0 0 1 872.198 939.558 cm
 []0 d 0 J 0.398 w 0 0 m 0 27.098 l S
 Q
 q
-1 0 0 1 306.627 870.417 cm
+1 0 0 1 306.627 939.358 cm
 []0 d 0 J 0.398 w 0 0 m 565.57 0 l S
 Q
 q
-1 0 0 1 306.627 816.021 cm
+1 0 0 1 306.627 884.962 cm
 []0 d 0 J 0.398 w 0 0 m 0 54.197 l S
 Q
 BT
-/F44 10.9091 Tf 312.605 860.733 Td [(mld_sub_restr_)]TJ 0 -13.549 Td [(SUB_RESTR)]TJ
+/F44 10.9091 Tf 312.605 929.675 Td [(mld_sub_restr_)]TJ 0 -13.549 Td [(SUB_RESTR)]TJ
 ET
 q
-1 0 0 1 417.795 816.021 cm
+1 0 0 1 417.795 884.962 cm
 []0 d 0 J 0.398 w 0 0 m 0 54.197 l S
 Q
 BT
-/F44 10.9091 Tf 423.773 860.733 Td [(character\050len=*\051)]TJ
+/F44 10.9091 Tf 423.773 929.675 Td [(character\050len=*\051)]TJ
 ET
 q
-1 0 0 1 536.768 816.021 cm
+1 0 0 1 536.768 884.962 cm
 []0 d 0 J 0.398 w 0 0 m 0 54.197 l S
 Q
 BT
-/F44 10.9091 Tf 542.746 860.733 Td [('HALO')]TJ 0 -13.549 Td [('NONE')]TJ
+/F44 10.9091 Tf 542.746 929.675 Td [('HALO')]TJ 0 -13.549 Td [('NONE')]TJ
 ET
 q
-1 0 0 1 639.432 816.021 cm
+1 0 0 1 639.432 884.962 cm
 []0 d 0 J 0.398 w 0 0 m 0 54.197 l S
 Q
 BT
-/F44 10.9091 Tf 645.41 860.733 Td [('HALO')]TJ
+/F44 10.9091 Tf 645.41 929.675 Td [('HALO')]TJ
 ET
 q
-1 0 0 1 718.51 816.021 cm
+1 0 0 1 718.51 884.962 cm
 []0 d 0 J 0.398 w 0 0 m 0 54.197 l S
 Q
 BT
-/F15 10.9091 Tf 724.488 860.733 Td [(T)28(yp)-28(e)-501(of)-501(restriction)-501(op)-28(erator:)]TJ/F44 10.9091 Tf 0 -13.549 Td [('HALO')]TJ/F15 10.9091 Tf 42.331 0 Td [(for)-731(taking)-730(in)27(to)-730(ac-)]TJ -42.331 -13.549 Td [(coun)28(t)-354(the)-354(o)28(v)28(erlap,)]TJ/F44 10.9091 Tf 90.45 0 Td [('NONE')]TJ/F15 10.9091 Tf 38.222 0 Td [(for)]TJ -128.672 -13.549 Td [(neglecting)-333(it.)]TJ
+/F15 10.9091 Tf 724.488 929.675 Td [(T)28(yp)-28(e)-501(of)-501(restriction)-501(op)-28(erator:)]TJ/F44 10.9091 Tf 0 -13.549 Td [('HALO')]TJ/F15 10.9091 Tf 42.331 0 Td [(for)-731(taking)-730(in)27(to)-730(ac-)]TJ -42.331 -13.55 Td [(coun)28(t)-354(the)-354(o)28(v)28(erlap,)]TJ/F44 10.9091 Tf 90.45 0 Td [('NONE')]TJ/F15 10.9091 Tf 38.222 0 Td [(for)]TJ -128.672 -13.549 Td [(neglecting)-333(it.)]TJ
 ET
 q
-1 0 0 1 872.198 816.021 cm
+1 0 0 1 872.198 884.962 cm
 []0 d 0 J 0.398 w 0 0 m 0 54.197 l S
 Q
 q
-1 0 0 1 306.627 815.822 cm
+1 0 0 1 306.627 884.763 cm
 []0 d 0 J 0.398 w 0 0 m 565.57 0 l S
 Q
 q
-1 0 0 1 306.627 761.426 cm
+1 0 0 1 306.627 830.367 cm
 []0 d 0 J 0.398 w 0 0 m 0 54.197 l S
 Q
 BT
-/F44 10.9091 Tf 312.605 806.138 Td [(mld_sub_prol_)]TJ 0 -13.549 Td [(SUB_PROL)]TJ
+/F44 10.9091 Tf 312.605 875.08 Td [(mld_sub_prol_)]TJ 0 -13.55 Td [(SUB_PROL)]TJ
 ET
 q
-1 0 0 1 417.795 761.426 cm
+1 0 0 1 417.795 830.367 cm
 []0 d 0 J 0.398 w 0 0 m 0 54.197 l S
 Q
 BT
-/F44 10.9091 Tf 423.773 806.138 Td [(character\050len=*\051)]TJ
+/F44 10.9091 Tf 423.773 875.08 Td [(character\050len=*\051)]TJ
 ET
 q
-1 0 0 1 536.768 761.426 cm
+1 0 0 1 536.768 830.367 cm
 []0 d 0 J 0.398 w 0 0 m 0 54.197 l S
 Q
 BT
-/F44 10.9091 Tf 542.746 806.138 Td [('SUM')]TJ 0 -13.549 Td [('NONE')]TJ
+/F44 10.9091 Tf 542.746 875.08 Td [('SUM')]TJ 0 -13.55 Td [('NONE')]TJ
 ET
 q
-1 0 0 1 639.432 761.426 cm
+1 0 0 1 639.432 830.367 cm
 []0 d 0 J 0.398 w 0 0 m 0 54.197 l S
 Q
 BT
-/F44 10.9091 Tf 645.41 806.138 Td [('NONE')]TJ
+/F44 10.9091 Tf 645.41 875.08 Td [('NONE')]TJ
 ET
 q
-1 0 0 1 718.51 761.426 cm
+1 0 0 1 718.51 830.367 cm
 []0 d 0 J 0.398 w 0 0 m 0 54.197 l S
 Q
 BT
-/F15 10.9091 Tf 724.488 806.138 Td [(T)28(yp)-28(e)-560(of)-561(prolon)1(gation)-561(op)-27(e)-1(r)1(a-)]TJ 0 -13.549 Td [(tor:)]TJ/F44 10.9091 Tf 21.445 0 Td [('SUM')]TJ/F15 10.9091 Tf 31.465 0 Td [(for)-259(adding)-260(th)1(e)-260(con-)]TJ -52.91 -13.549 Td [(tributions)-622(from)-623(t)1(he)-623(o)28(v)28(erlap,)]TJ/F44 10.9091 Tf 0 -13.55 Td [('NONE')]TJ/F15 10.9091 Tf 37.999 0 Td [(for)-333(neglecting)-334(them.)]TJ
+/F15 10.9091 Tf 724.488 875.08 Td [(T)28(yp)-28(e)-560(of)-561(prolon)1(gation)-561(op)-27(e)-1(r)1(a-)]TJ 0 -13.55 Td [(tor:)]TJ/F44 10.9091 Tf 21.445 0 Td [('SUM')]TJ/F15 10.9091 Tf 31.465 0 Td [(for)-259(adding)-260(th)1(e)-260(con-)]TJ -52.91 -13.549 Td [(tributions)-622(from)-623(t)1(he)-623(o)28(v)28(erlap,)]TJ/F44 10.9091 Tf 0 -13.549 Td [('NONE')]TJ/F15 10.9091 Tf 37.999 0 Td [(for)-333(neglecting)-334(them.)]TJ
 ET
 q
-1 0 0 1 872.198 761.426 cm
+1 0 0 1 872.198 830.367 cm
 []0 d 0 J 0.398 w 0 0 m 0 54.197 l S
 Q
 q
-1 0 0 1 306.627 761.226 cm
+1 0 0 1 306.627 830.168 cm
 []0 d 0 J 0.398 w 0 0 m 565.57 0 l S
 Q
 q
-1 0 0 1 306.627 693.281 cm
-[]0 d 0 J 0.398 w 0 0 m 0 67.746 l S
+1 0 0 1 306.627 735.124 cm
+[]0 d 0 J 0.398 w 0 0 m 0 94.844 l S
 Q
 BT
-/F44 10.9091 Tf 312.605 751.543 Td [(mld_sub_solve_)]TJ 0 -13.549 Td [(SUB_SOLVE)]TJ
+/F44 10.9091 Tf 312.605 820.484 Td [(mld_sub_solve_)]TJ 0 -13.549 Td [(SUB_SOLVE)]TJ
 ET
 q
-1 0 0 1 417.795 693.281 cm
-[]0 d 0 J 0.398 w 0 0 m 0 67.746 l S
+1 0 0 1 417.795 735.124 cm
+[]0 d 0 J 0.398 w 0 0 m 0 94.844 l S
 Q
 BT
-/F44 10.9091 Tf 423.773 751.543 Td [(character\050len=*\051)]TJ
+/F44 10.9091 Tf 423.773 820.484 Td [(character\050len=*\051)]TJ
 ET
 q
-1 0 0 1 536.768 693.281 cm
-[]0 d 0 J 0.398 w 0 0 m 0 67.746 l S
+1 0 0 1 536.768 735.124 cm
+[]0 d 0 J 0.398 w 0 0 m 0 94.844 l S
 Q
 BT
-/F44 10.9091 Tf 542.746 751.543 Td [('ILU')]TJ 0 -13.549 Td [('MILU')]TJ 0 -13.55 Td [('ILUT')]TJ 0 -13.549 Td [('UMF')]TJ 0 -13.549 Td [('SLU')]TJ
+/F44 10.9091 Tf 542.746 820.484 Td [('DIAG')]TJ 0 -13.549 Td [('GS')]TJ 0 -13.549 Td [('ILU')]TJ 0 -13.549 Td [('MILU')]TJ 0 -13.55 Td [('ILUT')]TJ 0 -13.549 Td [('UMF')]TJ 0 -13.549 Td [('SLU')]TJ
 ET
 q
-1 0 0 1 639.432 693.281 cm
-[]0 d 0 J 0.398 w 0 0 m 0 67.746 l S
+1 0 0 1 639.432 735.124 cm
+[]0 d 0 J 0.398 w 0 0 m 0 94.844 l S
 Q
 BT
-/F44 10.9091 Tf 645.41 751.543 Td [('ILU')]TJ
+/F44 10.9091 Tf 645.41 820.484 Td [('ILU')]TJ
 ET
 q
-1 0 0 1 718.51 693.281 cm
-[]0 d 0 J 0.398 w 0 0 m 0 67.746 l S
+1 0 0 1 718.51 735.124 cm
+[]0 d 0 J 0.398 w 0 0 m 0 94.844 l S
 Q
 BT
-/F15 10.9091 Tf 724.488 751.543 Td [(Prede\014ned)-1775(lo)-28(cal)-1775(solv)28(er:)]TJ 0 -13.549 Td [(ILU\050)]TJ/F22 10.9091 Tf 23.182 0 Td [(p)]TJ/F15 10.9091 Tf 5.488 0 Td [(\051,)-693(MILU\050)]TJ/F22 10.9091 Tf 48.012 0 Td [(p)]TJ/F15 10.9091 Tf 5.489 0 Td [(\051,)-693(ILU\050)]TJ/F22 10.9091 Tf 38.012 0 Td [(p;)-167(t)]TJ/F15 10.9091 Tf 14.276 0 Td [(\051,)]TJ -134.459 -13.55 Td [(LU)-931(from)-932(UMFP)83(A)28(CK,)-932(LU)]TJ 0 -13.549 Td [(from)-460(Sup)-28(erLU)-460(\050plus)-461(tr)1(iangu-)]TJ 0 -13.549 Td [(lar)-333(solv)28(e)-1(\051.)]TJ
+/F15 10.9091 Tf 724.488 820.484 Td [(Prede\014ned)-651(lo)-28(cal)-652(solv)28(er:)-1081(di-)]TJ 0 -13.549 Td [(agonal)-541(s)-1(calin)1(g,)-594(Gauss-Seidel,)]TJ 0 -13.549 Td [(ILU\050)]TJ/F22 10.9091 Tf 23.182 0 Td [(p)]TJ/F15 10.9091 Tf 5.488 0 Td [(\051,)-693(MILU\050)]TJ/F22 10.9091 Tf 48.012 0 Td [(p)]TJ/F15 10.9091 Tf 5.489 0 Td [(\051,)-693(ILU\050)]TJ/F22 10.9091 Tf 38.012 0 Td [(p;)-167(t)]TJ/F15 10.9091 Tf 14.276 0 Td [(\051,)]TJ -134.459 -13.549 Td [(LU)-931(from)-932(UMFP)83(A)28(CK,)-932(LU)]TJ 0 -13.55 Td [(from)-460(Sup)-28(erLU)-460(\050plus)-461(tr)1(iangu-)]TJ 0 -13.549 Td [(lar)-333(solv)28(e)-1(\051.)]TJ
 ET
 q
-1 0 0 1 872.198 693.281 cm
-[]0 d 0 J 0.398 w 0 0 m 0 67.746 l S
+1 0 0 1 872.198 735.124 cm
+[]0 d 0 J 0.398 w 0 0 m 0 94.844 l S
 Q
 q
-1 0 0 1 306.627 693.082 cm
+1 0 0 1 306.627 734.925 cm
 []0 d 0 J 0.398 w 0 0 m 565.57 0 l S
 Q
 q
-1 0 0 1 306.627 665.784 cm
+1 0 0 1 306.627 707.627 cm
 []0 d 0 J 0.398 w 0 0 m 0 27.098 l S
 Q
 BT
-/F44 10.9091 Tf 312.605 683.398 Td [(mld_sub_fillin_)]TJ 0 -13.549 Td [(SUB_FILLIN)]TJ
+/F44 10.9091 Tf 312.605 725.241 Td [(mld_sub_fillin_)]TJ 0 -13.549 Td [(SUB_FILLIN)]TJ
 ET
 q
-1 0 0 1 417.795 665.784 cm
+1 0 0 1 417.795 707.627 cm
 []0 d 0 J 0.398 w 0 0 m 0 27.098 l S
 Q
 BT
-/F44 10.9091 Tf 423.773 683.398 Td [(integer)]TJ
+/F44 10.9091 Tf 423.773 725.241 Td [(integer)]TJ
 ET
 q
-1 0 0 1 536.768 665.784 cm
+1 0 0 1 536.768 707.627 cm
 []0 d 0 J 0.398 w 0 0 m 0 27.098 l S
 Q
 BT
-/F15 10.9091 Tf 542.746 683.398 Td [(An)28(y)-333(in)27(t.)-333(n)28(um.)]TJ/F25 10.9091 Tf 70.606 0 Td [(\025)]TJ/F15 10.9091 Tf 11.515 0 Td [(0)]TJ
+/F15 10.9091 Tf 542.746 725.241 Td [(An)28(y)-333(in)27(t.)-333(n)28(um.)]TJ/F25 10.9091 Tf 70.606 0 Td [(\025)]TJ/F15 10.9091 Tf 11.515 0 Td [(0)]TJ
 ET
 q
-1 0 0 1 639.432 665.784 cm
+1 0 0 1 639.432 707.627 cm
 []0 d 0 J 0.398 w 0 0 m 0 27.098 l S
 Q
 BT
-/F15 10.9091 Tf 645.41 683.398 Td [(0)]TJ
+/F15 10.9091 Tf 645.41 725.241 Td [(0)]TJ
 ET
 q
-1 0 0 1 718.51 665.784 cm
+1 0 0 1 718.51 707.627 cm
 []0 d 0 J 0.398 w 0 0 m 0 27.098 l S
 Q
 BT
-/F15 10.9091 Tf 724.488 683.398 Td [(Fill-in)-550(lev)27(el)]TJ/F22 10.9091 Tf 62.164 0 Td [(p)]TJ/F15 10.9091 Tf 11.495 0 Td [(of)-551(th)1(e)-551(incom-)]TJ -73.659 -13.549 Td [(plete)-333(LU)-334(factorizations.)]TJ
+/F15 10.9091 Tf 724.488 725.241 Td [(Fill-in)-550(lev)27(el)]TJ/F22 10.9091 Tf 62.164 0 Td [(p)]TJ/F15 10.9091 Tf 11.495 0 Td [(of)-551(th)1(e)-551(incom-)]TJ -73.659 -13.549 Td [(plete)-333(LU)-334(factorizations.)]TJ
 ET
 q
-1 0 0 1 872.198 665.784 cm
+1 0 0 1 872.198 707.627 cm
 []0 d 0 J 0.398 w 0 0 m 0 27.098 l S
 Q
 q
-1 0 0 1 306.627 665.585 cm
+1 0 0 1 306.627 707.428 cm
 []0 d 0 J 0.398 w 0 0 m 565.57 0 l S
 Q
 q
-1 0 0 1 306.627 638.287 cm
+1 0 0 1 306.627 680.131 cm
 []0 d 0 J 0.398 w 0 0 m 0 27.098 l S
 Q
 BT
-/F44 10.9091 Tf 312.605 655.901 Td [(mld_sub_iluthrs_)]TJ 0 -13.549 Td [(SUB_ILUTHRS)]TJ
+/F44 10.9091 Tf 312.605 697.745 Td [(mld_sub_iluthrs_)]TJ 0 -13.55 Td [(SUB_ILUTHRS)]TJ
 ET
 q
-1 0 0 1 417.795 638.287 cm
+1 0 0 1 417.795 680.131 cm
 []0 d 0 J 0.398 w 0 0 m 0 27.098 l S
 Q
 BT
-/F44 10.9091 Tf 423.773 655.901 Td [(real\050)]TJ/F18 10.9091 Tf 28.636 0 Td [(kind)]TJ
+/F44 10.9091 Tf 423.773 697.745 Td [(real\050)]TJ/F18 10.9091 Tf 28.636 0 Td [(kind)]TJ
 ET
 q
-1 0 0 1 473.15 656.101 cm
+1 0 0 1 473.15 697.944 cm
 []0 d 0 J 0.398 w 0 0 m 3.345 0 l S
 Q
 BT
-/F18 10.9091 Tf 476.496 655.901 Td [(p)51(ar)51(ameter)]TJ/F44 10.9091 Tf 48.568 0 Td [(\051)]TJ
+/F18 10.9091 Tf 476.496 697.745 Td [(p)51(ar)51(ameter)]TJ/F44 10.9091 Tf 48.568 0 Td [(\051)]TJ
 ET
 q
-1 0 0 1 536.768 638.287 cm
+1 0 0 1 536.768 680.131 cm
 []0 d 0 J 0.398 w 0 0 m 0 27.098 l S
 Q
 BT
-/F15 10.9091 Tf 542.746 655.901 Td [(An)28(y)-333(re)-1(al)-333(n)28(um.)]TJ/F25 10.9091 Tf 72.152 0 Td [(\025)]TJ/F15 10.9091 Tf 11.515 0 Td [(0)]TJ
+/F15 10.9091 Tf 542.746 697.745 Td [(An)28(y)-333(re)-1(al)-333(n)28(um.)]TJ/F25 10.9091 Tf 72.152 0 Td [(\025)]TJ/F15 10.9091 Tf 11.515 0 Td [(0)]TJ
 ET
 q
-1 0 0 1 639.432 638.287 cm
+1 0 0 1 639.432 680.131 cm
 []0 d 0 J 0.398 w 0 0 m 0 27.098 l S
 Q
 BT
-/F15 10.9091 Tf 645.41 655.901 Td [(0)]TJ
+/F15 10.9091 Tf 645.41 697.745 Td [(0)]TJ
 ET
 q
-1 0 0 1 718.51 638.287 cm
+1 0 0 1 718.51 680.131 cm
 []0 d 0 J 0.398 w 0 0 m 0 27.098 l S
 Q
 BT
-/F15 10.9091 Tf 724.488 655.901 Td [(Drop)-1063(tol)1(e)-1(r)1(ance)]TJ/F22 10.9091 Tf 90.366 0 Td [(t)]TJ/F15 10.9091 Tf 15.532 0 Td [(in)-1063(t)1(he)]TJ -105.898 -13.549 Td [(ILU\050)]TJ/F22 10.9091 Tf 23.182 0 Td [(p;)-167(t)]TJ/F15 10.9091 Tf 14.276 0 Td [(\051)-333(factorization.)]TJ
+/F15 10.9091 Tf 724.488 697.745 Td [(Drop)-1063(tol)1(e)-1(r)1(ance)]TJ/F22 10.9091 Tf 90.366 0 Td [(t)]TJ/F15 10.9091 Tf 15.532 0 Td [(in)-1063(t)1(he)]TJ -105.898 -13.55 Td [(ILU\050)]TJ/F22 10.9091 Tf 23.182 0 Td [(p;)-167(t)]TJ/F15 10.9091 Tf 14.276 0 Td [(\051)-333(factorization.)]TJ
 ET
 q
-1 0 0 1 872.198 638.287 cm
+1 0 0 1 872.198 680.131 cm
 []0 d 0 J 0.398 w 0 0 m 0 27.098 l S
 Q
 q
-1 0 0 1 306.627 638.088 cm
+1 0 0 1 306.627 679.931 cm
 []0 d 0 J 0.398 w 0 0 m 565.57 0 l S
 Q
 q
-1 0 0 1 306.627 556.594 cm
+1 0 0 1 306.627 598.437 cm
 []0 d 0 J 0.398 w 0 0 m 0 81.295 l S
 Q
 BT
-/F44 10.9091 Tf 312.605 628.405 Td [(mld_sub_ren_)]TJ 0 -13.55 Td [(SUB_REN)]TJ
+/F44 10.9091 Tf 312.605 670.248 Td [(mld_sub_ren_)]TJ 0 -13.55 Td [(SUB_REN)]TJ
 ET
 q
-1 0 0 1 417.795 556.594 cm
+1 0 0 1 417.795 598.437 cm
 []0 d 0 J 0.398 w 0 0 m 0 81.295 l S
 Q
 BT
-/F44 10.9091 Tf 423.773 628.405 Td [(character\050len=*\051)]TJ
+/F44 10.9091 Tf 423.773 670.248 Td [(character\050len=*\051)]TJ
 ET
 q
-1 0 0 1 536.768 556.594 cm
+1 0 0 1 536.768 598.437 cm
 []0 d 0 J 0.398 w 0 0 m 0 81.295 l S
 Q
 BT
-/F44 10.9091 Tf 542.746 628.405 Td [('RENUM)]TJ
+/F44 10.9091 Tf 542.746 670.248 Td [('RENUM)]TJ
 ET
 q
-1 0 0 1 577.796 628.604 cm
+1 0 0 1 577.796 670.447 cm
 []0 d 0 J 0.398 w 0 0 m 3.436 0 l S
 Q
 BT
-/F44 10.9091 Tf 581.233 628.405 Td [(NONE')]TJ -38.487 -13.55 Td [('RENUM)]TJ
+/F44 10.9091 Tf 581.233 670.248 Td [(NONE')]TJ -38.487 -13.55 Td [('RENUM)]TJ
 ET
 q
-1 0 0 1 577.796 615.055 cm
+1 0 0 1 577.796 656.898 cm
 []0 d 0 J 0.398 w 0 0 m 3.436 0 l S
 Q
 BT
-/F44 10.9091 Tf 581.233 614.855 Td [(GLOBAL')]TJ
+/F44 10.9091 Tf 581.233 656.698 Td [(GLOBAL')]TJ
 ET
 q
-1 0 0 1 639.432 556.594 cm
+1 0 0 1 639.432 598.437 cm
 []0 d 0 J 0.398 w 0 0 m 0 81.295 l S
 Q
 BT
-/F44 10.9091 Tf 645.41 628.405 Td [('RENUM)]TJ
+/F44 10.9091 Tf 645.41 670.248 Td [('RENUM)]TJ
 ET
 q
-1 0 0 1 680.46 628.604 cm
+1 0 0 1 680.46 670.447 cm
 []0 d 0 J 0.398 w 0 0 m 3.436 0 l S
 Q
 BT
-/F44 10.9091 Tf 683.897 628.405 Td [(NONE')]TJ
+/F44 10.9091 Tf 683.897 670.248 Td [(NONE')]TJ
 ET
 q
-1 0 0 1 718.51 556.594 cm
+1 0 0 1 718.51 598.437 cm
 []0 d 0 J 0.398 w 0 0 m 0 81.295 l S
 Q
 BT
-/F15 10.9091 Tf 724.488 628.405 Td [(Ro)28(w)-257(and)-256(column)-256(reordering)-256(of)]TJ 0 -13.55 Td [(the)-423(lo)-28(cal)-422(s)-1(u)1(bmatrice)-1(s:)-623(no)-423(re-)]TJ 0 -13.549 Td [(ordering,)-248(reordering)-227(according)]TJ 0 -13.549 Td [(to)-253(the)-254(global)-253(n)28(um)27(b)-27(ering)-254(of)-253(the)]TJ 0 -13.549 Td [(ro)28(ws)-229(and)-229(columns)-229(of)-229(the)-229(whole)]TJ 0 -13.549 Td [(matrix.)]TJ
+/F15 10.9091 Tf 724.488 670.248 Td [(Ro)28(w)-257(and)-256(column)-256(reordering)-256(of)]TJ 0 -13.55 Td [(the)-423(lo)-28(cal)-422(s)-1(u)1(bmatrice)-1(s:)-623(no)-423(re-)]TJ 0 -13.549 Td [(ordering,)-248(reordering)-227(according)]TJ 0 -13.549 Td [(to)-253(the)-254(global)-253(n)28(um)27(b)-27(ering)-254(of)-253(the)]TJ 0 -13.549 Td [(ro)28(ws)-229(and)-229(columns)-229(of)-229(the)-229(whole)]TJ 0 -13.549 Td [(matrix.)]TJ
 ET
 q
-1 0 0 1 872.198 556.594 cm
+1 0 0 1 872.198 598.437 cm
 []0 d 0 J 0.398 w 0 0 m 0 81.295 l S
 Q
 q
-1 0 0 1 306.627 556.395 cm
+1 0 0 1 306.627 557.789 cm
+[]0 d 0 J 0.398 w 0 0 m 0 40.648 l S
+Q
+BT
+/F44 10.9091 Tf 312.605 588.952 Td [(mld_solver_sweeps_)]TJ 0 -27.098 Td [(SOLVER_SWEEPS)]TJ
+ET
+q
+1 0 0 1 417.795 557.789 cm
+[]0 d 0 J 0.398 w 0 0 m 0 40.648 l S
+Q
+BT
+/F44 10.9091 Tf 423.773 588.952 Td [(integer)]TJ
+ET
+q
+1 0 0 1 536.768 557.789 cm
+[]0 d 0 J 0.398 w 0 0 m 0 40.648 l S
+Q
+BT
+/F15 10.9091 Tf 542.746 588.952 Td [(An)28(y)-333(in)27(t.)-333(n)28(um.)]TJ/F25 10.9091 Tf 70.606 0 Td [(\025)]TJ/F15 10.9091 Tf 11.515 0 Td [(1)]TJ
+ET
+q
+1 0 0 1 639.432 557.789 cm
+[]0 d 0 J 0.398 w 0 0 m 0 40.648 l S
+Q
+BT
+/F15 10.9091 Tf 645.41 588.952 Td [(1)]TJ
+ET
+q
+1 0 0 1 718.51 557.789 cm
+[]0 d 0 J 0.398 w 0 0 m 0 40.648 l S
+Q
+BT
+/F15 10.9091 Tf 724.488 588.952 Td [(Num)28(b)-28(er)-674(of)-673(sw)27(eeps)-673(for)-674(iter-)]TJ 0 -13.549 Td [(ativ)28(e)-669(lo)-28(cal)-668(solv)27(er)-668(\050curren)28(tly)]TJ 0 -13.549 Td [(only)-333(Gauss-Seidel\051.)]TJ
+ET
+q
+1 0 0 1 872.198 557.789 cm
+[]0 d 0 J 0.398 w 0 0 m 0 40.648 l S
+Q
+q
+1 0 0 1 306.627 557.59 cm
+[]0 d 0 J 0.398 w 0 0 m 565.57 0 l S
+Q
+q
+1 0 0 1 306.627 489.645 cm
+[]0 d 0 J 0.398 w 0 0 m 0 67.746 l S
+Q
+BT
+/F44 10.9091 Tf 312.605 547.906 Td [(mld_solver_eps_)]TJ 0 -13.549 Td [(SOLVER_EPS)]TJ
+ET
+q
+1 0 0 1 417.795 489.645 cm
+[]0 d 0 J 0.398 w 0 0 m 0 67.746 l S
+Q
+BT
+/F44 10.9091 Tf 423.773 547.906 Td [(real)]TJ
+ET
+q
+1 0 0 1 536.768 489.645 cm
+[]0 d 0 J 0.398 w 0 0 m 0 67.746 l S
+Q
+BT
+/F15 10.9091 Tf 542.746 547.906 Td [(An)28(y)-333(re)-1(al)-333(n)28(um.)]TJ
+ET
+q
+1 0 0 1 639.432 489.645 cm
+[]0 d 0 J 0.398 w 0 0 m 0 67.746 l S
+Q
+BT
+/F15 10.9091 Tf 645.41 547.906 Td [(0)]TJ
+ET
+q
+1 0 0 1 718.51 489.645 cm
+[]0 d 0 J 0.398 w 0 0 m 0 67.746 l S
+Q
+BT
+/F15 10.9091 Tf 724.488 547.906 Td [(Stopping)-707(tolerance)-708(for)-707(iter-)]TJ 0 -13.549 Td [(ativ)28(e)-669(lo)-28(cal)-668(solv)27(er)-668(\050curren)28(tly)]TJ 0 -13.549 Td [(only)-629(Gauss-Seidel\051;)-778(if)]TJ/F25 10.9091 Tf 116.349 0 Td [(\024)]TJ/F15 10.9091 Tf 16.898 0 Td [(0,)]TJ -133.247 -13.549 Td [(then)-1253(p)-28(erform)-1253(presp)-27(eci\014ed)]TJ 0 -13.549 Td [(n)28(um)28(b)-28(er)-333(of)-334(iterations.)]TJ
+ET
+q
+1 0 0 1 872.198 489.645 cm
+[]0 d 0 J 0.398 w 0 0 m 0 67.746 l S
+Q
+q
+1 0 0 1 306.627 489.446 cm
+[]0 d 0 J 0.398 w 0 0 m 565.57 0 l S
+Q
+q
+1 0 0 1 306.627 487.453 cm
 []0 d 0 J 0.398 w 0 0 m 565.57 0 l S
 Q
 0 g 0 G
 BT
-/F15 10.9091 Tf 405.352 525.705 Td [(T)83(able)-333(3:)-444(P)27(arameters)-333(de\014ning)-333(the)-334(on)1(e)-1(-lev)28(el)-333(preconditioner)-333(used)-334(as)-333(smo)-28(other.)]TJ
+/F15 10.9091 Tf 405.352 456.764 Td [(T)83(able)-333(3:)-444(P)27(arameters)-333(de\014ning)-333(the)-334(on)1(e)-1(-lev)28(el)-333(preconditioner)-333(used)-334(as)-333(smo)-28(other.)]TJ
 0 g 0 G
 0 g 0 G
 ET
@@ -5357,7 +5443,7 @@ stream
 >>
 % 360 0 obj
 <<
-/D [434 0 R /XYZ 366.766 4378.384 null]
+/D [434 0 R /XYZ 366.767 4378.384 null]
 >>
 % 433 0 obj
 <<
@@ -5378,7 +5464,7 @@ stream
 >>
 % 442 0 obj
 <<
-/D [439 0 R /XYZ -3963.913 559.035 null]
+/D [439 0 R /XYZ -4032.854 559.035 null]
 >>
 % 438 0 obj
 <<
@@ -5408,7 +5494,7 @@ stream
 >>
 % 448 0 obj
 <<
-/D [445 0 R /XYZ 530.237 4405.505 null]
+/D [445 0 R /XYZ 530.198 4405.505 null]
 >>
 % 444 0 obj
 <<
@@ -7234,10 +7320,10 @@ endstream
 endobj
 540 0 obj
 <<
-/Length1 1647
-/Length2 7851
+/Length1 1669
+/Length2 7977
 /Length3 0
-/Length 9498      
+/Length 9646      
 >>
 stream
 %!PS-AdobeFont-1.0: CMSY10 003.002
@@ -7257,7 +7343,7 @@ FontDirectory/CMSY10 known{/CMSY10 findfont dup/UniqueID known{dup
 11 dict begin
 /FontType 1 def
 /FontMatrix [0.001 0 0 0.001 0 0 ]readonly def
-/FontName /WYYAGH+CMSY10 def
+/FontName /PLWDAY+CMSY10 def
 /FontBBox {-29 -960 1116 775 }readonly def
 /PaintType 0 def
 /FontInfo 9 dict dup begin
@@ -7281,6 +7367,7 @@ dup 103 /braceright put
 dup 15 /bullet put
 dup 50 /element put
 dup 21 /greaterequal put
+dup 20 /lessequal put
 dup 0 /minus put
 dup 54 /negationslash put
 dup 26 /propersubset put
@@ -7304,17 +7391,13 @@ A5
 zúåT¸ÞúË0ÙÝb#ö›
[›zyÏ’Õ$e1¬ÛkÕ5—dRAì`(Å…@{<7B>=yW>¿ðˆä>ú‘¯çy@æÅÿ¹Uáú;”œ±=ÃH;Çc{–ð<yBb7©l!&x<>ºð@©‡Ôï•-d¡Û½uäê\ÂÍæ"˜÷'uîÊ|!>Û3û0´<30>ü£^Õ®”)Õ™-{·O<C2B7>¹°OÖú2R¨Â´xÈèëˆQÅôwBgБâ¼ýSý»ž›"תl›bðöz;¹±˜IyÕ\E×ÁËê‰pr¨kM·@@Âe&À 3•˜lÃʉtS÷æÄ%<1D>í)è/pRÇøÞÙÀª"2µu
 ˜Ee’5ø+ÄÙ°s§Zò'±¹§õ˜àçbˆQ²ôð°SlŸjÅö <20>zvÞ€â~%.Kz\N¡Š„çüã÷Eê,ãeV2Ò=Mo)-<2D>JcÇVà¤:ÒUÙ2Ѹ>É1ˆ×È“"¿QWº ƒOýõÑ‘4*Ozûº.c¸ójc”«Y&üßÜØØÏ<13>ìê=ë7”Y™´öøÈðê,A5nÊIßFŠ•/œ?üÌŠ™á
]as.´N/ú	ç,Mceû•wW³¢0Ò#¸]ÊÙªÕü3wœ]ŠœŽ2“«À‘Öf˨øÅn'ä_¼	MuH½
 OëM‹lC7Þh®	?p|`e÷õ‘š(“ºÓ›·|*è{dgî²HÅPNDY®LèQ«´Ÿ›íðŠô‡œê£sRÑd
=©ÁÒ…ñ^\#ZHdŽ¤àõÑ;
-¶D¯Å+Z†”ÊQbË7QÑÚÁ±N¹Ì§ÝP<C39D>#?ã8 ©Õ“F,ùÃÒ(¼÷韓ðÀv>Ô%¢Uú§ô9Í«ÙÕN¤ÞDÔÄâÁÿOžQkç‰O´L€?<3F>™¿{™(-ÝaÇ£úüUs§êÑ™N·#·^Ã&÷Ô†Ÿ[üMJ¢¥Bð,± ÂVËðÌásP÷eQ­O»±ÇB<C387>Æ,s“%v6ù2å3}.D|`,;(`
©tQÉ)l!<00>TÀr+ìuì=Ù’[ŽÛÍéÏÒ¼4ß`hÇ-?àß2òD#<23>eã+®%-¹t‘U>kx²ß£**åó߾؂ÆÌXšý?Ç`ªîÑ]!&dè)Bgzó·á‹ùúìF<C3AC>mÊw¼GVÂM\nŸIç{Œý`Ýõ"S$OÙ05<30>A»G<C2BB>XVÐÛ<1C>NêIQkþƒι,àâyŠmoÊ-D$S8€dLÊ-OwãQu7<75>@¡g53vNµ<6á’¥¶)9}˜†'¿u~~ïVf<*ß‹¬8¢ôó<C3B4>Hº/2ÎÁ—ñT8°¬ñ¦;wØ‹Z¸§ôsxVÀ‘Âe¤í€vïƒ6qÕ/Ûðœ4w-7ÝkLPŠÍËCCxÝíó#u÷º‰ZÝ+‡™<E280A1>ƒ…à²û6#˜	b4R?›Ú‚ì$SÇm ã:ï|íÏ›<C38F>%¥Ló¿Øy
ÊØj#‘‘Â&+°I*p`r¸"O M6òº¯Ó–‚l«GàMèuT껵•ÜÒOâKlbÄ“
-_EG‹´SÞòËõÎù‡ùÒ/þ,ä+‡iÕø„ŽÆ¡ò÷Ûj£ÌFüŠ<¸0U [z)dÂG`‘õNÆdçÂZ=çåJ;Qâƒý<02>ìUNîß%¿‡ü¯xõg¢®À?Å—˜k
-ZQ³Dx/0~oÙ4<C399>Ï¡\Cª<>~ž¦qDâðè99¯ðÐpe:\
-ÕNÍ~E}Á¯Ö1Pº®ß…\Êmsy£ÄË;*™™[£@òk~
#è<11>MEø˜†hÁ!ÒEªi*X,H<>ÓË%ªµÃÏäásX…æQ kèPœtÿëJ‡YlÉÓ}tü¶@«ÝF›hæêá0îbhÈ'@¼<1F>Ìl„i:WæM†ì•ÀÄ<C380>Ý<EFBFBD>>$»MxŽÑí)ipÎné[_G×z:%^¯
›Â	<09>¥AÀÏT‡$,ú˜w°dí<64>koÖMÈéü9ÀΘ¿.}ÆÇöïDL›R<1A>g¨³¤u]”añ£g–{ã‹#Zd1[<´Ð^+hk
-
-ø~
çæ<C3A7>ȈQN£QýÓHÜ•‰4¬wã B2_³ÑF)ÙûY?Þ÷Çï,7OÎ¥7…neÓ½àÚ@Sd¥ÌâÇëká¾JF÷C"§_™Úl/²VÅñ5“–ÌÞåè0D£T·‡€<E280A1>µÒ¦Wöòå2(Š\¦Œíúí&tõôZÊQ<~Žß¾\8~öˆJ§ó(ä
NˆcñÄØ2™/’X¨t§Ž0`(Û–’ž³~À“4½U¥y1.¬Ø¬ uÈ"ìXÍBÊ«õiÙ;kþ×ÝQ)Cdº6¬Ø¢Þ—þ|Ýø
¨‹™È*Gs,÷ÄœJgÜçñЃڤã SÜÞ=WXòdeòn­«FB©a¯†H4+)'L€9ŠefD*Ô“¢¯õ¢!ÑæUýW‹ç6]‡doz9[ó;¦ÚͯñŸ9>*¾Ñ Cx»eð/ü•ö¸7a°Aú
ß›á€Á°þó¦•Wc1þñzsã<>l3)Ìn:Kê	‘ä²Ï³ã<>w†9j)½yý.ón¨>çû7V!àj«•Òþ+©(zîrØÅ<C398>>ܦÏâ÷š&¹bþbVòç
 bú6r`£<>ÉF¢R¹!†ë¢"H›îÃ-<2D>¡¦*
-oŽ‚Á3á1ò¸9Á(ÓÀ<C393>¥5[¦rHJÏ<Óò&8ßxs12PIÕ­Âä.2.µ‚¬”í[Èê­ÂD&ä<XÅö˜€
ê{¾ü‹¤
-m‰LcPA9¨Ýš éyÓSÏR†½¼Òý\»d°[ë.~hf`¶Œø`
 ¸	ÏâdÊ0ï »9šãl÷h9†‚!¤’«ã{Wñ†	6´ cåÔ¶ç%І»á¯wH2Øõ`•òK—JÖp‹|Ú)üß•j·°žóZÞJl6««¯¸°ú›Ô“˜Ð1BÐ<42>°"/3¥¹ìè¯ü—>)Th‘¦ÏG£ñaŒÜwL¬d<R>…çYýÃój+#®&u˜žo+ª§?ùvßÂ|“2s`(´v;¦å‹Æñúæýý-#£ô|“úL³oUX¾NÂ¥‡$î*ã𫦙Ôqá„8øLÁõozIKÔT<06>'ãyÎŒ=h¦æ2*eø9%y
ai½‹Ò)È€_Xe½û
-_<EFBFBD>po­Á{Ķ…Û?|Ën¦¿òîžM°2isB{Z‘ѽT$¼®P>LàE²ÐUzˆÓž.’T™È
ç¦ÝN¡ÒIA¡©´Ü”4ôþ¢ížŸ©Ë8xYsJª¿òš°L‹D´¶}¨+2í•+C~}`°˜æØ⨱ˆÊc7öºïô7¿³Ù=ök&—YÎ…I•Bø
r_±äZÌ4b¸c#³à\ëK<>A¨—Q „7<Ý9Ðjû(t°æ7
-Cï”Qœ¦4ÞFÅ@VfUÀ¿œëQÀו
˺8ÀÛ‘á—¡ÃFó	ÈD]²÷e’'sùdM¸ØRœúj?&_z¥¢ès²Ï¦ÿ&á#0Ê°w¿
-Äõ
±K8ë7‘œm‹
+¶D¯Å+Z†”ÊQbË7QÑÚÁ±N¹Ì§ÝP<C39D>#?ã8 ©Õ“F,ùÃÒ(¼÷韓ðÀv>Ô%¢Uú§ô9Í«ÙÕN¤ÞDÔÄâÁÿOžQkç‰O´L€?<3F>™¿{™(-ÝaÇ£úüUs§êÑ™N·#·^Ã&÷Ô†Ÿ[üMJ¢¥Bð,± ÂVËðÌásP÷eQ­O»±ÇB<C387>Æ,s“%v6ù2å3}.D|`,;(`
©tQÉ)l!<00>TÀr+ìuì=Ù’[ŽÛÍéÏÒ¼4ß`hÇ-?àß2òD#<23>eã+®%-¹t‘U>kx²ß£**åó߾؂ÆÌXšý?Ç`ªîÑ]!&dè)Bgzó·á‹ùúìF<C3AC>mÊw¼GVÂM\nŸIç{Œý`Ýõ"S$OÙ05<30>A»G<C2BB>XVÐÛ<1C>NêIPtžj‘É‚Éj(k_•ß˜ÆÀ€7Ü
|gè¡•f?û:5~fD×o-û:Ánè–|çè‚£è:×%œŽxcVýô-8:QîÔoΔaaBÖôÊǾn:ûRÈY‹sOÂì5"ÓwèܧKú%„}ÔÓÿ+oøgßûpãzïQDzœ˜lk$šà›á4½Ç“é?<=ô¨ˆ ýó¤*gãáš8œùžíš«§¼ÕÛíÝáVGõB
+²ãþB"Vùí•C‘€#º•æÐt^í?EZkóEçô °u5¥Î®úÊá¯ÖL!ò!e¯¶ÿ^¿ þlU>‰Î›Òj’ï›…à!–ÄZˆð™Ãdd!UJâ<4A>%	H汋ÕÅžtÂù;MÌZ¡Hè&”Y[œržÅ{°pJ×ÍYɃe$¡@ºñSÞä“ÿ<>ë,‰°›ôÝXê4ÓÉyâÙZ W¬¹aƱu;xß¾e¶²Zë"<22>hÊL%!ÐÁOìÂ.‰\âac᲌D$µôÅ)ÝØx¼URçÚž„<>èè­8Ùc,ºÍ
+Îò¿îõv—I4¾õïc vŽ¾
býaË<61>þЮҘPº:®³·]•"*RùX£Â¾Kr`(Xª½ålV'¼ÇvÝ÷#PÍ¡·=jÑty’åƒÎ9sT’ŽUGû¿¹q®{Ç®ÈÓÇš‰-˜î"X§
+6ñTÂwx~~.|kLó%Šc 0¸ª(ק_rú›¿[K:Fìg.+¡H çRˆF[í9¦2q¿}Y÷Z?/Ð<Èäiƒ9áhvI•?lÝè¢@x2#ã„SjÖ¬ñ¤ÄÙM´GêG!kVØ40bù'_ˆ¤|®P*7–`Î*?y
•ÄÃOB„#ùÿŸÐeã£SDâ<qúž>jâL<C3A2>
+m½2
ä¨,Ÿ
3S+Þ›ã¹E(ÇC¦ëÂþ<6o‰BÊíB—O'z`¾2wL9ÔagV±ÓÔ÷%<25>2Uflu™Ÿ²¬A%÷{dlÁ¿C\9ŒÌ<>fôù¿ŒÉ­ûáŠw<dÓ©us´·’~…ÈaŽÊT^­”SAs<¦¡×Í+[—\éÁA?b÷}ÑzËAfrŒn¦¤)Ÿ»°Ò/Õæèq±q³1‰Ÿ¹¬zÄNá0êLüü…\TÀ;,«
+A
bçÚbP¯rßìÖÓHcMX ªHYK©Ëÿ‚Ñ~X
屈<RMi¿­þŠfº¹œÞ¹<C39E>D;2Bpu§ U–~µ•MÐ^¯A{SŒ;·S¢æd›ü-gtúÜ\ž/-‰tþõ–ƒ©!ù8·‡üÃ!žðï"*òÜt`b5ÜœRÓž
[]ñæ>½;¦;OŸÆï
+8"£2<C2A3>5H¹^\5ZÒ]+¸Ç?诺î¾4>gå”2ÞNòìàó½'DIcM¡Ô™_èðÈvs|¬’Ø¿·YÆÜ6@F€—¿mÞØ6&R³lâ7\Ù<>cº:D%Âôq_ñMIK–÷ÜÎÎ=—õl|3øºñ
^ÜKN¿0É!@jÕÎã1TóÑË%nú¼§b$ßÔË™Î`®)ÿƒðõ´²Œ£cÜ­pkÒmâ2t&ù!<21>í¢qÍ7m&[u7£×`&	K‘¡%<25>­º“¹{³SÈz#VÙëc7ŸVâk°ºüªWŒøKó¾åH5óeÊàŽ›l?U¤ù<1E>ü}êç}€#„<M‚¨¸Rµ Ô.¾7¶ã©õ^
f8ú<38>
ˆì—×Ó̱VÚæ;«ŽàT*ÚaOú½=˃n›ÐÕjÓô{;ZT{—¾;ÝJùzk¢÷+m)¶£pž¾ÉÞ5­ò`{V[}nC#'KI†«TFõÙSß*xH T>HD‹Ñó^»½b‡¶’>ïpÕ,YÖaÁƒ0I¯>w½-æ’@ü°œ»Jê¿üÊö4®zù
q~)\En[’ú{<7B><>I¦R,²ë¹A)‡Ô³/µnnÉ}ÎŽÂRÏ£"DªÑH2¢$™ÒÞù'9ÊOC5Ò¥X¨xøCÞ“¥(Ó´<11>“È•¿/ÁK8ü.¤FÝ•+CýƒÒÃ7ð„ªkÇõ=«ÓHÐaK<Ï&”`ÃQo/y:éf‘}íS²
ÁÎ3Ðn$¸Z¿¨8×Ĭ2,ˆÂp†i°W‚êZtÜ%'Cgå<67>QÛ¢Lx-Ü	0YiÖírG“UÈò+¥×fî­Æ^ÃQâÕ¯Báke„Œ!ÄÆ
ãYűH§S³šÏ^YX sŒÿHå¡âŠxñ½f¥°?j%ûi‰6êéù©´õgaÁ…¾cö±d²ap¿E°[‚ÄnÞ&>
ÄÛã$_'#mæ¡>KÒàÄÅäÏ~¨ú@LL×â}°£y·™=䉲:ÿ0“Ù¬{›tŠ7ª¡„ÖìñÖ)Z·{]Ú½‰âK¦ÔQÙN±|ì»’…Ön³Þ¹<>0g0ç‰lžÚøv“TL3þ­Ð
]é
 endstream
 endobj
 542 0 obj
@@ -7972,19 +8055,19 @@ endobj
 /Type /ObjStm
 /N 100
 /First 940
-/Length 22123     
+/Length 22133     
 >>
 stream
 497 0 499 20 500 246 501 300 502 991 503 1550 504 1648 505 1680 506 2304 507 2694
 508 3325 509 3941 510 4605 511 5215 512 5722 513 6361 514 6840 515 7192 517 7828 519 8180
 521 8538 523 8882 525 9203 527 9527 529 9750 531 9973 533 10227 535 10791 537 11129 539 11349
-541 11579 543 11920 545 12140 547 12369 549 12782 551 13058 155 13627 108 13766 121 13905 279 14045
-173 14183 278 14322 321 14460 275 14598 120 14736 110 14874 320 15012 277 15148 174 15284 322 15422
-276 15557 171 15692 109 15831 200 15970 111 16109 175 16221 236 16333 308 16445 393 16557 437 16669
-472 16781 552 16893 553 17006 554 17078 555 17143 100 17211 96 17281 92 17363 88 17444 84 17513
-80 17595 76 17677 72 17759 68 17841 64 17910 60 18028 56 18084 52 18202 48 18271 44 18340
-40 18458 36 18527 32 18609 28 18691 24 18773 20 18842 16 18960 12 19041 8 19121 4 19186
-556 19278 557 19446 558 19607 559 19771 560 19999 561 20207 562 20444 563 20642 564 20849 565 21018
+541 11579 543 11930 545 12150 547 12379 549 12792 551 13068 155 13637 108 13776 121 13915 279 14055
+173 14193 278 14332 321 14470 275 14608 120 14746 110 14884 320 15022 277 15158 174 15294 322 15432
+276 15567 171 15702 109 15841 200 15980 111 16119 175 16231 236 16343 308 16455 393 16567 437 16679
+472 16791 552 16903 553 17016 554 17088 555 17153 100 17221 96 17291 92 17373 88 17454 84 17523
+80 17605 76 17687 72 17769 68 17851 64 17920 60 18038 56 18094 52 18212 48 18281 44 18350
+40 18468 36 18537 32 18619 28 18701 24 18783 20 18852 16 18970 12 19051 8 19131 4 19196
+556 19288 557 19456 558 19617 559 19781 560 20009 561 20217 562 20454 563 20652 564 20859 565 21028
 % 497 0 obj
 [892.9]
 % 499 0 obj
@@ -8204,7 +8287,7 @@ stream
 % 541 0 obj
 <<
 /Type /FontDescriptor
-/FontName /WYYAGH+CMSY10
+/FontName /PLWDAY+CMSY10
 /Flags 4
 /FontBBox [-29 -960 1116 775]
 /Ascent 750
@@ -8213,7 +8296,7 @@ stream
 /ItalicAngle -14
 /StemV 40
 /XHeight 431
-/CharSet (/Rfractur/backslash/bar/braceleft/braceright/bullet/element/greaterequal/minus/negationslash/propersubset/propersuperset/union)
+/CharSet (/Rfractur/backslash/bar/braceleft/braceright/bullet/element/greaterequal/lessequal/minus/negationslash/propersubset/propersuperset/union)
 /FontFile 540 0 R
 >>
 % 543 0 obj
@@ -8415,7 +8498,7 @@ stream
 <<
 /Type /Font
 /Subtype /Type1
-/BaseFont /WYYAGH+CMSY10
+/BaseFont /PLWDAY+CMSY10
 /FontDescriptor 541 0 R
 /FirstChar 0
 /LastChar 110
@@ -8806,8 +8889,8 @@ endobj
 585 0 obj
 <<
  /Title (MultiLevel Domain Decomposition Parallel Preconditioners Package based on PSBLAS, V. 2.0) /Subject (MultiLevel Domain Decomposition Parallel Preconditioners Package) /Keywords (Parallel Numerical Software, Algebraic Multilevel Preconditioners, Sparse Iterative Solvers, PSBLAS, MPI) /Creator (pdfLaTeX) /Producer ($Id: userguide.tex 2008-04-08 Pasqua D'Ambra, Daniela di Serafino, Salvatore Filippone$) /Author()/Title()/Subject()/Creator(LaTeX with hyperref package)/Producer(pdfTeX-1.40.15)/Keywords()
-/CreationDate (D:20160229132824Z)
-/ModDate (D:20160229132824Z)
+/CreationDate (D:20160229135950Z)
+/ModDate (D:20160229135950Z)
 /Trapped /False
 /PTEX.Fullbanner (This is pdfTeX, Version 3.14159265-2.6-1.40.15 (TeX Live 2014) kpathsea version 6.2.0)
 >>
@@ -8926,19 +9009,19 @@ endobj
 /W [1 3 1]
 /Root 584 0 R
 /Info 585 0 R
-/ID [<252E47806D66F0DAB10F5905A2354801> <252E47806D66F0DAB10F5905A2354801>]
+/ID [<3F2A25442D85405FC76ACEFD34DF1970> <3F2A25442D85405FC76ACEFD34DF1970>]
 /Length 2935      
 >>
 stream
 ÿ
Ny6
òY
=
òX±
òW±
òV±"
òU	
 ±#
òT±$
òS
±+
òR±6
òQ±7
òP±D
òO±T
òN
5
òM
5'
òL
58
òK
5W
òJ 
5]
òI!"
¤
òH#$
¤'
òG%&
¤,
òF'(
¤2
òE)*
¤7
òD+,
¤=
òC-.
¤A
òB/0
¤E
òA1237
45
ò%
ò4
ò-
ò6:8
×9>;
3<
ò,
ò&A?
ÿ@CDEFGHIJKLMNOPQRSTUVWXYZ]B
Š[\
ò$`^
-í_c±±
±±±±b±a
.x
ò3±
ò(
ò0
ò7
¤^

¬
¤F
¤b
¤Y
¤X
¤W
¤_±
 ±±±
±±±±±	
q_±±±±
ƒ‹±
ò5±±±±±±±± ±%±
ŠY±!
¤P
¤Q
¤[
¤Z
¤N±(±)±'±,±&
¤T±*
¤O
¤\±/±-
¸­±.±1±3±0
Çv±2
ò8±8±4
ÕÚ±5±:±<±9
è±;±>±?±@±A±B±H±E±=
î²±C
¤K
¤L±I±J±K±L±M±N±O±P±Q±R±Z±[±G±X±F

±S
ò+
ò2
ò/
ò)
ò'±U±V±W
¤`
¤a
¤G
¤H
¤I
¤M
¤]±\±]±_±Y

I±^
¤J
¤V±b±c
5
5

5±a
5±`

[
5
ò9
Ñò
5

5
5
5	
5
-
5
5
5

àT
5
ò.
ò*
ò1
5
5
5
5
5
5
5
5
5
	*
5
5
5
5
5
5
5
5 
5!
5"
5#
5$
5%
51
52
53
5.
5
µ
5&
5(
5)
5*
5+
5,
5-
57
¤
¤
54
55
5;
5<
50
59
5/
;{
56
5K
5=
5>
5?
5@
5A
5B
5C
5D
5E
5G
5:
_6
5F
5O
5P
5Q
5I
5L
5H
z<
5J
ò:
5R
5M
‡ç
5N
5U
5T
5X
5S
”
5V
5Z
5[
5^
5Y
«®
5\
5a
5b
5c
¤
¤

5`
¤
5_
¸
¤
9”
¤
¤
¤
¤
¤
ø
¤	
¤
-
¤
¤
¤

¤
¤
¤
=
¤
ò;
¤
¤
&
¤
¤
¤
¤
¤
Mx
¤
¤
¤ 
¤
‚…
¤
¤#
¤$
¤%
¤"
¤(
¤!
¬‡
¤&
¤*
¤-
¤)
¸õ
¤+
¤0
¤/
¤3
¤.
ܹ
¤1
ò<
¤5
¤8
¤4
ã
¤6
¤;
¤:
¤>
¤9
ì,
¤<
¤B
¤?
ót
¤@
¤R
¤C
W
¤D
¤c
¤S
°
¤T
¤U
ò
<1C>
ò

ò
ò
ò
ò
ò
ò
ò
ò	
ò
-
ò
ò
ò

ò
ò
ò
ò
`Z
ò
ªØ
ò
îc
ò
/
ò
R
ò
ŽÇ
ò
®
ò
Îg
ò
ûm
ò
b6
ò
šÁ
ò
¸(
ò
ÖÝ
ò
üZ
ò
(
ò 
6¬
ò!
<07>O
ò"
Ä"
ò#
ò=
ò>
ò?
ò@
òZ
ò[
ò\
ò]
ò^
ò_
ò`
òa
òb
òc7
v47
77777777	7
-777
7777
sd
<08>N
+
5
5
5

àT
5
ò.
ò*
ò1
5
5
5
5
5
5
5
5
5
	*
5
5
5
5
5
5
5
5 
5!
5"
5#
5$
5%
51
52
53
5.
5
µ
5&
5(
5)
5*
5+
5,
5-
57
¤
¤
54
55
5;
5<
50
59
5/
;{
56
5K
5=
5>
5?
5@
5A
5B
5C
5D
5E
5G
5:
_6
5F
5O
5P
5Q
5I
5L
5H
z<
5J
ò:
5R
5M
‡ç
5N
5U
5T
5X
5S
”
5V
5Z
5[
5^
5Y
«®
5\
5a
5b
5c
¤
¤

5`
¤
5_
¸
¤
C$
¤
¤
¤
¤
¤
ø
¤	
¤
+
¤
¤
¤

¤
¤
¤
=
¤
ò;
¤
¤
&H
¤
¤
¤
¤
¤
W
¤
¤
¤ 
¤
Œ
¤
¤#
¤$
¤%
¤"
¤(
¤!
¶
¤&
¤*
¤-
¤)
Â…
¤+
¤0
¤/
¤3
¤.
æI
¤1
ò<
¤5
¤8
¤4
ì¡
¤6
¤;
¤:
¤>
¤9
õ¼
¤<
¤B
¤?
ý
¤@
¤R
¤C
ç
¤D
¤c
¤S
'@
¤T
¤U
ò
&Á
ò

ò
ò
ò
ò
ò
ò
ò
ò	
ò
+
ò
ò
ò

ò
ò
ò
ò
iê
ò
´h
ò
÷ó
ò
9
ò
[—
ò
˜W
ò
·œ
ò
×÷
ò
ý
ò
kÆ
ò
¤Q
ò
Á¸
ò
àm
ò
~
ò
#L
ò 
@Ð
ò!
šs
ò"
ÎF
ò#
ò=
ò>
ò?
ò@
òZ
ò[
ò\
ò]
ò^
ò_
ò`
òa
òb
òc7
€b7
77777777	7
+777
7777
}’
‹|
 endstream
 endobj
 startxref
-557390
+559996
 %%EOF

docs/src/userinterface.tex

@@ -192,11 +192,11 @@ according to their needs.
                          & Basic multi-level framework: additive or multiplicative
                            among the levels (always additive inside a level).         \\ \hline 
 \verb|mld_smoother_type_| \break \verb|SMOOTHER_TYPE|  & \verb|character(len=*)|
-                         & \texttt{'DIAG'} \ \ \ \texttt{'BJAC'} \ \ \ \texttt{'AS'}
+                         & \texttt{'JACOBI'} \ \ \ \texttt{'BJAC'} \ \ \ \texttt{'AS'}
                          & \texttt{'AS'}
                          & Basic predefined one-level preconditioner
-                         (i.e.\ smoother): diagonal, 
-                           block Jacobi, AS. \\ \hline
+                         (i.e.\ smoother): Jacobi (with diagonal local
+                           solver), block Jacobi, AS. \\ \hline
 \verb|mld_smoother_pos_| \break \verb|SMOOTHER_POS| & \verb|character(len=*)|
                          & \texttt{'PRE'} \ \ \ \texttt{'POST'} \ \ \ \texttt{'TWOSIDE'}
                          & \texttt{'TWOSIDE'}
@@ -233,12 +233,16 @@ according to their needs.
                            \texttt{'SUM'} for adding the contributions from the overlap, \texttt{'NONE'}
                            for neglecting them.   \\ \hline
 \verb|mld_sub_solve_| \break \verb|SUB_SOLVE|    & \verb|character(len=*)|
-                         & \texttt{'ILU'} \hspace{2.5cm} \texttt{'MILU'} \hspace{2.5cm} \texttt{'ILUT'} 
+                         & \texttt{'DIAG'} \hspace{2.5cm}
+                           \texttt{'GS'} \hspace{2.5cm}
+                           \texttt{'ILU'} \hspace{2.5cm} 
+                           \texttt{'MILU'} \hspace{2.5cm} \texttt{'ILUT'} 
                            \hspace{2.5cm} \texttt{'UMF'} \hspace{2.5cm} \texttt{'SLU'}
                          & \texttt{'ILU'}
-                         & Predefined local solver: ILU($p$),
-                         MILU($p$), ILU($p,t$), LU from UMFPACK, LU  from SuperLU 
-                           (plus triangular solve). \\ \hline  
+                         & Predefined local solver: diagonal scaling,
+                           Gauss-Seidel, ILU($p$),  MILU($p$),
+                           ILU($p,t$), LU from UMFPACK, LU  from
+                           SuperLU (plus triangular solve). \\ \hline  
 \verb|mld_sub_fillin_| \break \verb|SUB_FILLIN|  & \verb|integer|
                          & Any~int.~num.~$\ge 0$
                          & 0
@@ -253,6 +257,17 @@ according to their needs.
                          & Row and column reordering of the local submatrices: no reordering,
                            reordering according to the global numbering of the rows and columns of
                            the whole matrix. \\
+\verb|mld_solver_sweeps_| \break \verb|SOLVER_SWEEPS|  & \verb|integer|
+                         & Any~int.~num.~$\ge 1$
+                         & 1
+                         & Number of sweeps for iterative local solver
+                           (currently only Gauss-Seidel). \\ \hline
+\verb|mld_solver_eps_| \break \verb|SOLVER_EPS|  & \verb|real|
+                         & Any~real~num.
+                         & 0
+                         & Stopping tolerance for iterative local solver
+                           (currently only Gauss-Seidel); if $\le0$, then
+                           perform prespecified number of iterations. \\ \hline
 \hline
 \end{tabular}
 \end{center}

mlprec/mld_base_prec_type.F90

@@ -379,7 +379,7 @@ contains
       val = mld_slu_
     case('SLUDIST')
       val = mld_sludist_
-    case('DSCALE')
+    case('DIAG')
       val = mld_diag_scale_
     case('ADD')
       val = mld_add_ml_
@@ -409,8 +409,6 @@ contains
       val = mld_twoside_smooth_
     case('NOPREC')
       val = mld_noprec_
-! !$    case('DIAG')
-! !$      val = mld_diag_
     case('BJAC')
       val = mld_bjac_
     case('JAC','JACOBI')

tests/pdegen/Makefile

@@ -11,7 +11,7 @@ FINCLUDES=$(FMFLAG). $(FMFLAG)$(MLDINCDIR)  $(FMFLAG)$(PSBINCDIR) $(FIFLAG).
 
 EXEDIR=./runs
 
-all: spde3d ppde3d spde2d ppde2d
+all: spde3d ppde3d spde2d ppde2d ppde3d_gs
 
 ppde3d: ppde3d.o  data_input.o
 	$(F90LINK) ppde3d.o  data_input.o -o ppde3d $(MLD_LIB) $(PSBLAS_LIB) $(LDLIBS)
@@ -30,8 +30,11 @@ spde2d: spde2d.o   data_input.o
 	$(F90LINK) spde2d.o  data_input.o -o spde2d $(MLD_LIB) $(PSBLAS_LIB) $(LDLIBS)
 	/bin/mv spde2d $(EXEDIR)
 
+ppde3d_gs: ppde3d_gs.o  data_input.o
+	$(F90LINK) ppde3d_gs.o  data_input.o -o ppde3d_gs $(MLD_LIB) $(PSBLAS_LIB) $(LDLIBS)
+	/bin/mv ppde3d_gs $(EXEDIR)
 
-ppde3d.o spde3d.o ppde2d.o spde2d.o: data_input.o 
+ppde3d_gs.o ppde3d.o spde3d.o ppde2d.o spde2d.o: data_input.o 
 
 check: all
 	cd runs && ./ppde2d <ppde.inp && ./spde2d<ppde.inp