GCC Code Coverage Report

Directory:	./
File:	src/pep/impls/peputils.c
Date:	2026-02-22 03:58:10

	Exec	Total	Coverage
Lines:	55	55	100.0%
Functions:	3	3	100.0%
Branches:	94	140	67.1%

  
      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).
    
         - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
    
      */
    
      /*
    
         Common subroutines for several PEP solvers
    
      */
    
      #include <slepc/private/pepimpl.h>    /*I "slepcpep.h" I*/
    
      #include <slepcblaslapack.h>
    
      /*
    
        Computes T_j=phy_idx(T) for a matrix T.
    
          Tp (if j>0) and Tpp (if j>1) are the evaluation
    
          of phy_(j-1) and phy(j-2)respectively.
    
      */
    
      451
      PetscErrorCode PEPEvaluateBasisMat(PEP pep,PetscInt k,PetscScalar *T,PetscInt ldt,PetscInt idx,PetscScalar *Tpp,PetscInt ldtpp,PetscScalar *Tp,PetscInt ldtp,PetscScalar *Tj,PetscInt ldtj)
    
      {
    
      451
        PetscInt       i;
    
      451
        PetscReal      *ca,*cb,*cg;
    
      451
        PetscScalar    g,a;
    
      451
        PetscBLASInt   k_,ldt_,ldtj_,ldtp_;
    
        1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.

      451
        PetscFunctionBegin;
    
        2/2✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.

      451
        if (idx==0) {
    
        2/2✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.

      418
          for (i=0;i<k;i++) {
    
        4/6✓ Branch 0 taken 5 times.
✓ Branch 1 taken 3 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      273
            PetscCall(PetscArrayzero(Tj+i*ldtj,k));
    
      273
            Tj[i+i*ldtj] = 1.0;
    
          }
    
        } else {
    
        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.

      306
          PetscCall(PetscBLASIntCast(ldt,&ldt_));
    
        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.

      306
          PetscCall(PetscBLASIntCast(ldtj,&ldtj_));
    
        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.

      306
          PetscCall(PetscBLASIntCast(ldtp,&ldtp_));
    
        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.

      306
          PetscCall(PetscBLASIntCast(k,&k_));
    
      306
          ca = pep->pbc; cb = pep->pbc+pep->nmat; cg = pep->pbc+2*pep->nmat;
    
        2/2✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.

      868
          for (i=0;i<k;i++) T[i*ldt+i] -= cb[idx-1];
    
        2/2✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.

      306
          if (idx>1) {
    
        7/8✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 8 times.
✓ Branch 3 taken 6 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
✓ Branch 6 taken 2 times.
✓ Branch 7 taken 2 times.

      450
            for (i=0;i<k;i++) PetscCall(PetscArraycpy(Tj+i*ldtj,Tpp+i*ldtpp,k));
    
          }
    
      306
          a = 1/ca[idx-1];
    
        2/2✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.

      306
          g = (idx==1)?0.0:-cg[idx-1]/ca[idx-1];
    
        10/20✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.
✓ Branch 8 taken 2 times.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✓ Branch 11 taken 2 times.
✓ Branch 12 taken 2 times.
✗ Branch 13 not taken.
✓ Branch 14 taken 2 times.
✗ Branch 15 not taken.
✗ Branch 16 not taken.
✓ Branch 17 taken 2 times.
✗ Branch 18 not taken.
✗ Branch 19 not taken.

      306
          PetscCallBLAS("BLASgemm",BLASgemm_("N","N",&k_,&k_,&k_,&a,T,&ldt_,Tp,&ldtp_,&g,Tj,&ldtj_));
    
        2/2✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.

      868
          for (i=0;i<k;i++) T[i*ldt+i] += cb[idx-1];
    
        }
    
        6/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 taken 2 times.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✓ Branch 11 taken 2 times.

      112
        PetscFunctionReturn(PETSC_SUCCESS);
    
      }
    
      /*
    
         PEPEvaluateBasis - evaluate the polynomial basis on a given parameter sigma
    
      */
    
      40661
      PetscErrorCode PEPEvaluateBasis(PEP pep,PetscScalar sigma,PetscScalar isigma,PetscScalar *vals,PetscScalar *ivals)
    
      {
    
      40661
        PetscInt   nmat=pep->nmat,k;
    
      40661
        PetscReal  *a=pep->pbc,*b=pep->pbc+nmat,*g=pep->pbc+2*nmat;
    
        1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.

      40661
        PetscFunctionBegin;
    
        4/4✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 8 times.
✓ Branch 3 taken 8 times.

      165166
        if (ivals) for (k=0;k<nmat;k++) ivals[k] = 0.0;
    
      40661
        vals[0] = 1.0;
    
      40661
        vals[1] = (sigma-b[0])/a[0];
    
      #if !defined(PETSC_USE_COMPLEX)
    
        2/2✓ Branch 0 taken 4 times.
✓ Branch 1 taken 4 times.

      17581
        if (ivals) ivals[1] = isigma/a[0];
    
      #endif
    
        2/2✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.

      101766
        for (k=2;k<nmat;k++) {
    
      61105
          vals[k] = ((sigma-b[k-1])*vals[k-1]-g[k-1]*vals[k-2])/a[k-1];
    
        2/2✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.

      61105
          if (ivals) vals[k] -= isigma*ivals[k-1]/a[k-1];
    
      #if !defined(PETSC_USE_COMPLEX)
    
      22395
          if (ivals) ivals[k] = ((sigma-b[k-1])*ivals[k-1]+isigma*vals[k-1]-g[k-1]*ivals[k-2])/a[k-1];
    
      #endif
    
        }
    
        6/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 taken 2 times.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✓ Branch 11 taken 2 times.

      40661
        PetscFunctionReturn(PETSC_SUCCESS);
    
      }
    
      /*
    
         PEPEvaluateBasisDerivative - evaluate the derivative of the polynomial basis on a given parameter sigma
    
      */
    
      4252
      PetscErrorCode PEPEvaluateBasisDerivative(PEP pep,PetscScalar sigma,PetscScalar isigma,PetscScalar *vals,PetscScalar *ivals)
    
      {
    
      4252
        PetscInt       nmat=pep->nmat,k;
    
      4252
        PetscReal      *a=pep->pbc,*b=pep->pbc+nmat,*g=pep->pbc+2*nmat;
    
        1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.

      4252
        PetscFunctionBegin;
    
        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.

      4252
        PetscCall(PEPEvaluateBasis(pep,sigma,isigma,vals,ivals));
    
        2/2✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.

      16026
        for (k=nmat-1;k>0;k--) {
    
      11774
          vals[k] = vals[k-1];
    
        2/2✓ Branch 0 taken 8 times.
✓ Branch 1 taken 4 times.

      11774
          if (ivals) ivals[k] = ivals[k-1];
    
        }
    
      4252
        vals[0] = 0.0;
    
      4252
        vals[1] = vals[1]/a[0];
    
      #if !defined(PETSC_USE_COMPLEX)
    
        1/2✓ Branch 0 taken 4 times.
✗ Branch 1 not taken.

      1129
        if (ivals) ivals[1] = ivals[1]/a[0];
    
      #endif
    
        2/2✓ Branch 0 taken 8 times.
✓ Branch 1 taken 8 times.

      11774
        for (k=2;k<nmat;k++) {
    
      7522
          vals[k] += (sigma-b[k-1])*vals[k-1]-g[k-1]*vals[k-2];
    
        2/2✓ Branch 0 taken 8 times.
✓ Branch 1 taken 4 times.

      7522
          if (ivals) vals[k] -= isigma*ivals[k-1];
    
      7522
          vals[k] /= a[k-1];
    
      #if !defined(PETSC_USE_COMPLEX)
    
        1/2✓ Branch 0 taken 4 times.
✗ Branch 1 not taken.

      1239
          if (ivals) {
    
      1239
            ivals[k] += (sigma-b[k-1])*ivals[k-1]+isigma*vals[k-1]-g[k-1]*ivals[k-2];
    
      1239
            ivals[k] /= a[k-1];
    
          }
    
      #endif
    
        }
    
        6/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 taken 2 times.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✓ Branch 11 taken 2 times.

      1062
        PetscFunctionReturn(PETSC_SUCCESS);
    
      }

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			Common subroutines for several PEP solvers
12			*/
13
14			#include <slepc/private/pepimpl.h> /I "slepcpep.h" I/
15			#include <slepcblaslapack.h>
16
17			/*
18			Computes T_j=phy_idx(T) for a matrix T.
19			Tp (if j>0) and Tpp (if j>1) are the evaluation
20			of phy_(j-1) and phy(j-2)respectively.
21			*/
22		451	PetscErrorCode PEPEvaluateBasisMat(PEP pep,PetscInt k,PetscScalar T,PetscInt ldt,PetscInt idx,PetscScalar Tpp,PetscInt ldtpp,PetscScalar Tp,PetscInt ldtp,PetscScalar Tj,PetscInt ldtj)
23			{
24		451	PetscInt i;
25		451	PetscReal ca,cb,*cg;
26		451	PetscScalar g,a;
27		451	PetscBLASInt k_,ldt_,ldtj_,ldtp_;
28
29	1/2 ✓ Branch 0 taken 2 times. ✗ Branch 1 not taken.	451	PetscFunctionBegin;
30	2/2 ✓ Branch 0 taken 8 times. ✓ Branch 1 taken 8 times.	451	if (idx==0) {
31	2/2 ✓ Branch 0 taken 8 times. ✓ Branch 1 taken 8 times.	418	for (i=0;i<k;i++) {
32	4/6 ✓ Branch 0 taken 5 times. ✓ Branch 1 taken 3 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	273	PetscCall(PetscArrayzero(Tj+i*ldtj,k));
33		273	Tj[i+i*ldtj] = 1.0;
34			}
35			} else {
36	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.	306	PetscCall(PetscBLASIntCast(ldt,&ldt_));
37	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 6 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	306	PetscCall(PetscBLASIntCast(ldtj,&ldtj_));
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.	306	PetscCall(PetscBLASIntCast(ldtp,&ldtp_));
39	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.	306	PetscCall(PetscBLASIntCast(k,&k_));
40		306	ca = pep->pbc; cb = pep->pbc+pep->nmat; cg = pep->pbc+2*pep->nmat;
41	2/2 ✓ Branch 0 taken 8 times. ✓ Branch 1 taken 8 times.	868	for (i=0;i<k;i++) T[i*ldt+i] -= cb[idx-1];
42	2/2 ✓ Branch 0 taken 8 times. ✓ Branch 1 taken 8 times.	306	if (idx>1) {
43	7/8 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 6 times. ✓ Branch 2 taken 8 times. ✓ Branch 3 taken 6 times. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times. ✓ Branch 6 taken 2 times. ✓ Branch 7 taken 2 times.	450	for (i=0;i<k;i++) PetscCall(PetscArraycpy(Tj+ildtj,Tpp+ildtpp,k));
44			}
45		306	a = 1/ca[idx-1];
46	2/2 ✓ Branch 0 taken 8 times. ✓ Branch 1 taken 8 times.	306	g = (idx==1)?0.0:-cg[idx-1]/ca[idx-1];
47	10/20 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 6 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✓ Branch 4 taken 2 times. ✗ Branch 5 not taken. ✗ Branch 6 not taken. ✓ Branch 7 taken 2 times. ✓ Branch 8 taken 2 times. ✗ Branch 9 not taken. ✗ Branch 10 not taken. ✓ Branch 11 taken 2 times. ✓ Branch 12 taken 2 times. ✗ Branch 13 not taken. ✓ Branch 14 taken 2 times. ✗ Branch 15 not taken. ✗ Branch 16 not taken. ✓ Branch 17 taken 2 times. ✗ Branch 18 not taken. ✗ Branch 19 not taken.	306	PetscCallBLAS("BLASgemm",BLASgemm_("N","N",&k_,&k_,&k_,&a,T,&ldt_,Tp,&ldtp_,&g,Tj,&ldtj_));
48	2/2 ✓ Branch 0 taken 8 times. ✓ Branch 1 taken 8 times.	868	for (i=0;i<k;i++) T[i*ldt+i] += cb[idx-1];
49			}
50	6/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 taken 2 times. ✗ Branch 9 not taken. ✗ Branch 10 not taken. ✓ Branch 11 taken 2 times.	112	PetscFunctionReturn(PETSC_SUCCESS);
51			}
52
53			/*
54			PEPEvaluateBasis - evaluate the polynomial basis on a given parameter sigma
55			*/
56		40661	PetscErrorCode PEPEvaluateBasis(PEP pep,PetscScalar sigma,PetscScalar isigma,PetscScalar vals,PetscScalar ivals)
57			{
58		40661	PetscInt nmat=pep->nmat,k;
59		40661	PetscReal a=pep->pbc,b=pep->pbc+nmat,g=pep->pbc+2nmat;
60
61	1/2 ✓ Branch 0 taken 2 times. ✗ Branch 1 not taken.	40661	PetscFunctionBegin;
62	4/4 ✓ Branch 0 taken 8 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 8 times. ✓ Branch 3 taken 8 times.	165166	if (ivals) for (k=0;k<nmat;k++) ivals[k] = 0.0;
63		40661	vals[0] = 1.0;
64		40661	vals[1] = (sigma-b[0])/a[0];
65			#if !defined(PETSC_USE_COMPLEX)
66	2/2 ✓ Branch 0 taken 4 times. ✓ Branch 1 taken 4 times.	17581	if (ivals) ivals[1] = isigma/a[0];
67			#endif
68	2/2 ✓ Branch 0 taken 8 times. ✓ Branch 1 taken 8 times.	101766	for (k=2;k<nmat;k++) {
69		61105	vals[k] = ((sigma-b[k-1])vals[k-1]-g[k-1]vals[k-2])/a[k-1];
70	2/2 ✓ Branch 0 taken 8 times. ✓ Branch 1 taken 8 times.	61105	if (ivals) vals[k] -= isigma*ivals[k-1]/a[k-1];
71			#if !defined(PETSC_USE_COMPLEX)
72		22395	if (ivals) ivals[k] = ((sigma-b[k-1])ivals[k-1]+isigmavals[k-1]-g[k-1]*ivals[k-2])/a[k-1];
73			#endif
74			}
75	6/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 taken 2 times. ✗ Branch 9 not taken. ✗ Branch 10 not taken. ✓ Branch 11 taken 2 times.	40661	PetscFunctionReturn(PETSC_SUCCESS);
76			}
77
78			/*
79			PEPEvaluateBasisDerivative - evaluate the derivative of the polynomial basis on a given parameter sigma
80			*/
81		4252	PetscErrorCode PEPEvaluateBasisDerivative(PEP pep,PetscScalar sigma,PetscScalar isigma,PetscScalar vals,PetscScalar ivals)
82			{
83		4252	PetscInt nmat=pep->nmat,k;
84		4252	PetscReal a=pep->pbc,b=pep->pbc+nmat,g=pep->pbc+2nmat;
85
86	1/2 ✓ Branch 0 taken 2 times. ✗ Branch 1 not taken.	4252	PetscFunctionBegin;
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.	4252	PetscCall(PEPEvaluateBasis(pep,sigma,isigma,vals,ivals));
88	2/2 ✓ Branch 0 taken 8 times. ✓ Branch 1 taken 8 times.	16026	for (k=nmat-1;k>0;k--) {
89		11774	vals[k] = vals[k-1];
90	2/2 ✓ Branch 0 taken 8 times. ✓ Branch 1 taken 4 times.	11774	if (ivals) ivals[k] = ivals[k-1];
91			}
92		4252	vals[0] = 0.0;
93		4252	vals[1] = vals[1]/a[0];
94			#if !defined(PETSC_USE_COMPLEX)
95	1/2 ✓ Branch 0 taken 4 times. ✗ Branch 1 not taken.	1129	if (ivals) ivals[1] = ivals[1]/a[0];
96			#endif
97	2/2 ✓ Branch 0 taken 8 times. ✓ Branch 1 taken 8 times.	11774	for (k=2;k<nmat;k++) {
98		7522	vals[k] += (sigma-b[k-1])vals[k-1]-g[k-1]vals[k-2];
99	2/2 ✓ Branch 0 taken 8 times. ✓ Branch 1 taken 4 times.	7522	if (ivals) vals[k] -= isigma*ivals[k-1];
100		7522	vals[k] /= a[k-1];
101			#if !defined(PETSC_USE_COMPLEX)
102	1/2 ✓ Branch 0 taken 4 times. ✗ Branch 1 not taken.	1239	if (ivals) {
103		1239	ivals[k] += (sigma-b[k-1])ivals[k-1]+isigmavals[k-1]-g[k-1]*ivals[k-2];
104		1239	ivals[k] /= a[k-1];
105			}
106			#endif
107			}
108	6/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 taken 2 times. ✗ Branch 9 not taken. ✗ Branch 10 not taken. ✓ Branch 11 taken 2 times.	1062	PetscFunctionReturn(PETSC_SUCCESS);
109			}
110

Function (Line)	Call count	Block coverage
PEPEvaluateBasis (line 56)	called 40661 times	100.0%
PEPEvaluateBasisDerivative (line 81)	called 4252 times	94.0%
PEPEvaluateBasisMat (line 22)	called 451 times	83.0%