New default value for itrace=-1, new behaviour for printouts. Docs.

new-parstruct
Salvatore Filippone 6 years ago
parent aee1aaa72e
commit d728186bcf

Binary file not shown.

Before

Width:  |  Height:  |  Size: 218 B

After

Width:  |  Height:  |  Size: 228 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 227 B

After

Width:  |  Height:  |  Size: 259 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 184 B

After

Width:  |  Height:  |  Size: 194 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 253 B

After

Width:  |  Height:  |  Size: 257 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 252 B

After

Width:  |  Height:  |  Size: 276 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 2.9 KiB

After

Width:  |  Height:  |  Size: 3.1 KiB

Binary file not shown.

Before

Width:  |  Height:  |  Size: 229 B

After

Width:  |  Height:  |  Size: 232 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 268 B

After

Width:  |  Height:  |  Size: 312 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 244 B

After

Width:  |  Height:  |  Size: 277 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 207 B

After

Width:  |  Height:  |  Size: 215 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 522 B

After

Width:  |  Height:  |  Size: 583 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 664 B

After

Width:  |  Height:  |  Size: 732 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 257 B

After

Width:  |  Height:  |  Size: 268 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 233 B

After

Width:  |  Height:  |  Size: 240 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 223 B

After

Width:  |  Height:  |  Size: 230 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 707 B

After

Width:  |  Height:  |  Size: 758 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 846 B

After

Width:  |  Height:  |  Size: 867 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 1.0 KiB

After

Width:  |  Height:  |  Size: 1.1 KiB

Binary file not shown.

Before

Width:  |  Height:  |  Size: 1.2 KiB

After

Width:  |  Height:  |  Size: 1.3 KiB

Binary file not shown.

Before

Width:  |  Height:  |  Size: 918 B

After

Width:  |  Height:  |  Size: 1.0 KiB

Binary file not shown.

Before

Width:  |  Height:  |  Size: 987 B

After

Width:  |  Height:  |  Size: 1.1 KiB

Binary file not shown.

Before

Width:  |  Height:  |  Size: 998 B

After

Width:  |  Height:  |  Size: 1.2 KiB

Binary file not shown.

Before

Width:  |  Height:  |  Size: 1012 B

After

Width:  |  Height:  |  Size: 1.1 KiB

Binary file not shown.

Before

Width:  |  Height:  |  Size: 187 B

After

Width:  |  Height:  |  Size: 196 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 797 B

After

Width:  |  Height:  |  Size: 915 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 601 B

After

Width:  |  Height:  |  Size: 659 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 436 B

After

Width:  |  Height:  |  Size: 449 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 444 B

After

Width:  |  Height:  |  Size: 478 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 393 B

After

Width:  |  Height:  |  Size: 410 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 335 B

After

Width:  |  Height:  |  Size: 354 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 353 B

After

Width:  |  Height:  |  Size: 208 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 291 B

After

Width:  |  Height:  |  Size: 395 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 178 B

After

Width:  |  Height:  |  Size: 184 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 185 B

After

Width:  |  Height:  |  Size: 201 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 436 B

After

Width:  |  Height:  |  Size: 482 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 258 B

After

Width:  |  Height:  |  Size: 267 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 2.4 KiB

After

Width:  |  Height:  |  Size: 4.1 KiB

Binary file not shown.

Before

Width:  |  Height:  |  Size: 239 B

After

Width:  |  Height:  |  Size: 245 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 894 B

After

Width:  |  Height:  |  Size: 1.1 KiB

Binary file not shown.

Before

Width:  |  Height:  |  Size: 292 B

After

Width:  |  Height:  |  Size: 311 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 2.2 KiB

After

Width:  |  Height:  |  Size: 3.9 KiB

Binary file not shown.

Before

Width:  |  Height:  |  Size: 713 B

After

Width:  |  Height:  |  Size: 799 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 308 B

After

Width:  |  Height:  |  Size: 311 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 2.2 KiB

After

Width:  |  Height:  |  Size: 3.9 KiB

Binary file not shown.

Before

Width:  |  Height:  |  Size: 783 B

After

Width:  |  Height:  |  Size: 909 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 315 B

After

Width:  |  Height:  |  Size: 320 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 2.2 KiB

After

Width:  |  Height:  |  Size: 3.9 KiB

Binary file not shown.

Before

Width:  |  Height:  |  Size: 575 B

After

Width:  |  Height:  |  Size: 655 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 416 B

After

Width:  |  Height:  |  Size: 476 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 459 B

After

Width:  |  Height:  |  Size: 498 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 549 B

After

Width:  |  Height:  |  Size: 572 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 389 B

After

Width:  |  Height:  |  Size: 415 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 412 B

After

Width:  |  Height:  |  Size: 421 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 1.3 KiB

After

Width:  |  Height:  |  Size: 1.9 KiB

Binary file not shown.

Before

Width:  |  Height:  |  Size: 2.3 KiB

After

Width:  |  Height:  |  Size: 3.7 KiB

Binary file not shown.

Before

Width:  |  Height:  |  Size: 2.2 KiB

After

Width:  |  Height:  |  Size: 3.0 KiB

Binary file not shown.

Before

Width:  |  Height:  |  Size: 1.6 KiB

After

Width:  |  Height:  |  Size: 2.3 KiB

Binary file not shown.

Before

Width:  |  Height:  |  Size: 4.7 KiB

After

Width:  |  Height:  |  Size: 5.0 KiB

Binary file not shown.

Before

Width:  |  Height:  |  Size: 4.6 KiB

After

Width:  |  Height:  |  Size: 5.3 KiB

Binary file not shown.

Before

Width:  |  Height:  |  Size: 333 B

After

Width:  |  Height:  |  Size: 366 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 1.2 KiB

After

Width:  |  Height:  |  Size: 2.4 KiB

Binary file not shown.

Before

Width:  |  Height:  |  Size: 222 B

After

Width:  |  Height:  |  Size: 231 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 714 B

After

Width:  |  Height:  |  Size: 800 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 298 B

After

Width:  |  Height:  |  Size: 369 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 1.2 KiB

After

Width:  |  Height:  |  Size: 1.4 KiB

Binary file not shown.

Before

Width:  |  Height:  |  Size: 180 B

After

Width:  |  Height:  |  Size: 186 B

Binary file not shown.

Before

Width:  |  Height:  |  Size: 211 B

After

Width:  |  Height:  |  Size: 219 B

@ -59,7 +59,7 @@ call psb_get_rank(rank, icontxt, id)
<P>
This subroutine returns the MPI rank of the PSBLAS process <SPAN CLASS="MATH"><IMG
WIDTH="18" HEIGHT="15" ALIGN="BOTTOM" BORDER="0"
WIDTH="18" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
SRC="img133.png"
ALT="$id$"></SPAN>
<DL>
@ -106,7 +106,7 @@ Specified as: an integer value. <!-- MATH
</DD>
<DT><STRONG>rank</STRONG></DT>
<DD>The MPI rank associated with the PSBLAS process <SPAN CLASS="MATH"><IMG
WIDTH="18" HEIGHT="15" ALIGN="BOTTOM" BORDER="0"
WIDTH="18" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
SRC="img133.png"
ALT="$id$"></SPAN>.
<BR>

