<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta content="text/html;charset=ISO-8859-1" http-equiv="Content-Type">
</head>
<body bgcolor="#ffffff" text="#000000">
<br>
<br>
Walter Landry wrote:<br>
<blockquote
cite="mid:20090507.155710.131509470607391773.walter@geodynamics.org"
type="cite">
<pre wrap="">
The "tolerance" variable in the Stokes_SLE_UzawaSolver controls both
the accuracy of the overall Uzawa iteration as well as the individual
(linear or direct) solves within that solution. </pre>
</blockquote>
If by "tolerance" you mean the "linearTolerance" variable then yes I
understand this.<br>
<blockquote
cite="mid:20090507.155710.131509470607391773.walter@geodynamics.org"
type="cite">
<pre wrap="">You can not specify
the accuracy of the direct solver because there is no iteration. It
directly inverts the matrix.
</pre>
</blockquote>
Yes, I understand this too.<br>
<blockquote
cite="mid:20090507.155710.131509470607391773.walter@geodynamics.org"
type="cite">
<pre wrap="">
Most of the time, setting the tolerance for linear solves to be very
small should give you the same result as a direct solve. </pre>
</blockquote>
I disagree as I have never seen evidence for this behavior. No matter
what linearTolerence is set to I ALWAYS require more nonlinear
iterations with the iterative solver than the direct solver. What I
was asking was whether it would be a good idea to control the tolerance
of that specific part within the Uzawa iteration the iterative solver
solves, separately than "linearTolerance" (which, as I understand
controls the tolerance of the whole Uzawa iteration). Basically, I'm
wondering if there is a (simple) way the code could be modified to
drive specifically the iterative solver towards more accurate
solutions. (i'm beginning to realize this may be a rhetorical question)<br>
<br>
What step within the Uzawa scheme is the direct/iterative solver
solving?<br>
<blockquote
cite="mid:20090507.155710.131509470607391773.walter@geodynamics.org"
type="cite">
<pre wrap="">However, if
the matrix is singular or nearly so, roundoff error will give you the
wrong answer.
Multigrid will probably have the same problems. Better
preconditioning might help.
</pre>
</blockquote>
<blockquote
cite="mid:20090507.155710.131509470607391773.walter@geodynamics.org"
type="cite">
<pre wrap="">
Cheers,
Walter Landry
<a class="moz-txt-link-abbreviated" href="mailto:walter@geodynamics.org">walter@geodynamics.org</a>
_______________________________________________
CIG-LONG mailing list
<a class="moz-txt-link-abbreviated" href="mailto:CIG-LONG@geodynamics.org">CIG-LONG@geodynamics.org</a>
<a class="moz-txt-link-freetext" href="http://geodynamics.org/cgi-bin/mailman/listinfo/cig-long">http://geodynamics.org/cgi-bin/mailman/listinfo/cig-long</a>
</pre>
</blockquote>
<br>
<div class="moz-signature">-- <br>
<meta http-equiv="Content-Type" content="text/html; ">
<meta name="ProgId" content="Word.Document">
<meta name="Generator" content="Microsoft Word 11">
<meta name="Originator" content="Microsoft Word 11">
<link rel="File-List" href="Signature_files/filelist.xml">
<title>Signature</title>
<!--[if gte mso 9]><xml>
<o:DocumentProperties>
<o:Author>Garrett</o:Author>
<o:LastAuthor>Garrett</o:LastAuthor>
<o:Revision>3</o:Revision>
<o:TotalTime>4</o:TotalTime>
<o:Created>2008-12-02T02:21:00Z</o:Created>
<o:LastSaved>2008-12-02T02:21:00Z</o:LastSaved>
<o:Pages>1</o:Pages>
<o:Words>42</o:Words>
<o:Characters>243</o:Characters>
<o:Company>SOEST</o:Company>
<o:Lines>2</o:Lines>
<o:Paragraphs>1</o:Paragraphs>
<o:CharactersWithSpaces>284</o:CharactersWithSpaces>
<o:Version>11.5606</o:Version>
</o:DocumentProperties>
</xml><![endif]--><!--[if gte mso 9]><xml>
<w:WordDocument>
<w:Zoom>75</w:Zoom>
<w:SpellingState>Clean</w:SpellingState>
<w:GrammarState>Clean</w:GrammarState>
<w:DrawingGridHorizontalSpacing>9.35 pt</w:DrawingGridHorizontalSpacing>
<w:DisplayVerticalDrawingGridEvery>2</w:DisplayVerticalDrawingGridEvery>
<w:ValidateAgainstSchemas/>
<w:SaveIfXMLInvalid>false</w:SaveIfXMLInvalid>
<w:IgnoreMixedContent>false</w:Ign
oreMixedContent>
<w:AlwaysShowPlaceholderText>false</w:AlwaysShowPlaceholderText>
<w:Compatibility>
<w:SelectEntireFieldWithStartOrEnd/>
<w:UseWord2002TableStyleRules/>
</w:Compatibility>
<w:BrowserLevel>MicrosoftInternetExplorer4</w:BrowserLevel>
</w:WordDocument>
</xml><![endif]--><!--[if gte mso 9]><xml>
<w:LatentStyles DefLockedState="false" LatentStyleCount="156">
</w:LatentStyles>
</xml><![endif]-->
<style>
<!--
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
        {mso-style-parent:"";
        margin:0in;
        margin-bottom:.0001pt;
        mso-pagination:widow-orphan;
        font-size:12.0pt;
        mso-bidi-font-size:10.0pt;
        font-family:"Times New Roman";
        mso-fareast-font-family:"Times New Roman";
        color:black;
        mso-bidi-font-weight:bold;}
@page Section1
        {size:8.5in 11.0in;
        margin:1.0in 1.0in 1.0in 1.0in;
        mso-header-margin:.5in;
        mso-footer-margin:.5in;
        mso-paper-source:0;}
div.Section1
        {page:Section1;}
-->
</style><!--[if gte mso 10]>
<style>
/* Style Definitions */
table.MsoNormalTable
        {mso-style-name:"Table Normal";
        mso-tstyle-rowband-size:0;
        mso-tstyle-colband-size:0;
        mso-style-noshow:yes;
        mso-style-parent:"";
        mso-padding-alt:0in 5.4pt 0in 5.4pt;
        mso-para-margin:0in;
        mso-para-margin-bottom:.0001pt;
        mso-pagination:widow-orphan;
        font-size:10.0pt;
        font-family:"Times New Roman";
        mso-ansi-language:#0400;
        mso-fareast-language:#0400;
        mso-bidi-language:#0400;}
</style>
<![endif]-->
<div class="Section1">
<p class="MsoNormal"><span style="font-size: 9pt; font-family: Arial;">=====================================<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size: 9pt; font-family: Arial;">Prof.
Garrett Ito<span style=""> </span><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size: 9pt; font-family: Arial;">SOEST,
Geology & Geophysics<span style=""> </span><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size: 9pt; font-family: Arial;">1680
East-West Rd., POST 810<span style=""> </span><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size: 9pt; font-family: Arial;">Univ.
of
Hawaii<span style=""> </span><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size: 9pt; font-family: Arial;">Honolulu,
HI
96822 USA<span style=""> </span><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size: 9pt; font-family: Arial;"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size: 9pt; font-family: Arial;">808-956-9717<span
style=""> </span><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size: 9pt; font-family: Arial;"><a class="moz-txt-link-abbreviated" href="mailto:gito@hawaii.edu">gito@hawaii.edu</a>
<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size: 9pt; font-family: Arial;"><a class="moz-txt-link-abbreviated" href="http://www.soest.hawaii.edu/GG/FACULTY/ITO">www.soest.hawaii.edu/GG/FACULTY/ITO</a><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size: 9pt; font-family: Arial;">======================================<o:p></o:p></span></p>
</div>
</div>
</body>
</html>