<div dir="ltr">It does not look like you are in fact using the branch with this fix. The branch code looks like:<div><br></div><div>if (n) <span style="font-family:arial,sans-serif;font-size:12.7272720336914px"> </span><span style="font-family:arial,sans-serif;font-size:12.7272720336914px">PetscValidScalarPointer(r,3);</span></div><div><span style="font-family:arial,sans-serif;font-size:12.7272720336914px"><br></span></div><div><span style="font-family:arial,sans-serif;font-size:12.7272720336914px">*and*</span></div><div><span style="font-family:arial,sans-serif;font-size:12.7272720336914px"><br></span></div><div>if (n) <span style="font-size:12.7272720336914px;font-family:arial,sans-serif"> </span><span style="font-size:12.7272720336914px;font-family:arial,sans-serif">PetscValidScalarPointer(r,4);</span><span style="font-family:arial,sans-serif;font-size:12.7272720336914px"><br></span></div><div><span style="font-size:12.7272720336914px;font-family:arial,sans-serif"><br></span></div><div><span style="font-size:12.7272720336914px;font-family:arial,sans-serif">So if you now put the check for argument 4 inside of your if statement then you will have reimplemented my branch, or at least most of it, and I think it might work or at least not give these errors. I don't think we ever fixed the issue that Cheby does will not have an eigen estimate now and so you will probably see terrible performance. I'm not sure if came up with a fix for this.</span></div><div><span style="font-size:12.7272720336914px;font-family:arial,sans-serif"><br></span></div><div><span style="font-family:arial,sans-serif;font-size:12.7272720336914px"><br></span></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Nov 10, 2014 at 9:18 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:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Super late update, sorry.<br>
<br>
Not sure<br>
if I understood you correctly:<br>
if(n) {<br>
PetscValidScalarPointer(r,3);<br>
ierr = PetscPrintf(PETSC_COMM_WORLD, " ---> n = %d nonzero if.\n",<br>
n);CHKERRQ(ierr);<br>
} else {<br>
ierr = PetscPrintf(PETSC_COMM_WORLD, " ---> n = 0 else.\n");CHKERRQ(ierr);<br>
}<br>
<br>
I get (btw this is also the case with a single process):<br>
./test1 -pc_type mg -pc_mg_levels 3<br>
<br>
---> n = 0 else.<br>
[0]PETSC ERROR: --------------------- Error Message<br>
--------------------------------------------------------------<br>
[0]PETSC ERROR: Null argument, when expecting valid pointer<br>
[0]PETSC ERROR: Null Pointer: Parameter # 4<br>
[0]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>
[0]PETSC ERROR: Petsc Development GIT revision: v3.5.2-925-g233c798 GIT Date:<br>
<span class="">2014-11-03 17:44:15 -0700<br>
</span>[0]PETSC ERROR: ./test1 on a arch-linux2-cxx-opt named Besikovitch-III by<br>
filippo Mon Nov 10 15:15:47 2014<br>
[0]PETSC ERROR: Configure options PETSC_ARCH=arch-linux2-cxx-opt<br>
[0]PETSC ERROR: #1 KSPComputeEigenvalues() line 125 in<br>
/home/filippo/Workspace/petsc-git/ksp-zero-eig/src/ksp/ksp/interface/itfunc.c<br>
[0]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>
[0]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>
[0]PETSC ERROR: #4 KSPSolve() line 465 in /home/filippo/Workspace/petsc-<br>
git/ksp-zero-eig/src/ksp/ksp/interface/itfunc.c<br>
[0]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>
[0]PETSC ERROR: #6 PCApply_MG() line 337 in /home/filippo/Workspace/petsc-<br>
git/ksp-zero-eig/src/ksp/pc/impls/mg/mg.c<br>
[0]PETSC ERROR: #7 PCApply() line 440 in /home/filippo/Workspace/petsc-<br>
git/ksp-zero-eig/src/ksp/pc/interface/precon.c<br>
[0]PETSC ERROR: #8 KSP_PCApply() line 230 in /home/filippo/Workspace/petsc-<br>
git/ksp-zero-eig/include/petsc-private/kspimpl.h<br>
[0]PETSC ERROR: #9 KSPInitialResidual() line 63 in<br>
/home/filippo/Workspace/petsc-git/ksp-zero-eig/src/ksp/ksp/interface/itres.c<br>
[0]PETSC ERROR: #10 KSPSolve_GMRES() line 234 in<br>
/home/filippo/Workspace/petsc-git/ksp-zero-eig/src/ksp/ksp/impls/gmres/gmres.c<br>
[0]PETSC ERROR: #11 KSPSolve() line 465 in /home/filippo/Workspace/petsc-<br>
git/ksp-zero-eig/src/ksp/ksp/interface/itfunc.c<br>
[0]PETSC ERROR: #12 main() line 149 in<br>
/home/filippo/Workspace/C++/MLMCTest/test1.cpp<br>
[0]PETSC ERROR: ----------------End of Error Message -------send entire error<br>
message to petsc-maint@mcs.anl.gov----------<br>
<div class="HOEnZb"><div class="h5"><br>
On Thursday 06 November 2014 11.27:59 you wrote:<br>
> if (n) PetscValidScalarPointer(r,3);<br>
> and put a print statement. n==0 in your case when it fails.</div></div></blockquote></div><br></div>