@ -143,17 +143,17 @@ Kind, rank and size must agree on all processes.
(local) NRM2 operations at the same time.
</LI>
<LI>Denoting by <SPAN CLASS="MATH"><IMG
WIDTH="31" HEIGHT="30" ALIGN="MIDDLE" BORDER="0"
WIDTH="32" HEIGHT="30" ALIGN="MIDDLE" BORDER="0"
SRC="img137.png"
ALT="$dat_i$"></SPAN> the value of the variable <SPAN CLASS="MATH"><IMG
WIDTH="26" HEIGHT="15" ALIGN="BOTTOM" BORDER="0"
WIDTH="26" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
SRC="img138.png"
ALT="$dat$"></SPAN> on process
<SPAN CLASS="MATH"><IMG
WIDTH="9" HEIGHT="17" ALIGN="BOTTOM" BORDER="0"
SRC="img4.png"
ALT="$i$"></SPAN>, the output <SPAN CLASS="MATH"><IMG
WIDTH="26" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="27" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img29.png"
ALT="$res$"></SPAN> is equivalent to the computation of
<BR><P></P>

@ -89,7 +89,7 @@ Intent: <B>in</B>.
<BR>
Specified as: an integer, real or complex variable, which may be a
scalar, or a rank 1 or 2 array, or a character or logical scalar. Type, kind and rank must agree on sender and receiver process; if <SPAN CLASS="MATH"><IMG
WIDTH="18" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="18" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img140.png"
ALT="$m$"></SPAN> is
not specified, size must agree as well.
@ -107,7 +107,7 @@ Specified as: an integer value <!-- MATH
$0<= dst <= np-1$
-->
<SPAN CLASS="MATH"><IMG
WIDTH="145" HEIGHT="30" ALIGN="MIDDLE" BORDER="0"
WIDTH="146" HEIGHT="30" ALIGN="MIDDLE" BORDER="0"
SRC="img141.png"
ALT="$0&lt;= dst &lt;= np-1$"></SPAN>.
<BR></DD>
@ -124,12 +124,12 @@ Specified as: an integer value <!-- MATH
$0<= m <= size(dat,1)$
-->
<SPAN CLASS="MATH"><IMG
WIDTH="171" HEIGHT="31" ALIGN="MIDDLE" BORDER="0"
WIDTH="171" HEIGHT="32" ALIGN="MIDDLE" BORDER="0"
SRC="img142.png"
ALT="$0&lt;= m &lt;= size(dat,1)$"></SPAN>.
<BR>
When <SPAN CLASS="MATH"><IMG
WIDTH="26" HEIGHT="15" ALIGN="BOTTOM" BORDER="0"
WIDTH="26" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
SRC="img138.png"
ALT="$dat$"></SPAN> is a rank 2 array, specifies the number of rows to be sent
independently of the leading dimension <SPAN CLASS="MATH"><IMG
@ -153,7 +153,7 @@ same value on sending and receiving processes.
<OL>
<LI>This subroutine implies a synchronization, but only between the
calling process and the destination process <SPAN CLASS="MATH"><IMG
WIDTH="25" HEIGHT="15" ALIGN="BOTTOM" BORDER="0"
WIDTH="25" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
SRC="img144.png"
ALT="$dst$"></SPAN>.
</LI>

@ -107,12 +107,12 @@ Specified as: an integer value <!-- MATH
$0<= m <= size(dat,1)$
-->
<SPAN CLASS="MATH"><IMG
WIDTH="171" HEIGHT="31" ALIGN="MIDDLE" BORDER="0"
WIDTH="171" HEIGHT="32" ALIGN="MIDDLE" BORDER="0"
SRC="img142.png"
ALT="$0&lt;= m &lt;= size(dat,1)$"></SPAN>.
<BR>
When <SPAN CLASS="MATH"><IMG
WIDTH="26" HEIGHT="15" ALIGN="BOTTOM" BORDER="0"
WIDTH="26" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
SRC="img138.png"
ALT="$dat$"></SPAN> is a rank 2 array, specifies the number of rows to be sent
independently of the leading dimension <SPAN CLASS="MATH"><IMG
@ -139,7 +139,7 @@ Intent: <B>inout</B>.
<BR>
Specified as: an integer, real or complex variable, which may be a
scalar, or a rank 1 or 2 array, or a character or logical scalar. Type, kind and rank must agree on sender and receiver process; if <SPAN CLASS="MATH"><IMG
WIDTH="18" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="18" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img140.png"
ALT="$m$"></SPAN> is
not specified, size must agree as well.
@ -152,7 +152,7 @@ not specified, size must agree as well.
<OL>
<LI>This subroutine implies a synchronization, but only between the
calling process and the source process <SPAN CLASS="MATH"><IMG
WIDTH="26" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="26" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img146.png"
ALT="$src$"></SPAN>.
</LI>

@ -212,7 +212,7 @@ ifstarsubroutinesubroutinepsb_errorPrints the error stack content and aborts
<P>
<BR>
<IMG
WIDTH="214" HEIGHT="16" ALIGN="BOTTOM" BORDER="0"
WIDTH="213" HEIGHT="16" ALIGN="BOTTOM" BORDER="0"
SRC="img151.png"
ALT="\begin{lstlisting}
call psb_error(icontxt)
@ -282,7 +282,7 @@ ifstarsubroutinesubroutinepsb_set_erractionSet the type of action to be
<P>
<BR>
<IMG
WIDTH="290" HEIGHT="16" ALIGN="BOTTOM" BORDER="0"
WIDTH="289" HEIGHT="16" ALIGN="BOTTOM" BORDER="0"
SRC="img153.png"
ALT="\begin{lstlisting}
call psb_set_erraction(err_act)

@ -56,7 +56,7 @@
<P>
<BR>
<IMG
WIDTH="462" HEIGHT="16" ALIGN="BOTTOM" BORDER="0"
WIDTH="461" HEIGHT="16" ALIGN="BOTTOM" BORDER="0"
SRC="img154.png"
ALT="\begin{lstlisting}
call hb_read(a, iret, iunit, filename, b, mtitle)

@ -56,7 +56,7 @@ hb_write -- Write a sparse matrix to a file
<P>
<BR>
<IMG
WIDTH="538" HEIGHT="16" ALIGN="BOTTOM" BORDER="0"
WIDTH="537" HEIGHT="16" ALIGN="BOTTOM" BORDER="0"
SRC="img155.png"
ALT="\begin{lstlisting}
call hb_write(a, iret, iunit, filename, key, rhs, mtitle)

@ -56,7 +56,7 @@ mm_mat_read -- Read a sparse matrix from a
<P>
<BR>
<IMG
WIDTH="396" HEIGHT="16" ALIGN="BOTTOM" BORDER="0"
WIDTH="394" HEIGHT="16" ALIGN="BOTTOM" BORDER="0"
SRC="img156.png"
ALT="\begin{lstlisting}
call mm_mat_read(a, iret, iunit, filename)

@ -56,7 +56,7 @@ mm_array_read -- Read a dense array from a
<P>
<BR>
<IMG
WIDTH="414" HEIGHT="16" ALIGN="BOTTOM" BORDER="0"
WIDTH="413" HEIGHT="16" ALIGN="BOTTOM" BORDER="0"
SRC="img157.png"
ALT="\begin{lstlisting}
call mm_array_read(b, iret, iunit, filename)

