config/pac.m4
 configure
 docs/html/WARNINGS
 docs/html/images.log
 docs/html/images.pl
 docs/html/images.tex
 docs/html/img100.png
 docs/html/img101.png
 docs/html/img92.png
 docs/html/img93.png
 docs/html/img94.png
 docs/html/img95.png
 docs/html/img96.png
 docs/html/img97.png
 docs/html/img98.png
 docs/html/img99.png
 docs/html/index.html
 docs/html/internals.pl
 docs/html/labels.pl
 docs/html/node1.html
 docs/html/node10.html
 docs/html/node11.html
 docs/html/node12.html
 docs/html/node13.html
 docs/html/node14.html
 docs/html/node15.html
 docs/html/node16.html
 docs/html/node17.html
 docs/html/node18.html
 docs/html/node19.html
 docs/html/node2.html
 docs/html/node20.html
 docs/html/node21.html
 docs/html/node22.html
 docs/html/node23.html
 docs/html/node24.html
 docs/html/node25.html
 docs/html/node26.html
 docs/html/node27.html
 docs/html/node28.html
 docs/html/node29.html
 docs/html/node3.html
 docs/html/node30.html
 docs/html/node31.html
 docs/html/node4.html
 docs/html/node5.html
 docs/html/node6.html
 docs/html/node7.html
 docs/html/node8.html
 docs/html/node9.html
 docs/html/userhtml.html
 docs/pdf/Makefile
 docs/pdf/abstract.tex
 docs/pdf/background.tex
 docs/pdf/bibliography.tex
 docs/pdf/building.tex
 docs/pdf/conventions.tex
 docs/pdf/distribution.tex
 docs/pdf/errors.tex
 docs/pdf/gettingstarted.tex
 docs/pdf/overview.tex
 docs/pdf/title.tex
 docs/pdf/userguide.tex
 docs/pdf/userhtml.tex
 docs/pdf/userinterface.tex

