GCC Code Coverage Report

Directory:	./
File:	src/nep/tutorials/nlevp/gun.c
Date:	2025-11-19 04:19:03

	Exec	Total	Coverage
Lines:	56	60	93.3%
Functions:	1	1	100.0%
Branches:	156	260	60.0%

  
      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).
    
         - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
    
      */
    
      /*
    
         This example implements one of the problems found at
    
             NLEVP: A Collection of Nonlinear Eigenvalue Problems,
    
             The University of Manchester.
    
         The details of the collection can be found at:
    
             [1] T. Betcke et al., "NLEVP: A Collection of Nonlinear Eigenvalue
    
                 Problems", ACM Trans. Math. Software 39(2), Article 7, 2013.
    
         The gun problem arises from model of a radio-frequency gun cavity, with
    
         the complex nonlinear function
    
         T(lambda) = K-lambda*M+i*lambda^(1/2)*W1+i*(lambda-108.8774^2)^(1/2)*W2
    
         Data files can be downloaded from https://gitlab.com/slepc/datafiles
    
      */
    
      static char help[] = "Radio-frequency gun cavity.\n\n"
    
        "The command line options are:\n"
    
        "-K <filename1> -M <filename2> -W1 <filename3> -W2 <filename4>, where filename1,..,filename4 are files containing the matrices in PETSc binary form defining the GUN problem.\n\n";
    
      #include <slepcnep.h>
    
      #define NMAT 4
    
      #define SIGMA 108.8774
    
      3
      int main(int argc,char **argv)
    
      {
    
      3
        Mat            A[NMAT];         /* problem matrices */
    
      3
        FN             f[NMAT];         /* functions to define the nonlinear operator */
    
      3
        FN             ff[2];           /* auxiliary functions to define the nonlinear operator */
    
      3
        NEP            nep;             /* nonlinear eigensolver context */
    
      3
        PetscBool      terse,flg;
    
      3
        const char*    string[NMAT]={"-K","-M","-W1","-W2"};
    
      3
        char           filename[PETSC_MAX_PATH_LEN];
    
      3
        PetscScalar    numer[2],sigma;
    
      3
        PetscInt       i;
    
      3
        PetscViewer    viewer;
    
        1/2✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.

      3
        PetscFunctionBeginUser;
    
        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(SlepcInitialize(&argc,&argv,NULL,help));
    
        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,"GUN problem\n\n"));
    
      #if !defined(PETSC_USE_COMPLEX)
    
        SETERRQ(PETSC_COMM_WORLD,PETSC_ERR_SUP,"This example requires complex scalars!");
    
      #endif
    
        /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
    
                             Load the problem matrices
    
           - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
    
        2/2✓ Branch 0 taken 3 times.
