fdurastante/amg4psblas
1
0
Fork 0
mirror of https://github.com/sfilippone/amg4psblas.git synced 2026-10-06 22:55:12 +00:00
Code Issues Packages Projects Releases Wiki Activity

Doc updates.

Browse Source
This commit is contained in:
Salvatore Filippone
2021-03-12 13:14:57 +01:00
parent 44acd20726
commit eb6ba73525
21 changed files with 787 additions and 768 deletions
Download Patch File Download Diff File Expand all files Collapse all files
docs/html/userhtmlsu6.html
+35 -36
View File
@@ -10,7 +10,7 @@
<link rel="stylesheet" type="text/css" href="userhtml.css">
</head><body
>
<!--l. 107--><div class="crosslinks"><p class="noindent"><span
<!--l. 114--><div class="crosslinks"><p class="noindent"><span
class="cmr-12">[</span><a
href="#tailuserhtmlsu6.html"><span
class="cmr-12">tail</span></a><span
@@ -22,7 +22,7 @@ class="cmr-12">] </span></p></div>
class="cmr-12">4.1 </span></span> <a
id="x14-130004.1"></a><span
class="cmr-12">Examples</span></h4>
<!--l. 109--><p class="noindent" ><span
<!--l. 116--><p class="noindent" ><span
class="cmr-12">The code reported in Figure</span><span
class="cmr-12">&#x00A0;</span><a
href="#x14-13001r1"><span
@@ -63,46 +63,45 @@ class="cmtt-12">psb_krylov_mod</span> <span
class="cmr-12">must be used by the example</span>
<span
class="cmr-12">program.</span>
<!--l. 119--><p class="indent" > <span
class="cmr-12">The part of the code concerning the reading and assembling of the sparse matrix</span>
<!--l. 126--><p class="indent" > <span
class="cmr-12">The part of the code dealing with reading and assembling the sparse matrix and the</span>
<span
class="cmr-12">and the right-hand side vector, performed through the PSBLAS routines for sparse</span>
class="cmr-12">right-hand side vector and the deallocation of the relevant data structured, performed</span>
<span
class="cmr-12">matrix and vector management, is not reported here for brevity; the statements</span>
class="cmr-12">through the PSBLAS routines for sparse matrix and vector management, is not</span>
<span
class="cmr-12">concerning the deallocation of the PSBLAS data structure are neglected too. The</span>
class="cmr-12">reported here for the sake of conciseness. The complete code can be found in the</span>
<span
class="cmr-12">complete code can be found in the example program file </span><span class="obeylines-h"><span class="verb"><span
class="cmr-12">example program file </span><span class="obeylines-h"><span class="verb"><span
class="cmtt-12">amg_dexample_ml.f90</span></span></span><span
class="cmr-12">, in the</span>
<span
class="cmr-12">directory </span><span class="obeylines-h"><span class="verb"><span
class="cmr-12">, in the directory </span><span class="obeylines-h"><span class="verb"><span
class="cmtt-12">examples/fileread</span></span></span> <span
class="cmr-12">of the AMG4PSBLAS implementation (see</span>
class="cmr-12">of</span>
<span
class="cmr-12">Section</span><span
class="cmr-12">the AMG4PSBLAS implementation (see Section</span><span
class="cmr-12">&#x00A0;</span><a
href="userhtmlsu5.html#x12-110003.5"><span
class="cmr-12">3.5</span><!--tex4ht:ref: sec:ex_and_test --></a><span
class="cmr-12">). A sample test problem along with the relevant input data is available in</span>
<span class="obeylines-h"><span class="verb"><span
class="cmtt-12">examples/fileread/runs</span></span></span><span
class="cmr-12">. For details on the use of the PSBLAS routines, see the</span>
class="cmr-12">). A sample test problem along</span>
<span
class="cmr-12">PSBLAS User&#8217;s Guide</span><span
class="cmr-12">with the relevant input data is available in </span><span class="obeylines-h"><span class="verb"><span
class="cmtt-12">examples/fileread/runs</span></span></span><span
class="cmr-12">. For details on</span>
<span
class="cmr-12">the use of the PSBLAS routines, see the PSBLAS User&#8217;s Guide</span><span
class="cmr-12">&#x00A0;</span><span class="cite"><span
class="cmr-12">[</span><a
href="userhtmlli4.html#XPSBLASGUIDE"><span
class="cmr-12">17</span></a><span
class="cmr-12">]</span></span><span
class="cmr-12">.</span>
<!--l. 131--><p class="indent" > <span
<!--l. 138--><p class="indent" > <span
class="cmr-12">The setup and application of the default multilevel preconditioner for the real single</span>
<span
class="cmr-12">precision and the complex, single and double precision, versions are obtained</span>
<span
class="cmr-12">precision and the complex, single and double precision, versions are obtained</span>
<span
class="cmr-12">with straightforward modifications of the previous example (see Section</span><span
class="cmr-12">&#x00A0;</span><a
@@ -117,7 +116,7 @@ class="cmr-12">.</span>
<!--l. 137--><p class="indent" > <a
<!--l. 144--><p class="indent" > <a
id="x14-13001r1"></a><hr class="float"><div class="float"
>
@@ -125,7 +124,7 @@ class="cmr-12">.</span>
<div class="center"
>
<!--l. 138--><p class="noindent" >
<!--l. 145--><p class="noindent" >
@@ -181,7 +180,7 @@ class="cmr-12">.</span>
&#x00A0;&#x00A0;call&#x00A0;psb_exit(ctxt)
&#x00A0;&#x00A0;stop
</pre>
<!--l. 248--><p class="nopar" >
<!--l. 255--><p class="nopar" >
</div>
@@ -194,7 +193,7 @@ class="content">setup and application of the default multilevel preconditioner (
</div><hr class="endfloat" />
<!--l. 257--><p class="indent" > <span
<!--l. 264--><p class="indent" > <span
class="cmr-12">Different versions of the multilevel preconditioner can be obtained by changing the</span>
<span
class="cmr-12">default values of the preconditioner parameters. The code reported in Figure</span><span
@@ -250,7 +249,7 @@ class="cmr-12">too. </span><span
class="cmbx-12">DA CORREGGERE NEL CODICE ESEMPIO</span>
<span
class="cmbx-12">3</span>
<!--l. 281--><p class="indent" > <span
<!--l. 288--><p class="indent" > <span
class="cmr-12">Finally, Figure</span><span
class="cmr-12">&#x00A0;</span><a
href="#x14-13004r4"><span
@@ -265,7 +264,7 @@ class="cmr-12">nonsymmetric. The corresponding example program is available in t
<span class="obeylines-h"><span class="verb"><span
class="cmtt-12">amg_dexample_1lev.f90</span></span></span><span
class="cmr-12">.</span>
<!--l. 288--><p class="indent" > <span
<!--l. 295--><p class="indent" > <span
class="cmr-12">For all the previous preconditioners, example programs where the sparse matrix</span>
<span
class="cmr-12">and the right-hand side are generated by discretizing a PDE with Dirichlet</span>
@@ -276,7 +275,7 @@ class="cmr-12">.</span>
<!--l. 291--><p class="indent" > <a
<!--l. 298--><p class="indent" > <a
id="x14-13002r2"></a><hr class="float"><div class="float"
>
@@ -284,7 +283,7 @@ class="cmr-12">.</span>
<div class="center"
>
<!--l. 305--><p class="noindent" >
<!--l. 312--><p class="noindent" >
<div class="minipage"><pre class="verbatim" id="verbatim-6">
...&#x00A0;...
!&#x00A0;build&#x00A0;a&#x00A0;V-cycle&#x00A0;preconditioner&#x00A0;with&#x00A0;1&#x00A0;block-Jacobi&#x00A0;sweep&#x00A0;(with
@@ -298,7 +297,7 @@ class="cmr-12">.</span>
&#x00A0;&#x00A0;call&#x00A0;P%smoothers_build(A,desc_A,info)
...&#x00A0;...
</pre>
<!--l. 320--><p class="nopar" ></div></div>
<!--l. 327--><p class="nopar" ></div></div>
<br /><div class="caption"
><span class="id">Listing 2: </span><span
class="content">setup of a multilevel preconditioner based on the default decoupled coarsening</span></div><!--tex4ht:label?: x14-13002r2 -->
@@ -309,7 +308,7 @@ class="content">setup of a multilevel preconditioner based on the default decoup
<!--l. 327--><p class="indent" > <a
<!--l. 334--><p class="indent" > <a
id="x14-13003r3"></a><hr class="float"><div class="float"
>
@@ -317,7 +316,7 @@ class="content">setup of a multilevel preconditioner based on the default decoup
<div class="center"
>
<!--l. 344--><p class="noindent" >
<!--l. 351--><p class="noindent" >
<div class="minipage"><pre class="verbatim" id="verbatim-7">
...&#x00A0;...
!&#x00A0;build&#x00A0;a&#x00A0;W-cycle&#x00A0;preconditioner&#x00A0;with&#x00A0;2&#x00A0;hybrid&#x00A0;Gauss-Seidel&#x00A0;sweeps
@@ -334,7 +333,7 @@ class="content">setup of a multilevel preconditioner based on the default decoup
&#x00A0;&#x00A0;call&#x00A0;P%smoothers_build(A,desc_A,info)
...&#x00A0;...
</pre>
<!--l. 362--><p class="nopar" ></div></div>
<!--l. 369--><p class="nopar" ></div></div>
<br /> <div class="caption"
><span class="id">Listing 3: </span><span
class="content">setup of a multilevel preconditioner based on the coupled coarsening based
@@ -346,7 +345,7 @@ on weighted matching</span></div><!--tex4ht:label?: x14-13003r3 -->
<!--l. 369--><p class="indent" > <a
<!--l. 376--><p class="indent" > <a
id="x14-13004r4"></a><hr class="float"><div class="float"
>
@@ -354,7 +353,7 @@ on weighted matching</span></div><!--tex4ht:label?: x14-13003r3 -->
<div class="center"
>
<!--l. 386--><p class="noindent" >
<!--l. 393--><p class="noindent" >
<div class="minipage"><pre class="verbatim" id="verbatim-8">
...&#x00A0;...
!&#x00A0;set&#x00A0;RAS&#x00A0;with&#x00A0;overlap&#x00A0;2&#x00A0;and&#x00A0;ILU(0)&#x00A0;on&#x00A0;the&#x00A0;local&#x00A0;blocks
@@ -365,7 +364,7 @@ on weighted matching</span></div><!--tex4ht:label?: x14-13003r3 -->
!&#x00A0;solve&#x00A0;Ax=b&#x00A0;with&#x00A0;preconditioned&#x00A0;BiCGSTAB
&#x00A0;&#x00A0;call&#x00A0;psb_krylov(&#8217;BICGSTAB&#8217;,A,P,b,x,tol,desc_A,info)
</pre>
<!--l. 398--><p class="nopar" ></div></div>
<!--l. 405--><p class="nopar" ></div></div>
<br /> <div class="caption"
><span class="id">Listing 4: </span><span
class="content">setup of a one-level Schwarz preconditioner.</span></div><!--tex4ht:label?: x14-13004r4 -->