mld2p4-2:

docs/html/index.html
 docs/html/node12.html
 docs/html/node14.html
 docs/html/node15.html
 docs/html/node18.html
 docs/html/node26.html
 docs/html/node5.html
 docs/html/node6.html
 docs/html/node8.html
 docs/html/userhtml.html
 docs/src/Makefile
 docs/src/building.tex
 docs/src/userguide.tex
 docs/src/userhtml.tex
 mlprec/impl/mld_c_ilu_solver_impl.f90
 mlprec/impl/mld_d_ilu_solver_impl.f90
 mlprec/impl/mld_s_ilu_solver_impl.f90
 mlprec/impl/mld_z_ilu_solver_impl.f90


Silly bug in ilu%apply( trans='C')

Doc fixes
stopcriterion
Salvatore Filippone 13 years ago
parent e2e8b05e33
commit d632c2ad03

@ -61,9 +61,9 @@ University of Rome ``Tor Vergata'', Italy
<BR>
<BR>
<BR>
Software version: 1.1
Software version: 2.0-beta
<BR>
March 13, 2009
March 20,2012
<BR>

@ -685,7 +685,7 @@ denoted by <IMG
SRC="img67.png"
ALT="$R_l$">.
<DIV ALIGN="CENTER"><A NAME="fig:mlhpost_alg"></A><A NAME="512"></A>
<DIV ALIGN="CENTER"><A NAME="fig:mlhpost_alg"></A><A NAME="509"></A>
<TABLE>
<CAPTION ALIGN="BOTTOM"><STRONG>Figure 1:</STRONG>
Application of the multi-level hybrid post-smoothed preconditioner.</CAPTION>

@ -140,7 +140,7 @@ compilers.
<P>
<BR><P></P>
<DIV ALIGN="CENTER"><A NAME="929"></A>
<DIV ALIGN="CENTER"><A NAME="926"></A>
<TABLE>
<CAPTION><STRONG>Table 1:</STRONG>
Preconditioner types, corresponding strings and default choices.

@ -86,7 +86,7 @@ the corresponding Fortran 95 codes are available in <code>examples/fileread/</co
<P>
<DIV ALIGN="CENTER"><A NAME="fig:ex_default"></A><A NAME="932"></A>
<DIV ALIGN="CENTER"><A NAME="fig:ex_default"></A><A NAME="929"></A>
<TABLE>
<CAPTION ALIGN="BOTTOM"><STRONG>Figure 2:</STRONG>
Setup and application of the default multi-level Schwarz preconditioner.
@ -203,7 +203,7 @@ and linked to the MLD2P4 library.
<BR>
<P>
<DIV ALIGN="CENTER"><A NAME="fig:ex_3lh"></A><A NAME="934"></A>
<DIV ALIGN="CENTER"><A NAME="fig:ex_3lh"></A><A NAME="931"></A>
<TABLE>
<CAPTION ALIGN="BOTTOM"><STRONG>Figure 3:</STRONG>
Setup of a hybrid three-level Schwarz preconditioner.</CAPTION>
@ -235,7 +235,7 @@ Setup of a hybrid three-level Schwarz preconditioner.</CAPTION>
<P>
<DIV ALIGN="CENTER"><A NAME="fig:ex_3la"></A><A NAME="936"></A>
<DIV ALIGN="CENTER"><A NAME="fig:ex_3la"></A><A NAME="933"></A>
<TABLE>
<CAPTION ALIGN="BOTTOM"><STRONG>Figure 4:</STRONG>
Setup of an additive three-level Schwarz preconditioner.</CAPTION>
@ -267,7 +267,7 @@ Setup of an additive three-level Schwarz preconditioner.</CAPTION>
<P>
<DIV ALIGN="CENTER"><A NAME="fig:ex_1l"></A><A NAME="938"></A>
<DIV ALIGN="CENTER"><A NAME="fig:ex_1l"></A><A NAME="935"></A>
<TABLE>
<CAPTION ALIGN="BOTTOM"><STRONG>Figure 5:</STRONG>
Setup of a one-level Schwarz preconditioner.</CAPTION>

