<div dir="ltr">Where is your matrix? It might be ending up with a very bad pivot. If the problem can be reproduced, it should be reported to the SuperLU_DIST developers to fix. (Note that we do not see this with other matrices.) You can also try MUMPS.</div>
<div class="gmail_extra"><br><br><div class="gmail_quote">On Sun, Dec 23, 2012 at 6:48 PM, Sanjay Govindjee <span dir="ltr"><<a href="mailto:s_g@berkeley.edu" target="_blank">s_g@berkeley.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">
    I wanted to use SuperLU Dist to perform a direct solve but seem to
    be encountering<br>
    a problem.  I was wonder if this is a know issue and if there is a
    solution for it.<br>
    <br>
    The problem is easily observed using ex6.c in
    src/ksp/ksp/examples/tests.<br>
    <br>
    Out of the box: make runex6 produces a residual error of O(1e-11),
    all is well.<br>
    <br>
    I then changed the run to run on two processors and add the flag<br>
    -pc_factor_mat_solver_package spooles  this produces a residual
    error of O(1e-11), all is still well.<br>
    <br>
    I then switch over to -pc_factor_mat_solver_package superlu_dist and
    the<br>
    residual error comes back as 22.6637!  Something seems very wrong.<br>
    <br>
    My build is perfectly vanilla:<br>
    <br>
    export PETSC_DIR=/Users/sg/petsc-3.3-p5/<br>
    export PETSC_ARCH=intel<br>
    <br>
    ./configure --with-cc=icc --with-fc=ifort  \
-download-{spooles,parmetis,superlu_dist,prometheus,mpich,ml,hypre,metis}<br>
    <br>
    make PETSC_DIR=/Users/sg/petsc-3.3-p5/ PETSC_ARCH=intel all<br>
    make PETSC_DIR=/Users/sg/petsc-3.3-p5/ PETSC_ARCH=intel test<span class="HOEnZb"><font color="#888888"><br>
    <br>
    -sanjay<br>
    
  </font></span></div>

</blockquote></div><br></div>