@ -56,7 +56,7 @@ mm_mat_write -- Write a sparse matrix to a
<P>
<BR>
<IMG
WIDTH="482" HEIGHT="16" ALIGN="BOTTOM" BORDER="0"
WIDTH="480" HEIGHT="16" ALIGN="BOTTOM" BORDER="0"
SRC="img158.png"
ALT="\begin{lstlisting}
call mm_mat_write(a, mtitle, iret, iunit, filename)

@ -72,7 +72,7 @@ err = \frac{\|r_i\|}{(\|A\|\|x_i\|+\|b\|)} < eps
-->
<IMG
WIDTH="205" HEIGHT="44" BORDER="0"
WIDTH="204" HEIGHT="44" BORDER="0"
SRC="img163.png"
ALT="\begin{displaymath}err = \frac{\Vert r_i\Vert}{(\Vert A\Vert\Vert x_i\Vert+\Vert b\Vert)} &lt; eps \end{displaymath}">
</DIV>
@ -110,7 +110,7 @@ err = \frac{\|r_i\|}{\|r_0\|_2} < eps
-->
<IMG
WIDTH="129" HEIGHT="44" BORDER="0"
WIDTH="128" HEIGHT="44" BORDER="0"
SRC="img165.png"
ALT="\begin{displaymath}err = \frac{\Vert r_i\Vert}{\Vert r_0\Vert _2} &lt; eps \end{displaymath}">
</DIV>
@ -288,13 +288,22 @@ Specified as: an integer variable <SPAN CLASS="MATH"><IMG
convergence every <SPAN CLASS="MATH"><IMG
WIDTH="45" HEIGHT="17" ALIGN="BOTTOM" BORDER="0"
SRC="img172.png"
ALT="$itrace$"></SPAN> iterations.
ALT="$itrace$"></SPAN> iterations. If <SPAN CLASS="MATH"><IMG
WIDTH="29" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img173.png"
ALT="$=0$"></SPAN> print a message in
case of convergence failure.
<BR>
Scope: <B>global</B>
<BR>
Type: <B>optional</B>
<BR>
Intent: <B>in</B>.
<BR>
Default: <SPAN CLASS="MATH"><IMG
WIDTH="87" HEIGHT="29" ALIGN="MIDDLE" BORDER="0"
SRC="img174.png"
ALT="$itrace = -1$"></SPAN>.
<BR></DD>
<DT><STRONG>irst</STRONG></DT>
<DD>An integer specifying the restart parameter.
@ -307,7 +316,7 @@ Intent: <B>in</B>.
<BR>
Values: <SPAN CLASS="MATH"><IMG
WIDTH="60" HEIGHT="29" ALIGN="MIDDLE" BORDER="0"
SRC="img173.png"
SRC="img175.png"
ALT="$irst&gt;0$"></SPAN>. This is employed for the BiCGSTABL or RGMRES
methods, otherwise it is ignored.
@ -367,8 +376,8 @@ Returned as: a real number.
SRC="img1.png"
ALT="$A$"></SPAN>; only
available with the <SPAN CLASS="MATH"><IMG
WIDTH="29" HEIGHT="15" ALIGN="BOTTOM" BORDER="0"
SRC="img174.png"
WIDTH="29" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
SRC="img176.png"
ALT="$CG$"></SPAN> method on real data.
<BR>
Scope: <B>global</B>

@ -132,7 +132,7 @@ subdomain by <SPAN CLASS="MATH"><IMG
WIDTH="14" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
SRC="img7.png"
ALT="$\cal I$"></SPAN>, <SPAN CLASS="MATH"><IMG
WIDTH="15" HEIGHT="16" ALIGN="BOTTOM" BORDER="0"
WIDTH="15" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
SRC="img8.png"
ALT="$\cal B$"></SPAN> and <SPAN CLASS="MATH"><IMG
WIDTH="18" HEIGHT="15" ALIGN="BOTTOM" BORDER="0"
@ -170,7 +170,7 @@ Point classfication.</CAPTION>
<TR><TD>
<DIV ALIGN="CENTER">
<IMG
WIDTH="318" HEIGHT="285" ALIGN="BOTTOM" BORDER="0"
WIDTH="318" HEIGHT="283" ALIGN="BOTTOM" BORDER="0"
SRC="img13.png"
ALT="\includegraphics[scale=0.65]{figures/points.eps}">
<!-- MATH

@ -89,15 +89,15 @@ Type: <B>optional</B>; default: entire vector.
<DT><STRONG>Function value</STRONG></DT>
<DD>An allocatable array holding a copy of the dense
vector contents. If the argument <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img16.png"
ALT="$n$"></SPAN> is specified, the size of the
returned array equals the minimum between <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img16.png"
ALT="$n$"></SPAN> and the internal size
of the vector, or 0 if <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img16.png"
ALT="$n$"></SPAN> is negative; otherwise, the size of the
array is the same as the internal size of the vector.

@ -86,13 +86,13 @@ Data types</CAPTION>
<DIV ALIGN="CENTER">
<TABLE CELLPADDING=3 BORDER="1" ALIGN="CENTER">
<TR><TD ALIGN="LEFT"><SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>, <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="30" ALIGN="MIDDLE" BORDER="0"
SRC="img21.png"
ALT="$y$"></SPAN>, <SPAN CLASS="MATH"><IMG
WIDTH="14" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="14" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img22.png"
ALT="$\alpha$"></SPAN>, <SPAN CLASS="MATH"><IMG
WIDTH="14" HEIGHT="30" ALIGN="MIDDLE" BORDER="0"
@ -129,7 +129,7 @@ Data types</CAPTION>
</DD>
<DT><STRONG>alpha</STRONG></DT>
<DD>the scalar <SPAN CLASS="MATH"><IMG
WIDTH="14" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="14" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img22.png"
ALT="$\alpha$"></SPAN>.
<BR>
@ -145,7 +145,7 @@ type indicated in Table&nbsp;<A HREF="#tab:f90axpby">1</A>.
<DT><STRONG>x</STRONG></DT>
<DD>the local portion of global dense matrix
<SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>.
<BR>
@ -158,7 +158,7 @@ Intent: <B>in</B>.
Specified as: a rank one or two array or an object of type vdata<TT>psb_T_vect_type</TT>
containing numbers of type
specified in Table&nbsp;<A HREF="#tab:f90axpby">1</A>. The rank of <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN> must be the same of <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="30" ALIGN="MIDDLE" BORDER="0"
@ -197,7 +197,7 @@ indicated in Table&nbsp;<A HREF="#tab:f90axpby">1</A>. The rank of <SPAN CLASS=
WIDTH="13" HEIGHT="30" ALIGN="MIDDLE" BORDER="0"
SRC="img21.png"
ALT="$y$"></SPAN> must be the same of <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>.
</DD>