@ -133,7 +133,7 @@ refer to Section&nbsp;<A HREF="node11.html#sec:background">4</A>.
<P>
<BR><P></P>
<DIV ALIGN="CENTER"><A NAME="1266"></A>
<DIV ALIGN="CENTER"><A NAME="1263"></A>
<TABLE>
<CAPTION><STRONG>Table 2:</STRONG>
Parameters defining the type of multi-level preconditioner.
@ -177,7 +177,7 @@ Parameters defining the type of multi-level preconditioner.
<P>
<BR><P></P>
<DIV ALIGN="CENTER"><A NAME="1268"></A>
<DIV ALIGN="CENTER"><A NAME="1265"></A>
<TABLE>
<CAPTION><STRONG>Table 3:</STRONG>
Parameters defining the one-level preconditioner used as smoother.
@ -277,7 +277,7 @@ Parameters defining the one-level preconditioner used as smoother.
<P>
<BR><P></P>
<DIV ALIGN="CENTER"><A NAME="1270"></A>
<DIV ALIGN="CENTER"><A NAME="1267"></A>
<TABLE>
<CAPTION><STRONG>Table 4:</STRONG>
Parameters defining the aggregation algorithm.
@ -380,7 +380,7 @@ Parameters defining the aggregation algorithm.
<P>
<BR><P></P>
<DIV ALIGN="CENTER"><A NAME="1273"></A>
<DIV ALIGN="CENTER"><A NAME="1270"></A>
<TABLE>
<CAPTION><STRONG>Table 5:</STRONG>
Parameters defining the coarse-space correction at the coarsest

@ -63,7 +63,7 @@ Mathematics Department, Macquarie University, Sydney.
The command line arguments were: <BR>
<STRONG>latex2html</STRONG> <TT>-local_icons -noaddress -dir ../../html userhtml.tex</TT>
<P>
The translation was initiated by Salvatore Filippone on 2011-03-25
The translation was initiated by Salvatore Filippone on 2012-03-20
<BR><HR>
</BODY>

@ -63,14 +63,13 @@ sources employed to generate the script, but usually this is not needed
to build the software.
<P>
MLD2P4 is implemented almost entirely in Fortran&nbsp;95, with some
MLD2P4 is implemented almost entirely in Fortran&nbsp;2003, with some
interfaces to external libraries in C; the Fortran compiler
must support the Fortran&nbsp;95 standard plus the extension TR15581, which
enhances the usability of <code>ALLOCATABLE</code> variables. Most modern
Fortran compilers support this language level. In particular, this is
supported by the GNU Fortran compiler as of version 4.2.0; however we
recommend to use the latest available release (4.3.1 at the time of
this writing).
must support the Fortran&nbsp;2003 standard plus the extension <code>MOLD=</code>
feature, which
enhances the usability of <code>ALLOCATE</code>.
Many compiles do this; in particular, this is
supported by the GNU Fortran compiler as of version 4.6.0
The software defines data types and interfaces for
real and complex data, in both single and double precision.

@ -80,13 +80,6 @@ The following base libraries are needed:
HREF="node25.html#MPI2">17</A>,<A
HREF="node25.html#MPI1">22</A>] A version of MPI is available on most
high-performance computing systems; only version 1.1 is required.
</DD>
<DT><STRONG>BLACS</STRONG></DT>
<DD>[<A
HREF="node25.html#BLACS">13</A>] The Basic Linear Algebra Communication Subprograms
are available in source form from <code>http://www.netlib.org/blacs</code>;
some vendors include them in their parallel computing
support libraries.
</DD>
<DT><STRONG>PSBLAS</STRONG></DT>
@ -94,17 +87,11 @@ The following base libraries are needed:
HREF="node25.html#PSBLASGUIDE">15</A>,<A
HREF="node25.html#psblas_00">16</A>] Parallel Sparse BLAS is
available from
<BR><code>http://www.ce.uniroma2.it/psblas</code>; version 2.3.1
(or later) is required. Indeed, all the prerequisites
<BR><code>http://www.ce.uniroma2.it/psblas</code>; version
3.0 (or later) is required. Indeed, all the prerequisites
listed so far are also prerequisites of PSBLAS.
To build the MLD2P4 library it is necessary to get access to
the source PSBLAS directory employed to build the version under use; after
the MLD2P4 build process completes, only the compiled form of the
PSBLAS library is necessary to build user applications.
</DD>
</DL>
<P>
Please note that the four previous libraries must have Fortran
interfaces compatible with MLD2P4;
usually this means that they should all be built with the same

