GCC Code Coverage Report

Directory:	./
File:	src/sys/classes/ds/tests/test12.c
Date:	2025-10-03 04:28:47

	Exec	Total	Coverage
Lines:	117	121	96.7%
Functions:	1	1	100.0%
Branches:	287	466	61.6%

  
      Line
      Branch
      Exec
      Source
    
      /*
    
         - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
    
         SLEPc - Scalable Library for Eigenvalue Problem Computations
    
         Copyright (c) 2002-, Universitat Politecnica de Valencia, Spain
    
         This file is part of SLEPc.
    
         SLEPc is distributed under a 2-clause BSD license (see LICENSE).
    
         - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
    
      */
    
      static char help[] = "Test DSNEP.\n\n";
    
      #include <slepcds.h>
    
      15
      int main(int argc,char **argv)
    
      {
    
      15
        DS             ds;
    
      15
        FN             f1,f2,f3,funs[3],qfun;
    
      15
        SlepcSC        sc;
    
      15
        PetscScalar    *Id,*A,*B,*wr,*wi,*X,*W,coeffs[2],auxr,alpha;
    
      15
        PetscReal      tol,tau=0.001,radius=10,h,a=20,xi,re,im,nrm,aux;
    
      15
        PetscInt       i,j,ii,jj,k,n=10,ld,nev,nfun;
    
      15
        PetscViewer    viewer;
    
      15
        PetscBool      verbose;
    
      15
        RG             rg;
    
      15
        DSMatType      mat[3]={DS_MAT_E0,DS_MAT_E1,DS_MAT_E2};
    
      #if !defined(PETSC_USE_COMPLEX)
    
      5
        PetscScalar    auxi;
    
      #endif
    
        1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.

      15
        PetscFunctionBeginUser;
    
        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.

      15
        PetscCall(SlepcInitialize(&argc,&argv,NULL,help));
    
        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.

      15
        PetscCall(PetscOptionsGetInt(NULL,NULL,"-n",&n,NULL));
    
        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.

      15
        PetscCall(PetscOptionsGetReal(NULL,NULL,"-tau",&tau,NULL));
    
        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.

      15
        PetscCall(PetscPrintf(PETSC_COMM_WORLD,"Solve a Dense System of type NEP - dimension %" PetscInt_FMT ", tau=%g.\n",n,(double)tau));
    
        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.

      15
        PetscCall(PetscOptionsHasName(NULL,NULL,"-verbose",&verbose));
    
        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.

      15
        PetscCall(PetscOptionsGetReal(NULL,NULL,"-radius",&radius,NULL));
    
        /* Create DS object */
    
        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.

      15
        PetscCall(DSCreate(PETSC_COMM_WORLD,&ds));
    
        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.

      15
        PetscCall(DSSetType(ds,DSNEP));
    
      15
        tol  = 1000*n*PETSC_MACHINE_EPSILON;
    
        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.

      15
        PetscCall(DSNEPSetRefine(ds,tol,PETSC_DECIDE));
    
        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.

      15
        PetscCall(DSSetFromOptions(ds));
    
        /* Set functions (prior to DSAllocate) */
    
        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.

      15
        PetscCall(FNCreate(PETSC_COMM_WORLD,&f1));
    
        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.

      15
        PetscCall(FNSetType(f1,FNRATIONAL));
    
      15
        coeffs[0] = -1.0; coeffs[1] = 0.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.

      15
        PetscCall(FNRationalSetNumerator(f1,2,coeffs));
    
        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.

      15
        PetscCall(FNCreate(PETSC_COMM_WORLD,&f2));
    
        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.

      15
        PetscCall(FNSetType(f2,FNRATIONAL));
    
      15
        coeffs[0] = 1.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.

      15
        PetscCall(FNRationalSetNumerator(f2,1,coeffs));
    
        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.

      15
        PetscCall(FNCreate(PETSC_COMM_WORLD,&f3));
    
        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.

      15
        PetscCall(FNSetType(f3,FNEXP));
    
        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.

      15
        PetscCall(FNSetScale(f3,-tau,1.0));
    
      15
        funs[0] = f1;
    
      15
        funs[1] = f2;
    
      15
        funs[2] = f3;
    
        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.

      15
        PetscCall(DSNEPSetFN(ds,3,funs));
    
        /* Set dimensions */
    
      15
        ld = n+2;  /* test leading dimension larger than n */
    
        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.

      15
        PetscCall(DSAllocate(ds,ld));
    
        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.

      15
        PetscCall(DSSetDimensions(ds,n,0,0));
    
        /* Set region (used only in method=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.

      15
        PetscCall(RGCreate(PETSC_COMM_WORLD,&rg));
    
        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.

      15
        PetscCall(RGSetType(rg,RGELLIPSE));
    
        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.

      15
        PetscCall(RGEllipseSetParameters(rg,0.0,radius,1.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.

      15
        PetscCall(DSNEPSetRG(ds,rg));
    
        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.

      15
        PetscCall(RGDestroy(&rg));
    
        /* Set up viewer */
    
        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.

      15
        PetscCall(PetscViewerASCIIGetStdout(PETSC_COMM_WORLD,&viewer));
    
        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.

      15
        PetscCall(PetscViewerPushFormat(viewer,PETSC_VIEWER_ASCII_INFO_DETAIL));
    
        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.

      15
        PetscCall(DSView(ds,viewer));
    
        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.

      15
        PetscCall(PetscViewerPopFormat(viewer));
    
        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.

      15
        if (verbose) PetscCall(PetscViewerPushFormat(viewer,PETSC_VIEWER_ASCII_MATLAB));
    
        /* Show info about functions */
    
        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.

      15
        PetscCall(DSNEPGetNumFN(ds,&nfun));
    
        2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.

      60
        for (i=0;i<nfun;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.

      45
          PetscCall(PetscPrintf(PETSC_COMM_WORLD,"Function %" PetscInt_FMT ":\n",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.

      45
          PetscCall(DSNEPGetFN(ds,i,&qfun));
    
        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.

      45
          PetscCall(FNView(qfun,NULL));
    
        }
    
        /* Fill matrices */
    
        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.

      15
        PetscCall(DSGetArray(ds,DS_MAT_E0,&Id));
    
        2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.

      165
        for (i=0;i<n;i++) Id[i+i*ld]=1.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.

      15
        PetscCall(DSRestoreArray(ds,DS_MAT_E0,&Id));
    
      15
        h = PETSC_PI/(PetscReal)(n+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.

      15
        PetscCall(DSGetArray(ds,DS_MAT_E1,&A));
    
        2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.

      165
        for (i=0;i<n;i++) A[i+i*ld]=-2.0/(h*h)+a;
    
        2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.

      150
        for (i=1;i<n;i++) {
    
      135
          A[i+(i-1)*ld]=1.0/(h*h);
    
      135
          A[(i-1)+i*ld]=1.0/(h*h);
    
        }
    
        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.

      15
        PetscCall(DSRestoreArray(ds,DS_MAT_E1,&A));
    
        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.

      15
        PetscCall(DSGetArray(ds,DS_MAT_E2,&B));
    
        2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.

      165
        for (i=0;i<n;i++) {
    
      150
          xi = (i+1)*h;
    
      150
          B[i+i*ld] = -4.1+xi*(1.0-PetscExpReal(xi-PETSC_PI));
    
        }
    
        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.

      15
        PetscCall(DSRestoreArray(ds,DS_MAT_E2,&B));
    
        1/2✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.

      15
        if (verbose) {
    
      ✗
          PetscCall(PetscPrintf(PETSC_COMM_WORLD,"Initial - - - - - - - - -\n"));
    
      ✗
          PetscCall(DSView(ds,viewer));
    
        }
    
        /* Solve */
    
        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.

      15
        PetscCall(PetscCalloc2(n,&wr,n,&wi));
    
        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.

      15
        PetscCall(DSGetSlepcSC(ds,&sc));
    
      15
        sc->comparison    = SlepcCompareLargestMagnitude;
    
      15
        sc->comparisonctx = NULL;
    
      15
        sc->map           = NULL;
    
      15
        sc->mapobj        = NULL;
    
        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.

      15
        PetscCall(DSSolve(ds,wr,wi));
    
        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.

      15
        PetscCall(DSSort(ds,wr,wi,NULL,NULL,NULL));
    
        1/2✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.

      15
        if (verbose) {
    
      ✗
          PetscCall(PetscPrintf(PETSC_COMM_WORLD,"After solve - - - - - - - - -\n"));
    
      ✗
          PetscCall(DSView(ds,viewer));
    
        }
    
        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.

      15
        PetscCall(DSGetDimensions(ds,NULL,NULL,NULL,&nev));
    
        /* Print computed eigenvalues */
    
        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.

      15
        PetscCall(PetscMalloc1(ld*ld,&W));
    
        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.

      15
        PetscCall(DSVectors(ds,DS_MAT_X,NULL,NULL));
    
        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.

      15
        PetscCall(DSGetArray(ds,DS_MAT_X,&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.

      15
        PetscCall(PetscPrintf(PETSC_COMM_WORLD,"Computed eigenvalues =\n"));
    
        2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.

      40
        for (i=0;i<nev;i++) {
    
      #if defined(PETSC_USE_COMPLEX)
    
      20
          re = PetscRealPart(wr[i]);
    
      20
          im = PetscImaginaryPart(wr[i]);
    
      #else
    
      5
          re = wr[i];
    
      5
          im = wi[i];
    
      #endif
    
          /* Residual */
    
        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.

      25
          PetscCall(PetscArrayzero(W,ld*ld));
    
        2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.

      100
          for (k=0;k<nfun;k++) {
    
        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.

      75
            PetscCall(FNEvaluateFunction(funs[k],wr[i],&alpha));
    
        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.

      75
            PetscCall(DSGetArray(ds,mat[k],&A));
    
        4/4✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ Branch 2 taken 10 times.
✓ Branch 3 taken 10 times.

      8325
            for (jj=0;jj<n;jj++) for (ii=0;ii<n;ii++) W[jj*ld+ii] += alpha*A[jj*ld+ii];
    
        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.

      75
            PetscCall(DSRestoreArray(ds,mat[k],&A));
    
          }
    
          nrm = 0.0;
    
        2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.

      275
          for (k=0;k<n;k++) {
    
            auxr = 0.0;
    
      #if !defined(PETSC_USE_COMPLEX)
    
            auxi = 0.0;
    
      #endif
    
        2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.

      2750
            for (j=0;j<n;j++) {
    
      2500
              auxr += W[k+j*ld]*X[i*ld+j];
    
      #if !defined(PETSC_USE_COMPLEX)
    
        1/2✗ Branch 0 not taken.
✓ Branch 1 taken 5 times.

      500
              if (PetscAbs(wi[j])!=0.0) auxi += W[k+j*ld]*X[(i+1)*ld+j];
    
      #endif
    
            }
    
      250
            aux = SlepcAbsEigenvalue(auxr,auxi);
    
      250
            nrm += aux*aux;
    
          }
    
      25
          nrm = PetscSqrtReal(nrm);
    
        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.

      25
          if (nrm/SlepcAbsEigenvalue(wr[i],wi[i])>tol) PetscCall(PetscPrintf(PETSC_COMM_WORLD,"Warning: the residual norm of the %" PetscInt_FMT "-th computed eigenpair %g\n",i,(double)nrm));
    
        7/10✓ Branch 0 taken 10 times.
✓ Branch 1 taken 5 times.
✓ Branch 2 taken 10 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 2 times.
✓ Branch 5 taken 8 times.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.

      25
          if (PetscAbs(im)<1e-10) PetscCall(PetscViewerASCIIPrintf(viewer,"  %.5f\n",(double)re));
    
        0/6✗ Branch 0 not taken.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.

      25
          else PetscCall(PetscViewerASCIIPrintf(viewer,"  %.5f%+.5fi\n",(double)re,(double)im));
    
        }
    
        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.

      15
        PetscCall(PetscFree(W));
    
        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.

      15
        PetscCall(DSRestoreArray(ds,DS_MAT_X,&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.

      15
        PetscCall(DSRestoreArray(ds,DS_MAT_W,&W));
    
        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.

      15
        PetscCall(PetscFree2(wr,wi));
    
        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.

      15
        PetscCall(FNDestroy(&f1));
    
        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.

      15
        PetscCall(FNDestroy(&f2));
    
        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.

      15
        PetscCall(FNDestroy(&f3));
    
        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.

      15
        PetscCall(DSDestroy(&ds));
    
        3/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 2 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.

      15
        PetscCall(SlepcFinalize());
    
        return 0;
    
      }
    
      /*TEST
    
         testset:
    
            test:
    
               filter: grep -v "solving the problem"
    
               suffix: 1
    
            test:
    
               suffix: 2
    
               args: -ds_method 1 -radius 10 -ds_nep_refine_its 1
    
               filter: grep -v "solving the problem" | sed -e "s/[+-]0\.0*i//g" | sed -e "s/37411/37410/" | sed -e "s/tolerance [0-9]\.[0-9]*e[+-]\([0-9]*\)/tolerance removed/" | sed -e "s/tolerance [0-9]\.\([0-9]*\)/tolerance removed/"
    
               requires: complex
    
      TEST*/

Line	Branch	Exec	Source
1			/*
2			- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
3			SLEPc - Scalable Library for Eigenvalue Problem Computations
4			Copyright (c) 2002-, Universitat Politecnica de Valencia, Spain
5
6			This file is part of SLEPc.
7			SLEPc is distributed under a 2-clause BSD license (see LICENSE).
8			- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
9			*/
10
11			static char help[] = "Test DSNEP.\n\n";
12
13			#include <slepcds.h>
14
15		15	int main(int argc,char **argv)
16			{
17		15	DS ds;
18		15	FN f1,f2,f3,funs[3],qfun;
19		15	SlepcSC sc;
20		15	PetscScalar Id,A,B,wr,wi,X,*W,coeffs[2],auxr,alpha;
21		15	PetscReal tol,tau=0.001,radius=10,h,a=20,xi,re,im,nrm,aux;
22		15	PetscInt i,j,ii,jj,k,n=10,ld,nev,nfun;
23		15	PetscViewer viewer;
24		15	PetscBool verbose;
25		15	RG rg;
26		15	DSMatType mat[3]={DS_MAT_E0,DS_MAT_E1,DS_MAT_E2};
27			#if !defined(PETSC_USE_COMPLEX)
28		5	PetscScalar auxi;
29			#endif
30
31	1/2 ✓ Branch 0 taken 2 times. ✗ Branch 1 not taken.	15	PetscFunctionBeginUser;
32	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.	15	PetscCall(SlepcInitialize(&argc,&argv,NULL,help));
33	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.	15	PetscCall(PetscOptionsGetInt(NULL,NULL,"-n",&n,NULL));
34	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.	15	PetscCall(PetscOptionsGetReal(NULL,NULL,"-tau",&tau,NULL));
35	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.	15	PetscCall(PetscPrintf(PETSC_COMM_WORLD,"Solve a Dense System of type NEP - dimension %" PetscInt_FMT ", tau=%g.\n",n,(double)tau));
36	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.	15	PetscCall(PetscOptionsHasName(NULL,NULL,"-verbose",&verbose));
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.	15	PetscCall(PetscOptionsGetReal(NULL,NULL,"-radius",&radius,NULL));
38
39			/* Create DS object */
40	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.	15	PetscCall(DSCreate(PETSC_COMM_WORLD,&ds));
41	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.	15	PetscCall(DSSetType(ds,DSNEP));
42		15	tol = 1000nPETSC_MACHINE_EPSILON;
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.	15	PetscCall(DSNEPSetRefine(ds,tol,PETSC_DECIDE));
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.	15	PetscCall(DSSetFromOptions(ds));
45
46			/* Set functions (prior to DSAllocate) */
47	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.	15	PetscCall(FNCreate(PETSC_COMM_WORLD,&f1));
48	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.	15	PetscCall(FNSetType(f1,FNRATIONAL));
49		15	coeffs[0] = -1.0; coeffs[1] = 0.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.	15	PetscCall(FNRationalSetNumerator(f1,2,coeffs));
51
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.	15	PetscCall(FNCreate(PETSC_COMM_WORLD,&f2));
53	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.	15	PetscCall(FNSetType(f2,FNRATIONAL));
54		15	coeffs[0] = 1.0;
55	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.	15	PetscCall(FNRationalSetNumerator(f2,1,coeffs));
56
57	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.	15	PetscCall(FNCreate(PETSC_COMM_WORLD,&f3));
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.	15	PetscCall(FNSetType(f3,FNEXP));
59	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.	15	PetscCall(FNSetScale(f3,-tau,1.0));
60
61		15	funs[0] = f1;
62		15	funs[1] = f2;
63		15	funs[2] = f3;
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.	15	PetscCall(DSNEPSetFN(ds,3,funs));
65
66			/* Set dimensions */
67		15	ld = n+2; /* test leading dimension larger than n */
68	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.	15	PetscCall(DSAllocate(ds,ld));
69	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.	15	PetscCall(DSSetDimensions(ds,n,0,0));
70
71			/* Set region (used only in method=1) */
72	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.	15	PetscCall(RGCreate(PETSC_COMM_WORLD,&rg));
73	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.	15	PetscCall(RGSetType(rg,RGELLIPSE));
74	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.	15	PetscCall(RGEllipseSetParameters(rg,0.0,radius,1.0));
75	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.	15	PetscCall(DSNEPSetRG(ds,rg));
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.	15	PetscCall(RGDestroy(&rg));
77
78			/* Set up viewer */
79	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.	15	PetscCall(PetscViewerASCIIGetStdout(PETSC_COMM_WORLD,&viewer));
80	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.	15	PetscCall(PetscViewerPushFormat(viewer,PETSC_VIEWER_ASCII_INFO_DETAIL));
81	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.	15	PetscCall(DSView(ds,viewer));
82	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.	15	PetscCall(PetscViewerPopFormat(viewer));
83	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.	15	if (verbose) PetscCall(PetscViewerPushFormat(viewer,PETSC_VIEWER_ASCII_MATLAB));
84
85			/* Show info about functions */
86	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.	15	PetscCall(DSNEPGetNumFN(ds,&nfun));
87	2/2 ✓ Branch 0 taken 10 times. ✓ Branch 1 taken 10 times.	60	for (i=0;i<nfun;i++) {
88	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.	45	PetscCall(PetscPrintf(PETSC_COMM_WORLD,"Function %" PetscInt_FMT ":\n",i));
89	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.	45	PetscCall(DSNEPGetFN(ds,i,&qfun));
90	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.	45	PetscCall(FNView(qfun,NULL));
91			}
92
93			/* Fill matrices */
94	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.	15	PetscCall(DSGetArray(ds,DS_MAT_E0,&Id));
95	2/2 ✓ Branch 0 taken 10 times. ✓ Branch 1 taken 10 times.	165	for (i=0;i<n;i++) Id[i+i*ld]=1.0;
96	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.	15	PetscCall(DSRestoreArray(ds,DS_MAT_E0,&Id));
97		15	h = PETSC_PI/(PetscReal)(n+1);
98	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.	15	PetscCall(DSGetArray(ds,DS_MAT_E1,&A));
99	2/2 ✓ Branch 0 taken 10 times. ✓ Branch 1 taken 10 times.	165	for (i=0;i<n;i++) A[i+ild]=-2.0/(hh)+a;
100	2/2 ✓ Branch 0 taken 10 times. ✓ Branch 1 taken 10 times.	150	for (i=1;i<n;i++) {
101		135	A[i+(i-1)ld]=1.0/(hh);
102		135	A[(i-1)+ild]=1.0/(hh);
103			}
104	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.	15	PetscCall(DSRestoreArray(ds,DS_MAT_E1,&A));
105	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.	15	PetscCall(DSGetArray(ds,DS_MAT_E2,&B));
106	2/2 ✓ Branch 0 taken 10 times. ✓ Branch 1 taken 10 times.	165	for (i=0;i<n;i++) {
107		150	xi = (i+1)*h;
108		150	B[i+ild] = -4.1+xi(1.0-PetscExpReal(xi-PETSC_PI));
109			}
110	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.	15	PetscCall(DSRestoreArray(ds,DS_MAT_E2,&B));
111
112	1/2 ✗ Branch 0 not taken. ✓ Branch 1 taken 10 times.	15	if (verbose) {
113		✗	PetscCall(PetscPrintf(PETSC_COMM_WORLD,"Initial - - - - - - - - -\n"));
114		✗	PetscCall(DSView(ds,viewer));
115			}
116
117			/* Solve */
118	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.	15	PetscCall(PetscCalloc2(n,&wr,n,&wi));
119	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.	15	PetscCall(DSGetSlepcSC(ds,&sc));
120		15	sc->comparison = SlepcCompareLargestMagnitude;
121		15	sc->comparisonctx = NULL;
122		15	sc->map = NULL;
123		15	sc->mapobj = NULL;
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.	15	PetscCall(DSSolve(ds,wr,wi));
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.	15	PetscCall(DSSort(ds,wr,wi,NULL,NULL,NULL));
126
127	1/2 ✗ Branch 0 not taken. ✓ Branch 1 taken 10 times.	15	if (verbose) {
128		✗	PetscCall(PetscPrintf(PETSC_COMM_WORLD,"After solve - - - - - - - - -\n"));
129		✗	PetscCall(DSView(ds,viewer));
130			}
131	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.	15	PetscCall(DSGetDimensions(ds,NULL,NULL,NULL,&nev));
132
133			/* Print computed eigenvalues */
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.	15	PetscCall(PetscMalloc1(ld*ld,&W));
135	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.	15	PetscCall(DSVectors(ds,DS_MAT_X,NULL,NULL));
136	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.	15	PetscCall(DSGetArray(ds,DS_MAT_X,&X));
137	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.	15	PetscCall(PetscPrintf(PETSC_COMM_WORLD,"Computed eigenvalues =\n"));
138	2/2 ✓ Branch 0 taken 10 times. ✓ Branch 1 taken 10 times.	40	for (i=0;i<nev;i++) {
139			#if defined(PETSC_USE_COMPLEX)
140		20	re = PetscRealPart(wr[i]);
141		20	im = PetscImaginaryPart(wr[i]);
142			#else
143		5	re = wr[i];
144		5	im = wi[i];
145			#endif
146			/* Residual */
147	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.	25	PetscCall(PetscArrayzero(W,ld*ld));
148	2/2 ✓ Branch 0 taken 10 times. ✓ Branch 1 taken 10 times.	100	for (k=0;k<nfun;k++) {
149	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.	75	PetscCall(FNEvaluateFunction(funs[k],wr[i],&alpha));
150	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.	75	PetscCall(DSGetArray(ds,mat[k],&A));
151	4/4 ✓ Branch 0 taken 10 times. ✓ Branch 1 taken 10 times. ✓ Branch 2 taken 10 times. ✓ Branch 3 taken 10 times.	8325	for (jj=0;jj<n;jj++) for (ii=0;ii<n;ii++) W[jjld+ii] += alphaA[jj*ld+ii];
152	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.	75	PetscCall(DSRestoreArray(ds,mat[k],&A));
153			}
154			nrm = 0.0;
155	2/2 ✓ Branch 0 taken 10 times. ✓ Branch 1 taken 10 times.	275	for (k=0;k<n;k++) {
156			auxr = 0.0;
157			#if !defined(PETSC_USE_COMPLEX)
158			auxi = 0.0;
159			#endif
160	2/2 ✓ Branch 0 taken 10 times. ✓ Branch 1 taken 10 times.	2750	for (j=0;j<n;j++) {
161		2500	auxr += W[k+jld]X[i*ld+j];
162			#if !defined(PETSC_USE_COMPLEX)
163	1/2 ✗ Branch 0 not taken. ✓ Branch 1 taken 5 times.	500	if (PetscAbs(wi[j])!=0.0) auxi += W[k+jld]X[(i+1)*ld+j];
164			#endif
165			}
166		250	aux = SlepcAbsEigenvalue(auxr,auxi);
167		250	nrm += aux*aux;
168			}
169		25	nrm = PetscSqrtReal(nrm);
170	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.	25	if (nrm/SlepcAbsEigenvalue(wr[i],wi[i])>tol) PetscCall(PetscPrintf(PETSC_COMM_WORLD,"Warning: the residual norm of the %" PetscInt_FMT "-th computed eigenpair %g\n",i,(double)nrm));
171	7/10 ✓ Branch 0 taken 10 times. ✓ Branch 1 taken 5 times. ✓ Branch 2 taken 10 times. ✗ Branch 3 not taken. ✓ Branch 4 taken 2 times. ✓ Branch 5 taken 8 times. ✓ Branch 6 taken 2 times. ✗ Branch 7 not taken. ✗ Branch 8 not taken. ✓ Branch 9 taken 2 times.	25	if (PetscAbs(im)<1e-10) PetscCall(PetscViewerASCIIPrintf(viewer," %.5f\n",(double)re));
172	0/6 ✗ Branch 0 not taken. ✗ Branch 1 not taken. ✗ Branch 2 not taken. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✗ Branch 5 not taken.	25	else PetscCall(PetscViewerASCIIPrintf(viewer," %.5f%+.5fi\n",(double)re,(double)im));
173			}
174	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.	15	PetscCall(PetscFree(W));
175	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.	15	PetscCall(DSRestoreArray(ds,DS_MAT_X,&X));
176	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.	15	PetscCall(DSRestoreArray(ds,DS_MAT_W,&W));
177	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.	15	PetscCall(PetscFree2(wr,wi));
178	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.	15	PetscCall(FNDestroy(&f1));
179	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.	15	PetscCall(FNDestroy(&f2));
180	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.	15	PetscCall(FNDestroy(&f3));
181	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.	15	PetscCall(DSDestroy(&ds));
182	3/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 2 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✗ Branch 5 not taken.	15	PetscCall(SlepcFinalize());
183			return 0;
184			}
185
186			/*TEST
187
188			testset:
189			test:
190			filter: grep -v "solving the problem"
191			suffix: 1
192			test:
193			suffix: 2
194			args: -ds_method 1 -radius 10 -ds_nep_refine_its 1
195			filter: grep -v "solving the problem" \| sed -e "s/[+-]0\.0i//g" \| sed -e "s/37411/37410/" \| sed -e "s/tolerance [0-9]\.[0-9]e[+-]\([0-9]\)/tolerance removed/" \| sed -e "s/tolerance [0-9]\.\([0-9]\)/tolerance removed/"
196			requires: complex
197
198			TEST*/
199