✓ Branch 1 taken 3 times.

      15
        for (i=0;i<NMAT;i++) {
    
        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.

      12
          PetscCall(PetscOptionsGetString(NULL,NULL,string[i],filename,sizeof(filename),&flg));
    
        1/4✗ Branch 0 not taken.
✓ Branch 1 taken 3 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.

      12
          PetscCheck(flg,PETSC_COMM_WORLD,PETSC_ERR_USER_INPUT,"Must indicate a filename with the %s option",string[i]);
    
        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.

      12
          PetscCall(PetscViewerBinaryOpen(PETSC_COMM_WORLD,filename,FILE_MODE_READ,&viewer));
    
        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.

      12
          PetscCall(MatCreate(PETSC_COMM_WORLD,&A[i]));
    
        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.

      12
          PetscCall(MatSetFromOptions(A[i]));
    
        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.

      12
          PetscCall(MatLoad(A[i],viewer));
    
        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.

      12
          PetscCall(PetscViewerDestroy(&viewer));
    
        }
    
        /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
    
                             Create the problem functions
    
           - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
    
        /* f1=1 */
    
        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(FNCreate(PETSC_COMM_WORLD,&f[0]));
    
        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(FNSetType(f[0],FNRATIONAL));
    
      3
        numer[0] = 1.0;
    
        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(FNRationalSetNumerator(f[0],1,numer));
    
        /* f2=-lambda */
    
        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(FNCreate(PETSC_COMM_WORLD,&f[1]));
    
        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(FNSetType(f[1],FNRATIONAL));
    
      3
        numer[0] = -1.0; numer[1] = 0.0;
    
        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(FNRationalSetNumerator(f[1],2,numer));
    
        /* f3=i*sqrt(lambda) */
    
        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(FNCreate(PETSC_COMM_WORLD,&f[2]));
    
        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(FNSetType(f[2],FNSQRT));
    
        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(FNSetScale(f[2],1.0,PETSC_i));
    
        /* f4=i*sqrt(lambda-sigma^2) */
    
      3
        sigma = SIGMA*SIGMA;
    
        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(FNCreate(PETSC_COMM_WORLD,&ff[0]));
    
        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(FNSetType(ff[0],FNSQRT));
    
        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(FNCreate(PETSC_COMM_WORLD,&ff[1]));
    
        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(FNSetType(ff[1],FNRATIONAL));
    
      3
        numer[0] = 1.0; numer[1] = -sigma;
    
        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(FNRationalSetNumerator(ff[1],2,numer));
    
        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(FNCreate(PETSC_COMM_WORLD,&f[3]));
    
        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(FNSetType(f[3],FNCOMBINE));
    
        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(FNCombineSetChildren(f[3],FN_COMBINE_COMPOSE,ff[1],ff[0]));
    
        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(FNSetScale(f[3],1.0,PETSC_i));
    
        /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
    
                      Create the eigensolver and solve the problem
    
           - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
    
        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(NEPCreate(PETSC_COMM_WORLD,&nep));
    
        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(NEPSetSplitOperator(nep,4,A,f,UNKNOWN_NONZERO_PATTERN));
    
        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(NEPSetFromOptions(nep));
    
        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(NEPSolve(nep));
    
        /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
    
                          Display solution and clean up
    
           - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
    
        /* show detailed info unless -terse option is given by user */
    
        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(PetscOptionsHasName(NULL,NULL,"-terse",&terse));
    
        5/8✓ Branch 0 taken 3 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 1 times.
✓ Branch 3 taken 2 times.
✓ Branch 4 taken 1 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 1 times.

      3
        if (terse) PetscCall(NEPErrorView(nep,NEP_ERROR_RELATIVE,NULL));
    
        else {
    
      ✗
          PetscCall(PetscViewerPushFormat(PETSC_VIEWER_STDOUT_WORLD,PETSC_VIEWER_ASCII_INFO_DETAIL));
    
      ✗
          PetscCall(NEPConvergedReasonView(nep,PETSC_VIEWER_STDOUT_WORLD));
    
      ✗
          PetscCall(NEPErrorView(nep,NEP_ERROR_RELATIVE,PETSC_VIEWER_STDOUT_WORLD));
    
      ✗
          PetscCall(PetscViewerPopFormat(PETSC_VIEWER_STDOUT_WORLD));
    
        }
    
        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(NEPDestroy(&nep));
    
        2/2✓ Branch 0 taken 3 times.
✓ Branch 1 taken 3 times.

      15
        for (i=0;i<NMAT;i++) {
    
        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.

      12
          PetscCall(MatDestroy(&A[i]));
    
        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.

      12
          PetscCall(FNDestroy(&f[i]));
    
        }
    
        7/8✓ Branch 0 taken 1 times.