@ -62,11 +62,11 @@ in the main directory to generate the necessary makefile(s).
<P>
As a minimal example consider the following:
<PRE>
./configure --with-psblas=/home/user/PSBLAS/psblas-2.3
./configure --with-psblas=PSB-INSTALL-DIR
</PRE>
which assumes that the various MPI compilers and support libraries are
available in the standard directories on the system, and specifies
only the PSBLAS build directory (note that the latter directory must
only the PSBLAS install directory (note that the latter directory must
be specified with an <EM>absolute</EM> path).
The full set of options may be looked at by issuing the command
<code>./configure --help</code>, which produces:
@ -180,13 +180,14 @@ Thus, a sample build with libraries in installation
directories specifics to the GNU 4.3 compiler suite might be as
follows, specifying only the UMFPACK external package:
<PRE>
./configure --with-psblas=/home/user/psblas-2.3/ \
--with-libs="-L/usr/local/BLAS/gnu43 -L/usr/local/BLACS/gnu43" \
--with-blacs=-lmpiblacs --with-umfpackdir=/usr/local/UMFPACK/gnu43
./configure --with-psblas=/home/user/psblas-3.0/ \
--with-libs="-L/usr/local/BLAS/gnu46 -L/usr/local/BLACS/gnu46" \
--with-umfpackdir=/usr/local/UMFPACK/gnu46
</PRE>
Once the configure script has completed execution, it will have
generated the file <code>Make.inc</code> which will then be used by all
Makefiles in the directory tree.
Makefiles in the directory tree; this file will be copied in the
install directory under the name <code>Make.inc.MLD2P4</code>.
<P>
To build the library the user will now enter

@ -61,9 +61,9 @@ University of Rome ``Tor Vergata'', Italy
<BR>
<BR>
<BR>
Software version: 1.1
Software version: 2.0-beta
<BR>
March 13, 2009
March 20,2012
<BR>

@ -139,7 +139,7 @@ PDF = $(join $(BASEFILE),.pdf)
PS = $(join $(BASEFILE),.ps)
GXS = $(join $(BASEFILE),.gxs)
GLX = $(join $(BASEFILE),.glx)
TARGETPDF= ../mld2p4-1.1-guide.pdf
TARGETPDF= ../mld2p4-2.0-guide.pdf
BASEHTML = $(patsubst %.tex,%,$(HTMLFILE))
HTML = $(join $(HTMLFILE),.html)
HTMLDIR = ../html

@ -7,14 +7,13 @@ script. The distribution also includes the autoconf and automake
sources employed to generate the script, but usually this is not needed
to build the software.
MLD2P4 is implemented almost entirely in Fortran~95, with some
MLD2P4 is implemented almost entirely in Fortran~2003, with some
interfaces to external libraries in C; the Fortran compiler
must support the Fortran~95 standard plus the extension TR15581, which
enhances the usability of \verb|ALLOCATABLE| variables. Most modern
Fortran compilers support this language level. In particular, this is
supported by the GNU Fortran compiler as of version 4.2.0; however we
recommend to use the latest available release (4.3.1 at the time of
this writing).
must support the Fortran~2003 standard plus the extension \verb|MOLD=|
feature, which
enhances the usability of \verb|ALLOCATE|.
Many compiles do this; in particular, this is
supported by the GNU Fortran compiler as of version 4.6.0
The software defines data types and interfaces for
real and complex data, in both single and double precision.
@ -36,20 +35,11 @@ The following base libraries are needed:
libraries.
\item[MPI] \cite{MPI2,MPI1} A version of MPI is available on most
high-performance computing systems; only version 1.1 is required.
\item[BLACS] \cite{BLACS} The Basic Linear Algebra Communication Subprograms
are available in source form from \verb|http://www.netlib.org/blacs|;
some vendors include them in their parallel computing
support libraries.
\item[PSBLAS] \cite{PSBLASGUIDE,psblas_00} Parallel Sparse BLAS is
available from \\ \verb|http://www.ce.uniroma2.it/psblas|; version 2.3.1
(or later) is required. Indeed, all the prerequisites
available from \\ \verb|http://www.ce.uniroma2.it/psblas|; version
3.0 (or later) is required. Indeed, all the prerequisites
listed so far are also prerequisites of PSBLAS.
To build the MLD2P4 library it is necessary to get access to
the source PSBLAS directory employed to build the version under use; after
the MLD2P4 build process completes, only the compiled form of the
PSBLAS library is necessary to build user applications.
\end{description}
Please note that the four previous libraries must have Fortran
interfaces compatible with MLD2P4;
usually this means that they should all be built with the same
@ -87,11 +77,11 @@ in the main directory to generate the necessary makefile(s).
As a minimal example consider the following:
\begin{verbatim}
./configure --with-psblas=/home/user/PSBLAS/psblas-2.3
./configure --with-psblas=PSB-INSTALL-DIR
\end{verbatim}
which assumes that the various MPI compilers and support libraries are
available in the standard directories on the system, and specifies
only the PSBLAS build directory (note that the latter directory must
only the PSBLAS install directory (note that the latter directory must
be specified with an {\em absolute} path).
The full set of options may be looked at by issuing the command
\verb|./configure --help|, which produces:
@ -205,13 +195,14 @@ Thus, a sample build with libraries in installation
directories specifics to the GNU 4.3 compiler suite might be as
follows, specifying only the UMFPACK external package:
\begin{verbatim}
./configure --with-psblas=/home/user/psblas-2.3/ \
--with-libs="-L/usr/local/BLAS/gnu43 -L/usr/local/BLACS/gnu43" \
--with-blacs=-lmpiblacs --with-umfpackdir=/usr/local/UMFPACK/gnu43
./configure --with-psblas=/home/user/psblas-3.0/ \
--with-libs="-L/usr/local/BLAS/gnu46 -L/usr/local/BLACS/gnu46" \
--with-umfpackdir=/usr/local/UMFPACK/gnu46
\end{verbatim}
Once the configure script has completed execution, it will have
generated the file \verb|Make.inc| which will then be used by all
Makefiles in the directory tree.
Makefiles in the directory tree; this file will be copied in the
install directory under the name \verb|Make.inc.MLD2P4|.
To build the library the user will now enter
\begin{verbatim}

@ -27,7 +27,7 @@
\pdfcompresslevel=0 %-- 0 = none, 9 = best
\pdfinfo{ %-- Info dictionary of PDF output /Author (PD, DdS, SF)
/Title (MultiLevel Domain Decomposition Parallel Preconditioners Package
based on PSBLAS, V. 1.1)
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)
@ -125,9 +125,9 @@ based on PSBLAS}
\vspace{\stretch{1}}
\noindent\hspace*{\centeroffset}\makebox[0pt][l]{\begin{minipage}{\textwidth}
\flushright
\large Software version: 1.1\\
\large Software version: 2.0-beta\\
%\today
\large March 13, 2009
\large March 20, 2012
\end{minipage}}
%\addtolength{\textwidth}{\centeroffset}
\vspace{\stretch{2}}

@ -102,9 +102,9 @@ based on PSBLAS}\\[3ex]
University of Rome ``Tor Vergata'', Italy\\[2ex]
%\\[10ex]
%\today
Software version: 1.1\\
Software version: 2.0-beta\\
%\today
March 13, 2009
March 20,2012
\clearpage
\ \\
\thispagestyle{empty}