@ -54,7 +54,7 @@ psb_gedot -- Dot Product</A>
<P>
This function computes dot product between two vectors <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN> and
<SPAN CLASS="MATH"><IMG
@ -63,7 +63,7 @@ This function computes dot product between two vectors <SPAN CLASS="MATH"><IMG
ALT="$y$"></SPAN>.
<BR>
If <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN> and <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="30" ALIGN="MIDDLE" BORDER="0"
@ -86,7 +86,7 @@ dot \leftarrow x^T y
<BR CLEAR="ALL">
<P></P>
Else if <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN> and <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="30" ALIGN="MIDDLE" BORDER="0"
@ -101,7 +101,7 @@ dot \leftarrow x^H y
-->
<IMG
WIDTH="75" HEIGHT="27" BORDER="0"
WIDTH="76" HEIGHT="27" BORDER="0"
SRC="img25.png"
ALT="\begin{displaymath}dot \leftarrow x^H y\end{displaymath}">
</DIV>
@ -121,10 +121,10 @@ Data types</CAPTION>
<DIV ALIGN="CENTER">
<TABLE CELLPADDING=3 BORDER="1" ALIGN="CENTER">
<TR><TD ALIGN="LEFT"><SPAN CLASS="MATH"><IMG
WIDTH="26" HEIGHT="15" ALIGN="BOTTOM" BORDER="0"
WIDTH="25" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
SRC="img26.png"
ALT="$dot$"></SPAN>, <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>, <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="30" ALIGN="MIDDLE" BORDER="0"
@ -162,7 +162,7 @@ Data types</CAPTION>
<DT><STRONG>x</STRONG></DT>
<DD>the local portion of global dense matrix
<SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>.
<BR>
@ -175,7 +175,7 @@ Intent: <B>in</B>.
Specified as: a rank one or two array or an object of type vdata<TT>psb_T_vect_type</TT>
containing numbers of type specified in
Table&nbsp;<A HREF="#tab:f90dot">2</A>. The rank of <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN> must be the same of <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="30" ALIGN="MIDDLE" BORDER="0"
@ -201,7 +201,7 @@ Table&nbsp;<A HREF="#tab:f90dot">2</A>. The rank of <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="30" ALIGN="MIDDLE" BORDER="0"
SRC="img21.png"
ALT="$y$"></SPAN> must be the same of <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>.
</DD>
@ -236,7 +236,7 @@ Default: <code>global=.true.</code>
</DD>
<DT><STRONG>Function value</STRONG></DT>
<DD>is the dot product of vectors <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN> and <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="30" ALIGN="MIDDLE" BORDER="0"
@ -272,7 +272,7 @@ An integer value; 0 means no error has been detected.
by using the following scheme:
<BR>
<IMG
WIDTH="500" HEIGHT="73" ALIGN="BOTTOM" BORDER="0"
WIDTH="499" HEIGHT="73" ALIGN="BOTTOM" BORDER="0"
SRC="img27.png"
ALT="\begin{lstlisting}
vres(1) = psb_gedot(x1,y1,desc_a,info,global=.false.)

@ -55,7 +55,7 @@ psb_gedots -- Generalized Dot Product</A>
<P>
This subroutine computes a series of dot products among the columns of
two dense matrices <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN> and <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="30" ALIGN="MIDDLE" BORDER="0"
@ -78,17 +78,17 @@ res(i) \leftarrow x(:,i)^T y(:,i)
<P></P>
If the matrices are complex, then the
usual convention applies, i.e. the conjugate transpose of <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN> is
used. If <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN> and <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="30" ALIGN="MIDDLE" BORDER="0"
SRC="img21.png"
ALT="$y$"></SPAN> are of rank one, then <SPAN CLASS="MATH"><IMG
WIDTH="26" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="27" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img29.png"
ALT="$res$"></SPAN> is a scalar, else it
is a rank one array.
@ -106,10 +106,10 @@ Data types</CAPTION>
<DIV ALIGN="CENTER">
<TABLE CELLPADDING=3 BORDER="1" ALIGN="CENTER">
<TR><TD ALIGN="LEFT"><SPAN CLASS="MATH"><IMG
WIDTH="26" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="27" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img29.png"
ALT="$res$"></SPAN>, <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>, <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="30" ALIGN="MIDDLE" BORDER="0"
@ -147,7 +147,7 @@ Data types</CAPTION>
<DT><STRONG>x</STRONG></DT>
<DD>the local portion of global dense matrix
<SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>.
<BR>
@ -160,7 +160,7 @@ Intent: <B>in</B>.
Specified as: a rank one or two array or an object of type vdata<TT>psb_T_vect_type</TT>
containing numbers of type specified in
Table&nbsp;<A HREF="#tab:f90mdot">3</A>. The rank of <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN> must be the same of <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="30" ALIGN="MIDDLE" BORDER="0"
@ -186,7 +186,7 @@ Table&nbsp;<A HREF="#tab:f90mdot">3</A>. The rank of <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="30" ALIGN="MIDDLE" BORDER="0"
SRC="img21.png"
ALT="$y$"></SPAN> must be the same of <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>.
</DD>
@ -206,7 +206,7 @@ Specified as: an object of type descdata<TT>psb_desc_type</TT>.
</DD>
<DT><STRONG>res</STRONG></DT>
<DD>is the dot product of vectors <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN> and <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="30" ALIGN="MIDDLE" BORDER="0"

@ -55,12 +55,12 @@ psb_normi -- Infinity-Norm of Vector</A>
<P>
This function computes
the infinity-norm of a vector <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>.
<BR>
If <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN> is a real vector
it computes infinity norm as:
@ -80,7 +80,7 @@ amax \leftarrow \max_i |x_i|
<BR CLEAR="ALL">
<P></P>
else if <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN> is a complex vector then it computes the infinity-norm as:
<BR><P></P>
@ -92,7 +92,7 @@ amax \leftarrow \max_i {(|re(x_i)| + |im(x_i)|)}
-->
<IMG
WIDTH="233" HEIGHT="36" BORDER="0"
WIDTH="232" HEIGHT="36" BORDER="0"
SRC="img31.png"
ALT="\begin{displaymath}amax \leftarrow \max_i {(\vert re(x_i)\vert + \vert im(x_i)\vert)}\end{displaymath}">
</DIV>
@ -115,11 +115,11 @@ Data types</CAPTION>
<DIV ALIGN="CENTER">
<TABLE CELLPADDING=3 BORDER="1" ALIGN="CENTER">
<TR><TD ALIGN="LEFT"><SPAN CLASS="MATH"><IMG
WIDTH="44" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="43" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img32.png"
ALT="$amax$"></SPAN></TD>
<TD ALIGN="LEFT"><SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN></TD>
<TD ALIGN="LEFT"><B>Function</B></TD>
@ -158,7 +158,7 @@ Data types</CAPTION>
<DT><STRONG>x</STRONG></DT>
<DD>the local portion of global dense matrix
<SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>.
@ -205,7 +205,7 @@ Default: <code>global=.true.</code>
</DD>
<DT><STRONG>Function value</STRONG></DT>
<DD>is the infinity norm of vector <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>.
<BR>
@ -238,7 +238,7 @@ An integer value; 0 means no error has been detected.
by using the following scheme:
<BR>
<IMG
WIDTH="481" HEIGHT="73" ALIGN="BOTTOM" BORDER="0"
WIDTH="480" HEIGHT="73" ALIGN="BOTTOM" BORDER="0"
SRC="img33.png"
ALT="\begin{lstlisting}
vres(1) = psb_geamax(x1,desc_a,info,global=.false.)

