GCC Code Coverage Report

Directory:	./
File:	src/eps/tests/test2.c
Date:	2025-11-19 04:19:03

	Exec	Total	Coverage
Lines:	54	54	100.0%
Functions:	1	1	100.0%
Branches:	182	268	67.9%

  
      Line
      Branch
      Exec
      Source
    
      /*
    
         - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
    
         SLEPc - Scalable Library for Eigenvalue Problem Computations
    
         Copyright (c) 2002-, Universitat Politecnica de Valencia, Spain
    
         This file is part of SLEPc.
    
         SLEPc is distributed under a 2-clause BSD license (see LICENSE).
    
         - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
    
      */
    
      static char help[] = "Tests multiple calls to EPSSolve with the same matrix.\n\n";
    
      #include <slepceps.h>
    
      212
      int main(int argc,char **argv)
    
      {
    
      212
        Mat            A;           /* problem matrix */
    
      212
        EPS            eps;         /* eigenproblem solver context */
    
      212
        ST             st;
    
      212
        PetscReal      tol=PetscMax(1000*PETSC_MACHINE_EPSILON,1e-9);
    
      212
        PetscInt       n=30,i,Istart,Iend;
    
      212
        PetscBool      flg;
    
      212
        EPSLanczosReorthogType reorth;
    
        1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.

      212
        PetscFunctionBeginUser;
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      212
        PetscCall(SlepcInitialize(&argc,&argv,NULL,help));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      212
        PetscCall(PetscOptionsGetInt(NULL,NULL,"-n",&n,NULL));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      212
        PetscCall(PetscPrintf(PETSC_COMM_WORLD,"\n1-D Laplacian Eigenproblem, n=%" PetscInt_FMT "\n\n",n));
    
        /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
    
           Compute the operator matrix that defines the eigensystem, Ax=kx
    
           - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      212
        PetscCall(MatCreate(PETSC_COMM_WORLD,&A));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      212
        PetscCall(MatSetSizes(A,PETSC_DECIDE,PETSC_DECIDE,n,n));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      212
        PetscCall(MatSetFromOptions(A));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      212
        PetscCall(MatGetOwnershipRange(A,&Istart,&Iend));
    
        2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.

      5752
        for (i=Istart;i<Iend;i++) {
    
        6/8✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 8 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.

      5540
          if (i>0) PetscCall(MatSetValue(A,i,i-1,-1.0,INSERT_VALUES));
    
        6/8✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 8 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.

      5540
          if (i<n-1) PetscCall(MatSetValue(A,i,i+1,-1.0,INSERT_VALUES));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      5540
          PetscCall(MatSetValue(A,i,i,2.0,INSERT_VALUES));
    
        }
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      212
        PetscCall(MatAssemblyBegin(A,MAT_FINAL_ASSEMBLY));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      212
        PetscCall(MatAssemblyEnd(A,MAT_FINAL_ASSEMBLY));
    
        /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
    
                              Create the eigensolver
    
           - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      212
        PetscCall(EPSCreate(PETSC_COMM_WORLD,&eps));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      212
        PetscCall(EPSSetOperators(eps,A,NULL));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      212
        PetscCall(EPSSetProblemType(eps,EPS_HEP));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      212
        PetscCall(EPSSetTolerances(eps,tol,PETSC_CURRENT));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      212
        PetscCall(EPSSetFromOptions(eps));
    
        /* illustrate how to extract parameters from specific solver types */
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      212
        PetscCall(PetscObjectTypeCompare((PetscObject)eps,EPSLANCZOS,&flg));
    
        2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.

      212
        if (flg) {
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      50
          PetscCall(EPSLanczosGetReorthog(eps,&reorth));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      50
          PetscCall(PetscPrintf(PETSC_COMM_WORLD,"Reorthogonalization type used in Lanczos: %s\n",EPSLanczosReorthogTypes[reorth]));
    
        }
    
        /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
    
                          Solve for largest eigenvalues
    
           - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      212
        PetscCall(EPSSetWhichEigenpairs(eps,EPS_LARGEST_REAL));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      212
        PetscCall(EPSSolve(eps));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      212
        PetscCall(PetscPrintf(PETSC_COMM_WORLD," - - - Largest eigenvalues - - -\n"));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      212
        PetscCall(EPSErrorView(eps,EPS_ERROR_RELATIVE,NULL));
    
        /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
    
                          Solve for smallest eigenvalues
    
           - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      212
        PetscCall(EPSSetWhichEigenpairs(eps,EPS_SMALLEST_REAL));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      212
        PetscCall(EPSSolve(eps));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      212
        PetscCall(PetscPrintf(PETSC_COMM_WORLD," - - - Smallest eigenvalues - - -\n"));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      212
        PetscCall(EPSErrorView(eps,EPS_ERROR_RELATIVE,NULL));
    
        /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
    
                          Solve for interior eigenvalues (target=2.1)
    
           - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      212
        PetscCall(EPSSetWhichEigenpairs(eps,EPS_TARGET_MAGNITUDE));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      212
        PetscCall(EPSSetTarget(eps,2.1));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      212
        PetscCall(PetscObjectTypeCompare((PetscObject)eps,EPSLANCZOS,&flg));
    
        2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.

      212
        if (flg) {
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      50
          PetscCall(EPSGetST(eps,&st));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      50
          PetscCall(STSetType(st,STSINVERT));
    
        } else {
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      162
          PetscCall(PetscObjectTypeCompare((PetscObject)eps,EPSKRYLOVSCHUR,&flg));
    
        6/8✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 8 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.

      162
          if (!flg) PetscCall(PetscObjectTypeCompare((PetscObject)eps,EPSARNOLDI,&flg));
    
        6/8✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 8 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.

      162
          if (flg) PetscCall(EPSSetExtraction(eps,EPS_HARMONIC));
    
        }
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      212
        PetscCall(EPSSolve(eps));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      212
        PetscCall(PetscPrintf(PETSC_COMM_WORLD," - - - Interior eigenvalues - - -\n"));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      212
        PetscCall(EPSErrorView(eps,EPS_ERROR_RELATIVE,NULL));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      212
        PetscCall(EPSDestroy(&eps));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      212
        PetscCall(MatDestroy(&A));
    
        3/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 2 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.

      212
        PetscCall(SlepcFinalize());
    
        return 0;
    
      }
    
      /*TEST
    
         testset:
    
            args: -eps_nev 4
    
            requires: !single
    
            output_file: output/test2_1.out
    
            test:
    
               suffix: 1
    
               args: -eps_type {{arnoldi gd jd lapack}}
    
            test:
    
               suffix: 1_gd2
    
               args: -eps_type gd -eps_gd_double_expansion
    
               timeoutfactor: 2
    
            test:
    
               suffix: 1_krylovschur
    
               args: -eps_type krylovschur -eps_krylovschur_locking {{0 1}}
    
            test:
    
               suffix: 1_scalapack
    
               requires: scalapack !__float128
    
               args: -eps_type scalapack
    
            test:
    
               suffix: 1_elpa
    
               requires: elpa
    
               args: -eps_type elpa
    
            test:
    
               suffix: 1_elemental
    
               requires: elemental
    
               args: -eps_type elemental
    
         testset:
    
            args: -eps_type lanczos -eps_nev 4
    
            requires: !single
    
            filter: grep -v "Lanczos"
    
            output_file: output/test2_1.out
    
            test:
    
               suffix: 2
    
               args: -eps_lanczos_reorthog {{local full periodic partial}}
    
            test:
    
               suffix: 2_selective
    
               args: -eps_lanczos_reorthog selective
    
               requires: !defined(PETSCTEST_VALGRIND)
    
         testset:
    
            args: -n 32 -eps_nev 4
    
            requires: !single
    
            output_file: output/test2_3.out
    
            test:
    
               nsize: 2
    
               suffix: 3
    
               args: -eps_type {{krylovschur lapack}}
    
            test:
    
               nsize: 2
    
               suffix: 3_gd
    
               args: -eps_type gd -eps_gd_krylov_start
    
               timeoutfactor: 2
    
            test:
    
               suffix: 3_jd
    
               args: -eps_type jd -eps_jd_krylov_start -eps_ncv 18
    
         testset:
    
            args: -eps_nev 4 -eps_type {{krylovschur arnoldi gd jd}}
    
            requires: !single
    
            output_file: output/test2_1.out
    
            test:
    
               suffix: 4_cuda
    
               requires: cuda
    
               args: -mat_type aijcusparse
    
            test:
    
               suffix: 4_hip
    
               requires: hip
    
               args: -mat_type aijhipsparse
    
      TEST*/

