GCC Code Coverage Report

Directory:	./
File:	src/sys/classes/fn/tests/test12.c
Date:	2025-11-19 04:19:03

	Exec	Total	Coverage
Lines:	51	57	89.5%
Functions:	1	1	100.0%
Branches:	124	236	52.5%

  
      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 matrix function evaluation via diagonalization.\n\n";
    
      #include <slepcfn.h>
    
      30
      int main(int argc,char **argv)
    
      {
    
      30
        FN             fn;
    
      30
        Mat            A,F,G;
    
      30
        PetscInt       i,j,n=10;
    
      30
        PetscReal      nrm;
    
      30
        PetscScalar    *As,alpha,beta;
    
      30
        PetscViewer    viewer;
    
      30
        PetscBool      verbose;
    
        1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.

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

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

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

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

      30
        PetscCall(PetscPrintf(PETSC_COMM_WORLD,"Matrix function of symmetric/Hermitian matrix, n=%" PetscInt_FMT ".\n",n));
    
        /* Create function 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.

      30
        PetscCall(FNCreate(PETSC_COMM_WORLD,&fn));
    
        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
        PetscCall(FNSetType(fn,FNEXP));   /* default to exponential */
    
      #if defined(PETSC_USE_COMPLEX)
    
        1/2✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.

      15
        alpha = PetscCMPLX(0.3,0.8);
    
      15
        beta  = PetscCMPLX(1.1,-0.1);
    
      #else
    
      15
        alpha = 0.3;
    
      15
        beta  = 1.1;
    
      #endif
    
        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
        PetscCall(FNSetScale(fn,alpha,beta));
    
        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
        PetscCall(FNSetFromOptions(fn));
    
        /* 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.

      30
        PetscCall(PetscViewerASCIIGetStdout(PETSC_COMM_WORLD,&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.

      30
        if (verbose) PetscCall(PetscViewerPushFormat(viewer,PETSC_VIEWER_ASCII_MATLAB));
    
        /* Create a symmetric/Hermitian Toeplitz matrix */
    
        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
        PetscCall(MatCreateSeqDense(PETSC_COMM_SELF,n,n,NULL,&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.

      30
        PetscCall(PetscObjectSetName((PetscObject)A,"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.

      30
        PetscCall(MatDenseGetArray(A,&As));
    
        2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.

      330
        for (i=0;i<n;i++) As[i+i*n]=2.0;
    
        2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.

      90
        for (j=1;j<3;j++) {
    
        2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.

      570
          for (i=0;i<n-j;i++) {
    
      #if defined(PETSC_USE_COMPLEX)
    
      255
            As[i+(i+j)*n]=PetscCMPLX(1.0,0.1); As[(i+j)+i*n]=PetscCMPLX(1.0,-0.1);
    
      #else
    
      255
            As[i+(i+j)*n]=0.5; As[(i+j)+i*n]=0.5;
    
      #endif
    
          }
    
        }
    
        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
        PetscCall(MatDenseRestoreArray(A,&As));
    
        1/2✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.

      30
        if (verbose) {
    
      ✗
          PetscCall(PetscPrintf(PETSC_COMM_WORLD,"Matrix A - - - - - - - -\n"));
    
      ✗
          PetscCall(MatView(A,viewer));
    
        }
    
        /* compute matrix function */
    
        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
        PetscCall(MatCreateSeqDense(PETSC_COMM_SELF,n,n,NULL,&F));
    
        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
        PetscCall(PetscObjectSetName((PetscObject)F,"F"));
    
        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
        PetscCall(FNEvaluateFunctionMat(fn,A,F));
    
        1/2✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.

      30
        if (verbose) {
    
      ✗
          PetscCall(PetscPrintf(PETSC_COMM_WORLD,"Computed f(A) - - - - - - -\n"));
    
      ✗
          PetscCall(MatView(F,viewer));
    
        }
    
        /* Repeat with MAT_HERMITIAN flag set */
    
        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
        PetscCall(MatSetOption(A,MAT_HERMITIAN,PETSC_TRUE));
    
        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
        PetscCall(MatCreateSeqDense(PETSC_COMM_SELF,n,n,NULL,&G));
    
        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
        PetscCall(PetscObjectSetName((PetscObject)G,"G"));
    
        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
        PetscCall(FNEvaluateFunctionMat(fn,A,G));
    
        1/2✗ Branch 0 not taken.
✓ Branch 1 taken 10 times.

      30
        if (verbose) {
    
      ✗
          PetscCall(PetscPrintf(PETSC_COMM_WORLD,"Computed f(A) symm - - - - - - -\n"));
    
      ✗
          PetscCall(MatView(G,viewer));
    
        }
    
        /* compare the two results */
    
        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
        PetscCall(MatAXPY(F,-1.0,G,SAME_NONZERO_PATTERN));
    
        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
        PetscCall(MatNorm(F,NORM_FROBENIUS,&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.

      30
        if (nrm>100*PETSC_MACHINE_EPSILON) PetscCall(PetscPrintf(PETSC_COMM_WORLD,"Warning: the norm of F-G is %g\n",(double)nrm));
    
        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(PetscPrintf(PETSC_COMM_WORLD,"Computed results match.\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.

      30
        PetscCall(MatDestroy(&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.

      30
        PetscCall(MatDestroy(&F));
    
        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
        PetscCall(MatDestroy(&G));
    
        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
        PetscCall(FNDestroy(&fn));
    
        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.

      30
        PetscCall(SlepcFinalize());
    
        return 0;
    
      }
    
      /*TEST
    
         test:
    
            suffix: 1
    
            nsize: 1
    
            args: -fn_type {{exp sqrt}shared output}
    
            output_file: output/test12_1.out
    
         test:
    
            suffix: 1_rational
    
            nsize: 1
    
            args: -fn_type rational -fn_rational_numerator 2,-1.5 -fn_rational_denominator 1,0.8
    
            output_file: output/test12_1.out
    
      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 matrix function evaluation via diagonalization.\n\n";
12
13			#include <slepcfn.h>
14
15		30	int main(int argc,char **argv)
16			{
17		30	FN fn;
18		30	Mat A,F,G;
19		30	PetscInt i,j,n=10;
20		30	PetscReal nrm;
21		30	PetscScalar *As,alpha,beta;
22		30	PetscViewer viewer;
23		30	PetscBool verbose;
24
25	1/2 ✓ Branch 0 taken 2 times. ✗ Branch 1 not taken.	30	PetscFunctionBeginUser;
26	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	PetscCall(SlepcInitialize(&argc,&argv,NULL,help));
27	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	PetscCall(PetscOptionsGetInt(NULL,NULL,"-n",&n,NULL));
28	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	PetscCall(PetscOptionsHasName(NULL,NULL,"-verbose",&verbose));
29	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	PetscCall(PetscPrintf(PETSC_COMM_WORLD,"Matrix function of symmetric/Hermitian matrix, n=%" PetscInt_FMT ".\n",n));
30
31			/* Create function object */
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.	30	PetscCall(FNCreate(PETSC_COMM_WORLD,&fn));
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.	30	PetscCall(FNSetType(fn,FNEXP)); /* default to exponential */
34			#if defined(PETSC_USE_COMPLEX)
35	1/2 ✓ Branch 0 taken 1 times. ✗ Branch 1 not taken.	15	alpha = PetscCMPLX(0.3,0.8);
36		15	beta = PetscCMPLX(1.1,-0.1);
37			#else
38		15	alpha = 0.3;
39		15	beta = 1.1;
40			#endif
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.	30	PetscCall(FNSetScale(fn,alpha,beta));
42	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	PetscCall(FNSetFromOptions(fn));
43
44			/* Set up viewer */
45	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	PetscCall(PetscViewerASCIIGetStdout(PETSC_COMM_WORLD,&viewer));
46	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.	30	if (verbose) PetscCall(PetscViewerPushFormat(viewer,PETSC_VIEWER_ASCII_MATLAB));
47
48			/* Create a symmetric/Hermitian Toeplitz matrix */
49	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	PetscCall(MatCreateSeqDense(PETSC_COMM_SELF,n,n,NULL,&A));
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.	30	PetscCall(PetscObjectSetName((PetscObject)A,"A"));
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.	30	PetscCall(MatDenseGetArray(A,&As));
52	2/2 ✓ Branch 0 taken 10 times. ✓ Branch 1 taken 10 times.	330	for (i=0;i<n;i++) As[i+i*n]=2.0;
53	2/2 ✓ Branch 0 taken 10 times. ✓ Branch 1 taken 10 times.	90	for (j=1;j<3;j++) {
54	2/2 ✓ Branch 0 taken 10 times. ✓ Branch 1 taken 10 times.	570	for (i=0;i<n-j;i++) {
55			#if defined(PETSC_USE_COMPLEX)
56		255	As[i+(i+j)n]=PetscCMPLX(1.0,0.1); As[(i+j)+in]=PetscCMPLX(1.0,-0.1);
57			#else
58		255	As[i+(i+j)n]=0.5; As[(i+j)+in]=0.5;
59			#endif
60			}
61			}
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.	30	PetscCall(MatDenseRestoreArray(A,&As));
63	1/2 ✗ Branch 0 not taken. ✓ Branch 1 taken 10 times.	30	if (verbose) {
64		✗	PetscCall(PetscPrintf(PETSC_COMM_WORLD,"Matrix A - - - - - - - -\n"));
65		✗	PetscCall(MatView(A,viewer));
66			}
67
68			/* compute matrix function */
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.	30	PetscCall(MatCreateSeqDense(PETSC_COMM_SELF,n,n,NULL,&F));
70	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	PetscCall(PetscObjectSetName((PetscObject)F,"F"));
71	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	PetscCall(FNEvaluateFunctionMat(fn,A,F));
72	1/2 ✗ Branch 0 not taken. ✓ Branch 1 taken 10 times.	30	if (verbose) {
73		✗	PetscCall(PetscPrintf(PETSC_COMM_WORLD,"Computed f(A) - - - - - - -\n"));
74		✗	PetscCall(MatView(F,viewer));
75			}
76
77			/* Repeat with MAT_HERMITIAN flag set */
78	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	PetscCall(MatSetOption(A,MAT_HERMITIAN,PETSC_TRUE));
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.	30	PetscCall(MatCreateSeqDense(PETSC_COMM_SELF,n,n,NULL,&G));
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.	30	PetscCall(PetscObjectSetName((PetscObject)G,"G"));
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.	30	PetscCall(FNEvaluateFunctionMat(fn,A,G));
82	1/2 ✗ Branch 0 not taken. ✓ Branch 1 taken 10 times.	30	if (verbose) {
83		✗	PetscCall(PetscPrintf(PETSC_COMM_WORLD,"Computed f(A) symm - - - - - - -\n"));
84		✗	PetscCall(MatView(G,viewer));
85			}
86
87			/* compare the two results */
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.	30	PetscCall(MatAXPY(F,-1.0,G,SAME_NONZERO_PATTERN));
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.	30	PetscCall(MatNorm(F,NORM_FROBENIUS,&nrm));
90	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.	30	if (nrm>100*PETSC_MACHINE_EPSILON) PetscCall(PetscPrintf(PETSC_COMM_WORLD,"Warning: the norm of F-G is %g\n",(double)nrm));
91	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(PetscPrintf(PETSC_COMM_WORLD,"Computed results match.\n"));
92
93	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	PetscCall(MatDestroy(&A));
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.	30	PetscCall(MatDestroy(&F));
95	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	PetscCall(MatDestroy(&G));
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.	30	PetscCall(FNDestroy(&fn));
97	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.	30	PetscCall(SlepcFinalize());
98			return 0;
99			}
100
101			/*TEST
102
103			test:
104			suffix: 1
105			nsize: 1
106			args: -fn_type {{exp sqrt}shared output}
107			output_file: output/test12_1.out
108
109			test:
110			suffix: 1_rational
111			nsize: 1
112			args: -fn_type rational -fn_rational_numerator 2,-1.5 -fn_rational_denominator 1,0.8
113			output_file: output/test12_1.out
114
115			TEST*/
116