@ -55,7 +55,7 @@ psb_geamaxs -- Generalized Infinity Norm</A>
<P>
This subroutine computes a series of infinity norms on the columns of
a dense matrix <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>:
<BR><P></P>
@ -67,7 +67,7 @@ res(i) \leftarrow \max_k |x(k,i)|
-->
<IMG
WIDTH="148" HEIGHT="36" BORDER="0"
WIDTH="147" HEIGHT="36" BORDER="0"
SRC="img34.png"
ALT="\begin{displaymath}res(i) \leftarrow \max_k \vert x(k,i)\vert \end{displaymath}">
</DIV>
@ -89,11 +89,11 @@ Data types</CAPTION>
<DIV ALIGN="CENTER">
<TABLE CELLPADDING=3 BORDER="1" ALIGN="CENTER">
<TR><TD ALIGN="LEFT"><SPAN CLASS="MATH"><IMG
WIDTH="26" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="27" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img29.png"
ALT="$res$"></SPAN></TD>
<TD ALIGN="LEFT"><SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN></TD>
<TD ALIGN="LEFT"><B>Subroutine</B></TD>
@ -132,7 +132,7 @@ Data types</CAPTION>
<DT><STRONG>x</STRONG></DT>
<DD>the local portion of global dense matrix
<SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>.
<BR>
@ -162,7 +162,7 @@ Specified as: an object of type descdata<TT>psb_desc_type</TT>.
</DD>
<DT><STRONG>res</STRONG></DT>
<DD>is the infinity norm of the columns of <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>.
<BR>

@ -54,12 +54,12 @@ psb_norm1 -- 1-Norm of Vector</A>
<P>
This function computes the 1-norm of a vector <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>.
<BR>
If <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN> is a real vector
it computes 1-norm as:
@ -79,7 +79,7 @@ asum \leftarrow \|x_i\|
<BR CLEAR="ALL">
<P></P>
else if <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN> is a complex vector then it computes 1-norm as:
<BR><P></P>
@ -114,11 +114,11 @@ Data types</CAPTION>
<DIV ALIGN="CENTER">
<TABLE CELLPADDING=3 BORDER="1" ALIGN="CENTER">
<TR><TD ALIGN="LEFT"><SPAN CLASS="MATH"><IMG
WIDTH="43" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="43" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img37.png"
ALT="$asum$"></SPAN></TD>
<TD ALIGN="LEFT"><SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN></TD>
<TD ALIGN="LEFT"><B>Function</B></TD>
@ -157,7 +157,7 @@ Data types</CAPTION>
<DT><STRONG>x</STRONG></DT>
<DD>the local portion of global dense matrix
<SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>.
@ -203,7 +203,7 @@ Default: <code>global=.true.</code>
</DD>
<DT><STRONG>Function value</STRONG></DT>
<DD>is the 1-norm of vector <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>.
<BR>
@ -236,7 +236,7 @@ An integer value; 0 means no error has been detected.
by using the following scheme:
<BR>
<IMG
WIDTH="481" HEIGHT="73" ALIGN="BOTTOM" BORDER="0"
WIDTH="480" HEIGHT="73" ALIGN="BOTTOM" BORDER="0"
SRC="img38.png"
ALT="\begin{lstlisting}
vres(1) = psb_geasum(x1,desc_a,info,global=.false.)

@ -55,7 +55,7 @@ psb_geasums -- Generalized 1-Norm of Vector</A>
<P>
This subroutine computes a series of 1-norms on the columns of
a dense matrix <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>:
<BR><P></P>
@ -67,19 +67,19 @@ res(i) \leftarrow \max_k |x(k,i)|
-->
<IMG
WIDTH="148" HEIGHT="36" BORDER="0"
WIDTH="147" HEIGHT="36" BORDER="0"
SRC="img34.png"
ALT="\begin{displaymath}res(i) \leftarrow \max_k \vert x(k,i)\vert \end{displaymath}">
</DIV>
<BR CLEAR="ALL">
<P></P>
This function computes the 1-norm of a vector <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>.
<BR>
If <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN> is a real vector
it computes 1-norm as:
@ -99,7 +99,7 @@ res(i) \leftarrow \|x_i\|
<BR CLEAR="ALL">
<P></P>
else if <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN> is a complex vector then it computes 1-norm as:
<BR><P></P>
@ -111,7 +111,7 @@ res(i) \leftarrow \|re(x)\|_1 + \|im(x)\|_1
-->
<IMG
WIDTH="207" HEIGHT="29" BORDER="0"
WIDTH="206" HEIGHT="28" BORDER="0"
SRC="img40.png"
ALT="\begin{displaymath}res(i) \leftarrow \Vert re(x)\Vert _1 + \Vert im(x)\Vert _1\end{displaymath}">
</DIV>
@ -133,11 +133,11 @@ Data types</CAPTION>
<DIV ALIGN="CENTER">
<TABLE CELLPADDING=3 BORDER="1" ALIGN="CENTER">
<TR><TD ALIGN="LEFT"><SPAN CLASS="MATH"><IMG
WIDTH="26" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="27" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img29.png"
ALT="$res$"></SPAN></TD>
<TD ALIGN="LEFT"><SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN></TD>
<TD ALIGN="LEFT"><B>Subroutine</B></TD>
@ -176,7 +176,7 @@ Data types</CAPTION>
<DT><STRONG>x</STRONG></DT>
<DD>the local portion of global dense matrix
<SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>.
@ -209,7 +209,7 @@ Specified as: an object of type descdata<TT>psb_desc_type</TT>.
</DD>
<DT><STRONG>res</STRONG></DT>
<DD>contains the 1-norm of (the columns of) <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>.
<BR>

@ -61,7 +61,7 @@ space to which there corresponds an index space and a matrix sparsity
pattern. As an example, consider a cell-centered finite-volume
discretization of the Navier-Stokes equations on a simulation domain;
the index space <SPAN CLASS="MATH"><IMG
WIDTH="45" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="46" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img15.png"
ALT="$1\dots n$"></SPAN> is isomorphic to the set of cell centers,
whereas the pattern of the associated linear system matrix is
@ -72,7 +72,7 @@ by the discretization stencil.
Thus the first order of business is to establish an index space, and
this is done with a call to <code>psb_cdall</code> in which we specify the
size of the index space <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img16.png"
ALT="$n$"></SPAN> and the allocation of the elements of the
index space to the various processes making up the MPI (virtual)
@ -81,7 +81,7 @@ parallel machine.
<P>
The index space is partitioned among processes, and this creates a
mapping from the &ldquo;global&rdquo; numbering <SPAN CLASS="MATH"><IMG
WIDTH="45" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="46" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img15.png"
ALT="$1\dots n$"></SPAN> to a numbering
&ldquo;local&rdquo; to each process; each process <SPAN CLASS="MATH"><IMG
@ -96,7 +96,7 @@ mapping from the &ldquo;global&rdquo; numbering <SPAN CLASS="MATH"><IMG
SRC="img17.png"
ALT="$1\dots n_{\hbox{row}_i}$"></SPAN>, each element of which corresponds to a certain
element of <SPAN CLASS="MATH"><IMG
WIDTH="45" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="46" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img15.png"
ALT="$1\dots n$"></SPAN>. The user does not set explicitly this mapping;
when the application needs to indicate to which element of the index
@ -106,7 +106,7 @@ library will translate into the appropriate &ldquo;local&rdquo; numbering.
<P>
For a given index space <SPAN CLASS="MATH"><IMG
WIDTH="45" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="46" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img15.png"
ALT="$1\dots n$"></SPAN> there are many possible associated
topologies, i.e. many different discretization stencils; thus the