✓ Branch 1 taken 2 times.
✓ Branch 2 taken 3 times.
✓ Branch 3 taken 2 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
✓ Branch 6 taken 1 times.
✓ Branch 7 taken 1 times.

      9
        for (i=0;i<2;i++) PetscCall(FNDestroy(&ff[i]));
    
        2/6✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.

      3
        PetscCall(SlepcFinalize());
    
        return 0;
    
      }
    
      /*TEST
    
         build:
    
            requires: complex
    
         test:
    
            suffix: 1
    
            args: -K ${DATAFILESPATH}/matrices/complex/gun_K.petsc -M ${DATAFILESPATH}/matrices/complex/gun_M.petsc -W1 ${DATAFILESPATH}/matrices/complex/gun_W1.petsc -W2 ${DATAFILESPATH}/matrices/complex/gun_W2.petsc -nep_type nleigs -rg_type polygon -rg_polygon_vertices 12500-1i,120500-1i,120500+30000i,70000+30000i -nep_target 65000 -nep_nev 24 -terse
    
            requires: double complex datafilespath !defined(PETSC_USE_64BIT_INDICES) !defined(PETSCTEST_VALGRIND)
    
            timeoutfactor: 10
    
      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			This example implements one of the problems found at
12			NLEVP: A Collection of Nonlinear Eigenvalue Problems,
13			The University of Manchester.
14			The details of the collection can be found at:
15			[1] T. Betcke et al., "NLEVP: A Collection of Nonlinear Eigenvalue
16			Problems", ACM Trans. Math. Software 39(2), Article 7, 2013.
17
18			The gun problem arises from model of a radio-frequency gun cavity, with
19			the complex nonlinear function
20			T(lambda) = K-lambdaM+ilambda^(1/2)W1+i(lambda-108.8774^2)^(1/2)*W2
21
22			Data files can be downloaded from https://gitlab.com/slepc/datafiles
23			*/
24
25			static char help[] = "Radio-frequency gun cavity.\n\n"
26			"The command line options are:\n"
27			"-K <filename1> -M <filename2> -W1 <filename3> -W2 <filename4>, where filename1,..,filename4 are files containing the matrices in PETSc binary form defining the GUN problem.\n\n";
28
29			#include <slepcnep.h>
30
31			#define NMAT 4
32			#define SIGMA 108.8774
33
34		3	int main(int argc,char **argv)
35			{
36		3	Mat A[NMAT]; /* problem matrices */
37		3	FN f[NMAT]; /* functions to define the nonlinear operator */
38		3	FN ff[2]; /* auxiliary functions to define the nonlinear operator */
39		3	NEP nep; /* nonlinear eigensolver context */
40		3	PetscBool terse,flg;
41		3	const char* string[NMAT]={"-K","-M","-W1","-W2"};
42		3	char filename[PETSC_MAX_PATH_LEN];
43		3	PetscScalar numer[2],sigma;
44		3	PetscInt i;
45		3	PetscViewer viewer;
46
47	1/2 ✓ Branch 0 taken 1 times. ✗ Branch 1 not taken.	3	PetscFunctionBeginUser;
48	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(SlepcInitialize(&argc,&argv,NULL,help));
49
50	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,"GUN problem\n\n"));
51			#if !defined(PETSC_USE_COMPLEX)
52			SETERRQ(PETSC_COMM_WORLD,PETSC_ERR_SUP,"This example requires complex scalars!");
53			#endif
54
55			/* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
56			Load the problem matrices
57			- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
58
59	2/2 ✓ Branch 0 taken 3 times. ✓ Branch 1 taken 3 times.	15	for (i=0;i<NMAT;i++) {
60	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.	12	PetscCall(PetscOptionsGetString(NULL,NULL,string[i],filename,sizeof(filename),&flg));
61	1/4 ✗ Branch 0 not taken. ✓ Branch 1 taken 3 times. ✗ Branch 2 not taken. ✗ Branch 3 not taken.	12	PetscCheck(flg,PETSC_COMM_WORLD,PETSC_ERR_USER_INPUT,"Must indicate a filename with the %s option",string[i]);
62	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.	12	PetscCall(PetscViewerBinaryOpen(PETSC_COMM_WORLD,filename,FILE_MODE_READ,&viewer));
63	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.	12	PetscCall(MatCreate(PETSC_COMM_WORLD,&A[i]));
64	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.	12	PetscCall(MatSetFromOptions(A[i]));
65	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.	12	PetscCall(MatLoad(A[i],viewer));
66	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.	12	PetscCall(PetscViewerDestroy(&viewer));
67			}
68
69			/* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
70			Create the problem functions
71			- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
72
73			/* f1=1 */
74	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(FNCreate(PETSC_COMM_WORLD,&f[0]));
75	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(FNSetType(f[0],FNRATIONAL));
76		3	numer[0] = 1.0;
77	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(FNRationalSetNumerator(f[0],1,numer));
78
79			/* f2=-lambda */
80	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(FNCreate(PETSC_COMM_WORLD,&f[1]));
81	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(FNSetType(f[1],FNRATIONAL));
82		3	numer[0] = -1.0; numer[1] = 0.0;
83	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(FNRationalSetNumerator(f[1],2,numer));
84
85			/* f3=isqrt(lambda) /
86	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(FNCreate(PETSC_COMM_WORLD,&f[2]));
87	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(FNSetType(f[2],FNSQRT));
88	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(FNSetScale(f[2],1.0,PETSC_i));
89
90			/* f4=isqrt(lambda-sigma^2) /
91		3	sigma = SIGMA*SIGMA;
92	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(FNCreate(PETSC_COMM_WORLD,&ff[0]));
93	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(FNSetType(ff[0],FNSQRT));
94	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(FNCreate(PETSC_COMM_WORLD,&ff[1]));
95	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(FNSetType(ff[1],FNRATIONAL));
96		3	numer[0] = 1.0; numer[1] = -sigma;
97	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(FNRationalSetNumerator(ff[1],2,numer));
98	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(FNCreate(PETSC_COMM_WORLD,&f[3]));
99	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(FNSetType(f[3],FNCOMBINE));
100	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(FNCombineSetChildren(f[3],FN_COMBINE_COMPOSE,ff[1],ff[0]));
101	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(FNSetScale(f[3],1.0,PETSC_i));
102
103			/* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
104			Create the eigensolver and solve the problem
105			- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
106
107	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(NEPCreate(PETSC_COMM_WORLD,&nep));
108	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(NEPSetSplitOperator(nep,4,A,f,UNKNOWN_NONZERO_PATTERN));
109	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(NEPSetFromOptions(nep));
110
111	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(NEPSolve(nep));
112
113			/* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
114			Display solution and clean up
115			- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
116
117			/* show detailed info unless -terse option is given by user */
118	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(PetscOptionsHasName(NULL,NULL,"-terse",&terse));
119	5/8 ✓ Branch 0 taken 3 times. ✗ Branch 1 not taken. ✓ Branch 2 taken 1 times. ✓ Branch 3 taken 2 times. ✓ Branch 4 taken 1 times. ✗ Branch 5 not taken. ✗ Branch 6 not taken. ✓ Branch 7 taken 1 times.	3	if (terse) PetscCall(NEPErrorView(nep,NEP_ERROR_RELATIVE,NULL));
120			else {
121		✗	PetscCall(PetscViewerPushFormat(PETSC_VIEWER_STDOUT_WORLD,PETSC_VIEWER_ASCII_INFO_DETAIL));
122		✗	PetscCall(NEPConvergedReasonView(nep,PETSC_VIEWER_STDOUT_WORLD));
123		✗	PetscCall(NEPErrorView(nep,NEP_ERROR_RELATIVE,PETSC_VIEWER_STDOUT_WORLD));
124		✗	PetscCall(PetscViewerPopFormat(PETSC_VIEWER_STDOUT_WORLD));
125			}
126	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(NEPDestroy(&nep));
127	2/2 ✓ Branch 0 taken 3 times. ✓ Branch 1 taken 3 times.	15	for (i=0;i<NMAT;i++) {
128	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.	12	PetscCall(MatDestroy(&A[i]));
129	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.	12	PetscCall(FNDestroy(&f[i]));
130			}
131	7/8 ✓ Branch 0 taken 1 times. ✓ Branch 1 taken 2 times. ✓ Branch 2 taken 3 times. ✓ Branch 3 taken 2 times. ✗ Branch 4 not taken. ✓ Branch 5 taken 1 times. ✓ Branch 6 taken 1 times. ✓ Branch 7 taken 1 times.	9	for (i=0;i<2;i++) PetscCall(FNDestroy(&ff[i]));
132	2/6 ✓ Branch 0 taken 1 times. ✗ Branch 1 not taken. ✓ Branch 2 taken 1 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✗ Branch 5 not taken.	3	PetscCall(SlepcFinalize());
133			return 0;
134			}
135
136			/*TEST
137
138			build:
139			requires: complex
140
141			test:
142			suffix: 1
143			args: -K ${DATAFILESPATH}/matrices/complex/gun_K.petsc -M ${DATAFILESPATH}/matrices/complex/gun_M.petsc -W1 ${DATAFILESPATH}/matrices/complex/gun_W1.petsc -W2 ${DATAFILESPATH}/matrices/complex/gun_W2.petsc -nep_type nleigs -rg_type polygon -rg_polygon_vertices 12500-1i,120500-1i,120500+30000i,70000+30000i -nep_target 65000 -nep_nev 24 -terse
144			requires: double complex datafilespath !defined(PETSC_USE_64BIT_INDICES) !defined(PETSCTEST_VALGRIND)
145			timeoutfactor: 10
146
147			TEST*/
148