GCC Code Coverage Report

Directory:	./
File:	src/nep/interface/nepresolv.c
Date:	2025-10-04 04:19:13

	Exec	Total	Coverage
Lines:	65	65	100.0%
Functions:	3	3	100.0%
Branches:	185	350	52.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).
    
         - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
    
      */
    
      /*
    
         NEP routines related to resolvent T^{-1}(z) = sum_i (z-lambda_i)^{-1} x_i y_i'
    
      */
    
      #include <slepc/private/nepimpl.h>       /*I "slepcnep.h" I*/
    
      typedef struct {
    
        NEP              nep;
    
        RG               rg;
    
        PetscScalar      omega;
    
        PetscScalar      *nfactor;         /* normalization factors y_i'*T'(lambda_i)*x_i */
    
        PetscBool        *nfactor_avail;
    
        PetscScalar      *dots;            /* inner products y_i'*v */
    
        PetscBool        *dots_avail;
    
        PetscObjectId    vid;
    
        PetscObjectState vstate;
    
      } NEP_RESOLVENT_MATSHELL;
    
      40
      static PetscErrorCode MatMult_Resolvent(Mat M,Vec v,Vec r)
    
      {
    
      40
        NEP_RESOLVENT_MATSHELL *ctx;
    
      40
        NEP                    nep;
    
      40
        PetscInt               i,inside=1;
    
      40
        PetscScalar            alpha;
    
      40
        Vec                    x,y,z,w;
    
        1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.

      40
        PetscFunctionBegin;
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      40
        PetscCall(MatShellGetContext(M,&ctx));
    
      40
        nep = ctx->nep;
    
      40
        w = nep->work[0];
    
      40
        z = nep->work[1];
    
        4/4✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 10 times.
✓ Branch 3 taken 10 times.

      40
        if (((PetscObject)v)->id != ctx->vid || ((PetscObject)v)->state != ctx->vstate) {
    
        4/6✓ Branch 0 taken 5 times.
✓ Branch 1 taken 5 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      20
          PetscCall(PetscArrayzero(ctx->dots_avail,ctx->nep->nconv));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      20
          PetscCall(PetscObjectGetId((PetscObject)v,&ctx->vid));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      20
          PetscCall(VecGetState(v,&ctx->vstate));
    
        }
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      40
        PetscCall(VecSet(r,0.0));
    
        2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.

      360
        for (i=0;i<nep->nconv;i++) {
    
        1/8✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✗ Branch 7 not taken.

      320
          if (ctx->rg) PetscCall(RGCheckInside(ctx->rg,1,&nep->eigr[i],&nep->eigi[i],&inside));
    
        1/2✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.

      320
          if (inside>=0) {
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      320
            PetscCall(BVGetColumn(nep->V,i,&x));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      320
            PetscCall(BVGetColumn(nep->W,i,&y));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      320
            PetscCall(NEPApplyJacobian(nep,nep->eigr[i],x,z,w,NULL));
    
        2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.

      320
            if (!ctx->dots_avail[i]) {
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      160
              PetscCall(VecDot(v,y,&ctx->dots[i]));
    
      160
              ctx->dots_avail[i] = PETSC_TRUE;
    
            }
    
        2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.

      320
            if (!ctx->nfactor_avail[i]) {
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      80
              PetscCall(VecDot(w,y,&ctx->nfactor[i]));
    
      80
              ctx->nfactor_avail[i] = PETSC_TRUE;
    
            }
    
      320
            alpha = ctx->dots[i]/(ctx->nfactor[i]*(ctx->omega-nep->eigr[i]));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      320
            PetscCall(VecAXPY(r,alpha,x));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      320
            PetscCall(BVRestoreColumn(nep->V,i,&x));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      320
            PetscCall(BVRestoreColumn(nep->W,i,&y));
    
          }
    
        }
    
        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.

      8
        PetscFunctionReturn(PETSC_SUCCESS);
    
      }
    
      10
      static PetscErrorCode MatDestroy_Resolvent(Mat M)
    
      {
    
      10
        NEP_RESOLVENT_MATSHELL *ctx;
    
        1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.

      10
        PetscFunctionBegin;
    
        1/2✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.

      10
        if (M) {
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      10
          PetscCall(MatShellGetContext(M,&ctx));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      10
          PetscCall(PetscFree4(ctx->nfactor,ctx->nfactor_avail,ctx->dots,ctx->dots_avail));
    
        5/8✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 8 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.

      10
          PetscCall(PetscFree(ctx));
    
        }
    
        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.

      2
        PetscFunctionReturn(PETSC_SUCCESS);
    
      }
    
      /*@
    
         NEPApplyResolvent - Applies the resolvent T^{-1}(z) to a given vector.
    
         Collective
    
         Input Parameters:
    
      +  nep   - eigensolver context obtained from NEPCreate()
    
      .  rg    - optional region
    
      .  omega - value where the resolvent must be evaluated
    
      -  v     - input vector
    
         Output Parameter:
    
      .  r     - result vector
    
         Notes:
    
         The resolvent T^{-1}(z) = sum_i (z-lambda_i)^{-1}*x_i*y_i' is evaluated at
    
         z=omega and the matrix-vector multiplication r = T^{-1}(omega)*v is computed.
    
         Vectors x_i and y_i are right and left eigenvectors, respectively, normalized
    
         so that y_i'*T'(lambda_i)*x_i=1. The sum contains only eigenvectors that have
    
         been previously computed with NEPSolve(), and if a region rg is given then only
    
         those corresponding to eigenvalues inside the region are considered.
    
         Level: intermediate
    
      .seealso: NEPGetLeftEigenvector(), NEPSolve()
    
      @*/
    
      40
      PetscErrorCode NEPApplyResolvent(NEP nep,RG rg,PetscScalar omega,Vec v,Vec r)
    
      {
    
      40
        NEP_RESOLVENT_MATSHELL *ctx;
    
        1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.

      40
        PetscFunctionBegin;
    
        3/16✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
✗ Branch 10 not taken.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✗ Branch 15 not taken.

      40
        PetscValidHeaderSpecific(nep,NEP_CLASSID,1);
    
        30/68✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✓ 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 taken 2 times.
✗ Branch 11 not taken.
✓ Branch 12 taken 2 times.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✓ Branch 15 taken 2 times.
✗ Branch 16 not taken.
✓ Branch 17 taken 2 times.
✓ Branch 18 taken 2 times.
✗ Branch 19 not taken.
✗ Branch 20 not taken.
✓ Branch 21 taken 2 times.
✓ Branch 22 taken 2 times.
✗ Branch 23 not taken.
✗ Branch 24 not taken.
✓ Branch 25 taken 2 times.
✓ Branch 26 taken 2 times.
✗ Branch 27 not taken.
✗ Branch 28 not taken.
✓ Branch 29 taken 2 times.
✓ Branch 30 taken 2 times.
✗ Branch 31 not taken.
✗ Branch 32 not taken.
✓ Branch 33 taken 2 times.
✓ Branch 34 taken 2 times.
✗ Branch 35 not taken.
✓ Branch 36 taken 2 times.
✗ Branch 37 not taken.
✗ Branch 38 not taken.
✓ Branch 39 taken 2 times.
✗ Branch 40 not taken.
✗ Branch 41 not taken.
✗ Branch 42 not taken.
✓ Branch 43 taken 2 times.
✗ Branch 44 not taken.
✗ Branch 45 not taken.
✗ Branch 46 not taken.
✓ Branch 47 taken 2 times.
✓ Branch 48 taken 2 times.
✗ Branch 49 not taken.
✗ Branch 50 not taken.
✓ Branch 51 taken 2 times.
✓ Branch 52 taken 2 times.
✗ Branch 53 not taken.
✓ Branch 54 taken 2 times.
✗ Branch 55 not taken.
✗ Branch 56 not taken.
✓ Branch 57 taken 2 times.
✗ Branch 58 not taken.
✗ Branch 59 not taken.
✓ Branch 60 taken 2 times.
✗ Branch 61 not taken.
✓ Branch 62 taken 2 times.
✗ Branch 63 not taken.
✗ Branch 64 not taken.
✓ Branch 65 taken 2 times.
✗ Branch 66 not taken.
✗ Branch 67 not taken.

      40
        PetscValidLogicalCollectiveScalar(nep,omega,3);
    
        3/16✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
✗ Branch 10 not taken.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✗ Branch 15 not taken.

      40
        PetscValidHeaderSpecific(v,VEC_CLASSID,4);
    
        3/16✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
✗ Branch 10 not taken.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✗ Branch 13 not taken.
✗ Branch 14 not taken.
✗ Branch 15 not taken.

      40
        PetscValidHeaderSpecific(r,VEC_CLASSID,5);
    
        1/4✗ Branch 0 not taken.
✓ Branch 1 taken 2 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.

      40
        NEPCheckSolved(nep,1);
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      40
        PetscCall(PetscLogEventBegin(NEP_Resolvent,nep,0,0,0));
    
        2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.

      40
        if (!nep->resolvent) {
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      10
          PetscCall(PetscNew(&ctx));
    
      10
          ctx->nep = nep;
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      10
          PetscCall(PetscCalloc4(nep->nconv,&ctx->nfactor,nep->nconv,&ctx->nfactor_avail,nep->nconv,&ctx->dots,nep->nconv,&ctx->dots_avail));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      10
          PetscCall(MatCreateShell(PetscObjectComm((PetscObject)nep),nep->nloc,nep->nloc,nep->n,nep->n,ctx,&nep->resolvent));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      10
          PetscCall(MatShellSetOperation(nep->resolvent,MATOP_MULT,(PetscErrorCodeFn*)MatMult_Resolvent));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      10
          PetscCall(MatShellSetOperation(nep->resolvent,MATOP_DESTROY,(PetscErrorCodeFn*)MatDestroy_Resolvent));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      30
        } else PetscCall(MatShellGetContext(nep->resolvent,&ctx));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      40
        PetscCall(NEPComputeVectors(nep));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      40
        PetscCall(NEPSetWorkVecs(nep,2));
    
      40
        ctx->rg    = rg;
    
      40
        ctx->omega = omega;
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      40
        PetscCall(MatMult(nep->resolvent,v,r));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      40
        PetscCall(PetscLogEventEnd(NEP_Resolvent,nep,0,0,0));
    
        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.

      8
        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			NEP routines related to resolvent T^{-1}(z) = sum_i (z-lambda_i)^{-1} x_i y_i'
12			*/
13
14			#include <slepc/private/nepimpl.h> /I "slepcnep.h" I/
15
16			typedef struct {
17			NEP nep;
18			RG rg;
19			PetscScalar omega;
20			PetscScalar nfactor; / normalization factors y_i'T'(lambda_i)x_i */
21			PetscBool *nfactor_avail;
22			PetscScalar dots; / inner products y_i'v /
23			PetscBool *dots_avail;
24			PetscObjectId vid;
25			PetscObjectState vstate;
26			} NEP_RESOLVENT_MATSHELL;
27
28		40	static PetscErrorCode MatMult_Resolvent(Mat M,Vec v,Vec r)
29			{
30		40	NEP_RESOLVENT_MATSHELL *ctx;
31		40	NEP nep;
32		40	PetscInt i,inside=1;
33		40	PetscScalar alpha;
34		40	Vec x,y,z,w;
35
36	1/2 ✓ Branch 0 taken 2 times. ✗ Branch 1 not taken.	40	PetscFunctionBegin;
37	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	40	PetscCall(MatShellGetContext(M,&ctx));
38		40	nep = ctx->nep;
39		40	w = nep->work[0];
40		40	z = nep->work[1];
41	4/4 ✓ Branch 0 taken 10 times. ✓ Branch 1 taken 10 times. ✓ Branch 2 taken 10 times. ✓ Branch 3 taken 10 times.	40	if (((PetscObject)v)->id != ctx->vid \|\| ((PetscObject)v)->state != ctx->vstate) {
42	4/6 ✓ Branch 0 taken 5 times. ✓ Branch 1 taken 5 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	20	PetscCall(PetscArrayzero(ctx->dots_avail,ctx->nep->nconv));
43	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	20	PetscCall(PetscObjectGetId((PetscObject)v,&ctx->vid));
44	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	20	PetscCall(VecGetState(v,&ctx->vstate));
45			}
46	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	40	PetscCall(VecSet(r,0.0));
47	2/2 ✓ Branch 0 taken 10 times. ✓ Branch 1 taken 10 times.	360	for (i=0;i<nep->nconv;i++) {
48	1/8 ✗ Branch 0 not taken. ✓ Branch 1 taken 10 times. ✗ Branch 2 not taken. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✗ Branch 5 not taken. ✗ Branch 6 not taken. ✗ Branch 7 not taken.	320	if (ctx->rg) PetscCall(RGCheckInside(ctx->rg,1,&nep->eigr[i],&nep->eigi[i],&inside));
49	1/2 ✓ Branch 0 taken 10 times. ✗ Branch 1 not taken.	320	if (inside>=0) {
50	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	320	PetscCall(BVGetColumn(nep->V,i,&x));
51	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	320	PetscCall(BVGetColumn(nep->W,i,&y));
52	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	320	PetscCall(NEPApplyJacobian(nep,nep->eigr[i],x,z,w,NULL));
53	2/2 ✓ Branch 0 taken 10 times. ✓ Branch 1 taken 10 times.	320	if (!ctx->dots_avail[i]) {
54	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	160	PetscCall(VecDot(v,y,&ctx->dots[i]));
55		160	ctx->dots_avail[i] = PETSC_TRUE;
56			}
57	2/2 ✓ Branch 0 taken 10 times. ✓ Branch 1 taken 10 times.	320	if (!ctx->nfactor_avail[i]) {
58	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	80	PetscCall(VecDot(w,y,&ctx->nfactor[i]));
59		80	ctx->nfactor_avail[i] = PETSC_TRUE;
60			}
61		320	alpha = ctx->dots[i]/(ctx->nfactor[i]*(ctx->omega-nep->eigr[i]));
62	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	320	PetscCall(VecAXPY(r,alpha,x));
63	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	320	PetscCall(BVRestoreColumn(nep->V,i,&x));
64	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	320	PetscCall(BVRestoreColumn(nep->W,i,&y));
65			}
66			}
67	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.	8	PetscFunctionReturn(PETSC_SUCCESS);
68			}
69
70		10	static PetscErrorCode MatDestroy_Resolvent(Mat M)
71			{
72		10	NEP_RESOLVENT_MATSHELL *ctx;
73
74	1/2 ✓ Branch 0 taken 2 times. ✗ Branch 1 not taken.	10	PetscFunctionBegin;
75	1/2 ✓ Branch 0 taken 10 times. ✗ Branch 1 not taken.	10	if (M) {
76	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	10	PetscCall(MatShellGetContext(M,&ctx));
77	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	10	PetscCall(PetscFree4(ctx->nfactor,ctx->nfactor_avail,ctx->dots,ctx->dots_avail));
78	5/8 ✓ Branch 0 taken 10 times. ✗ Branch 1 not taken. ✓ Branch 2 taken 2 times. ✓ Branch 3 taken 8 times. ✓ Branch 4 taken 2 times. ✗ Branch 5 not taken. ✗ Branch 6 not taken. ✓ Branch 7 taken 2 times.	10	PetscCall(PetscFree(ctx));
79			}
80	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.	2	PetscFunctionReturn(PETSC_SUCCESS);
81			}
82
83			/*@
84			NEPApplyResolvent - Applies the resolvent T^{-1}(z) to a given vector.
85
86			Collective
87
88			Input Parameters:
89			+ nep - eigensolver context obtained from NEPCreate()
90			. rg - optional region
91			. omega - value where the resolvent must be evaluated
92			- v - input vector
93
94			Output Parameter:
95			. r - result vector
96
97			Notes:
98			The resolvent T^{-1}(z) = sum_i (z-lambda_i)^{-1}x_iy_i' is evaluated at
99			z=omega and the matrix-vector multiplication r = T^{-1}(omega)*v is computed.
100			Vectors x_i and y_i are right and left eigenvectors, respectively, normalized
101			so that y_i'T'(lambda_i)x_i=1. The sum contains only eigenvectors that have
102			been previously computed with NEPSolve(), and if a region rg is given then only
103			those corresponding to eigenvalues inside the region are considered.
104
105			Level: intermediate
106
107			.seealso: NEPGetLeftEigenvector(), NEPSolve()
108			@*/
109		40	PetscErrorCode NEPApplyResolvent(NEP nep,RG rg,PetscScalar omega,Vec v,Vec r)
110			{
111		40	NEP_RESOLVENT_MATSHELL *ctx;
112
113	1/2 ✓ Branch 0 taken 2 times. ✗ Branch 1 not taken.	40	PetscFunctionBegin;
114	3/16 ✗ Branch 0 not taken. ✓ Branch 1 taken 2 times. ✗ Branch 2 not taken. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times. ✗ Branch 6 not taken. ✗ Branch 7 not taken. ✗ Branch 8 not taken. ✓ Branch 9 taken 2 times. ✗ Branch 10 not taken. ✗ Branch 11 not taken. ✗ Branch 12 not taken. ✗ Branch 13 not taken. ✗ Branch 14 not taken. ✗ Branch 15 not taken.	40	PetscValidHeaderSpecific(nep,NEP_CLASSID,1);
115	30/68 ✗ Branch 0 not taken. ✓ Branch 1 taken 2 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✓ 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 taken 2 times. ✗ Branch 11 not taken. ✓ Branch 12 taken 2 times. ✗ Branch 13 not taken. ✗ Branch 14 not taken. ✓ Branch 15 taken 2 times. ✗ Branch 16 not taken. ✓ Branch 17 taken 2 times. ✓ Branch 18 taken 2 times. ✗ Branch 19 not taken. ✗ Branch 20 not taken. ✓ Branch 21 taken 2 times. ✓ Branch 22 taken 2 times. ✗ Branch 23 not taken. ✗ Branch 24 not taken. ✓ Branch 25 taken 2 times. ✓ Branch 26 taken 2 times. ✗ Branch 27 not taken. ✗ Branch 28 not taken. ✓ Branch 29 taken 2 times. ✓ Branch 30 taken 2 times. ✗ Branch 31 not taken. ✗ Branch 32 not taken. ✓ Branch 33 taken 2 times. ✓ Branch 34 taken 2 times. ✗ Branch 35 not taken. ✓ Branch 36 taken 2 times. ✗ Branch 37 not taken. ✗ Branch 38 not taken. ✓ Branch 39 taken 2 times. ✗ Branch 40 not taken. ✗ Branch 41 not taken. ✗ Branch 42 not taken. ✓ Branch 43 taken 2 times. ✗ Branch 44 not taken. ✗ Branch 45 not taken. ✗ Branch 46 not taken. ✓ Branch 47 taken 2 times. ✓ Branch 48 taken 2 times. ✗ Branch 49 not taken. ✗ Branch 50 not taken. ✓ Branch 51 taken 2 times. ✓ Branch 52 taken 2 times. ✗ Branch 53 not taken. ✓ Branch 54 taken 2 times. ✗ Branch 55 not taken. ✗ Branch 56 not taken. ✓ Branch 57 taken 2 times. ✗ Branch 58 not taken. ✗ Branch 59 not taken. ✓ Branch 60 taken 2 times. ✗ Branch 61 not taken. ✓ Branch 62 taken 2 times. ✗ Branch 63 not taken. ✗ Branch 64 not taken. ✓ Branch 65 taken 2 times. ✗ Branch 66 not taken. ✗ Branch 67 not taken.	40	PetscValidLogicalCollectiveScalar(nep,omega,3);
116	3/16 ✗ Branch 0 not taken. ✓ Branch 1 taken 2 times. ✗ Branch 2 not taken. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times. ✗ Branch 6 not taken. ✗ Branch 7 not taken. ✗ Branch 8 not taken. ✓ Branch 9 taken 2 times. ✗ Branch 10 not taken. ✗ Branch 11 not taken. ✗ Branch 12 not taken. ✗ Branch 13 not taken. ✗ Branch 14 not taken. ✗ Branch 15 not taken.	40	PetscValidHeaderSpecific(v,VEC_CLASSID,4);
117	3/16 ✗ Branch 0 not taken. ✓ Branch 1 taken 2 times. ✗ Branch 2 not taken. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times. ✗ Branch 6 not taken. ✗ Branch 7 not taken. ✗ Branch 8 not taken. ✓ Branch 9 taken 2 times. ✗ Branch 10 not taken. ✗ Branch 11 not taken. ✗ Branch 12 not taken. ✗ Branch 13 not taken. ✗ Branch 14 not taken. ✗ Branch 15 not taken.	40	PetscValidHeaderSpecific(r,VEC_CLASSID,5);
118	1/4 ✗ Branch 0 not taken. ✓ Branch 1 taken 2 times. ✗ Branch 2 not taken. ✗ Branch 3 not taken.	40	NEPCheckSolved(nep,1);
119
120	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	40	PetscCall(PetscLogEventBegin(NEP_Resolvent,nep,0,0,0));
121	2/2 ✓ Branch 0 taken 10 times. ✓ Branch 1 taken 10 times.	40	if (!nep->resolvent) {
122	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	10	PetscCall(PetscNew(&ctx));
123		10	ctx->nep = nep;
124	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	10	PetscCall(PetscCalloc4(nep->nconv,&ctx->nfactor,nep->nconv,&ctx->nfactor_avail,nep->nconv,&ctx->dots,nep->nconv,&ctx->dots_avail));
125	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	10	PetscCall(MatCreateShell(PetscObjectComm((PetscObject)nep),nep->nloc,nep->nloc,nep->n,nep->n,ctx,&nep->resolvent));
126	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	10	PetscCall(MatShellSetOperation(nep->resolvent,MATOP_MULT,(PetscErrorCodeFn*)MatMult_Resolvent));
127	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	10	PetscCall(MatShellSetOperation(nep->resolvent,MATOP_DESTROY,(PetscErrorCodeFn*)MatDestroy_Resolvent));
128	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	30	} else PetscCall(MatShellGetContext(nep->resolvent,&ctx));
129	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	40	PetscCall(NEPComputeVectors(nep));
130	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	40	PetscCall(NEPSetWorkVecs(nep,2));
131		40	ctx->rg = rg;
132		40	ctx->omega = omega;
133	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	40	PetscCall(MatMult(nep->resolvent,v,r));
134	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 8 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	40	PetscCall(PetscLogEventEnd(NEP_Resolvent,nep,0,0,0));
135	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.	8	PetscFunctionReturn(PETSC_SUCCESS);
136			}
137

Function (Line)	Call count	Block coverage
MatDestroy_Resolvent (line 70)	called 10 times	81.0%
MatMult_Resolvent (line 28)	called 40 times	79.0%
NEPApplyResolvent (line 109)	called 40 times	69.0%