@ -54,12 +54,12 @@ psb_norm2 -- 2-Norm of Vector</A>
<P>
This function computes the 2-norm of a vector <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>.
<BR>
If <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN> is a real vector
it computes 2-norm as:
@ -79,7 +79,7 @@ nrm2 \leftarrow \sqrt{x^T x}
<BR CLEAR="ALL">
<P></P>
else if <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN> is a complex vector then it computes 2-norm as:
<BR><P></P>
@ -108,11 +108,11 @@ Data types</CAPTION>
<DIV ALIGN="CENTER">
<TABLE CELLPADDING=3 BORDER="1" ALIGN="CENTER">
<TR><TD ALIGN="LEFT"><SPAN CLASS="MATH"><IMG
WIDTH="43" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="44" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img43.png"
ALT="$nrm2$"></SPAN></TD>
<TD ALIGN="LEFT"><SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN></TD>
<TD ALIGN="LEFT"><B>Function</B></TD>
@ -157,7 +157,7 @@ psb_norm2(x, desc_a, info [,global])
<DT><STRONG>x</STRONG></DT>
<DD>the local portion of global dense matrix
<SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>.
<BR>
@ -203,7 +203,7 @@ Default: <code>global=.true.</code>
</DD>
<DT><STRONG>Function Value</STRONG></DT>
<DD>is the 2-norm of vector <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>.
<BR>
@ -238,7 +238,7 @@ An integer value; 0 means no error has been detected.
by using the following scheme:
<BR>
<IMG
WIDTH="481" HEIGHT="73" ALIGN="BOTTOM" BORDER="0"
WIDTH="480" HEIGHT="73" ALIGN="BOTTOM" BORDER="0"
SRC="img44.png"
ALT="\begin{lstlisting}
vres(1) = psb_genrm2(x1,desc_a,info,global=.false.)

@ -55,7 +55,7 @@ psb_genrm2s -- Generalized 2-Norm of Vector</A>
<P>
This subroutine computes a series of 2-norms on the columns of
a dense matrix <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>:
<BR><P></P>
@ -67,7 +67,7 @@ res(i) \leftarrow \|x(:,i)\|_2
-->
<IMG
WIDTH="125" HEIGHT="29" BORDER="0"
WIDTH="124" HEIGHT="28" BORDER="0"
SRC="img45.png"
ALT="\begin{displaymath}res(i) \leftarrow \Vert x(:,i)\Vert _2 \end{displaymath}">
</DIV>
@ -89,11 +89,11 @@ Data types</CAPTION>
<DIV ALIGN="CENTER">
<TABLE CELLPADDING=3 BORDER="1" ALIGN="CENTER">
<TR><TD ALIGN="LEFT"><SPAN CLASS="MATH"><IMG
WIDTH="26" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="27" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img29.png"
ALT="$res$"></SPAN></TD>
<TD ALIGN="LEFT"><SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN></TD>
<TD ALIGN="LEFT"><B>Subroutine</B></TD>
@ -132,7 +132,7 @@ Data types</CAPTION>
<DT><STRONG>x</STRONG></DT>
<DD>the local portion of global dense matrix
<SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>.
@ -165,7 +165,7 @@ Specified as: an object of type descdata<TT>psb_desc_type</TT>.
</DD>
<DT><STRONG>res</STRONG></DT>
<DD>contains the 1-norm of (the columns of) <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>.
<BR>

@ -68,7 +68,7 @@ nrm1 \leftarrow \|A\|_1
-->
<IMG
WIDTH="98" HEIGHT="29" BORDER="0"
WIDTH="97" HEIGHT="28" BORDER="0"
SRC="img46.png"
ALT="\begin{displaymath}nrm1 \leftarrow \Vert A\Vert _1 \end{displaymath}">
</DIV>

@ -68,7 +68,7 @@ nrmi \leftarrow \|A\|_\infty
-->
<IMG
WIDTH="101" HEIGHT="29" BORDER="0"
WIDTH="101" HEIGHT="28" BORDER="0"
SRC="img47.png"
ALT="\begin{displaymath}nrmi \leftarrow \Vert A\Vert _\infty \end{displaymath}">
</DIV>

