<div class="gmail_quote">On Fri, Dec 3, 2010 at 18:00, Hong Zhang <span dir="ltr"><<a href="mailto:hzhang@mcs.anl.gov">hzhang@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div id=":43s">We replaced old numerical factorization and MatSolve with faster data structure,<br>
then cleaned up MatPivotCheck() (still using previous algorithm).<br>
Changes, such as<br>
changeset:   15909:fbc1ba94d8d9<br>
user:        hong@hong-zhangs-macbook-air.local<br>
date:        Wed May 12 22:17:30 2010 -0500<br>
files:       src/mat/impls/aij/seq/aijfact.c<br>
description:<br>
replace macros with inline function MatPivotCheck() for inplace lu (aij format)<br>
<br>
could accelerate performance.</div></blockquote></div><br><div>I know the data structures changed, but I see significantly different results (number of iterations, presence of negative pivots) for ICC(0) with SBAIJ(2).  The new results are better, but I had an explanation for poor convergence that involved negative pivots needing to be shifted.  This was the case at that time, but now I don't see negative pivots, and unless I'm misunderstanding what I'm seeing, no shift is being applied now.  If that is really the case, then I'm seeing more than a data structure change.  If a shift is still being applied, I should be able to find where that happens.</div>
<div><br></div><div>I currently see only positive determinants at sbaijfact.c:1122 (current petsc-dev).</div><div><br></div><div>    ierr = Kernel_A_gets_inverse_A_2(diag,shift);CHKERRQ(ierr);</div><div><br></div><div>If I hg update -d 2010-02-01 (changeset aef32a8bffe7), then on the corresponding line (sbaijfact.c:964), I see some determinants on the order of -1.e-9 (only when the problem starts to get hard).</div>
<div><br></div><div>I'm trying to understand what changed.</div><div><br></div><div>Jed</div>