<div dir="ltr">Humm, not sure what the problem is. I used your test code to debug this branch.<div><br></div><div>You might go to this line:</div><div><br></div><div>[1]PETSC ERROR: #1 KSPComputeEigenvalues() line 119 in<br>/home/filippo/Workspace/petsc-git/ksp-zero-eig/src/ksp/ksp/interface/itfunc.c<br></div><div>







<p class=""><span class="">if (n) PetscValidScalarPointer(r,3);</span></p><div class="gmail_extra">and put a print statement.  n==0 in your case when it fails.</div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Nov 5, 2014 at 2:37 AM, Filippo Leonardi <span dir="ltr"><<a href="mailto:filippo.leonardi@sam.math.ethz.ch" target="_blank">filippo.leonardi@sam.math.ethz.ch</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">@ Jed<br>
<br>
Can't figure out what solver is causing the problem actually, I tried<br>
evrerything: mg_levels, mg_coarse and nothing.<br>
<span class=""><br>
On Tuesday 04 November 2014 12.21:29 Mark Adams wrote:<br>
> On Tue, Nov 4, 2014 at 2:55 AM, Filippo Leonardi <<br>
><br>
> <a href="mailto:filippo.leonardi@sam.math.ethz.ch">filippo.leonardi@sam.math.ethz.ch</a>> wrote:<br>
> > Sorry for the late update.<br>
> ><br>
> > I tested the code with the branch:<br>
> > ksp-zero-eig<br>
><br>
> Did you make PETSc after the checking out the branch?<br>
<br>
</span>Yes.<br>
<span class=""><br>
> Please send us the<br>
> error message and run in debug mode if possible.<br>
<br>
</span>See below.<br>
<span class=""><br>
><br>
> > But the problem seems to still be there, also tried:<br>
> > -ksp_chebyshev_estimate_eigenvalues_random<br>
> > but does not change.<br>
> ><br>
> > Just tell me what test you need and I'll try to do them.<br>
> ><br>
> > Best,<br>
> > Filippo<br>
> ><br>
> > On Wednesday 29 October 2014 10.30:14 Mark Adams wrote:<br>
> > > > Yeah, that should work; an alternative is to defer the estimate until<br>
> > > > the next solve.<br>
> > ><br>
> > > Sounds better.  As I said one of my users has a zero RHS on the first<br>
> ><br>
> > solve.<br>
<br>
</span>( Despite ARCH being arch-linux2-cxx-opt  it's compiled in debug mode )<br>
<br>
mpirun -np 4 ./test1 -mg_levels_ksp_chebyshev_estimate_eigenvalues_random -<br>
ksp_chebyshev_estimate_eigenvalues_random -pc_type mg -pc_mg_levels 5 -<br>
draw_pause -1 -ksp_view<br>
<br>
[0]PETSC ERROR: [1]PETSC ERROR: --------------------- Error Message<br>
--------------------------------------------------------------<br>
[1]PETSC ERROR: Null argument, when expecting valid pointer<br>
[1]PETSC ERROR: Null Pointer: Parameter # 3<br>
[1]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html" target="_blank">http://www.mcs.anl.gov/petsc/documentation/faq.html</a> for<br>
trouble shooting.<br>
[1]PETSC ERROR: Petsc Development GIT revision: v3.5.2-925-g233c798  GIT Date:<br>
2014-11-03 17:44:15 -0700<br>
[1]PETSC ERROR: ./test1 on a arch-linux2-cxx-opt named Besikovitch-III by<br>
filippo Wed Nov  5 08:32:52 2014<br>
[1]PETSC ERROR: Configure options<br>
[1]PETSC ERROR: #1 KSPComputeEigenvalues() line 119 in<br>
/home/filippo/Workspace/petsc-git/ksp-zero-eig/src/ksp/ksp/interface/itfunc.c<br>
[1]PETSC ERROR: #2 KSPChebyshevComputeExtremeEigenvalues_Private() line 328 in<br>
/home/filippo/Workspace/petsc-git/ksp-zero-eig/src/ksp/ksp/impls/cheby/cheby.c<br>
[1]PETSC ERROR: #3 KSPSolve_Chebyshev() line 373 in<br>
/home/filippo/Workspace/petsc-git/ksp-zero-eig/src/ksp/ksp/impls/cheby/cheby.c<br>
[1]PETSC ERROR: #4 KSPSolve() line 460 in /home/filippo/Workspace/petsc-<br>
git/ksp-zero-eig/src/ksp/ksp/interface/itfunc.c<br>
[1]PETSC ERROR: #5 PCMGMCycle_Private() line 19 in<br>
/home/filippo/Workspace/petsc-git/ksp-zero-eig/src/ksp/pc/impls/mg/mg.c<br>
[1]PETSC ERROR: #6 PCApply_MG() line 337 in /home/filippo/Workspace/petsc-<br>
git/ksp-zero-eig/src/ksp/pc/impls/mg/mg.c     </blockquote></div><br></div></div></div>