@ -114,8 +114,6 @@ subroutine mld_c_ilu_solver_apply_vect(alpha,sv,x,beta,y,desc_data,trans,work,in
allocate(ww(n_col),aux(4*n_col),stat=info)
endif
call wv%bld(n_col,mold=x%v)
if (info /= psb_success_) then
info=psb_err_alloc_request_
call psb_errpush(info,name,i_err=(/5*n_col,0,0,0,0/),&
@ -147,7 +145,7 @@ subroutine mld_c_ilu_solver_apply_vect(alpha,sv,x,beta,y,desc_data,trans,work,in
call wv1%mlt(cone,sv%dv,wv,czero,info,conjgx=trans_)
if (info == psb_success_) call psb_spsm(alpha,sv%l,wv,beta,y,desc_data,info,&
if (info == psb_success_) call psb_spsm(alpha,sv%l,wv1,beta,y,desc_data,info,&
& trans=trans_,scale='U',choice=psb_none_,work=aux)
case default

@ -114,8 +114,6 @@ subroutine mld_d_ilu_solver_apply_vect(alpha,sv,x,beta,y,desc_data,trans,work,in
allocate(ww(n_col),aux(4*n_col),stat=info)
endif
call wv%bld(n_col,mold=x%v)
if (info /= psb_success_) then
info=psb_err_alloc_request_
call psb_errpush(info,name,i_err=(/5*n_col,0,0,0,0/),&
@ -147,7 +145,7 @@ subroutine mld_d_ilu_solver_apply_vect(alpha,sv,x,beta,y,desc_data,trans,work,in
call wv1%mlt(done,sv%dv,wv,dzero,info,conjgx=trans_)
if (info == psb_success_) call psb_spsm(alpha,sv%l,wv,beta,y,desc_data,info,&
if (info == psb_success_) call psb_spsm(alpha,sv%l,wv1,beta,y,desc_data,info,&
& trans=trans_,scale='U',choice=psb_none_,work=aux)
case default

@ -114,8 +114,6 @@ subroutine mld_s_ilu_solver_apply_vect(alpha,sv,x,beta,y,desc_data,trans,work,in
allocate(ww(n_col),aux(4*n_col),stat=info)
endif
call wv%bld(n_col,mold=x%v)
if (info /= psb_success_) then
info=psb_err_alloc_request_
call psb_errpush(info,name,i_err=(/5*n_col,0,0,0,0/),&
@ -147,7 +145,7 @@ subroutine mld_s_ilu_solver_apply_vect(alpha,sv,x,beta,y,desc_data,trans,work,in
call wv1%mlt(sone,sv%dv,wv,szero,info,conjgx=trans_)
if (info == psb_success_) call psb_spsm(alpha,sv%l,wv,beta,y,desc_data,info,&
if (info == psb_success_) call psb_spsm(alpha,sv%l,wv1,beta,y,desc_data,info,&
& trans=trans_,scale='U',choice=psb_none_,work=aux)
case default

@ -114,8 +114,6 @@ subroutine mld_z_ilu_solver_apply_vect(alpha,sv,x,beta,y,desc_data,trans,work,in
allocate(ww(n_col),aux(4*n_col),stat=info)
endif
call wv%bld(n_col,mold=x%v)
if (info /= psb_success_) then
info=psb_err_alloc_request_
call psb_errpush(info,name,i_err=(/5*n_col,0,0,0,0/),&
@ -147,7 +145,7 @@ subroutine mld_z_ilu_solver_apply_vect(alpha,sv,x,beta,y,desc_data,trans,work,in
call wv1%mlt(zone,sv%dv,wv,zzero,info,conjgx=trans_)
if (info == psb_success_) call psb_spsm(alpha,sv%l,wv,beta,y,desc_data,info,&
if (info == psb_success_) call psb_spsm(alpha,sv%l,wv1,beta,y,desc_data,info,&
& trans=trans_,scale='U',choice=psb_none_,work=aux)
case default

Loading…
Cancel
Save