Dear Satish,<div> thanks for your fast reply! I am compiling PETSc with the patch you provided.</div><div><br></div><div>I also have already found a workaround: I've skip the automatic download of <meta charset="utf-8">SuperLU_dist and I've specified manually the version 2.3 instead of the 2.4 (I had a backup of that package). Well, using that version of SuperLU_dist the problem does not appear.</div>
<div><br></div><div>I will keep you updated about the efficacy of the patch.</div><div><br></div><div>Cheers</div><div><br clear="all">--<br>Filippo SPIGA, MSc Computer Science<br><div>~ homepage: <a href="http://tinyurl.com/fspiga" target="_blank">http://tinyurl.com/fspiga</a> ~<br>
<br>«Nobody will drive us out of Cantor's paradise.»<br>
-- David Hilbert<br><br>*****<br>Disclaimer: "Please note this message and any attachment are CONFIDENTIAL an may be privileged or otherwise protected from disclosure. The contents are not to be disclosed to anyone other than the addressee. Unauthorized recipients are requested to preserve this confidentiality and to advise the sender immediately of any error in transmission."</div>
<br>
<br><br><div class="gmail_quote">On Wed, Dec 15, 2010 at 2:46 PM, Satish Balay <span dir="ltr"><<a href="mailto:balay@mcs.anl.gov" target="_blank">balay@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">
sorry missed the attachment.. its attached now.<br>
<font color="#888888"><br>
satish<br>
</font><div><div></div><div><br>
On Wed, 15 Dec 2010, Satish Balay wrote:<br>
<br>
> Can you try the attached patch - and see if this goes away?<br>
><br>
> patch -Np1 < superlu_dist.patch<br>
><br>
> Satish<br>
><br>
> On Wed, 15 Dec 2010, Filippo Spiga wrote:<br>
><br>
> > Dear all,<br>
> > Sorry for the subject of this email, I had no better ideas to introduce<br>
> > the topic (-: I spent a couple of hours yesterday night trying to figure out<br>
> > about a strange behavior of PETSc-3.1-p6 using SUPERLU_dist. This problem<br>
> > arises if I combine the usage of KSPSetOptionsPrefix()/PCSetOptionsPrefix()<br>
> > and "-ksp_view" option. This is the error that I obtain:<br>
> ><br>
> > KSP Object:(lop_)<br>
> > type: gmres<br>
> > GMRES: restart=30, using Classical (unmodified) Gram-Schmidt<br>
> > Orthogonalization with no iterative refinement<br>
> > GMRES: happy breakdown tolerance 1e-30<br>
> > maximum iterations=1000, initial guess is zero<br>
> > tolerances: relative=1e-14, absolute=1e-14, divergence=10000<br>
> > left preconditioning<br>
> > using PRECONDITIONED norm type for convergence test<br>
> > PC Object:(lop_)<br>
> > type: lu<br>
> > LU: out-of-place factorization<br>
> > tolerance for zero pivot 1e-12<br>
> > matrix ordering: nd<br>
> > factor fill ratio given 0, needed 0<br>
> > Factored matrix follows:<br>
> > Matrix Object:<br>
> > type=seqaij, rows=242, cols=242<br>
> > package used to perform factorization: superlu_dist<br>
> > total: nonzeros=0, allocated nonzeros=242<br>
> > SuperLU_DIST run parameters:<br>
> > Process grid nprow 1 x npcol 1<br>
> > Equilibrate matrix TRUE<br>
> > Matrix input mode 0<br>
> > Replace tiny pivots TRUE<br>
> > Use iterative refinement FALSE<br>
> > Processors in row 1 col partition 1<br>
> > Row permutation LargeDiag<br>
> > [0|09:06:54]: unknown: MatFactorInfo_SuperLU_DIST() line 705 in<br>
> > src/mat/impls/aij/mpi/superlu_dist/superlu_dist.c: Unknown column<br>
> > permutation<br>
> ><br>
> ><br>
> > Inside the code, I have these lines:<br>
> ><br>
> > ierr = KSPCreate(PETSC_COMM_WORLD,&ksp); CHKERRQ(ierr);<br>
> > ierr = KSPSetOptionsPrefix(ksp,"lop_"); CHKERRQ(ierr);<br>
> > ierr = KSPGetPC(ksp,&precond); CHKERRQ(ierr);<br>
> > ierr = PCSetOptionsPrefix(precond,"lop_"); CHKERRQ(ierr);<br>
> > ierr = KSPSetFromOptions(ksp); CHKERRQ(ierr);<br>
> > ierr = KSPSolve(ksp,b,*u); CHKERRQ(ierr);<br>
> ><br>
> > so all the options that I pass to PETSc (using "PetscOptionsSetValue()" or<br>
> > the command line) must have the "lop_" prefix. These are my options:<br>
> ><br>
> > -lop_ksp_rtol 1.e-10 -lop_ksp_atol 1.e-10 -lop_ksp_monitor<br>
> > -lop_ksp_converged_reason -lop_ksp_type gmres -lop_ksp_monitor_true_residual<br>
> > -lop_pc_type lu -lop_pc_factor_mat_solver_package superlu_dist -lop_ksp_view<br>
> ><br>
> ><br>
> > If I remove "-lop_ksp_view", the error that I reported disappears (and also<br>
> > disappear the final summary<br>
> > If I remove -lop_pc_type lu -lop_pc_factor_mat_solver_package superlu_dist"<br>
> > and I use HYPRE ("-lop_pc_type hypre") and I keep the option "-lop_ksp_view"<br>
> > as valid option, the problem does not appear.<br>
> > If I remove "-lop_ksp_view" from command line but after "KSPSolve()" I add<br>
> > "KSPView()", the problem persists.<br>
> > If I suppress the prefix ("-ksp_rtol 1.e-10 -ksp_atol ... -ksp_view"), the<br>
> > problem does not appear.<br>
> ><br>
> > PETSc-3.1-p3, PETSc-3.1-p4 and PETSc-3.1-p5 are not affected of this problem<br>
> > (other machines that I use have these versions installed, the code runs<br>
> > without problems). In the meanwhile I am recompiling PETSc-3.1-p6 on another<br>
> > machine.<br>
> ><br>
> > Do you have any suggestion to skip/solve this problem? (if it is a real<br>
> > problem, maybe I did something wrong...)<br>
> ><br>
> > Many thanks in advance,<br>
> > Cheers<br>
> ><br>
> > --<br>
> > Filippo SPIGA, MSc Computer Science<br>
> > ~ homepage: <a href="http://tinyurl.com/fspiga" target="_blank">http://tinyurl.com/fspiga</a> ~<br>
> ><br>
> > «Nobody will drive us out of Cantor's paradise.»<br>
> > -- David Hilbert<br>
> ><br>
> > *****<br>
> > Disclaimer: "Please note this message and any attachment are CONFIDENTIAL an<br>
> > may be privileged or otherwise protected from disclosure. The contents are<br>
> > not to be disclosed to anyone other than the addressee. Unauthorized<br>
> > recipients are requested to preserve this confidentiality and to advise the<br>
> > sender immediately of any error in transmission."<br>
> ><br>
><br>
</div></div></blockquote></div><br></div>