GCC Code Coverage Report

Directory:	./
File:	src/eps/tests/test25.c
Date:	2025-10-04 04:19:13

	Exec	Total	Coverage
Lines:	53	53	100.0%
Functions:	1	1	100.0%
Branches:	160	244	65.6%

  
      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[] = "Solves a GNHEP problem with contour integral.\n\n"
    
        "Based on ex7.\n"
    
        "The command line options are:\n"
    
        "  -f1 <filename> -f2 <filename>, PETSc binary files containing A and B.\n\n";
    
      #include <slepceps.h>
    
      10
      int main(int argc,char **argv)
    
      {
    
      10
        EPS               eps;
    
      10
        RG                rg;
    
      10
        Mat               A,B;
    
      10
        PetscBool         flg;
    
      10
        EPSCISSExtraction extr;
    
      10
        EPSCISSQuadRule   quad;
    
      10
        char              filename[PETSC_MAX_PATH_LEN];
    
      10
        PetscViewer       viewer;
    
        1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.

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

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

      10
        PetscCall(PetscPrintf(PETSC_COMM_WORLD,"\nGNHEP problem with contour integral\n\n"));
    
        /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
    
              Load the matrices that define the eigensystem, Ax=kBx
    
           - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      10
        PetscCall(PetscOptionsGetString(NULL,NULL,"-f1",filename,sizeof(filename),&flg));
    
        1/4✗ Branch 0 not taken.
✓ Branch 1 taken 8 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.

      10
        PetscCheck(flg,PETSC_COMM_WORLD,PETSC_ERR_USER_INPUT,"Must indicate a file name for matrix A with the -f1 option");
    
      #if defined(PETSC_USE_COMPLEX)
    
        4/6✓ Branch 0 taken 1 times.
✓ Branch 1 taken 2 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.

      3
        PetscCall(PetscPrintf(PETSC_COMM_WORLD," Reading COMPLEX matrices from binary files...\n"));
    
      #else
    
        4/6✓ Branch 0 taken 1 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.

      7
        PetscCall(PetscPrintf(PETSC_COMM_WORLD," Reading REAL matrices from binary files...\n"));
    
      #endif
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      10
        PetscCall(PetscViewerBinaryOpen(PETSC_COMM_WORLD,filename,FILE_MODE_READ,&viewer));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

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

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

      10
        PetscCall(MatLoad(A,viewer));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      10
        PetscCall(PetscViewerDestroy(&viewer));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      10
        PetscCall(PetscOptionsGetString(NULL,NULL,"-f2",filename,sizeof(filename),&flg));
    
        2/2✓ Branch 0 taken 5 times.
✓ Branch 1 taken 3 times.

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

      7
          PetscCall(PetscViewerBinaryOpen(PETSC_COMM_WORLD,filename,FILE_MODE_READ,&viewer));
    
        4/6✓ Branch 0 taken 1 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.

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

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

      7
          PetscCall(MatLoad(B,viewer));
    
        4/6✓ Branch 0 taken 1 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.

      7
          PetscCall(PetscViewerDestroy(&viewer));
    
        } else {
    
        4/6✓ Branch 0 taken 1 times.
✓ Branch 1 taken 2 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.

      3
          PetscCall(PetscPrintf(PETSC_COMM_WORLD," Matrix B was not provided, setting B=I\n\n"));
    
      3
          B = NULL;
    
        }
    
        /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
    
                      Create the eigensolver and solve the problem
    
           - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

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

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

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

      10
        PetscCall(EPSSetTolerances(eps,1e-9,PETSC_CURRENT));
    
        /* set CISS solver with various options */
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      10
        PetscCall(EPSSetType(eps,EPSCISS));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      10
        PetscCall(EPSCISSSetExtraction(eps,EPS_CISS_EXTRACTION_HANKEL));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      10
        PetscCall(EPSCISSSetQuadRule(eps,EPS_CISS_QUADRULE_CHEBYSHEV));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      10
        PetscCall(EPSCISSSetUseST(eps,PETSC_TRUE));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      10
        PetscCall(EPSGetRG(eps,&rg));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      10
        PetscCall(RGSetType(rg,RGINTERVAL));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      10
        PetscCall(RGIntervalSetEndpoints(rg,-3000.0,0.0,0.0,0.0));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

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

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

      10
        PetscCall(PetscObjectTypeCompare((PetscObject)eps,EPSCISS,&flg));
    
        1/2✓ Branch 0 taken 8 times.