Configure minro fix: require psblas 2.3
Doc fixes: dual version of title for pdf/html, fixed tables.
This commit is contained in:
Salvatore Filippone
2008-07-23 08:51:18 +00:00
parent b7445f740c
commit 3f302b360c
65 changed files with 2450 additions and 2042 deletions
+92 -87
View File
@@ -1,4 +1,4 @@
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2 Final//EN">
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<!--Converted with LaTeX2HTML 2002-2-1 (1.71)
original version by: Nikos Drakos, CBLU, University of Leeds
@@ -24,34 +24,35 @@ original version by: Nikos Drakos, CBLU, University of Leeds
</HEAD>
<BODY >
<!--Navigation Panel-->
<A NAME="tex2html211"
<DIV CLASS="navigation"><!--Navigation Panel-->
<A NAME="tex2html210"
HREF="node13.html">
<IMG WIDTH="37" HEIGHT="24" ALIGN="BOTTOM" BORDER="0" ALT="next"
SRC="file:/usr/share/latex2html/icons/next.png"></A>
<A NAME="tex2html207"
<A NAME="tex2html206"
HREF="node10.html">
<IMG WIDTH="26" HEIGHT="24" ALIGN="BOTTOM" BORDER="0" ALT="up"
SRC="file:/usr/share/latex2html/icons/up.png"></A>
<A NAME="tex2html203"
<A NAME="tex2html202"
HREF="node11.html">
<IMG WIDTH="63" HEIGHT="24" ALIGN="BOTTOM" BORDER="0" ALT="previous"
SRC="file:/usr/share/latex2html/icons/prev.png"></A>
<A NAME="tex2html209"
HREF="node1.html">
<A NAME="tex2html208"
HREF="node2.html">
<IMG WIDTH="65" HEIGHT="24" ALIGN="BOTTOM" BORDER="0" ALT="contents"
SRC="file:/usr/share/latex2html/icons/contents.png"></A>
<BR>
<B> Next:</B> <A NAME="tex2html212"
<B> Next:</B> <A NAME="tex2html211"
HREF="node13.html">Getting Started</A>
<B> Up:</B> <A NAME="tex2html208"
<B> Up:</B> <A NAME="tex2html207"
HREF="node10.html">Multi-level Domain Decomposition Background</A>
<B> Previous:</B> <A NAME="tex2html204"
<B> Previous:</B> <A NAME="tex2html203"
HREF="node11.html">Multi-level Schwarz Preconditioners</A>
&nbsp; <B> <A NAME="tex2html210"
HREF="node1.html">Contents</A></B>
<BR>
&nbsp; <B> <A NAME="tex2html209"
HREF="node2.html">Contents</A></B>
<BR>
<BR></DIV>
<!--End of Navigation Panel-->
<H2><A NAME="SECTION00062000000000000000"></A><A NAME="sec:aggregation"></A>
@@ -60,29 +61,29 @@ Smoothed Aggregation
</H2>
<P>
In order to define the restriction operator <IMG
In order to define the restriction operator <SPAN CLASS="MATH"><IMG
WIDTH="29" HEIGHT="32" ALIGN="MIDDLE" BORDER="0"
SRC="img67.png"
ALT="$R_C$">, which is used to compute
the coarse-level matrix <IMG
ALT="$R_C$"></SPAN>, which is used to compute
the coarse-level matrix <SPAN CLASS="MATH"><IMG
WIDTH="29" HEIGHT="32" ALIGN="MIDDLE" BORDER="0"
SRC="img43.png"
ALT="$A_C$">, MLD2P4 uses the <I>smoothed aggregation</I>
ALT="$A_C$"></SPAN>, MLD2P4 uses the <SPAN CLASS="textit">smoothed aggregation</SPAN>
algorithm described in [<A
HREF="node30.html#BREZINA_VANEK">1</A>,<A
HREF="node30.html#VANEK_MANDEL_BREZINA">19</A>].
HREF="node30.html#VANEK_MANDEL_BREZINA">24</A>].
The basic idea of this algorithm is to build a coarse set of vertices
<IMG
<SPAN CLASS="MATH"><IMG
WIDTH="32" HEIGHT="32" ALIGN="MIDDLE" BORDER="0"
SRC="img44.png"
ALT="$W_C$"> by suitably grouping the vertices of <IMG
ALT="$W_C$"></SPAN> by suitably grouping the vertices of <SPAN CLASS="MATH"><IMG
WIDTH="24" HEIGHT="15" ALIGN="BOTTOM" BORDER="0"
SRC="img10.png"
ALT="$W$"> into disjoint subsets
(aggregates), and to define the coarse-to-fine space transfer operator <IMG
ALT="$W$"></SPAN> into disjoint subsets
(aggregates), and to define the coarse-to-fine space transfer operator <SPAN CLASS="MATH"><IMG
WIDTH="29" HEIGHT="40" ALIGN="MIDDLE" BORDER="0"
SRC="img68.png"
ALT="$R_C^T$"> by
ALT="$R_C^T$"></SPAN> by
applying a suitable smoother to a simple piecewise constant
prolongation operator, to improve the quality of the coarse-space correction.
@@ -90,29 +91,29 @@ prolongation operator, to improve the quality of the coarse-space correction.
Three main steps can be identified in the smoothed aggregation procedure:
<OL>
<LI>coarsening of the vertex set <IMG
<LI>coarsening of the vertex set <SPAN CLASS="MATH"><IMG
WIDTH="24" HEIGHT="15" ALIGN="BOTTOM" BORDER="0"
SRC="img10.png"
ALT="$W$">, to obtain <IMG
ALT="$W$"></SPAN>, to obtain <SPAN CLASS="MATH"><IMG
WIDTH="32" HEIGHT="32" ALIGN="MIDDLE" BORDER="0"
SRC="img44.png"
ALT="$W_C$">;
ALT="$W_C$"></SPAN>;
</LI>
<LI>construction of the prolongator <IMG
<LI>construction of the prolongator <SPAN CLASS="MATH"><IMG
WIDTH="29" HEIGHT="40" ALIGN="MIDDLE" BORDER="0"
SRC="img68.png"
ALT="$R_C^T$">;
ALT="$R_C^T$"></SPAN>;
</LI>
<LI>application of <IMG
<LI>application of <SPAN CLASS="MATH"><IMG
WIDTH="29" HEIGHT="32" ALIGN="MIDDLE" BORDER="0"
SRC="img67.png"
ALT="$R_C$"> and <IMG
ALT="$R_C$"></SPAN> and <SPAN CLASS="MATH"><IMG
WIDTH="29" HEIGHT="40" ALIGN="MIDDLE" BORDER="0"
SRC="img68.png"
ALT="$R_C^T$"> to build <IMG
ALT="$R_C^T$"></SPAN> to build <SPAN CLASS="MATH"><IMG
WIDTH="29" HEIGHT="32" ALIGN="MIDDLE" BORDER="0"
SRC="img43.png"
ALT="$A_C$">.
ALT="$A_C$"></SPAN>.
</LI>
</OL>
@@ -120,20 +121,20 @@ Three main steps can be identified in the smoothed aggregation procedure:
To perform the coarsening step, we have implemented the aggregation algorithm sketched
in [<A
HREF="node30.html#apnum_07">4</A>]. According to [<A
HREF="node30.html#VANEK_MANDEL_BREZINA">19</A>], a modification of
HREF="node30.html#VANEK_MANDEL_BREZINA">24</A>], a modification of
this algorithm has been actually considered,
in which each aggregate <IMG
in which each aggregate <SPAN CLASS="MATH"><IMG
WIDTH="26" HEIGHT="32" ALIGN="MIDDLE" BORDER="0"
SRC="img69.png"
ALT="$N_r$"> is made of vertices of <IMG
ALT="$N_r$"></SPAN> is made of vertices of <SPAN CLASS="MATH"><IMG
WIDTH="24" HEIGHT="15" ALIGN="BOTTOM" BORDER="0"
SRC="img10.png"
ALT="$W$"> that are <I>strongly coupled</I>
to a certain root vertex <IMG
ALT="$W$"></SPAN> that are <SPAN CLASS="textit">strongly coupled</SPAN>
to a certain root vertex <SPAN CLASS="MATH"><IMG
WIDTH="53" HEIGHT="32" ALIGN="MIDDLE" BORDER="0"
SRC="img70.png"
ALT="$r \in W$">, i.e. <BR><P></P>
<DIV ALIGN="CENTER">
ALT="$r \in W$"></SPAN>, i.e. <BR><P></P>
<DIV ALIGN="CENTER" CLASS="mathdisplay">
<!-- MATH
\begin{displaymath}
N_r = \left\{s \in W: |a_{rs}| > \theta \sqrt{|a_{rr}a_{ss}|} \right\}
@@ -153,47 +154,47 @@ N_r = \left\{s \in W: |a_{rs}| > \theta \sqrt{|a_{rr}a_{ss}|} \right\}
for a given <!-- MATH
$\theta \in [0,1]$
-->
<IMG
<SPAN CLASS="MATH"><IMG
WIDTH="69" HEIGHT="36" ALIGN="MIDDLE" BORDER="0"
SRC="img72.png"
ALT="$\theta \in [0,1]$">.
Since this algorithm has a sequential nature, a <I>decoupled</I> version of
it has been chosen, where each processor <IMG
ALT="$\theta \in [0,1]$"></SPAN>.
Since this algorithm has a sequential nature, a <SPAN CLASS="textit">decoupled</SPAN> version of
it has been chosen, where each processor <SPAN CLASS="MATH"><IMG
WIDTH="10" HEIGHT="18" ALIGN="BOTTOM" BORDER="0"
SRC="img73.png"
ALT="$i$"> independently applies the algorithm to
the set of vertices <IMG
ALT="$i$"></SPAN> independently applies the algorithm to
the set of vertices <SPAN CLASS="MATH"><IMG
WIDTH="31" HEIGHT="39" ALIGN="MIDDLE" BORDER="0"
SRC="img74.png"
ALT="$W_i^0$"> assigned to it in the initial data distribution. This
ALT="$W_i^0$"></SPAN> assigned to it in the initial data distribution. This
version is embarrassingly parallel, since it does not require any data communication.
On the other hand, it may produce non-uniform aggregates near boundary vertices,
i.e. near vertices adjacent to vertices in other processors, and is strongly
dependent on the number of processors and on the initial partitioning of the matrix <IMG
dependent on the number of processors and on the initial partitioning of the matrix <SPAN CLASS="MATH"><IMG
WIDTH="18" HEIGHT="15" ALIGN="BOTTOM" BORDER="0"
SRC="img2.png"
ALT="$A$">.
ALT="$A$"></SPAN>.
Nevertheless, this algorithm has been chosen for the implementation in MLD2P4,
since it has been shown to produce good results in practice
[<A
HREF="node30.html#aaecc_07">3</A>,<A
HREF="node30.html#apnum_07">4</A>,<A
HREF="node30.html#TUMINARO_TONG">18</A>].
HREF="node30.html#TUMINARO_TONG">23</A>].
<P>
The prolongator <IMG
The prolongator <SPAN CLASS="MATH"><IMG
WIDTH="75" HEIGHT="40" ALIGN="MIDDLE" BORDER="0"
SRC="img75.png"
ALT="$P_C=R_C^T$"> is built starting from a <I>tentative prolongator</I>
ALT="$P_C=R_C^T$"></SPAN> is built starting from a <SPAN CLASS="textit">tentative prolongator</SPAN>
<!-- MATH
$P \in \Re^{n \times n_C}$
-->
<IMG
<SPAN CLASS="MATH"><IMG
WIDTH="90" HEIGHT="38" ALIGN="MIDDLE" BORDER="0"
SRC="img76.png"
ALT="$P \in \Re^{n \times n_C}$">, defined as
ALT="$P \in \Re^{n \times n_C}$"></SPAN>, defined as
<BR>
<DIV ALIGN="RIGHT">
<DIV ALIGN="RIGHT" CLASS="mathdisplay">
<!-- MATH
\begin{equation}
@@ -204,6 +205,7 @@ P=(p_{ij}), \quad p_{ij}=
\end{array} \right. .
\end{equation}
-->
<A NAME="eq:tent_prol"></A>
<TABLE WIDTH="100%" ALIGN="CENTER">
<TR VALIGN="MIDDLE"><TD ALIGN="CENTER" NOWRAP><A NAME="eq:tent_prol"></A><IMG
WIDTH="290" HEIGHT="52" BORDER="0"
@@ -216,32 +218,33 @@ P=(p_{ij}), \quad p_{ij}=
0 &amp; \quad \mbox{otherwise}
\end{array} \right. .
\end{displaymath}"></TD>
<TD WIDTH=10 ALIGN="RIGHT">
(2)</TD></TR>
<TD CLASS="eqno" WIDTH=10 ALIGN="RIGHT">
(<SPAN CLASS="arabic">2</SPAN>)</TD></TR>
</TABLE>
<BR CLEAR="ALL"></DIV><P></P>
<IMG
<SPAN CLASS="MATH"><IMG
WIDTH="27" HEIGHT="32" ALIGN="MIDDLE" BORDER="0"
SRC="img78.png"
ALT="$P_C$"> is obtained by
applying to <IMG
ALT="$P_C$"></SPAN> is obtained by
applying to <SPAN CLASS="MATH"><IMG
WIDTH="18" HEIGHT="15" ALIGN="BOTTOM" BORDER="0"
SRC="img79.png"
ALT="$P$"> a smoother <!-- MATH
ALT="$P$"></SPAN> a smoother <!-- MATH
$S \in \Re^{n \times n}$
-->
<IMG
<SPAN CLASS="MATH"><IMG
WIDTH="78" HEIGHT="38" ALIGN="MIDDLE" BORDER="0"
SRC="img80.png"
ALT="$S \in \Re^{n \times n}$">:
ALT="$S \in \Re^{n \times n}$"></SPAN>:
<BR>
<DIV ALIGN="RIGHT">
<DIV ALIGN="RIGHT" CLASS="mathdisplay">
<!-- MATH
\begin{equation}
P_C = S P,
\end{equation}
-->
<A NAME="eq:smoothed_prol"></A>
<TABLE WIDTH="100%" ALIGN="CENTER">
<TR VALIGN="MIDDLE"><TD ALIGN="CENTER" NOWRAP><A NAME="eq:smoothed_prol"></A><IMG
WIDTH="73" HEIGHT="30" BORDER="0"
@@ -249,27 +252,28 @@ P_C = S P,
ALT="\begin{displaymath}
P_C = S P,
\end{displaymath}"></TD>
<TD WIDTH=10 ALIGN="RIGHT">
(3)</TD></TR>
<TD CLASS="eqno" WIDTH=10 ALIGN="RIGHT">
(<SPAN CLASS="arabic">3</SPAN>)</TD></TR>
</TABLE>
<BR CLEAR="ALL"></DIV><P></P>
in order to remove oscillatory components from the range of the prolongator
and hence to improve the convergence properties of the multi-level
Schwarz method [<A
HREF="node30.html#BREZINA_VANEK">1</A>,<A
HREF="node30.html#StubenGMD69_99">17</A>].
A simple choice for <IMG
HREF="node30.html#StubenGMD69_99">22</A>].
A simple choice for <SPAN CLASS="MATH"><IMG
WIDTH="16" HEIGHT="15" ALIGN="BOTTOM" BORDER="0"
SRC="img82.png"
ALT="$S$"> is the damped Jacobi smoother:
ALT="$S$"></SPAN> is the damped Jacobi smoother:
<BR>
<DIV ALIGN="RIGHT">
<DIV ALIGN="RIGHT" CLASS="mathdisplay">
<!-- MATH
\begin{equation}
S = I - \omega D^{-1} A ,
\end{equation}
-->
<A NAME="eq:jac_smoother"></A>
<TABLE WIDTH="100%" ALIGN="CENTER">
<TR VALIGN="MIDDLE"><TD ALIGN="CENTER" NOWRAP><A NAME="eq:jac_smoother"></A><IMG
WIDTH="126" HEIGHT="30" BORDER="0"
@@ -277,52 +281,53 @@ S = I - \omega D^{-1} A ,
ALT="\begin{displaymath}
S = I - \omega D^{-1} A ,
\end{displaymath}"></TD>
<TD WIDTH=10 ALIGN="RIGHT">
(4)</TD></TR>
<TD CLASS="eqno" WIDTH=10 ALIGN="RIGHT">
(<SPAN CLASS="arabic">4</SPAN>)</TD></TR>
</TABLE>
<BR CLEAR="ALL"></DIV><P></P>
where the value of <IMG
where the value of <SPAN CLASS="MATH"><IMG
WIDTH="16" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
SRC="img84.png"
ALT="$\omega$"> can be chosen
using some estimate of the spectral radius of <IMG
ALT="$\omega$"></SPAN> can be chosen
using some estimate of the spectral radius of <SPAN CLASS="MATH"><IMG
WIDTH="50" HEIGHT="21" ALIGN="BOTTOM" BORDER="0"
SRC="img85.png"
ALT="$D^{-1}A$"> [<A
ALT="$D^{-1}A$"></SPAN> [<A
HREF="node30.html#BREZINA_VANEK">1</A>].
<P>
<HR>
<DIV CLASS="navigation"><HR>
<!--Navigation Panel-->
<A NAME="tex2html211"
<A NAME="tex2html210"
HREF="node13.html">
<IMG WIDTH="37" HEIGHT="24" ALIGN="BOTTOM" BORDER="0" ALT="next"
SRC="file:/usr/share/latex2html/icons/next.png"></A>
<A NAME="tex2html207"
<A NAME="tex2html206"
HREF="node10.html">
<IMG WIDTH="26" HEIGHT="24" ALIGN="BOTTOM" BORDER="0" ALT="up"
SRC="file:/usr/share/latex2html/icons/up.png"></A>
<A NAME="tex2html203"
<A NAME="tex2html202"
HREF="node11.html">
<IMG WIDTH="63" HEIGHT="24" ALIGN="BOTTOM" BORDER="0" ALT="previous"
SRC="file:/usr/share/latex2html/icons/prev.png"></A>
<A NAME="tex2html209"
HREF="node1.html">
<A NAME="tex2html208"
HREF="node2.html">
<IMG WIDTH="65" HEIGHT="24" ALIGN="BOTTOM" BORDER="0" ALT="contents"
SRC="file:/usr/share/latex2html/icons/contents.png"></A>
<BR>
<B> Next:</B> <A NAME="tex2html212"
<B> Next:</B> <A NAME="tex2html211"
HREF="node13.html">Getting Started</A>
<B> Up:</B> <A NAME="tex2html208"
<B> Up:</B> <A NAME="tex2html207"
HREF="node10.html">Multi-level Domain Decomposition Background</A>
<B> Previous:</B> <A NAME="tex2html204"
<B> Previous:</B> <A NAME="tex2html203"
HREF="node11.html">Multi-level Schwarz Preconditioners</A>
&nbsp; <B> <A NAME="tex2html210"
HREF="node1.html">Contents</A></B>
&nbsp; <B> <A NAME="tex2html209"
HREF="node2.html">Contents</A></B> </DIV>
<!--End of Navigation Panel-->
<ADDRESS>
Salvatore Filippone
2008-07-22
2008-07-23
</ADDRESS>
</BODY>
</HTML>