Line	Branch	Exec	Source
1			/*
2			- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
3			SLEPc - Scalable Library for Eigenvalue Problem Computations
4			Copyright (c) 2002-, Universitat Politecnica de Valencia, Spain
5
6			This file is part of SLEPc.
7			SLEPc is distributed under a 2-clause BSD license (see LICENSE).
8			- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
9			*/
10
11			static char help[] = "Tests multiple calls to EPSSolve with the same matrix.\n\n";
12
13			#include <slepceps.h>
14
15		212	int main(int argc,char **argv)
16			{
17		212	Mat A; /* problem matrix */
18		212	EPS eps; /* eigenproblem solver context */
19		212	ST st;
20		212	PetscReal tol=PetscMax(1000*PETSC_MACHINE_EPSILON,1e-9);
21		212	PetscInt n=30,i,Istart,Iend;
22		212	PetscBool flg;
23		212	EPSLanczosReorthogType reorth;
24
25	1/2 ✓ Branch 0 taken 2 times. ✗ Branch 1 not taken.	212	PetscFunctionBeginUser;
26	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	212	PetscCall(SlepcInitialize(&argc,&argv,NULL,help));
27
28	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	212	PetscCall(PetscOptionsGetInt(NULL,NULL,"-n",&n,NULL));
29	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	212	PetscCall(PetscPrintf(PETSC_COMM_WORLD,"\n1-D Laplacian Eigenproblem, n=%" PetscInt_FMT "\n\n",n));
30
31			/* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
32			Compute the operator matrix that defines the eigensystem, Ax=kx
33			- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
34
35	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	212	PetscCall(MatCreate(PETSC_COMM_WORLD,&A));
36	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	212	PetscCall(MatSetSizes(A,PETSC_DECIDE,PETSC_DECIDE,n,n));
37	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	212	PetscCall(MatSetFromOptions(A));
38
39	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	212	PetscCall(MatGetOwnershipRange(A,&Istart,&Iend));
40	2/2 ✓ Branch 0 taken 10 times. ✓ Branch 1 taken 10 times.	5752	for (i=Istart;i<Iend;i++) {
41	6/8 ✓ Branch 0 taken 10 times. ✓ Branch 1 taken 10 times. ✓ Branch 2 taken 2 times. ✓ Branch 3 taken 8 times. ✓ Branch 4 taken 2 times. ✗ Branch 5 not taken. ✗ Branch 6 not taken. ✓ Branch 7 taken 2 times.	5540	if (i>0) PetscCall(MatSetValue(A,i,i-1,-1.0,INSERT_VALUES));
42	6/8 ✓ Branch 0 taken 10 times. ✓ Branch 1 taken 10 times. ✓ Branch 2 taken 2 times. ✓ Branch 3 taken 8 times. ✓ Branch 4 taken 2 times. ✗ Branch 5 not taken. ✗ Branch 6 not taken. ✓ Branch 7 taken 2 times.	5540	if (i<n-1) PetscCall(MatSetValue(A,i,i+1,-1.0,INSERT_VALUES));
43	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	5540	PetscCall(MatSetValue(A,i,i,2.0,INSERT_VALUES));
44			}
45	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	212	PetscCall(MatAssemblyBegin(A,MAT_FINAL_ASSEMBLY));
46	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	212	PetscCall(MatAssemblyEnd(A,MAT_FINAL_ASSEMBLY));
47
48			/* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
49			Create the eigensolver
50			- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
51	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	212	PetscCall(EPSCreate(PETSC_COMM_WORLD,&eps));
52	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	212	PetscCall(EPSSetOperators(eps,A,NULL));
53	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	212	PetscCall(EPSSetProblemType(eps,EPS_HEP));
54	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	212	PetscCall(EPSSetTolerances(eps,tol,PETSC_CURRENT));
55	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	212	PetscCall(EPSSetFromOptions(eps));
56
57			/* illustrate how to extract parameters from specific solver types */
58	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	212	PetscCall(PetscObjectTypeCompare((PetscObject)eps,EPSLANCZOS,&flg));
59	2/2 ✓ Branch 0 taken 10 times. ✓ Branch 1 taken 10 times.	212	if (flg) {
60	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	50	PetscCall(EPSLanczosGetReorthog(eps,&reorth));
61	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	50	PetscCall(PetscPrintf(PETSC_COMM_WORLD,"Reorthogonalization type used in Lanczos: %s\n",EPSLanczosReorthogTypes[reorth]));
62			}
63
64			/* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
65			Solve for largest eigenvalues
66			- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
67	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	212	PetscCall(EPSSetWhichEigenpairs(eps,EPS_LARGEST_REAL));
68	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	212	PetscCall(EPSSolve(eps));
69	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	212	PetscCall(PetscPrintf(PETSC_COMM_WORLD," - - - Largest eigenvalues - - -\n"));
70	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	212	PetscCall(EPSErrorView(eps,EPS_ERROR_RELATIVE,NULL));
71
72			/* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
73			Solve for smallest eigenvalues
74			- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
75	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	212	PetscCall(EPSSetWhichEigenpairs(eps,EPS_SMALLEST_REAL));
76	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	212	PetscCall(EPSSolve(eps));
77	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	212	PetscCall(PetscPrintf(PETSC_COMM_WORLD," - - - Smallest eigenvalues - - -\n"));
78	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	212	PetscCall(EPSErrorView(eps,EPS_ERROR_RELATIVE,NULL));
79
80			/* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
81			Solve for interior eigenvalues (target=2.1)
82			- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
83	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	212	PetscCall(EPSSetWhichEigenpairs(eps,EPS_TARGET_MAGNITUDE));
84	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	212	PetscCall(EPSSetTarget(eps,2.1));
85	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	212	PetscCall(PetscObjectTypeCompare((PetscObject)eps,EPSLANCZOS,&flg));
86	2/2 ✓ Branch 0 taken 10 times. ✓ Branch 1 taken 10 times.	212	if (flg) {
87	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	50	PetscCall(EPSGetST(eps,&st));
88	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	50	PetscCall(STSetType(st,STSINVERT));
89			} else {
90	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	162	PetscCall(PetscObjectTypeCompare((PetscObject)eps,EPSKRYLOVSCHUR,&flg));
91	6/8 ✓ Branch 0 taken 10 times. ✓ Branch 1 taken 10 times. ✓ Branch 2 taken 2 times. ✓ Branch 3 taken 8 times. ✓ Branch 4 taken 2 times. ✗ Branch 5 not taken. ✗ Branch 6 not taken. ✓ Branch 7 taken 2 times.	162	if (!flg) PetscCall(PetscObjectTypeCompare((PetscObject)eps,EPSARNOLDI,&flg));
92	6/8 ✓ Branch 0 taken 10 times. ✓ Branch 1 taken 10 times. ✓ Branch 2 taken 2 times. ✓ Branch 3 taken 8 times. ✓ Branch 4 taken 2 times. ✗ Branch 5 not taken. ✗ Branch 6 not taken. ✓ Branch 7 taken 2 times.	162	if (flg) PetscCall(EPSSetExtraction(eps,EPS_HARMONIC));
93			}
94	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	212	PetscCall(EPSSolve(eps));
95	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	212	PetscCall(PetscPrintf(PETSC_COMM_WORLD," - - - Interior eigenvalues - - -\n"));
96	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	212	PetscCall(EPSErrorView(eps,EPS_ERROR_RELATIVE,NULL));
97
98	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	212	PetscCall(EPSDestroy(&eps));
99	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	212	PetscCall(MatDestroy(&A));
100	3/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 2 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✗ Branch 5 not taken.	212	PetscCall(SlepcFinalize());
101			return 0;
102			}
103
104			/*TEST
105
106			testset:
107			args: -eps_nev 4
108			requires: !single
109			output_file: output/test2_1.out
110			test:
111			suffix: 1
112			args: -eps_type {{arnoldi gd jd lapack}}
113			test:
114			suffix: 1_gd2
115			args: -eps_type gd -eps_gd_double_expansion
116			timeoutfactor: 2
117			test:
118			suffix: 1_krylovschur
119			args: -eps_type krylovschur -eps_krylovschur_locking {{0 1}}
120			test:
121			suffix: 1_scalapack
122			requires: scalapack !__float128
123			args: -eps_type scalapack
124			test:
125			suffix: 1_elpa
126			requires: elpa
127			args: -eps_type elpa
128			test:
129			suffix: 1_elemental
130			requires: elemental
131			args: -eps_type elemental
132
133			testset:
134			args: -eps_type lanczos -eps_nev 4
135			requires: !single
136			filter: grep -v "Lanczos"
137			output_file: output/test2_1.out
138			test:
139			suffix: 2
140			args: -eps_lanczos_reorthog {{local full periodic partial}}
141			test:
142			suffix: 2_selective
143			args: -eps_lanczos_reorthog selective
144			requires: !defined(PETSCTEST_VALGRIND)
145
146			testset:
147			args: -n 32 -eps_nev 4
148			requires: !single
149			output_file: output/test2_3.out
150			test:
151			nsize: 2
152			suffix: 3
153			args: -eps_type {{krylovschur lapack}}
154			test:
155			nsize: 2
156			suffix: 3_gd
157			args: -eps_type gd -eps_gd_krylov_start
158			timeoutfactor: 2
159			test:
160			suffix: 3_jd
161			args: -eps_type jd -eps_jd_krylov_start -eps_ncv 18
162
163			testset:
164			args: -eps_nev 4 -eps_type {{krylovschur arnoldi gd jd}}
165			requires: !single
166			output_file: output/test2_1.out
167			test:
168			suffix: 4_cuda
169			requires: cuda
170			args: -mat_type aijcusparse
171			test:
172			suffix: 4_hip
173			requires: hip
174			args: -mat_type aijhipsparse
175
176			TEST*/
177