✗ Branch 1 not taken.

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

      10
          PetscCall(EPSCISSGetExtraction(eps,&extr));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      10
          PetscCall(EPSCISSGetQuadRule(eps,&quad));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      10
          PetscCall(PetscPrintf(PETSC_COMM_WORLD," Solved with CISS using %s extraction with %s quadrature rule\n\n",EPSCISSExtractions[extr],EPSCISSQuadRules[quad]));
    
        }
    
        /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
    
                          Display solution and clean up
    
           - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

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

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

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

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

      10
        PetscCall(SlepcFinalize());
    
        return 0;
    
      }
    
      /*TEST
    
         testset:
    
            args: -f1 ${SLEPC_DIR}/share/slepc/datafiles/matrices/bfw62a.petsc -f2 ${SLEPC_DIR}/share/slepc/datafiles/matrices/bfw62b.petsc
    
            output_file: output/test25_1.out
    
            test:
    
               suffix: 1
    
               requires: double !complex !defined(PETSC_USE_64BIT_INDICES)
    
            test:
    
               suffix: 1_cuda
    
               args: -mat_type aijcusparse -st_pc_factor_mat_solver_type cusparse
    
               requires: cuda double !complex !defined(PETSC_USE_64BIT_INDICES)
    
            test:
    
               suffix: 1_hip
    
               args: -mat_type aijhipsparse -st_pc_factor_mat_solver_type hipsparse
    
               requires: hip double !complex !defined(PETSC_USE_64BIT_INDICES)
    
         testset:
    
            args: -f1 ${DATAFILESPATH}/matrices/complex/qc324.petsc -rg_type ellipse -rg_ellipse_center 1-0.09i -rg_ellipse_radius 0.15 -rg_ellipse_vscale 0.1
    
            output_file: output/test25_2.out
    
            test:
    
               suffix: 2
    
               requires: double complex datafilespath !defined(PETSC_USE_64BIT_INDICES)
    
            test:
    
               suffix: 2_cuda
    
               args: -mat_type aijcusparse -st_pc_factor_mat_solver_type cusparse
    
               requires: cuda double complex datafilespath !defined(PETSC_USE_64BIT_INDICES)
    
            test:
    
               suffix: 2_hip
    
               args: -mat_type aijhipsparse -st_pc_factor_mat_solver_type hipsparse
    
               requires: hip double complex datafilespath !defined(PETSC_USE_64BIT_INDICES)
    
      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[] = "Solves a GNHEP problem with contour integral.\n\n"
12			"Based on ex7.\n"
13			"The command line options are:\n"
14			" -f1 <filename> -f2 <filename>, PETSc binary files containing A and B.\n\n";
15
16			#include <slepceps.h>
17
18		10	int main(int argc,char **argv)
19			{
20		10	EPS eps;
21		10	RG rg;
22		10	Mat A,B;
23		10	PetscBool flg;
24		10	EPSCISSExtraction extr;
25		10	EPSCISSQuadRule quad;
26		10	char filename[PETSC_MAX_PATH_LEN];
27		10	PetscViewer viewer;
28
29	1/2 ✓ Branch 0 taken 2 times. ✗ Branch 1 not taken.	10	PetscFunctionBeginUser;
30	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 6 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	10	PetscCall(SlepcInitialize(&argc,&argv,NULL,help));
31	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 6 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	10	PetscCall(PetscPrintf(PETSC_COMM_WORLD,"\nGNHEP problem with contour integral\n\n"));
32
33			/* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
34			Load the matrices that define the eigensystem, Ax=kBx
35			- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
36
37	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 6 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	10	PetscCall(PetscOptionsGetString(NULL,NULL,"-f1",filename,sizeof(filename),&flg));
38	1/4 ✗ Branch 0 not taken. ✓ Branch 1 taken 8 times. ✗ Branch 2 not taken. ✗ Branch 3 not taken.	10	PetscCheck(flg,PETSC_COMM_WORLD,PETSC_ERR_USER_INPUT,"Must indicate a file name for matrix A with the -f1 option");
39
40			#if defined(PETSC_USE_COMPLEX)
41	4/6 ✓ Branch 0 taken 1 times. ✓ Branch 1 taken 2 times. ✓ Branch 2 taken 1 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 1 times.	3	PetscCall(PetscPrintf(PETSC_COMM_WORLD," Reading COMPLEX matrices from binary files...\n"));
42			#else
43	4/6 ✓ Branch 0 taken 1 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 1 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 1 times.	7	PetscCall(PetscPrintf(PETSC_COMM_WORLD," Reading REAL matrices from binary files...\n"));
44			#endif
45	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 6 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	10	PetscCall(PetscViewerBinaryOpen(PETSC_COMM_WORLD,filename,FILE_MODE_READ,&viewer));
46	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 6 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	10	PetscCall(MatCreate(PETSC_COMM_WORLD,&A));
47	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 6 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	10	PetscCall(MatSetFromOptions(A));
48	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 6 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	10	PetscCall(MatLoad(A,viewer));
49	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 6 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	10	PetscCall(PetscViewerDestroy(&viewer));
50
51	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 6 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	10	PetscCall(PetscOptionsGetString(NULL,NULL,"-f2",filename,sizeof(filename),&flg));
52	2/2 ✓ Branch 0 taken 5 times. ✓ Branch 1 taken 3 times.	10	if (flg) {
53	4/6 ✓ Branch 0 taken 1 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 1 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 1 times.	7	PetscCall(PetscViewerBinaryOpen(PETSC_COMM_WORLD,filename,FILE_MODE_READ,&viewer));
54	4/6 ✓ Branch 0 taken 1 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 1 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 1 times.	7	PetscCall(MatCreate(PETSC_COMM_WORLD,&B));
55	4/6 ✓ Branch 0 taken 1 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 1 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 1 times.	7	PetscCall(MatSetFromOptions(B));
56	4/6 ✓ Branch 0 taken 1 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 1 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 1 times.	7	PetscCall(MatLoad(B,viewer));
57	4/6 ✓ Branch 0 taken 1 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 1 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 1 times.	7	PetscCall(PetscViewerDestroy(&viewer));
58			} else {
59	4/6 ✓ Branch 0 taken 1 times. ✓ Branch 1 taken 2 times. ✓ Branch 2 taken 1 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 1 times.	3	PetscCall(PetscPrintf(PETSC_COMM_WORLD," Matrix B was not provided, setting B=I\n\n"));
60		3	B = NULL;
61			}
62
63			/* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
64			Create the eigensolver and solve the problem
65			- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
66
67	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 6 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	10	PetscCall(EPSCreate(PETSC_COMM_WORLD,&eps));
68	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 6 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	10	PetscCall(EPSSetOperators(eps,A,B));
69	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 6 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	10	PetscCall(EPSSetProblemType(eps,EPS_GNHEP));
70	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 6 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	10	PetscCall(EPSSetTolerances(eps,1e-9,PETSC_CURRENT));
71
72			/* set CISS solver with various options */
73	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 6 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	10	PetscCall(EPSSetType(eps,EPSCISS));
74	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 6 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	10	PetscCall(EPSCISSSetExtraction(eps,EPS_CISS_EXTRACTION_HANKEL));
75	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 6 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	10	PetscCall(EPSCISSSetQuadRule(eps,EPS_CISS_QUADRULE_CHEBYSHEV));
76	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 6 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	10	PetscCall(EPSCISSSetUseST(eps,PETSC_TRUE));
77	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 6 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	10	PetscCall(EPSGetRG(eps,&rg));
78	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 6 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	10	PetscCall(RGSetType(rg,RGINTERVAL));
79	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 6 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	10	PetscCall(RGIntervalSetEndpoints(rg,-3000.0,0.0,0.0,0.0));
80
81	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 6 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	10	PetscCall(EPSSetFromOptions(eps));
82
83	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 6 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	10	PetscCall(EPSSolve(eps));
84	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 6 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	10	PetscCall(PetscObjectTypeCompare((PetscObject)eps,EPSCISS,&flg));
85	1/2 ✓ Branch 0 taken 8 times. ✗ Branch 1 not taken.	10	if (flg) {
86	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 6 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	10	PetscCall(EPSCISSGetExtraction(eps,&extr));
87	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 6 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	10	PetscCall(EPSCISSGetQuadRule(eps,&quad));
88	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 6 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	10	PetscCall(PetscPrintf(PETSC_COMM_WORLD," Solved with CISS using %s extraction with %s quadrature rule\n\n",EPSCISSExtractions[extr],EPSCISSQuadRules[quad]));
89			}
90
91			/* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
92			Display solution and clean up
93			- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
94
95	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 6 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	10	PetscCall(EPSErrorView(eps,EPS_ERROR_BACKWARD,NULL));
96	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 6 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	10	PetscCall(EPSDestroy(&eps));
97	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 6 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	10	PetscCall(MatDestroy(&A));
98	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 6 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	10	PetscCall(MatDestroy(&B));
99	3/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 1 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✗ Branch 5 not taken.	10	PetscCall(SlepcFinalize());
100			return 0;
101			}
102
103			/*TEST
104
105			testset:
106			args: -f1 ${SLEPC_DIR}/share/slepc/datafiles/matrices/bfw62a.petsc -f2 ${SLEPC_DIR}/share/slepc/datafiles/matrices/bfw62b.petsc
107			output_file: output/test25_1.out
108			test:
109			suffix: 1
110			requires: double !complex !defined(PETSC_USE_64BIT_INDICES)
111			test:
112			suffix: 1_cuda
113			args: -mat_type aijcusparse -st_pc_factor_mat_solver_type cusparse
114			requires: cuda double !complex !defined(PETSC_USE_64BIT_INDICES)
115			test:
116			suffix: 1_hip
117			args: -mat_type aijhipsparse -st_pc_factor_mat_solver_type hipsparse
118			requires: hip double !complex !defined(PETSC_USE_64BIT_INDICES)
119
120			testset:
121			args: -f1 ${DATAFILESPATH}/matrices/complex/qc324.petsc -rg_type ellipse -rg_ellipse_center 1-0.09i -rg_ellipse_radius 0.15 -rg_ellipse_vscale 0.1
122			output_file: output/test25_2.out
123			test:
124			suffix: 2
125			requires: double complex datafilespath !defined(PETSC_USE_64BIT_INDICES)
126			test:
127			suffix: 2_cuda
128			args: -mat_type aijcusparse -st_pc_factor_mat_solver_type cusparse
129			requires: cuda double complex datafilespath !defined(PETSC_USE_64BIT_INDICES)
130			test:
131			suffix: 2_hip
132			args: -mat_type aijhipsparse -st_pc_factor_mat_solver_type hipsparse
133			requires: hip double complex datafilespath !defined(PETSC_USE_64BIT_INDICES)
134
135			TEST*/
136