GCC Code Coverage Report

Directory: ./
File: src/sys/classes/bv/tests/test6.c
Date: 2025-10-04 04:19:13
Exec Total Coverage
Lines: 59 60 98.3%
Functions: 1 1 100.0%
Branches: 169 276 61.2%
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 BV orthogonalization functions with constraints.\n\n";
12
13 #include <slepcbv.h>
14
15 176 int main(int argc,char **argv)
16 {
17 176 BV X;
18 176 Mat M;
19 176 Vec v,t,*C;
20 176 PetscInt i,j,n=20,k=8,nc=2,nloc;
21 176 PetscViewer view;
22 176 PetscBool verbose;
23 176 PetscReal norm;
24 176 PetscScalar alpha;
25
26
1/2
✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
 176 PetscFunctionBeginUser;
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.
 176 PetscCall(SlepcInitialize(&argc,&argv,NULL,help));
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.
 176 PetscCall(PetscOptionsGetInt(NULL,NULL,"-n",&n,NULL));
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.
 176 PetscCall(PetscOptionsGetInt(NULL,NULL,"-k",&k,NULL));
30
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.
 176 PetscCall(PetscOptionsGetInt(NULL,NULL,"-nc",&nc,NULL));
31
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.
 176 PetscCall(PetscOptionsHasName(NULL,NULL,"-verbose",&verbose));
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.
 176 PetscCall(PetscPrintf(PETSC_COMM_WORLD,"Test BV orthogonalization with %" PetscInt_FMT " columns + %" PetscInt_FMT " constraints, of length %" PetscInt_FMT ".\n",k,nc,n));
33
34 /* Create template vector */
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.
 176 PetscCall(VecCreate(PETSC_COMM_WORLD,&t));
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.
 176 PetscCall(VecSetSizes(t,PETSC_DECIDE,n));
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.
 176 PetscCall(VecSetFromOptions(t));
38
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.
 176 PetscCall(VecGetLocalSize(t,&nloc));
39
40 /* Create BV object X */
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.
 176 PetscCall(BVCreate(PETSC_COMM_WORLD,&X));
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.
 176 PetscCall(PetscObjectSetName((PetscObject)X,"X"));
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.
 176 PetscCall(BVSetSizes(X,nloc,n,k));
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.
 176 PetscCall(BVSetFromOptions(X));
45
46 /* Generate constraints and attach them to X */
47
1/2
✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
 176 if (nc>0) {
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.
 176 PetscCall(VecDuplicateVecs(t,nc,&C));
49
2/2
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
 528 for (j=0;j<nc;j++) {
50
7/8
✓ Branch 0 taken 2 times.
✓ Branch 1 taken 8 times.
✓ Branch 2 taken 10 times.
✓ Branch 3 taken 8 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
✓ Branch 6 taken 2 times.
✓ Branch 7 taken 2 times.
 880 for (i=0;i<=j;i++) PetscCall(VecSetValue(C[j],i,1.0,INSERT_VALUES));
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.
 352 PetscCall(VecAssemblyBegin(C[j]));
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.
 352 PetscCall(VecAssemblyEnd(C[j]));
53 }
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.
 176 PetscCall(BVInsertConstraints(X,&nc,C));
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.
 176 PetscCall(VecDestroyVecs(nc,&C));
56 }
57
58 /* Set up viewer */
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.
 176 PetscCall(PetscViewerASCIIGetStdout(PETSC_COMM_WORLD,&view));
60
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.
 176 if (verbose) PetscCall(PetscViewerPushFormat(view,PETSC_VIEWER_ASCII_MATLAB));
61
62 /* Fill X entries */
63
2/2
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
 1584 for (j=0;j<k;j++) {
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.
 1408 PetscCall(BVGetColumn(X,j,&v));
65
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.
 1408 PetscCall(VecSet(v,0.0));
66
2/2
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
 16896 for (i=0;i<=n/2;i++) {
67
1/2
✓ Branch 0 taken 10 times.
✗ Branch 1 not taken.
 15488 if (i+j<n) {
68 15488 alpha = (3.0*i+j-2)/(2*(i+j+1));
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.
 15488 PetscCall(VecSetValue(v,i+j,alpha,INSERT_VALUES));
70 }
71 }
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.
 1408 PetscCall(VecAssemblyBegin(v));
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.
 1408 PetscCall(VecAssemblyEnd(v));
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.
 1408 PetscCall(BVRestoreColumn(X,j,&v));
75 }
76
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.
 176 if (verbose) PetscCall(BVView(X,view));
77
78 /* Test BVOrthogonalizeColumn */
79
2/2
✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
 1584 for (j=0;j<k;j++) {
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.
 1408 PetscCall(BVOrthogonalizeColumn(X,j,NULL,&norm,NULL));
81 1408 alpha = 1.0/norm;
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.
 1408 PetscCall(BVScaleColumn(X,j,alpha));
83 }
84
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.
 176 if (verbose) PetscCall(BVView(X,view));
85
86 /* Check orthogonality */
87
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.
 176 PetscCall(MatCreateSeqDense(PETSC_COMM_SELF,k,k,NULL,&M));
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.
 176 PetscCall(BVDot(X,X,M));
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.
 176 PetscCall(MatShift(M,-1.0));
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.
 176 PetscCall(MatNorm(M,NORM_1,&norm));
91
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.
 176 if (norm<100*PETSC_MACHINE_EPSILON) PetscCall(PetscPrintf(PETSC_COMM_WORLD,"Level of orthogonality < 100*eps\n"));
92 ✗ else PetscCall(PetscPrintf(PETSC_COMM_WORLD,"Level of orthogonality: %g\n",(double)norm));
93
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.
 176 PetscCall(MatDestroy(&M));
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.
 176 PetscCall(BVDestroy(&X));
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.
 176 PetscCall(VecDestroy(&t));
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.
 176 PetscCall(SlepcFinalize());
98 return 0;
99 }
100
101 /*TEST
102
103 testset:
104 output_file: output/test6_1.out
105 test:
106 suffix: 1
107 args: -bv_type {{vecs contiguous svec mat}shared output}
108 test:
109 suffix: 1_cuda
110 args: -bv_type {{svec mat}} -vec_type cuda
111 requires: cuda
112 test:
113 suffix: 1_hip
114 args: -bv_type {{svec mat}} -vec_type hip
115 requires: hip
116 test:
117 suffix: 2
118 nsize: 2
119 args: -bv_type {{vecs contiguous svec mat}shared output}
120 test:
121 suffix: 3
122 args: -bv_type {{vecs contiguous svec mat}shared output} -bv_orthog_type mgs
123 test:
124 suffix: 3_cuda
125 args: -bv_type {{svec mat}} -vec_type cuda -bv_orthog_type mgs
126 requires: cuda
127 test:
128 suffix: 3_hip
129 args: -bv_type {{svec mat}} -vec_type hip -bv_orthog_type mgs
130 requires: hip
131
132 TEST*/
133