@ -88,7 +88,7 @@ y \leftarrow \alpha A^T x + \beta y
<A NAME="eq:f90spmm_tra"></A>
<TABLE WIDTH="100%" ALIGN="CENTER">
<TR VALIGN="MIDDLE"><TD ALIGN="CENTER" NOWRAP><A NAME="eq:f90spmm_tra"></A><IMG
WIDTH="111" HEIGHT="27" BORDER="0"
WIDTH="112" HEIGHT="27" BORDER="0"
SRC="img49.png"
ALT="\begin{displaymath}
y \leftarrow \alpha A^T x + \beta y
@ -122,7 +122,7 @@ y \leftarrow \alpha A^H x + \beta y
where:
<DL>
<DT><STRONG><SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN></STRONG></DT>
<DD>is the global dense matrix <SPAN CLASS="MATH"><IMG
@ -163,13 +163,13 @@ Data types</CAPTION>
WIDTH="16" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
SRC="img1.png"
ALT="$A$"></SPAN>, <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>, <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="30" ALIGN="MIDDLE" BORDER="0"
SRC="img21.png"
ALT="$y$"></SPAN>, <SPAN CLASS="MATH"><IMG
WIDTH="14" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="14" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img22.png"
ALT="$\alpha$"></SPAN>, <SPAN CLASS="MATH"><IMG
WIDTH="14" HEIGHT="30" ALIGN="MIDDLE" BORDER="0"
@ -213,7 +213,7 @@ call psb_spmm(alpha, a, x, beta, y,desc_a, info, &amp;
</DD>
<DT><STRONG>alpha</STRONG></DT>
<DD>the scalar <SPAN CLASS="MATH"><IMG
WIDTH="14" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="14" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img22.png"
ALT="$\alpha$"></SPAN>.
<BR>
@ -244,7 +244,7 @@ Specified as: an object of type spdata<TT>psb_Tspmat_type</TT>.
<DT><STRONG>x</STRONG></DT>
<DD>the local portion of global dense matrix
<SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>.
@ -258,7 +258,7 @@ Intent: <B>in</B>.
Specified as: a rank one or two array or an object of type vdata<TT>psb_T_vect_type</TT>
containing numbers of type specified in
Table&nbsp;<A HREF="#tab:f90spmm">12</A>. The rank of <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN> must be the same of <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="30" ALIGN="MIDDLE" BORDER="0"
@ -299,7 +299,7 @@ Table&nbsp;<A HREF="#tab:f90spmm">12</A>. The rank of <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="30" ALIGN="MIDDLE" BORDER="0"
SRC="img21.png"
ALT="$y$"></SPAN> must be the same of <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>.
</DD>
@ -336,7 +336,7 @@ Type: <B>optional</B>
Intent: <B>in</B>.
<BR>
Default: <SPAN CLASS="MATH"><IMG
WIDTH="78" HEIGHT="16" ALIGN="BOTTOM" BORDER="0"
WIDTH="79" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
SRC="img53.png"
ALT="$trans = N$"></SPAN>
<BR>
@ -354,7 +354,7 @@ Type: <B>optional</B>
Intent: <B>inout</B>.
<BR>
Specified as: a rank one array of the same type of <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN> and <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="30" ALIGN="MIDDLE" BORDER="0"

@ -86,7 +86,7 @@ y &amp;\leftarrow&amp; \alpha T^{-H} D x + \beta y\\
where:
<DL>
<DT><STRONG><SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN></STRONG></DT>
<DD>is the global dense matrix <SPAN CLASS="MATH"><IMG
@ -140,7 +140,7 @@ Data types</CAPTION>
WIDTH="16" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
SRC="img55.png"
ALT="$T$"></SPAN>, <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>, <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="30" ALIGN="MIDDLE" BORDER="0"
@ -149,7 +149,7 @@ Data types</CAPTION>
WIDTH="17" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
SRC="img56.png"
ALT="$D$"></SPAN>, <SPAN CLASS="MATH"><IMG
WIDTH="14" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="14" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img22.png"
ALT="$\alpha$"></SPAN>, <SPAN CLASS="MATH"><IMG
WIDTH="14" HEIGHT="30" ALIGN="MIDDLE" BORDER="0"
@ -186,7 +186,7 @@ Data types</CAPTION>
</DD>
<DT><STRONG>alpha</STRONG></DT>
<DD>the scalar <SPAN CLASS="MATH"><IMG
WIDTH="14" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="14" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img22.png"
ALT="$\alpha$"></SPAN>.
<BR>
@ -218,7 +218,7 @@ Specified as: an object type specified in
<DT><STRONG>x</STRONG></DT>
<DD>the local portion of global dense matrix
<SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>.
@ -232,7 +232,7 @@ Intent: <B>in</B>.
Specified as: a rank one or two array or an object of type vdata<TT>psb_T_vect_type</TT>
containing numbers of type specified in
Table&nbsp;<A HREF="#tab:f90spsm">13</A>. The rank of <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN> must be the same of <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="30" ALIGN="MIDDLE" BORDER="0"
@ -273,7 +273,7 @@ Table&nbsp;<A HREF="#tab:f90spsm">13</A>. The rank of <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="30" ALIGN="MIDDLE" BORDER="0"
SRC="img21.png"
ALT="$y$"></SPAN> must be the same of <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>.
</DD>
@ -308,7 +308,7 @@ Type: <B>optional</B>
Intent: <B>in</B>.
<BR>
Default: <SPAN CLASS="MATH"><IMG
WIDTH="78" HEIGHT="16" ALIGN="BOTTOM" BORDER="0"
WIDTH="79" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
SRC="img53.png"
ALT="$trans = N$"></SPAN>
<BR>
@ -334,7 +334,7 @@ Type: <B>optional</B>
Intent: <B>in</B>.
<BR>
Default: <SPAN CLASS="MATH"><IMG
WIDTH="76" HEIGHT="15" ALIGN="BOTTOM" BORDER="0"
WIDTH="76" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
SRC="img57.png"
ALT="$unitd = U$"></SPAN>
<BR>
@ -397,7 +397,7 @@ Type: <B>optional</B>
Intent: <B>inout</B>.
<BR>
Specified as: a rank one array of the same type of <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN> with the
TARGET attribute.

@ -75,7 +75,7 @@ x \leftarrow x
where:
<DL>
<DT><STRONG><SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN></STRONG></DT>
<DD>is a global dense submatrix.
@ -92,10 +92,10 @@ Data types</CAPTION>
<DIV ALIGN="CENTER">
<TABLE CELLPADDING=3 BORDER="1" ALIGN="CENTER">
<TR><TD ALIGN="LEFT"><SPAN CLASS="MATH"><IMG
WIDTH="14" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="14" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img22.png"
ALT="$\alpha$"></SPAN>, <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN></TD>
<TD ALIGN="LEFT"><B>Subroutine</B></TD>
@ -125,7 +125,7 @@ Data types</CAPTION>
<P>
<BR>
<IMG
WIDTH="395" HEIGHT="35" ALIGN="BOTTOM" BORDER="0"
WIDTH="394" HEIGHT="35" ALIGN="BOTTOM" BORDER="0"
SRC="img60.png"
ALT="\begin{lstlisting}
call psb_halo(x, desc_a, info)
@ -143,7 +143,7 @@ call psb_halo(x, desc_a, info, work, data)
</DD>
<DT><STRONG>x</STRONG></DT>
<DD>global dense matrix <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>.
<BR>
@ -178,7 +178,7 @@ Type: <B>optional</B>
Intent: <B>inout</B>.
<BR>
Specified as: a rank one array of the same type of <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>.
</DD>
@ -200,7 +200,7 @@ index list on which to base the data exchange.
</DD>
<DT><STRONG>x</STRONG></DT>
<DD>global dense result matrix <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>.
<BR>
@ -237,12 +237,12 @@ Sample discretization mesh.</CAPTION>
<TR><TD>
<DIV ALIGN="CENTER">
<IMG
WIDTH="358" HEIGHT="313" ALIGN="BOTTOM" BORDER="0"
WIDTH="357" HEIGHT="313" ALIGN="BOTTOM" BORDER="0"
SRC="img61.png"
ALT="\includegraphics[scale=0.45]{figures/try8x8.eps}">
<IMG
WIDTH="358" HEIGHT="313" ALIGN="BOTTOM" BORDER="0"
WIDTH="357" HEIGHT="313" ALIGN="BOTTOM" BORDER="0"
SRC="img62.png"
ALT="\includegraphics[scale=0.45]{figures/try8x8}">
@ -258,7 +258,7 @@ distribution is such that each process will own 32 entries in the
index space, with a halo made of 8 entries placed at local indices 33
through 40. If process 0 assigns an initial value of 1 to its entries
in the <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN> vector, and process 1 assigns a value of 2, then after a
call to <code>psb_halo</code> the contents of the local vectors will be the

@ -74,11 +74,11 @@ x \leftarrow Q x
where:
<DL>
<DT><STRONG><SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN></STRONG></DT>
<DD>is the global dense submatrix <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>
</DD>
@ -90,10 +90,10 @@ where:
operators <SPAN CLASS="MATH"><IMG
WIDTH="22" HEIGHT="30" ALIGN="MIDDLE" BORDER="0"
SRC="img65.png"
ALT="$ P_a$"></SPAN> and <SPAN CLASS="MATH"><IMG
ALT="$P_a$"></SPAN> and <SPAN CLASS="MATH"><IMG
WIDTH="26" HEIGHT="20" ALIGN="BOTTOM" BORDER="0"
SRC="img66.png"
ALT="$ P^{T}$"></SPAN>.
ALT="$P^{T}$"></SPAN>.
</DD>
</DL>
@ -107,7 +107,7 @@ Data types</CAPTION>
<DIV ALIGN="CENTER">
<TABLE CELLPADDING=3 BORDER="1" ALIGN="CENTER">
<TR><TD ALIGN="LEFT"><SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN></TD>
<TD ALIGN="LEFT"><B>Subroutine</B></TD>
@ -134,7 +134,7 @@ Data types</CAPTION>
<P>
<BR>
<IMG
WIDTH="577" HEIGHT="35" ALIGN="BOTTOM" BORDER="0"
WIDTH="576" HEIGHT="35" ALIGN="BOTTOM" BORDER="0"
SRC="img67.png"
ALT="\begin{lstlisting}
call psb_ovrl(x, desc_a, info)
@ -152,7 +152,7 @@ call psb_ovrl(x, desc_a, info, update=update_type, work=work)
</DD>
<DT><STRONG>x</STRONG></DT>
<DD>global dense matrix <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>.
<BR>
@ -222,7 +222,7 @@ Type: <B>optional</B>
Intent: <B>inout</B>.
<BR>
Specified as: a one dimensional array of the same type of <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>.
@ -233,7 +233,7 @@ Specified as: a one dimensional array of the same type of <SPAN CLASS="MATH"><IM
</DD>
<DT><STRONG>x</STRONG></DT>
<DD>global dense result matrix <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN>.
<BR>
@ -270,7 +270,7 @@ An integer value; 0 means no error has been detected.
<LI>The operator <SPAN CLASS="MATH"><IMG
WIDTH="26" HEIGHT="20" ALIGN="BOTTOM" BORDER="0"
SRC="img66.png"
ALT="$ P^{T}$"></SPAN> performs the reduction sum of overlap
ALT="$P^{T}$"></SPAN> performs the reduction sum of overlap
elements; it is a &ldquo;prolongation&rdquo; operator <SPAN CLASS="MATH"><IMG
WIDTH="26" HEIGHT="20" ALIGN="BOTTOM" BORDER="0"
SRC="img68.png"
@ -281,7 +281,7 @@ of data;
<LI>The operator <SPAN CLASS="MATH"><IMG
WIDTH="22" HEIGHT="30" ALIGN="MIDDLE" BORDER="0"
SRC="img65.png"
ALT="$ P_a$"></SPAN> performs a scaling on the overlap elements by
ALT="$P_a$"></SPAN> performs a scaling on the overlap elements by
the amount of replication; thus, when combined with the reduction
operator, it implements the average of replicated elements over all of
their instances.
@ -297,12 +297,12 @@ Sample discretization mesh.</CAPTION>
<TR><TD>
<DIV ALIGN="CENTER">
<IMG
WIDTH="533" HEIGHT="454" ALIGN="BOTTOM" BORDER="0"
WIDTH="533" HEIGHT="453" ALIGN="BOTTOM" BORDER="0"
SRC="img72.png"
ALT="\includegraphics[scale=0.65]{figures/try8x8_ov.eps}">
<IMG
WIDTH="534" HEIGHT="454" ALIGN="BOTTOM" BORDER="0"
WIDTH="533" HEIGHT="453" ALIGN="BOTTOM" BORDER="0"
SRC="img71.png"
ALT="\includegraphics[scale=0.65]{figures/try8x8_ov}">
@ -319,7 +319,7 @@ distribution is such that each process will own 40 entries in the
index space, with an overlap of 16 entries placed at local indices 25
through 40; the halo will run from local index 41 through local index 48.. If process 0 assigns an initial value of 1 to its entries
in the <SPAN CLASS="MATH"><IMG
WIDTH="13" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="13" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img20.png"
ALT="$x$"></SPAN> vector, and process 1 assigns a value of 2, then after a
call to <code>psb_ovrl</code> with <code>psb_avg_</code> and a call to

@ -98,7 +98,7 @@ process <SPAN CLASS="MATH"><IMG
ALT="$i$"></SPAN>.
</DD>
<DT><STRONG><SPAN CLASS="MATH"><IMG
WIDTH="49" HEIGHT="15" ALIGN="BOTTOM" BORDER="0"
WIDTH="49" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
SRC="img77.png"
ALT="$collect$"></SPAN></STRONG></DT>
<DD>is the collect function.
@ -145,7 +145,7 @@ Data types</CAPTION>
<P>
<BR>
<IMG
WIDTH="472" HEIGHT="35" ALIGN="BOTTOM" BORDER="0"
WIDTH="471" HEIGHT="35" ALIGN="BOTTOM" BORDER="0"
SRC="img79.png"
ALT="\begin{lstlisting}
call psb_gather(glob_x, loc_x, desc_a, info, root)

@ -65,7 +65,7 @@ loc\_x_i \leftarrow scatter(glob\_x)
-->
<IMG
WIDTH="168" HEIGHT="29" BORDER="0"
WIDTH="167" HEIGHT="28" BORDER="0"
SRC="img83.png"
ALT="\begin{displaymath}loc\_x_i \leftarrow scatter(glob\_x) \end{displaymath}">
</DIV>
@ -96,7 +96,7 @@ process <SPAN CLASS="MATH"><IMG
ALT="$i$"></SPAN>.
</DD>
<DT><STRONG><SPAN CLASS="MATH"><IMG
WIDTH="53" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="53" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img84.png"
ALT="$scatter$"></SPAN></STRONG></DT>
<DD>is the scatter function.
@ -143,7 +143,7 @@ Data types</CAPTION>
<P>
<BR>
<IMG
WIDTH="538" HEIGHT="16" ALIGN="BOTTOM" BORDER="0"
WIDTH="537" HEIGHT="16" ALIGN="BOTTOM" BORDER="0"
SRC="img85.png"
ALT="\begin{lstlisting}
call psb_scatter(glob_x, loc_x, desc_a, info, root, mold)

@ -122,7 +122,7 @@ Specified as: an integer value <SPAN CLASS="MATH"><IMG
WIDTH="27" HEIGHT="30" ALIGN="MIDDLE" BORDER="0"
SRC="img89.png"
ALT="$0,1$"></SPAN>, default <SPAN CLASS="MATH"><IMG
WIDTH="12" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="12" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img90.png"
ALT="$0$"></SPAN>.
@ -207,7 +207,7 @@ Specified as: a logical value, default: <code>.false.</code>
WIDTH="80" HEIGHT="32" ALIGN="MIDDLE" BORDER="0"
SRC="img92.png"
ALT="$lidx(1:nl)$"></SPAN> to be assigned to the global indices <SPAN CLASS="MATH"><IMG
WIDTH="17" HEIGHT="15" ALIGN="BOTTOM" BORDER="0"
WIDTH="18" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
SRC="img93.png"
ALT="$vl$"></SPAN>.
<BR>

@ -106,7 +106,7 @@ Type: <B>required</B>.
Intent: <B>in</B>.
<BR>
Specified as: an integer array of length <SPAN CLASS="MATH"><IMG
WIDTH="21" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="22" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img105.png"
ALT="$nz$"></SPAN>.
</DD>
@ -120,7 +120,7 @@ Type: <B>required</B>.
Intent: <B>in</B>.
<BR>
Specified as: an integer array of length <SPAN CLASS="MATH"><IMG
WIDTH="21" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="22" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img105.png"
ALT="$nz$"></SPAN>.
</DD>
@ -135,7 +135,7 @@ Type: <B>optional</B>.
Intent: <B>in</B>.
<BR>
Specified as: a logical array of length <SPAN CLASS="MATH"><IMG
WIDTH="21" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="22" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img105.png"
ALT="$nz$"></SPAN>, default <code>.true.</code>.
</DD>
@ -149,7 +149,7 @@ Type: <B>optional</B>.
Intent: <B>in</B>.
<BR>
Specified as: an integer array of length <SPAN CLASS="MATH"><IMG
WIDTH="21" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="22" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img105.png"
ALT="$nz$"></SPAN>.
</DD>
@ -192,7 +192,7 @@ Type: <B>optional</B>.
Intent: <B>out</B>.
<BR>
Specified as: an integer array of length <SPAN CLASS="MATH"><IMG
WIDTH="21" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="22" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img105.png"
ALT="$nz$"></SPAN>.
</DD>
@ -206,7 +206,7 @@ Type: <B>optional</B>.
Intent: <B>out</B>.
<BR>
Specified as: an integer array of length <SPAN CLASS="MATH"><IMG
WIDTH="21" HEIGHT="14" ALIGN="BOTTOM" BORDER="0"
WIDTH="22" HEIGHT="13" ALIGN="BOTTOM" BORDER="0"
SRC="img105.png"
ALT="$nz$"></SPAN>.

Some files were not shown because too many files have changed in this diff Show More

Loading…
Cancel
Save