GCC Code Coverage Report

Directory:	./
File:	src/eps/tests/test11.c
Date:	2026-02-22 03:58:10

	Exec	Total	Coverage
Lines:	81	81	100.0%
Functions:	3	3	100.0%
Branches:	191	288	66.3%

  
      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 the same problem as in ex5, but with a user-defined sorting criterion.\n\n"
    
        "It is a standard nonsymmetric eigenproblem with real eigenvalues and the rightmost eigenvalue is known to be 1.\n"
    
        "This example illustrates how the user can set a custom spectrum selection.\n\n"
    
        "The command line options are:\n"
    
        "  -m <m>, where <m> = number of grid subdivisions in each dimension.\n\n";
    
      #include <slepceps.h>
    
      /*
    
         User-defined routines
    
      */
    
      PetscErrorCode MyEigenSort(PetscScalar ar,PetscScalar ai,PetscScalar br,PetscScalar bi,PetscInt *r,void *ctx);
    
      PetscErrorCode MatMarkovModel(PetscInt m,Mat A);
    
      40
      int main(int argc,char **argv)
    
      {
    
      40
        Vec            v0;              /* initial vector */
    
      40
        Mat            A;               /* operator matrix */
    
      40
        EPS            eps;             /* eigenproblem solver context */
    
      40
        ST             st;              /* spectral transformation associated */
    
      40
        PetscReal      tol=PETSC_SMALL;
    
      40
        PetscScalar    target=0.5;
    
      40
        PetscInt       N,m=15,nev;
    
      40
        char           str[50];
    
        1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.

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

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

      40
        PetscCall(PetscOptionsGetInt(NULL,NULL,"-m",&m,NULL));
    
      40
        N = m*(m+1)/2;
    
        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.

      40
        PetscCall(PetscPrintf(PETSC_COMM_WORLD,"\nMarkov Model, N=%" PetscInt_FMT " (m=%" PetscInt_FMT ")\n",N,m));
    
        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.

      40
        PetscCall(PetscOptionsGetScalar(NULL,NULL,"-target",&target,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.

      40
        PetscCall(SlepcSNPrintfScalar(str,sizeof(str),target,PETSC_FALSE));
    
        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.

      40
        PetscCall(PetscPrintf(PETSC_COMM_WORLD,"Searching closest eigenvalues to the right of %s.\n\n",str));
    
        /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
    
           Compute the operator matrix that defines the eigensystem, Ax=kx
    
           - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
    
        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.

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

      40
        PetscCall(MatSetSizes(A,PETSC_DECIDE,PETSC_DECIDE,N,N));
    
        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.

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

      40
        PetscCall(MatMarkovModel(m,A));
    
        /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
    
                      Create the eigensolver and set various options
    
           - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
    
        /*
    
           Create eigensolver context
    
        */
    
        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.

      40
        PetscCall(EPSCreate(PETSC_COMM_WORLD,&eps));
    
        /*
    
           Set operators. In this case, it is a standard eigenvalue 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.

      40
        PetscCall(EPSSetOperators(eps,A,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.

      40
        PetscCall(EPSSetProblemType(eps,EPS_NHEP));
    
        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.

      40
        PetscCall(EPSSetTolerances(eps,tol,PETSC_CURRENT));
    
        /*
    
           Set the custom comparing routine in order to obtain the eigenvalues
    
           closest to the target on the right only
    
        */
    
        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.

      40
        PetscCall(EPSSetEigenvalueComparison(eps,MyEigenSort,&target));
    
        /*
    
           Set solver parameters at runtime
    
        */
    
        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.

      40
        PetscCall(EPSSetFromOptions(eps));
    
        /*
    
           Set the preconditioner based on A - target * I
    
        */
    
        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.

      40
        PetscCall(EPSGetST(eps,&st));
    
        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.

      40
        PetscCall(STSetShift(st,target));
    
        /*
    
           Set the initial vector. This is optional, if not done the initial
    
           vector is set to random values
    
        */
    
        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.

      40
        PetscCall(MatCreateVecs(A,&v0,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.

      40
        PetscCall(VecSet(v0,1.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.

      40
        PetscCall(EPSSetInitialSpace(eps,1,&v0));
    
        /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
    
                            Solve the eigensystem
    
           - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
    
        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.

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

      40
        PetscCall(EPSGetDimensions(eps,&nev,NULL,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.

      40
        PetscCall(PetscPrintf(PETSC_COMM_WORLD," Number of requested eigenvalues: %" PetscInt_FMT "\n",nev));
    
        /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
    
                          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.

      40
        PetscCall(EPSErrorView(eps,EPS_ERROR_RELATIVE,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.

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

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

      40
        PetscCall(VecDestroy(&v0));
    
        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.

      40
        PetscCall(SlepcFinalize());
    
        return 0;
    
      }
    
      40
      PetscErrorCode MatMarkovModel(PetscInt m,Mat A)
    
      {
    
      40
        const PetscReal cst = 0.5/(PetscReal)(m-1);
    
      40
        PetscReal       pd,pu;
    
      40
        PetscInt        Istart,Iend,i,j,jmax,ix=0;
    
        1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.

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

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

      640
        for (i=1;i<=m;i++) {
    
      600
          jmax = m-i+1;
    
        2/2✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.

      5400
          for (j=1;j<=jmax;j++) {
    
      4800
            ix = ix + 1;
    
        2/4✓ Branch 0 taken 8 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 8 times.

      4800
            if (ix-1<Istart || ix>Iend) continue;  /* compute only owned rows */
    
        2/2✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.

      4800
            if (j!=jmax) {
    
      4200
              pd = cst*(PetscReal)(i+j-1);
    
              /* north */
    
        6/8✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 6 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.

      4200
              if (i==1) PetscCall(MatSetValue(A,ix-1,ix,2*pd,INSERT_VALUES));
    
        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.

      3640
              else PetscCall(MatSetValue(A,ix-1,ix,pd,INSERT_VALUES));
    
              /* east */
    
        6/8✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 6 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.

      4200
              if (j==1) PetscCall(MatSetValue(A,ix-1,ix+jmax-1,2*pd,INSERT_VALUES));
    
        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.

      3640
              else PetscCall(MatSetValue(A,ix-1,ix+jmax-1,pd,INSERT_VALUES));
    
            }
    
            /* south */
    
      4800
            pu = 0.5 - cst*(PetscReal)(i+j-3);
    
        6/8✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 6 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.

      4800
            if (j>1) PetscCall(MatSetValue(A,ix-1,ix-2,pu,INSERT_VALUES));
    
            /* west */
    
        6/8✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 6 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.

      4800
            if (i>1) PetscCall(MatSetValue(A,ix-1,ix-jmax-2,pu,INSERT_VALUES));
    
          }
    
        }
    
        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.

      40
        PetscCall(MatAssemblyBegin(A,MAT_FINAL_ASSEMBLY));
    
        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.

      40
        PetscCall(MatAssemblyEnd(A,MAT_FINAL_ASSEMBLY));
    
        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.

      10
        PetscFunctionReturn(PETSC_SUCCESS);
    
      }
    
      /*
    
          Function for user-defined eigenvalue ordering criterion.
    
          Given two eigenvalues ar+i*ai and br+i*bi, the subroutine must choose
    
          one of them as the preferred one according to the criterion.
    
          In this example, the preferred value is the one closest to the target,
    
          but on the right side.
    
      */
    
      162795
      PetscErrorCode MyEigenSort(PetscScalar ar,PetscScalar ai,PetscScalar br,PetscScalar bi,PetscInt *r,void *ctx)
    
      {
    
      162795
        PetscScalar target = *(PetscScalar*)ctx;
    
      162795
        PetscReal   da,db;
    
      162795
        PetscBool   aisright,bisright;
    
        1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.

      162795
        PetscFunctionBeginUser;
    
        2/2✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.

      162795
        if (PetscRealPart(target) < PetscRealPart(ar)) aisright = PETSC_TRUE;
    
      97166
        else aisright = PETSC_FALSE;
    
        2/2✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.

      162795
        if (PetscRealPart(target) < PetscRealPart(br)) bisright = PETSC_TRUE;
    
      140341
        else bisright = PETSC_FALSE;
    
        2/2✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.

      162795
        if (aisright == bisright) {
    
          /* both are on the same side of the target */
    
      112548
          da = SlepcAbsEigenvalue(ar-target,ai);
    
      112548
          db = SlepcAbsEigenvalue(br-target,bi);
    
        2/2✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.

      112548
          if (da < db) *r = -1;
    
        2/2✓ Branch 0 taken 8 times.
✓ Branch 1 taken 4 times.

      30963
          else if (da > db) *r = 1;
    
      22
          else *r = 0;
    
        2/2✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.

      50247
        } else if (aisright && !bisright) *r = -1; /* 'a' is on the right */
    
      3536
        else *r = 1;  /* 'b' is on the right */
    
        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.

      162795
        PetscFunctionReturn(PETSC_SUCCESS);
    
      }
    
      /*TEST
    
         testset:
    
            args: -eps_nev 4
    
            requires: !single
    
            output_file: output/test11_1.out
    
            test:
    
               suffix: 1
    
               args: -eps_type {{krylovschur arnoldi lapack}} -st_type sinvert
    
            test:
    
               suffix: 1_ks_cayley
    
               args: -st_type cayley -st_cayley_antishift 1
    
         test:
    
            suffix: 2
    
            args: -target 0.77 -eps_type gd -eps_nev 4 -eps_tol 1e-7 -eps_gd_krylov_start -eps_gd_blocksize 3
    
            requires: double
    
            filter: sed -e "s/[+-]0\.0*i//g"
    
      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 the same problem as in ex5, but with a user-defined sorting criterion.\n\n"
12			"It is a standard nonsymmetric eigenproblem with real eigenvalues and the rightmost eigenvalue is known to be 1.\n"
13			"This example illustrates how the user can set a custom spectrum selection.\n\n"
14			"The command line options are:\n"
15			" -m <m>, where <m> = number of grid subdivisions in each dimension.\n\n";
16
17			#include <slepceps.h>
18
19			/*
20			User-defined routines
21			*/
22
23			PetscErrorCode MyEigenSort(PetscScalar ar,PetscScalar ai,PetscScalar br,PetscScalar bi,PetscInt r,void ctx);
24			PetscErrorCode MatMarkovModel(PetscInt m,Mat A);
25
26		40	int main(int argc,char **argv)
27			{
28		40	Vec v0; /* initial vector */
29		40	Mat A; /* operator matrix */
30		40	EPS eps; /* eigenproblem solver context */
31		40	ST st; /* spectral transformation associated */
32		40	PetscReal tol=PETSC_SMALL;
33		40	PetscScalar target=0.5;
34		40	PetscInt N,m=15,nev;
35		40	char str[50];
36
37	1/2 ✓ Branch 0 taken 2 times. ✗ Branch 1 not taken.	40	PetscFunctionBeginUser;
38	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.	40	PetscCall(SlepcInitialize(&argc,&argv,NULL,help));
39
40	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.	40	PetscCall(PetscOptionsGetInt(NULL,NULL,"-m",&m,NULL));
41		40	N = m*(m+1)/2;
42	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.	40	PetscCall(PetscPrintf(PETSC_COMM_WORLD,"\nMarkov Model, N=%" PetscInt_FMT " (m=%" PetscInt_FMT ")\n",N,m));
43	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.	40	PetscCall(PetscOptionsGetScalar(NULL,NULL,"-target",&target,NULL));
44	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.	40	PetscCall(SlepcSNPrintfScalar(str,sizeof(str),target,PETSC_FALSE));
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.	40	PetscCall(PetscPrintf(PETSC_COMM_WORLD,"Searching closest eigenvalues to the right of %s.\n\n",str));
46
47			/* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
48			Compute the operator matrix that defines the eigensystem, Ax=kx
49			- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
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.	40	PetscCall(MatCreate(PETSC_COMM_WORLD,&A));
52	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.	40	PetscCall(MatSetSizes(A,PETSC_DECIDE,PETSC_DECIDE,N,N));
53	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.	40	PetscCall(MatSetFromOptions(A));
54	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.	40	PetscCall(MatMarkovModel(m,A));
55
56			/* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
57			Create the eigensolver and set various options
58			- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
59
60			/*
61			Create eigensolver context
62			*/
63	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.	40	PetscCall(EPSCreate(PETSC_COMM_WORLD,&eps));
64
65			/*
66			Set operators. In this case, it is a standard eigenvalue problem
67			*/
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.	40	PetscCall(EPSSetOperators(eps,A,NULL));
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.	40	PetscCall(EPSSetProblemType(eps,EPS_NHEP));
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.	40	PetscCall(EPSSetTolerances(eps,tol,PETSC_CURRENT));
71
72			/*
73			Set the custom comparing routine in order to obtain the eigenvalues
74			closest to the target on the right only
75			*/
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.	40	PetscCall(EPSSetEigenvalueComparison(eps,MyEigenSort,&target));
77
78			/*
79			Set solver parameters at runtime
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.	40	PetscCall(EPSSetFromOptions(eps));
82
83			/*
84			Set the preconditioner based on A - target * I
85			*/
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.	40	PetscCall(EPSGetST(eps,&st));
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.	40	PetscCall(STSetShift(st,target));
88
89			/*
90			Set the initial vector. This is optional, if not done the initial
91			vector is set to random values
92			*/
93	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.	40	PetscCall(MatCreateVecs(A,&v0,NULL));
94	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.	40	PetscCall(VecSet(v0,1.0));
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.	40	PetscCall(EPSSetInitialSpace(eps,1,&v0));
96
97			/* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
98			Solve the eigensystem
99			- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
100
101	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.	40	PetscCall(EPSSolve(eps));
102	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.	40	PetscCall(EPSGetDimensions(eps,&nev,NULL,NULL));
103	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.	40	PetscCall(PetscPrintf(PETSC_COMM_WORLD," Number of requested eigenvalues: %" PetscInt_FMT "\n",nev));
104
105			/* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
106			Display solution and clean up
107			- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
108
109	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.	40	PetscCall(EPSErrorView(eps,EPS_ERROR_RELATIVE,NULL));
110	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.	40	PetscCall(EPSDestroy(&eps));
111	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.	40	PetscCall(MatDestroy(&A));
112	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.	40	PetscCall(VecDestroy(&v0));
113	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.	40	PetscCall(SlepcFinalize());
114			return 0;
115			}
116
117		40	PetscErrorCode MatMarkovModel(PetscInt m,Mat A)
118			{
119		40	const PetscReal cst = 0.5/(PetscReal)(m-1);
120		40	PetscReal pd,pu;
121		40	PetscInt Istart,Iend,i,j,jmax,ix=0;
122
123	1/2 ✓ Branch 0 taken 2 times. ✗ Branch 1 not taken.	40	PetscFunctionBeginUser;
124	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.	40	PetscCall(MatGetOwnershipRange(A,&Istart,&Iend));
125	2/2 ✓ Branch 0 taken 8 times. ✓ Branch 1 taken 8 times.	640	for (i=1;i<=m;i++) {
126		600	jmax = m-i+1;
127	2/2 ✓ Branch 0 taken 8 times. ✓ Branch 1 taken 8 times.	5400	for (j=1;j<=jmax;j++) {
128		4800	ix = ix + 1;
129	2/4 ✓ Branch 0 taken 8 times. ✗ Branch 1 not taken. ✗ Branch 2 not taken. ✓ Branch 3 taken 8 times.	4800	if (ix-1<Istart \|\| ix>Iend) continue; /* compute only owned rows */
130	2/2 ✓ Branch 0 taken 8 times. ✓ Branch 1 taken 8 times.	4800	if (j!=jmax) {
131		4200	pd = cst*(PetscReal)(i+j-1);
132			/* north */
133	6/8 ✓ Branch 0 taken 8 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✓ Branch 3 taken 6 times. ✓ Branch 4 taken 2 times. ✗ Branch 5 not taken. ✗ Branch 6 not taken. ✓ Branch 7 taken 2 times.	4200	if (i==1) PetscCall(MatSetValue(A,ix-1,ix,2*pd,INSERT_VALUES));
134	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.	3640	else PetscCall(MatSetValue(A,ix-1,ix,pd,INSERT_VALUES));
135			/* east */
136	6/8 ✓ Branch 0 taken 8 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✓ Branch 3 taken 6 times. ✓ Branch 4 taken 2 times. ✗ Branch 5 not taken. ✗ Branch 6 not taken. ✓ Branch 7 taken 2 times.	4200	if (j==1) PetscCall(MatSetValue(A,ix-1,ix+jmax-1,2*pd,INSERT_VALUES));
137	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.	3640	else PetscCall(MatSetValue(A,ix-1,ix+jmax-1,pd,INSERT_VALUES));
138			}
139			/* south */
140		4800	pu = 0.5 - cst*(PetscReal)(i+j-3);
141	6/8 ✓ Branch 0 taken 8 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✓ Branch 3 taken 6 times. ✓ Branch 4 taken 2 times. ✗ Branch 5 not taken. ✗ Branch 6 not taken. ✓ Branch 7 taken 2 times.	4800	if (j>1) PetscCall(MatSetValue(A,ix-1,ix-2,pu,INSERT_VALUES));
142			/* west */
143	6/8 ✓ Branch 0 taken 8 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✓ Branch 3 taken 6 times. ✓ Branch 4 taken 2 times. ✗ Branch 5 not taken. ✗ Branch 6 not taken. ✓ Branch 7 taken 2 times.	4800	if (i>1) PetscCall(MatSetValue(A,ix-1,ix-jmax-2,pu,INSERT_VALUES));
144			}
145			}
146	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.	40	PetscCall(MatAssemblyBegin(A,MAT_FINAL_ASSEMBLY));
147	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.	40	PetscCall(MatAssemblyEnd(A,MAT_FINAL_ASSEMBLY));
148	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.	10	PetscFunctionReturn(PETSC_SUCCESS);
149			}
150
151			/*
152			Function for user-defined eigenvalue ordering criterion.
153
154			Given two eigenvalues ar+iai and br+ibi, the subroutine must choose
155			one of them as the preferred one according to the criterion.
156			In this example, the preferred value is the one closest to the target,
157			but on the right side.
158			*/
159		162795	PetscErrorCode MyEigenSort(PetscScalar ar,PetscScalar ai,PetscScalar br,PetscScalar bi,PetscInt r,void ctx)
160			{
161		162795	PetscScalar target = (PetscScalar)ctx;
162		162795	PetscReal da,db;
163		162795	PetscBool aisright,bisright;
164
165	1/2 ✓ Branch 0 taken 2 times. ✗ Branch 1 not taken.	162795	PetscFunctionBeginUser;
166	2/2 ✓ Branch 0 taken 8 times. ✓ Branch 1 taken 8 times.	162795	if (PetscRealPart(target) < PetscRealPart(ar)) aisright = PETSC_TRUE;
167		97166	else aisright = PETSC_FALSE;
168	2/2 ✓ Branch 0 taken 8 times. ✓ Branch 1 taken 8 times.	162795	if (PetscRealPart(target) < PetscRealPart(br)) bisright = PETSC_TRUE;
169		140341	else bisright = PETSC_FALSE;
170	2/2 ✓ Branch 0 taken 8 times. ✓ Branch 1 taken 8 times.	162795	if (aisright == bisright) {
171			/* both are on the same side of the target */
172		112548	da = SlepcAbsEigenvalue(ar-target,ai);
173		112548	db = SlepcAbsEigenvalue(br-target,bi);
174	2/2 ✓ Branch 0 taken 8 times. ✓ Branch 1 taken 8 times.	112548	if (da < db) *r = -1;
175	2/2 ✓ Branch 0 taken 8 times. ✓ Branch 1 taken 4 times.	30963	else if (da > db) *r = 1;
176		22	else *r = 0;
177	2/2 ✓ Branch 0 taken 8 times. ✓ Branch 1 taken 8 times.	50247	} else if (aisright && !bisright) r = -1; / 'a' is on the right */
178		3536	else r = 1; / 'b' is on the right */
179	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.	162795	PetscFunctionReturn(PETSC_SUCCESS);
180			}
181
182			/*TEST
183
184			testset:
185			args: -eps_nev 4
186			requires: !single
187			output_file: output/test11_1.out
188			test:
189			suffix: 1
190			args: -eps_type {{krylovschur arnoldi lapack}} -st_type sinvert
191			test:
192			suffix: 1_ks_cayley
193			args: -st_type cayley -st_cayley_antishift 1
194
195			test:
196			suffix: 2
197			args: -target 0.77 -eps_type gd -eps_nev 4 -eps_tol 1e-7 -eps_gd_krylov_start -eps_gd_blocksize 3
198			requires: double
199			filter: sed -e "s/[+-]0\.0*i//g"
200
201			TEST*/
202

Function (Line)	Call count	Block coverage
MatMarkovModel (line 117)	called 40 times	82.0%
MyEigenSort (line 159)	called 162795 times	100.0%
main (line 26)	called 40 times	83.0%