GCC Code Coverage Report

Directory:	./
File:	src/eps/tests/test34.c
Date:	2025-12-10 04:20:18

	Exec	Total	Coverage
Lines:	51	51	100.0%
Functions:	1	1	100.0%
Branches:	158	236	66.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[] = "Test interface to external package PRIMME.\n\n"
    
        "This is based on ex3.c. The command line options are:\n"
    
        "  -n <n>, where <n> = number of grid subdivisions in x dimension.\n"
    
        "  -m <m>, where <m> = number of grid subdivisions in y dimension.\n\n";
    
      #include <slepceps.h>
    
      24
      int main(int argc,char **argv)
    
      {
    
      24
        Mat             A;           /* matrix */
    
      24
        EPS             eps;         /* eigenproblem solver context */
    
      24
        ST              st;          /* spectral transformation context */
    
      24
        KSP             ksp;
    
      24
        PC              pc;
    
      24
        PetscInt        N,n=35,m,Istart,Iend,II,i,j,bs;
    
      24
        PetscBool       flag;
    
      24
        EPSPRIMMEMethod meth;
    
        1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.

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

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

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

      24
        PetscCall(PetscOptionsGetInt(NULL,NULL,"-m",&m,&flag));
    
        1/2✓ Branch 0 taken 6 times.
✗ Branch 1 not taken.

      24
        if (!flag) m=n;
    
      24
        N = n*m;
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      24
        PetscCall(PetscPrintf(PETSC_COMM_WORLD,"\nStandard eigenproblem with PRIMME, N=%" PetscInt_FMT " (%" PetscInt_FMT "x%" PetscInt_FMT " grid)\n\n",N,n,m));
    
        /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
    
               Compute the matrices that define the eigensystem, Ax=kBx
    
           - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

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

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

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

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

      22074
        for (II=Istart;II<Iend;II++) {
    
      22050
          i = II/n; j = II-i*n;
    
        6/8✓ Branch 0 taken 6 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 4 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.

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

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

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

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

      22050
          PetscCall(MatSetValue(A,II,II,4.0,INSERT_VALUES));
    
        }
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

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

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

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

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

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

      24
        PetscCall(EPSSetType(eps,EPSPRIMME));
    
        /*
    
           Set several options
    
        */
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

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

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

      24
        PetscCall(STSetType(st,STPRECOND));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      24
        PetscCall(STGetKSP(st,&ksp));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      24
        PetscCall(KSPGetPC(ksp,&pc));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      24
        PetscCall(KSPSetType(ksp,KSPPREONLY));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      24
        PetscCall(PCSetType(pc,PCICC));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

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

      24
        PetscCall(EPSPRIMMESetMethod(eps,EPS_PRIMME_GD_OLSEN_PLUSK));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      24
        PetscCall(EPSSetFromOptions(eps));
    
        /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
    
                       Compute eigenvalues and display info
    
           - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

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

      24
        PetscCall(EPSPRIMMEGetBlockSize(eps,&bs));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      24
        PetscCall(EPSPRIMMEGetMethod(eps,&meth));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      24
        PetscCall(PetscPrintf(PETSC_COMM_WORLD," PRIMME: using block size %" PetscInt_FMT ", method %s\n",bs,EPSPRIMMEMethods[meth]));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

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

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

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

      24
        PetscCall(SlepcFinalize());
    
        return 0;
    
      }
    
      /*TEST
    
         build:
    
            requires: primme
    
         testset:
    
            args: -eps_nev 4
    
            requires: primme
    
            output_file: output/test34_1.out
    
            test:
    
               suffix: 1
    
            test:
    
               suffix: 2
    
               args: -st_pc_type bjacobi -eps_target 0.01 -eps_target_real -eps_refined
    
               nsize: 2
    
            test:
    
               suffix: 3
    
               args: -eps_smallest_magnitude -eps_harmonic
    
      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 interface to external package PRIMME.\n\n"
12			"This is based on ex3.c. The command line options are:\n"
13			" -n <n>, where <n> = number of grid subdivisions in x dimension.\n"
14			" -m <m>, where <m> = number of grid subdivisions in y dimension.\n\n";
15
16			#include <slepceps.h>
17
18		24	int main(int argc,char **argv)
19			{
20		24	Mat A; /* matrix */
21		24	EPS eps; /* eigenproblem solver context */
22		24	ST st; /* spectral transformation context */
23		24	KSP ksp;
24		24	PC pc;
25		24	PetscInt N,n=35,m,Istart,Iend,II,i,j,bs;
26		24	PetscBool flag;
27		24	EPSPRIMMEMethod meth;
28
29	1/2 ✓ Branch 0 taken 2 times. ✗ Branch 1 not taken.	24	PetscFunctionBeginUser;
30	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	24	PetscCall(SlepcInitialize(&argc,&argv,NULL,help));
31
32	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	24	PetscCall(PetscOptionsGetInt(NULL,NULL,"-n",&n,NULL));
33	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	24	PetscCall(PetscOptionsGetInt(NULL,NULL,"-m",&m,&flag));
34	1/2 ✓ Branch 0 taken 6 times. ✗ Branch 1 not taken.	24	if (!flag) m=n;
35		24	N = n*m;
36	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	24	PetscCall(PetscPrintf(PETSC_COMM_WORLD,"\nStandard eigenproblem with PRIMME, N=%" PetscInt_FMT " (%" PetscInt_FMT "x%" PetscInt_FMT " grid)\n\n",N,n,m));
37
38			/* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
39			Compute the matrices that define the eigensystem, Ax=kBx
40			- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
41
42	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	24	PetscCall(MatCreate(PETSC_COMM_WORLD,&A));
43	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	24	PetscCall(MatSetSizes(A,PETSC_DECIDE,PETSC_DECIDE,N,N));
44	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	24	PetscCall(MatSetFromOptions(A));
45
46	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	24	PetscCall(MatGetOwnershipRange(A,&Istart,&Iend));
47	2/2 ✓ Branch 0 taken 6 times. ✓ Branch 1 taken 6 times.	22074	for (II=Istart;II<Iend;II++) {
48		22050	i = II/n; j = II-i*n;
49	6/8 ✓ Branch 0 taken 6 times. ✓ Branch 1 taken 6 times. ✓ Branch 2 taken 2 times. ✓ Branch 3 taken 4 times. ✓ Branch 4 taken 2 times. ✗ Branch 5 not taken. ✗ Branch 6 not taken. ✓ Branch 7 taken 2 times.	22050	if (i>0) PetscCall(MatSetValue(A,II,II-n,-1.0,INSERT_VALUES));
50	6/8 ✓ Branch 0 taken 6 times. ✓ Branch 1 taken 6 times. ✓ Branch 2 taken 2 times. ✓ Branch 3 taken 4 times. ✓ Branch 4 taken 2 times. ✗ Branch 5 not taken. ✗ Branch 6 not taken. ✓ Branch 7 taken 2 times.	22050	if (i<m-1) PetscCall(MatSetValue(A,II,II+n,-1.0,INSERT_VALUES));
51	6/8 ✓ Branch 0 taken 6 times. ✓ Branch 1 taken 6 times. ✓ Branch 2 taken 2 times. ✓ Branch 3 taken 4 times. ✓ Branch 4 taken 2 times. ✗ Branch 5 not taken. ✗ Branch 6 not taken. ✓ Branch 7 taken 2 times.	22050	if (j>0) PetscCall(MatSetValue(A,II,II-1,-1.0,INSERT_VALUES));
52	6/8 ✓ Branch 0 taken 6 times. ✓ Branch 1 taken 6 times. ✓ Branch 2 taken 2 times. ✓ Branch 3 taken 4 times. ✓ Branch 4 taken 2 times. ✗ Branch 5 not taken. ✗ Branch 6 not taken. ✓ Branch 7 taken 2 times.	22050	if (j<n-1) PetscCall(MatSetValue(A,II,II+1,-1.0,INSERT_VALUES));
53	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	22050	PetscCall(MatSetValue(A,II,II,4.0,INSERT_VALUES));
54			}
55
56	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	24	PetscCall(MatAssemblyBegin(A,MAT_FINAL_ASSEMBLY));
57	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	24	PetscCall(MatAssemblyEnd(A,MAT_FINAL_ASSEMBLY));
58
59			/* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
60			Create the eigensolver and set various options
61			- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
62
63	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	24	PetscCall(EPSCreate(PETSC_COMM_WORLD,&eps));
64	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	24	PetscCall(EPSSetOperators(eps,A,NULL));
65	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	24	PetscCall(EPSSetProblemType(eps,EPS_HEP));
66	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	24	PetscCall(EPSSetType(eps,EPSPRIMME));
67
68			/*
69			Set several options
70			*/
71	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	24	PetscCall(EPSSetWhichEigenpairs(eps,EPS_SMALLEST_REAL));
72	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	24	PetscCall(EPSGetST(eps,&st));
73	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	24	PetscCall(STSetType(st,STPRECOND));
74	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	24	PetscCall(STGetKSP(st,&ksp));
75	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	24	PetscCall(KSPGetPC(ksp,&pc));
76	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	24	PetscCall(KSPSetType(ksp,KSPPREONLY));
77	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	24	PetscCall(PCSetType(pc,PCICC));
78
79	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	24	PetscCall(EPSPRIMMESetBlockSize(eps,4));
80	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	24	PetscCall(EPSPRIMMESetMethod(eps,EPS_PRIMME_GD_OLSEN_PLUSK));
81	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	24	PetscCall(EPSSetFromOptions(eps));
82
83			/* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
84			Compute eigenvalues and display info
85			- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
86
87	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	24	PetscCall(EPSSolve(eps));
88	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	24	PetscCall(EPSPRIMMEGetBlockSize(eps,&bs));
89	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	24	PetscCall(EPSPRIMMEGetMethod(eps,&meth));
90	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	24	PetscCall(PetscPrintf(PETSC_COMM_WORLD," PRIMME: using block size %" PetscInt_FMT ", method %s\n",bs,EPSPRIMMEMethods[meth]));
91
92	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	24	PetscCall(EPSErrorView(eps,EPS_ERROR_ABSOLUTE,NULL));
93
94	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	24	PetscCall(EPSDestroy(&eps));
95	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	24	PetscCall(MatDestroy(&A));
96	2/6 ✓ Branch 0 taken 2 times. ✗ Branch 1 not taken. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✗ Branch 5 not taken.	24	PetscCall(SlepcFinalize());
97			return 0;
98			}
99
100			/*TEST
101
102			build:
103			requires: primme
104
105			testset:
106			args: -eps_nev 4
107			requires: primme
108			output_file: output/test34_1.out
109			test:
110			suffix: 1
111			test:
112			suffix: 2
113			args: -st_pc_type bjacobi -eps_target 0.01 -eps_target_real -eps_refined
114			nsize: 2
115			test:
116			suffix: 3
117			args: -eps_smallest_magnitude -eps_harmonic
118
119			TEST*/
120