<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Tue, Oct 28, 2014 at 12:05 PM, Tabrez Ali <span dir="ltr"><<a href="mailto:stali@geology.wisc.edu" target="_blank">stali@geology.wisc.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
  
    
  
  <div bgcolor="#FFFFFF" text="#000000">
    <div>Matt<br>
      <br>
      The analytical (Okada) solution is wrong at the fault. PETSc is
      right. It is a strike slip fault and constraints are being
      enforced correctly.</div></div></blockquote><div><br></div><div>If the solution is right and the residual is large, then the function calculating the residual is wrong, since</div><div>the residual is supposed to define the solution.</div><div><br></div><div>   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div bgcolor="#FFFFFF" text="#000000"><div><span class="HOEnZb"><font color="#888888"><br>
      Tabrez</font></span><div><div class="h5"><br>
      <br>
      <br>
      On 10/28/2014 11:57 AM, Matthew Knepley wrote:<br>
    </div></div></div><div><div class="h5">
    <blockquote type="cite">
      <div dir="ltr">
        <div class="gmail_extra">
          <div class="gmail_quote">On Tue, Oct 28, 2014 at 11:56 AM,
            Tabrez Ali <span dir="ltr"><<a href="mailto:stali@geology.wisc.edu" target="_blank">stali@geology.wisc.edu</a>></span>
            wrote:<br>
            <blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Jed<br>
              <br>
              You're right and it is precisely what I am trying to
              understand as well. I am attaching the comparison to the
              analytical solution. As you can see the agreement is
              almost perfect, except at the discontinuity (y=0 in the
              analytical solution).<br>
              <br>
              So why am I getting the correct solution when according to
              true residuals I should not?</blockquote>
            <div><br>
            </div>
            <div>The discontinuty is generating the large residuals.
              Your saddle is not actually enforcing</div>
            <div>the discontinuty, so you get the bulk mostly right by
              solving the (1,1) block, but not</div>
            <div>the slip. At least it looks that way to me.</div>
            <div><br>
            </div>
            <div>   Matt</div>
            <div> </div>
            <blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span><font color="#888888"><br>
                  Tabrez</font></span>
              <div>
                <div><br>
                  <br>
                  On 10/28/2014 11:18 AM, Jed Brown wrote:<br>
                  <blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
                    Tabrez Ali <<a href="mailto:stali@geology.wisc.edu" target="_blank">stali@geology.wisc.edu</a>>
                    writes:<br>
                    <br>
                    <blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
                      Matt and Jed<br>
                      <br>
                      But it does appear to work well in practice (as in
                      I get the correct<br>
                      solution).<br>
                    </blockquote>
                    Your "correct" solution has a residual seven orders
                    of magnitude bigger<br>
                    than the initial guess.  Does that not give you
                    pause?<br>
                  </blockquote>
                  <br>
                </div>
              </div>
            </blockquote>
          </div>
          <br>
          <br clear="all">
          <div><br>
          </div>
          -- <br>
          What most experimenters take for granted before they begin
          their experiments is infinitely more interesting than any
          results to which their experiments lead.<br>
          -- Norbert Wiener
        </div>
      </div>
    </blockquote>
    <br>
  </div></div></div>

</blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener
</div></div>