[petsc-users] slepc NHEP error

[Barry Smith]

  Why do you have the CHKERRQ(ierr); commented out in your code? 

   Because of this you are getting mangled confusing error messages. 

   Put a ierr = in front of all calls and a CHKERRQ(ierr); after each call. 

   Then resend the new error message which will be much clearer.



> On Jun 14, 2017, at 2:58 PM, Kannan, Ramakrishnan <kannanr at ornl.gov> wrote:
> 
> Hello,
>  
> I am running NHEP across 16 MPI processors over 16 nodes in a matrix of global size of 1,000,000x1,000,000 with approximately global 16,000,000 non-zeros. Each node has approximately 1million non-zeros.
>  
> The following is my slepc code for EPS.
>  
> PetscInt nev;
> ierr = EPSCreate(PETSC_COMM_WORLD, &eps);  // CHKERRQ(ierr);
> ierr = EPSSetOperators(eps, A, NULL);  // CHKERRQ(ierr);
> ierr = EPSSetProblemType(eps, EPS_NHEP);  // CHKERRQ(ierr);
> EPSSetWhichEigenpairs(eps, EPS_LARGEST_REAL);
> EPSSetDimensions(eps, 100, PETSC_DEFAULT, PETSC_DEFAULT);
> PRINTROOT("calling epssolve");
> ierr = EPSSolve(eps);  // CHKERRQ(ierr);
> ierr = EPSGetType(eps, &type);  // CHKERRQ(ierr);
> ierr = PetscPrintf(PETSC_COMM_WORLD, " Solution method: %s\n\n", type);
> // CHKERRQ(ierr);
> ierr = EPSGetDimensions(eps, &nev, NULL, NULL);  // CHKERRQ(ierr);
> ierr = PetscPrintf(PETSC_COMM_WORLD,
>                        " Number of requested eigenvalues: %D\n",
>                        nev);  // CHKERRQ(ierr);
>  
> I am getting the following error. Attached is the entire error file for your reference. Please let me know what should I fix in this code.
>  
> 2]PETSC ERROR: Argument out of range
> [2]PETSC ERROR: Argument 2 out of range
> [2]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting.
> [2]PETSC ERROR: Petsc Release Version 3.7.4, Oct, 02, 2016
> [2]PETSC ERROR: ./miniapps on a sandybridge named nid00300 by d3s Wed Jun 14 15:32:00 2017
> [13]PETSC ERROR: #1 BVDotVec_BLAS_Private() line 272 in /autofs/nccs-svm1_home1/ramki/libraries/slepc-3.7.3/src/sys/classes/bv/interface/bvblas.c
> [13]PETSC ERROR: #2 BVDotVec_Svec() line 150 in /autofs/nccs-svm1_home1/ramki/libraries/slepc-3.7.3/src/sys/classes/bv/impls/svec/svec.c
> [13]PETSC ERROR: #3 BVDotVec() line 191 in /autofs/nccs-svm1_home1/ramki/libraries/slepc-3.7.3/src/sys/classes/bv/interface/bvglobal.c
> [13]PETSC ERROR: #4 BVOrthogonalizeCGS1() line 81 in /autofs/nccs-svm1_home1/ramki/libraries/slepc-3.7.3/src/sys/classes/bv/interface/bvorthog.c
> [13]PETSC ERROR: #5 BVOrthogonalizeCGS() line 214 in /autofs/nccs-svm1_home1/ramki/libraries/slepc-3.7.3/src/sys/classes/bv/interface/bvorthog.c
> [13]PETSC ERROR: #6 BVOrthogonalizeColumn() line 371 in /autofs/nccs-svm1_home1/ramki/libraries/slepc-3.7.3/src/sys/classes/bv/interface/bvorthog.c
> [13]PETSC ERROR: #7 EPSBasicArnoldi() line 59 in /autofs/nccs-svm1_home1/ramki/libraries/slepc-3.7.3/src/eps/impls/krylov/epskrylov.c
> [13]PETSC ERROR: #8 EPSSolve_KrylovSchur_Default() line 203 in /autofs/nccs-svm1_home1/ramki/libraries/slepc-3.7.3/src/eps/impls/krylov/krylovschur/krylovschur.c
> [13]PETSC ERROR: #9 EPSSolve() line 101 in /autofs/nccs-svm1_home1/ramki/libraries/slepc-3.7.3/src/eps/interface/epssolve.c
>  
> -- 
> Regards,
> Ramki
>  
> <test.pbs.e613713>