GCC Code Coverage Report

Directory:	./
File:	src/svd/tests/test2.c
Date:	2025-11-19 04:19:03

	Exec	Total	Coverage
Lines:	33	34	97.1%
Functions:	1	1	100.0%
Branches:	72	118	61.0%

  
      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 SVD with different builds with a matrix loaded from a file.\n\n"
    
        " (matrices available in PETSc's distribution).\n\n";
    
      #include <slepcsvd.h>
    
      48
      int main(int argc,char **argv)
    
      {
    
      48
        Mat            A;               /* operator matrix */
    
      48
        SVD            svd;             /* singular value problem solver context */
    
      48
        char           filename[PETSC_MAX_PATH_LEN],path[PETSC_MAX_PATH_LEN];
    
      48
        const char     *prefix,*scalar,*ints,*floats;
    
      48
        PetscReal      tol=PETSC_SMALL;
    
      48
        PetscViewer    viewer;
    
      48
        PetscBool      flg;
    
        1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.

      48
        PetscFunctionBeginUser;
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      48
        PetscCall(SlepcInitialize(&argc,&argv,NULL,help));
    
        /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
    
              Load the matrix for which the SVD must be computed
    
           - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
    
      #if defined(PETSC_USE_COMPLEX)
    
      24
        prefix = "nh";
    
      24
        scalar = "complex";
    
      #else
    
      24
        prefix = "ns";
    
      24
        scalar = "real";
    
      #endif
    
      #if defined(PETSC_USE_64BIT_INDICES)
    
        ints   = "int64";
    
      #else
    
      48
        ints   = "int32";
    
      #endif
    
      #if defined(PETSC_USE_REAL_DOUBLE)
    
      48
        floats = "float64";
    
      #elif defined(PETSC_USE_REAL_SINGLE)
    
        floats = "float32";
    
      #endif
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      48
        PetscCall(PetscOptionsGetString(NULL,NULL,"-path",path,sizeof(path),&flg));
    
        5/8✓ Branch 0 taken 6 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 4 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.

      48
        if (flg) PetscCall(PetscSNPrintf(filename,sizeof(filename),"%s/%s-%s-%s-%s",path,prefix,scalar,ints,floats));
    
      ✗
        else PetscCall(PetscSNPrintf(filename,sizeof(filename),"%s/share/petsc/datafiles/matrices/%s-%s-%s-%s",PETSC_DIR,prefix,scalar,ints,floats));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      48
        PetscCall(PetscPrintf(PETSC_COMM_WORLD,"\nReading matrix from binary file...\n\n"));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      48
        PetscCall(PetscViewerBinaryOpen(PETSC_COMM_WORLD,filename,FILE_MODE_READ,&viewer));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      48
        PetscCall(MatCreate(PETSC_COMM_WORLD,&A));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      48
        PetscCall(MatSetFromOptions(A));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      48
        PetscCall(MatLoad(A,viewer));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      48
        PetscCall(PetscViewerDestroy(&viewer));
    
        /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
    
                           Create the SVD solver
    
           - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      48
        PetscCall(SVDCreate(PETSC_COMM_WORLD,&svd));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      48
        PetscCall(SVDSetOperators(svd,A,NULL));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      48
        PetscCall(SVDSetTolerances(svd,tol,PETSC_CURRENT));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      48
        PetscCall(SVDSetFromOptions(svd));
    
        /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
    
                      Compute the singular triplets and display solution
    
           - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      48
        PetscCall(SVDSolve(svd));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      48
        PetscCall(SVDErrorView(svd,SVD_ERROR_RELATIVE,NULL));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      48
        PetscCall(SVDDestroy(&svd));
    
        4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.

      48
        PetscCall(MatDestroy(&A));
    
        2/6✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.

      48
        PetscCall(SlepcFinalize());
    
        return 0;
    
      }
    
      /*TEST
    
         build:
    
            requires: !__float128
    
         test:
    
            args: -svd_nsv 7 -svd_type {{lanczos trlanczos cross cyclic lapack randomized}} -path ${DATAFILESPATH}/matrices/petsc-small
    
            requires: !single datafilespath
    
         testset:
    
            args: -svd_nsv 7 -svd_mpd 11 -svd_type primme -path ${DATAFILESPATH}/matrices/petsc-small
    
            requires: primme !single datafilespath
    
            output_file: output/test2_1.out
    
            test:
    
               suffix: 1_primme
    
            test:
    
               suffix: 1_primme_args
    
               args: -svd_primme_blocksize 2 -svd_primme_method hybrid
    
      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 SVD with different builds with a matrix loaded from a file.\n\n"
12			" (matrices available in PETSc's distribution).\n\n";
13
14			#include <slepcsvd.h>
15
16		48	int main(int argc,char **argv)
17			{
18		48	Mat A; /* operator matrix */
19		48	SVD svd; /* singular value problem solver context */
20		48	char filename[PETSC_MAX_PATH_LEN],path[PETSC_MAX_PATH_LEN];
21		48	const char prefix,scalar,ints,floats;
22		48	PetscReal tol=PETSC_SMALL;
23		48	PetscViewer viewer;
24		48	PetscBool flg;
25
26	1/2 ✓ Branch 0 taken 2 times. ✗ Branch 1 not taken.	48	PetscFunctionBeginUser;
27	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	48	PetscCall(SlepcInitialize(&argc,&argv,NULL,help));
28
29			/* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
30			Load the matrix for which the SVD must be computed
31			- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
32			#if defined(PETSC_USE_COMPLEX)
33		24	prefix = "nh";
34		24	scalar = "complex";
35			#else
36		24	prefix = "ns";
37		24	scalar = "real";
38			#endif
39			#if defined(PETSC_USE_64BIT_INDICES)
40			ints = "int64";
41			#else
42		48	ints = "int32";
43			#endif
44			#if defined(PETSC_USE_REAL_DOUBLE)
45		48	floats = "float64";
46			#elif defined(PETSC_USE_REAL_SINGLE)
47			floats = "float32";
48			#endif
49
50	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	48	PetscCall(PetscOptionsGetString(NULL,NULL,"-path",path,sizeof(path),&flg));
51	5/8 ✓ Branch 0 taken 6 times. ✗ Branch 1 not taken. ✓ Branch 2 taken 2 times. ✓ Branch 3 taken 4 times. ✓ Branch 4 taken 2 times. ✗ Branch 5 not taken. ✗ Branch 6 not taken. ✓ Branch 7 taken 2 times.	48	if (flg) PetscCall(PetscSNPrintf(filename,sizeof(filename),"%s/%s-%s-%s-%s",path,prefix,scalar,ints,floats));
52		✗	else PetscCall(PetscSNPrintf(filename,sizeof(filename),"%s/share/petsc/datafiles/matrices/%s-%s-%s-%s",PETSC_DIR,prefix,scalar,ints,floats));
53	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	48	PetscCall(PetscPrintf(PETSC_COMM_WORLD,"\nReading matrix from binary file...\n\n"));
54	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	48	PetscCall(PetscViewerBinaryOpen(PETSC_COMM_WORLD,filename,FILE_MODE_READ,&viewer));
55	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	48	PetscCall(MatCreate(PETSC_COMM_WORLD,&A));
56	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	48	PetscCall(MatSetFromOptions(A));
57	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	48	PetscCall(MatLoad(A,viewer));
58	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	48	PetscCall(PetscViewerDestroy(&viewer));
59
60			/* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
61			Create the SVD solver
62			- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
63	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	48	PetscCall(SVDCreate(PETSC_COMM_WORLD,&svd));
64	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	48	PetscCall(SVDSetOperators(svd,A,NULL));
65	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	48	PetscCall(SVDSetTolerances(svd,tol,PETSC_CURRENT));
66	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	48	PetscCall(SVDSetFromOptions(svd));
67
68			/* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
69			Compute the singular triplets and display solution
70			- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
71	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	48	PetscCall(SVDSolve(svd));
72	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	48	PetscCall(SVDErrorView(svd,SVD_ERROR_RELATIVE,NULL));
73	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	48	PetscCall(SVDDestroy(&svd));
74	4/6 ✓ Branch 0 taken 2 times. ✓ Branch 1 taken 4 times. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✓ Branch 5 taken 2 times.	48	PetscCall(MatDestroy(&A));
75	2/6 ✓ Branch 0 taken 2 times. ✗ Branch 1 not taken. ✓ Branch 2 taken 2 times. ✗ Branch 3 not taken. ✗ Branch 4 not taken. ✗ Branch 5 not taken.	48	PetscCall(SlepcFinalize());
76			return 0;
77			}
78
79			/*TEST
80
81			build:
82			requires: !__float128
83
84			test:
85			args: -svd_nsv 7 -svd_type {{lanczos trlanczos cross cyclic lapack randomized}} -path ${DATAFILESPATH}/matrices/petsc-small
86			requires: !single datafilespath
87
88			testset:
89			args: -svd_nsv 7 -svd_mpd 11 -svd_type primme -path ${DATAFILESPATH}/matrices/petsc-small
90			requires: primme !single datafilespath
91			output_file: output/test2_1.out
92			test:
93			suffix: 1_primme
94			test:
95			suffix: 1_primme_args
96			args: -svd_primme_blocksize 2 -svd_primme_method hybrid
97
98			TEST*/
99