http://www.cfd-online.com/W/index.php?title=Introduction_and_need&feed=atom&action=historyIntroduction and need - Revision history2016-06-28T06:59:40ZRevision history for this page on the wikiMediaWiki 1.16.5http://www.cfd-online.com/W/index.php?title=Introduction_and_need&diff=4769&oldid=prevTsaad at 22:11, 13 December 20052005-12-13T22:11:27Z<p></p>
<table style="background-color: white; color:black;">
<col class='diff-marker' />
<col class='diff-content' />
<col class='diff-marker' />
<col class='diff-content' />
<tr valign='top'>
<td colspan='2' style="background-color: white; color:black;">← Older revision</td>
<td colspan='2' style="background-color: white; color:black;">Revision as of 22:11, 13 December 2005</td>
</tr><tr><td colspan="2" class="diff-lineno">Line 4:</td>
<td colspan="2" class="diff-lineno">Line 4:</td></tr>
<tr><td class='diff-marker'> </td><td style="background: #eee; color:black; font-size: smaller;"></td><td class='diff-marker'> </td><td style="background: #eee; color:black; font-size: smaller;"></td></tr>
<tr><td class='diff-marker'> </td><td style="background: #eee; color:black; font-size: smaller;"><div>The choice of implementors is an iterative solver. Direct solvers are computationally expensive, and inefficient for large sparse matrices. The most compelling reason for choosing iterative solvers is the inherent non-linearity of the coefficient matrix. As the coefficient matrix is updated at the end of each outer iteration, it is unjustifiable to spend the extra cost for a direct solution that will be iteratively driven to the final solution.</div></td><td class='diff-marker'> </td><td style="background: #eee; color:black; font-size: smaller;"><div>The choice of implementors is an iterative solver. Direct solvers are computationally expensive, and inefficient for large sparse matrices. The most compelling reason for choosing iterative solvers is the inherent non-linearity of the coefficient matrix. As the coefficient matrix is updated at the end of each outer iteration, it is unjustifiable to spend the extra cost for a direct solution that will be iteratively driven to the final solution.</div></td></tr>
<tr><td colspan="2"> </td><td class='diff-marker'>+</td><td style="background: #cfc; color:black; font-size: smaller;"><div><ins style="color: red; font-weight: bold; text-decoration: none;"></ins></div></td></tr>
<tr><td colspan="2"> </td><td class='diff-marker'>+</td><td style="background: #cfc; color:black; font-size: smaller;"><div><ins style="color: red; font-weight: bold; text-decoration: none;">=External links=</ins></div></td></tr>
<tr><td colspan="2"> </td><td class='diff-marker'>+</td><td style="background: #cfc; color:black; font-size: smaller;"><div><ins style="color: red; font-weight: bold; text-decoration: none;">[http://netlib2.cs.utk.edu/linalg/html_templates/Templates.html Templates for the solution of linear systems] - An excellent free book that presents a summary of most iterative methods used in the solution of linear systems of equations. It also includes a pseudocode for each method.</ins></div></td></tr>
</table>Tsaadhttp://www.cfd-online.com/W/index.php?title=Introduction_and_need&diff=4329&oldid=prevTsaad at 07:26, 5 December 20052005-12-05T07:26:04Z<p></p>
<p><b>New page</b></p><div>Almost any method used for the discretization of the Navier-Stokes equations (and most other physical phenomena) yield a linear system of equations that needs to be solved. The need for robust, efficient, and relatively fast linear solvers becomes a crucial issue. While some linear systems exhibit specific behavior (e.g. symmetric, positive definite, tridiagonal...), the general rule excludes this specialty. For example, matrices arising from the finite volume method when ussed on unstructured grids are unsymmetric, sparse, and sometimes singular. Therefore, a huge body of research dedicated for designing robust algorithms has been accumulated over the last two decades.<br />
<br />
Without loss of generality, most linear systems arising form the FV method are sparse (8 or 9 nonzero elements per row at most for a 3D solver). It is thus essential to plan ahead and choose a general solver that is able to handle the sparsity and skewness of the coefficient matrix.<br />
<br />
The choice of implementors is an iterative solver. Direct solvers are computationally expensive, and inefficient for large sparse matrices. The most compelling reason for choosing iterative solvers is the inherent non-linearity of the coefficient matrix. As the coefficient matrix is updated at the end of each outer iteration, it is unjustifiable to spend the extra cost for a direct solution that will be iteratively driven to the final solution.</div>Tsaad