GCC Code Coverage Report

Directory:	./
File:	src/eps/tests/test13.c
Date:	2026-01-13 04:10:43

	Exec	Total	Coverage
Lines:	54	55	98.2%
Functions:	3	3	100.0%
Branches:	156	250	62.4%

  
      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[] = "Test EPSSetArbitrarySelection.\n\n";
    
      #include <slepceps.h>
    
      141956
      PetscErrorCode MyArbitrarySelection(PetscScalar eigr,PetscScalar eigi,Vec xr,Vec xi,PetscScalar *rr,PetscScalar *ri,void *ctx)
    
      {
    
      141956
        Vec             xref = *(Vec*)ctx;
    
        1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.

      141956
        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.

      141956
        PetscCall(VecDot(xr,xref,rr));
    
      141956
        *rr = PetscAbsScalar(*rr);
    
        1/2✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.

      141956
        if (ri) *ri = 0.0;
    
        5/12✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 2 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
✗ Branch 10 not taken.
✗ Branch 11 not taken.

      27196
        PetscFunctionReturn(PETSC_SUCCESS);
    
      }
    
      80
      PetscErrorCode MyArbitrarySelectionDestroy(void **ctx)
    
      {
    
      80
        Vec             xref = *(Vec*)*ctx;
    
        1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.

      80
        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.

      80
        PetscCall(VecDestroy(&xref));
    
        5/12✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 2 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
✗ Branch 10 not taken.
✗ Branch 11 not taken.

      16
        PetscFunctionReturn(PETSC_SUCCESS);
    
      }
    
      80
      int main(int argc,char **argv)
    
      {
    
      80
        Mat            A;           /* problem matrices */
    
      80
        EPS            eps;         /* eigenproblem solver context */
    
      80
        PetscScalar    seigr,seigi;
    
      80
        PetscReal      tol=1000*PETSC_MACHINE_EPSILON;
    
      80
        Vec            sxr,sxi;
    
      80
        PetscInt       n=30,i,Istart,Iend,nconv;
    
        1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.

      80
        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.

      80
        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.

      80
        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.

      80
        PetscCall(PetscPrintf(PETSC_COMM_WORLD,"\nTridiagonal with zero diagonal, n=%" PetscInt_FMT "\n\n",n));
    
        /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
    
                 Create matrix tridiag([-1 0 -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.

      80
        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.

      80
        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.

      80
        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.

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

      2480
        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.

      2400
          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.

      2400
          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.

      80
        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.

      80
        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.

      80
        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.

      80
        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.

      80
        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.

      80
        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.

      80
        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.

      80
        PetscCall(EPSSetFromOptions(eps));
    
        /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
    
                      Solve eigenproblem and store some solution
    
           - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
    
        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.

      80
        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.

      80
        PetscCall(EPSGetConverged(eps,&nconv));
    
        1/2✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.

      80
        if (nconv>0) {
    
        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.

      80
          PetscCall(MatCreateVecs(A,&sxr,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.

      80
          PetscCall(MatCreateVecs(A,&sxi,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.

      80
          PetscCall(EPSGetEigenpair(eps,0,&seigr,&seigi,sxr,sxi));
    
        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.

      80
          PetscCall(VecDestroy(&sxi));
    
        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.

      80
          PetscCall(EPSErrorView(eps,EPS_ERROR_RELATIVE,NULL));
    
          /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
    
                       Solve eigenproblem using an arbitrary selection
    
             - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
    
        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.

      80
          PetscCall(EPSSetArbitrarySelection(eps,MyArbitrarySelection,&sxr));
    
        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.

      80
          PetscCall(EPSSetArbitrarySelectionContextDestroy(eps,MyArbitrarySelectionDestroy));
    
        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.

      80
          PetscCall(EPSSetWhichEigenpairs(eps,EPS_LARGEST_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.

      80
          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.

      80
          PetscCall(EPSErrorView(eps,EPS_ERROR_RELATIVE,NULL));
    
      ✗
        } else PetscCall(PetscPrintf(PETSC_COMM_WORLD,"Problem: no eigenpairs converged.\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.

      80
        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.

      80
        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.

      80
        PetscCall(SlepcFinalize());
    
        return 0;
    
      }
    
      /*TEST
    
         testset:
    
            args: -eps_max_it 5000 -st_pc_type jacobi
    
            output_file: output/test13_1.out
    
            filter: sed -e "s/-1.98975/-1.98974/"
    
            test:
    
               suffix: 1
    
               args: -eps_type {{krylovschur gd}}
    
            test:
    
               suffix: 1_jd
    
               args: -eps_type jd -st_ksp_type gmres
    
            test:
    
               suffix: 1_gd2
    
               args: -eps_type gd -eps_gd_double_expansion
    
            test:
    
               suffix: 2
    
               args: -eps_non_hermitian -eps_type {{krylovschur gd jd}}
    
            test:
    
               suffix: 2_gd2
    
               args: -eps_non_hermitian -eps_type gd -eps_gd_double_expansion
    
               timeoutfactor: 2
    
      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[] = "Test EPSSetArbitrarySelection.\n\n";
12
13			#include <slepceps.h>
14
15		141956	PetscErrorCode MyArbitrarySelection(PetscScalar eigr,PetscScalar eigi,Vec xr,Vec xi,PetscScalar rr,PetscScalar ri,void *ctx)
16			{
17		141956	Vec xref = (Vec)ctx;
18
19	1/2 ✓ Branch 0 taken 2 times. ✗ Branch 1 not taken.	141956	PetscFunctionBeginUser;
20	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.	141956	PetscCall(VecDot(xr,xref,rr));
21		141956	rr = PetscAbsScalar(rr);
22	1/2 ✓ Branch 0 taken 10 times. ✗ Branch 1 not taken.	141956	if (ri) *ri = 0.0;
23	5/12 ✓ Branch 0 taken 2 times. ✗ Branch 1 not taken. ✗ Branch 2 not taken. ✓ Branch 3 taken 2 times. ✓ Branch 4 taken 2 times. ✗ Branch 5 not taken. ✓ Branch 6 taken 2 times. ✗ Branch 7 not taken. ✗ Branch 8 not taken. ✓ Branch 9 taken 2 times. ✗ Branch 10 not taken. ✗ Branch 11 not taken.	27196	PetscFunctionReturn(PETSC_SUCCESS);
24			}
25
26		80	PetscErrorCode MyArbitrarySelectionDestroy(void **ctx)
27			{
28		80	Vec xref = (Vec)*ctx;
29
30	1/2 ✓ Branch 0 taken 2 times. ✗ Branch 1 not taken.	80	PetscFunctionBeginUser;
31	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.	80	PetscCall(VecDestroy(&xref));
32	5/12 ✓ Branch 0 taken 2 times. ✗ Branch 1 not taken. ✗ Branch 2 not taken. ✓ Branch 3 taken 2 times. ✓ Branch 4 taken 2 times. ✗ Branch 5 not taken. ✓ Branch 6 taken 2 times. ✗ Branch 7 not taken. ✗ Branch 8 not taken. ✓ Branch 9 taken 2 times. ✗ Branch 10 not taken. ✗ Branch 11 not taken.	16	PetscFunctionReturn(PETSC_SUCCESS);
33			}
34
35		80	int main(int argc,char **argv)
36			{
37		80	Mat A; /* problem matrices */
38		80	EPS eps; /* eigenproblem solver context */
39		80	PetscScalar seigr,seigi;
40		80	PetscReal tol=1000*PETSC_MACHINE_EPSILON;
41		80	Vec sxr,sxi;
42		80	PetscInt n=30,i,Istart,Iend,nconv;
43
44	1/2 ✓ Branch 0 taken 2 times. ✗ Branch 1 not taken.	80	PetscFunctionBeginUser;
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.	80	PetscCall(SlepcInitialize(&argc,&argv,NULL,help));
46
47	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.	80	PetscCall(PetscOptionsGetInt(NULL,NULL,"-n",&n,NULL));
48	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.	80	PetscCall(PetscPrintf(PETSC_COMM_WORLD,"\nTridiagonal with zero diagonal, n=%" PetscInt_FMT "\n\n",n));
49
50			/* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
51			Create matrix tridiag([-1 0 -1])
52			- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
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.	80	PetscCall(MatCreate(PETSC_COMM_WORLD,&A));
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.	80	PetscCall(MatSetSizes(A,PETSC_DECIDE,PETSC_DECIDE,n,n));
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.	80	PetscCall(MatSetFromOptions(A));
56
57	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.	80	PetscCall(MatGetOwnershipRange(A,&Istart,&Iend));
58	2/2 ✓ Branch 0 taken 10 times. ✓ Branch 1 taken 10 times.	2480	for (i=Istart;i<Iend;i++) {
59	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.	2400	if (i>0) PetscCall(MatSetValue(A,i,i-1,-1.0,INSERT_VALUES));
60	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.	2400	if (i<n-1) PetscCall(MatSetValue(A,i,i+1,-1.0,INSERT_VALUES));
61			}
62	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.	80	PetscCall(MatAssemblyBegin(A,MAT_FINAL_ASSEMBLY));
63	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.	80	PetscCall(MatAssemblyEnd(A,MAT_FINAL_ASSEMBLY));
64
65			/* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
66			Create the eigensolver
67			- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
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.	80	PetscCall(EPSCreate(PETSC_COMM_WORLD,&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.	80	PetscCall(EPSSetProblemType(eps,EPS_HEP));
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.	80	PetscCall(EPSSetTolerances(eps,tol,PETSC_CURRENT));
71	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.	80	PetscCall(EPSSetOperators(eps,A,NULL));
72	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.	80	PetscCall(EPSSetWhichEigenpairs(eps,EPS_SMALLEST_REAL));
73	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.	80	PetscCall(EPSSetFromOptions(eps));
74
75			/* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
76			Solve eigenproblem and store some solution
77			- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
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.	80	PetscCall(EPSSolve(eps));
79	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.	80	PetscCall(EPSGetConverged(eps,&nconv));
80	1/2 ✓ Branch 0 taken 10 times. ✗ Branch 1 not taken.	80	if (nconv>0) {
81	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.	80	PetscCall(MatCreateVecs(A,&sxr,NULL));
82	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.	80	PetscCall(MatCreateVecs(A,&sxi,NULL));
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.	80	PetscCall(EPSGetEigenpair(eps,0,&seigr,&seigi,sxr,sxi));
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.	80	PetscCall(VecDestroy(&sxi));
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.	80	PetscCall(EPSErrorView(eps,EPS_ERROR_RELATIVE,NULL));
86
87			/* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
88			Solve eigenproblem using an arbitrary selection
89			- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
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.	80	PetscCall(EPSSetArbitrarySelection(eps,MyArbitrarySelection,&sxr));
91	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.	80	PetscCall(EPSSetArbitrarySelectionContextDestroy(eps,MyArbitrarySelectionDestroy));
92	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.	80	PetscCall(EPSSetWhichEigenpairs(eps,EPS_LARGEST_MAGNITUDE));
93	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.	80	PetscCall(EPSSolve(eps));
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.	80	PetscCall(EPSErrorView(eps,EPS_ERROR_RELATIVE,NULL));
95		✗	} else PetscCall(PetscPrintf(PETSC_COMM_WORLD,"Problem: no eigenpairs converged.\n"));
96
97	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.	80	PetscCall(EPSDestroy(&eps));
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.	80	PetscCall(MatDestroy(&A));
99	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.	80	PetscCall(SlepcFinalize());
100			return 0;
101			}
102
103			/*TEST
104
105			testset:
106			args: -eps_max_it 5000 -st_pc_type jacobi
107			output_file: output/test13_1.out
108			filter: sed -e "s/-1.98975/-1.98974/"
109			test:
110			suffix: 1
111			args: -eps_type {{krylovschur gd}}
112			test:
113			suffix: 1_jd
114			args: -eps_type jd -st_ksp_type gmres
115			test:
116			suffix: 1_gd2
117			args: -eps_type gd -eps_gd_double_expansion
118			test:
119			suffix: 2
120			args: -eps_non_hermitian -eps_type {{krylovschur gd jd}}
121			test:
122			suffix: 2_gd2
123			args: -eps_non_hermitian -eps_type gd -eps_gd_double_expansion
124			timeoutfactor: 2
125
126			TEST*/
127

Function (Line)	Call count	Block coverage
MyArbitrarySelection (line 15)	called 141956 times	83.0%
MyArbitrarySelectionDestroy (line 26)	called 80 times	75.0%
main (line 35)	called 80 times	81.0%