[petsc-users] Problem Using EPSGetEigenpair

Jose E. Roman jroman at dsic.upv.es
Fri Feb 8 08:05:14 CST 2019


You are printing the variable as an integer, and it is not an integer, it's a PetscScalar. Use %g
Jose


> El 8 feb 2019, a las 15:02, Eda Oktay via petsc-users <petsc-users at mcs.anl.gov> escribió:
> 
> Hello,
> 
> I am trying to test EPSGetEigenpair function on ex11 in EPS. 
> 
> First, I executed ex11 and get 0.097887 as eigenvalue (k).
> 
> Then under the code 
> 
> ierr = PetscPrintf(PETSC_COMM_WORLD," Number of requested eigenvalues: %D\n",nev);CHKERRQ(ierr);
> 
> I added the following code:
> 
>   ierr = EPSGetEigenpair(eps,0,&kr,NULL,NULL,NULL);  
>   ierr = PetscPrintf(PETSC_COMM_WORLD," Requested eigenvalue: %D\n",kr);CHKERRQ(ierr);
> 
> After execution, I expected to find 0.097887 but the requested eigenvalue was found as -876957631.
> 
> Why EPSGetEigenpair couldn't give 0.097887, isn't  0.097887 the eigenvalue that EPSSetWhichEigenpairs(eps,EPS_SMALLEST_REAL) calculated?
> 
> Best regards,
> 
> Eda



More information about the petsc-users mailing list