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[] = "Solves a generalized eigensystem Ax=kBx with matrices loaded from a file.\n\n" | ||
12 | "The command line options are:\n" | ||
13 | " -f1 <filename> -f2 <filename>, PETSc binary files containing A and B.\n" | ||
14 | " -evecs <filename>, output file to save computed eigenvectors.\n" | ||
15 | " -ninitial <nini>, number of user-provided initial guesses.\n" | ||
16 | " -finitial <filename>, binary file containing <nini> vectors.\n" | ||
17 | " -nconstr <ncon>, number of user-provided constraints.\n" | ||
18 | " -fconstr <filename>, binary file containing <ncon> vectors.\n\n"; | ||
19 | |||
20 | #include <slepceps.h> | ||
21 | |||
22 | 13 | int main(int argc,char **argv) | |
23 | { | ||
24 | 13 | Mat A,B; /* matrices */ | |
25 | 13 | EPS eps; /* eigenproblem solver context */ | |
26 | 13 | ST st; | |
27 | 13 | KSP ksp; | |
28 | 13 | EPSType type; | |
29 | 13 | PetscReal tol; | |
30 | 13 | Vec xr,xi,*Iv,*Cv; | |
31 | 13 | PetscInt nev,maxit,i,its,lits,nconv,nini=0,ncon=0; | |
32 | 13 | char filename[PETSC_MAX_PATH_LEN]; | |
33 | 13 | PetscViewer viewer; | |
34 | 13 | PetscBool flg,evecs,ishermitian,terse; | |
35 | |||
36 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
13 | PetscFunctionBeginUser; |
37 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(SlepcInitialize(&argc,&argv,NULL,help)); |
38 | |||
39 | /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - | ||
40 | Load the matrices that define the eigensystem, Ax=kBx | ||
41 | - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */ | ||
42 | |||
43 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(PetscPrintf(PETSC_COMM_WORLD,"\nGeneralized eigenproblem stored in file.\n\n")); |
44 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(PetscOptionsGetString(NULL,NULL,"-f1",filename,sizeof(filename),&flg)); |
45 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 8 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
13 | PetscCheck(flg,PETSC_COMM_WORLD,PETSC_ERR_USER_INPUT,"Must indicate a file name for matrix A with the -f1 option"); |
46 | |||
47 | #if defined(PETSC_USE_COMPLEX) | ||
48 |
4/6✓ Branch 0 taken 1 times.
✓ Branch 1 taken 2 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
|
3 | PetscCall(PetscPrintf(PETSC_COMM_WORLD," Reading COMPLEX matrices from binary files...\n")); |
49 | #else | ||
50 |
4/6✓ Branch 0 taken 1 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 1 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 1 times.
|
10 | PetscCall(PetscPrintf(PETSC_COMM_WORLD," Reading REAL matrices from binary files...\n")); |
51 | #endif | ||
52 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(PetscViewerBinaryOpen(PETSC_COMM_WORLD,filename,FILE_MODE_READ,&viewer)); |
53 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(MatCreate(PETSC_COMM_WORLD,&A)); |
54 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(MatSetFromOptions(A)); |
55 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(MatLoad(A,viewer)); |
56 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(PetscViewerDestroy(&viewer)); |
57 | |||
58 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(PetscOptionsGetString(NULL,NULL,"-f2",filename,sizeof(filename),&flg)); |
59 |
1/2✓ Branch 0 taken 8 times.
✗ Branch 1 not taken.
|
13 | if (flg) { |
60 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(PetscViewerBinaryOpen(PETSC_COMM_WORLD,filename,FILE_MODE_READ,&viewer)); |
61 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(MatCreate(PETSC_COMM_WORLD,&B)); |
62 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(MatSetFromOptions(B)); |
63 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(MatLoad(B,viewer)); |
64 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(PetscViewerDestroy(&viewer)); |
65 | } else { | ||
66 | ✗ | PetscCall(PetscPrintf(PETSC_COMM_WORLD," Matrix B was not provided, setting B=I\n\n")); | |
67 | ✗ | B = NULL; | |
68 | } | ||
69 | |||
70 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(MatCreateVecs(A,NULL,&xr)); |
71 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(MatCreateVecs(A,NULL,&xi)); |
72 | |||
73 | /* | ||
74 | Read user constraints if available | ||
75 | */ | ||
76 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(PetscOptionsGetInt(NULL,NULL,"-nconstr",&ncon,&flg)); |
77 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 8 times.
|
13 | if (flg) { |
78 | ✗ | PetscCheck(ncon>0,PETSC_COMM_WORLD,PETSC_ERR_USER_INPUT,"The number of constraints must be >0"); | |
79 | ✗ | PetscCall(PetscOptionsGetString(NULL,NULL,"-fconstr",filename,sizeof(filename),&flg)); | |
80 | ✗ | PetscCheck(flg,PETSC_COMM_WORLD,PETSC_ERR_USER_INPUT,"Must specify the name of the file storing the constraints"); | |
81 | ✗ | PetscCall(PetscViewerBinaryOpen(PETSC_COMM_WORLD,filename,FILE_MODE_READ,&viewer)); | |
82 | ✗ | PetscCall(VecDuplicateVecs(xr,ncon,&Cv)); | |
83 | ✗ | for (i=0;i<ncon;i++) PetscCall(VecLoad(Cv[i],viewer)); | |
84 | ✗ | PetscCall(PetscViewerDestroy(&viewer)); | |
85 | } | ||
86 | |||
87 | /* | ||
88 | Read initial guesses if available | ||
89 | */ | ||
90 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(PetscOptionsGetInt(NULL,NULL,"-ninitial",&nini,&flg)); |
91 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 8 times.
|
13 | if (flg) { |
92 | ✗ | PetscCheck(nini>0,PETSC_COMM_WORLD,PETSC_ERR_USER_INPUT,"The number of initial vectors must be >0"); | |
93 | ✗ | PetscCall(PetscOptionsGetString(NULL,NULL,"-finitial",filename,sizeof(filename),&flg)); | |
94 | ✗ | PetscCheck(flg,PETSC_COMM_WORLD,PETSC_ERR_USER_INPUT,"Must specify the name of the file containing the initial vectors"); | |
95 | ✗ | PetscCall(PetscViewerBinaryOpen(PETSC_COMM_WORLD,filename,FILE_MODE_READ,&viewer)); | |
96 | ✗ | PetscCall(VecDuplicateVecs(xr,nini,&Iv)); | |
97 | ✗ | for (i=0;i<nini;i++) PetscCall(VecLoad(Iv[i],viewer)); | |
98 | ✗ | PetscCall(PetscViewerDestroy(&viewer)); | |
99 | } | ||
100 | |||
101 | /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - | ||
102 | Create the eigensolver and set various options | ||
103 | - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */ | ||
104 | |||
105 | /* | ||
106 | Create eigensolver context | ||
107 | */ | ||
108 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(EPSCreate(PETSC_COMM_WORLD,&eps)); |
109 | |||
110 | /* | ||
111 | Set operators. In this case, it is a generalized eigenvalue problem | ||
112 | */ | ||
113 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(EPSSetOperators(eps,A,B)); |
114 | |||
115 | /* | ||
116 | If the user provided initial guesses or constraints, pass them here | ||
117 | */ | ||
118 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(EPSSetInitialSpace(eps,nini,Iv)); |
119 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(EPSSetDeflationSpace(eps,ncon,Cv)); |
120 | |||
121 | /* | ||
122 | Set solver parameters at runtime | ||
123 | */ | ||
124 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(EPSSetFromOptions(eps)); |
125 | |||
126 | /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - | ||
127 | Solve the eigensystem | ||
128 | - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */ | ||
129 | |||
130 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(EPSSolve(eps)); |
131 | |||
132 | /* | ||
133 | Optional: Get some information from the solver and display it | ||
134 | */ | ||
135 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(EPSGetIterationNumber(eps,&its)); |
136 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(PetscPrintf(PETSC_COMM_WORLD," Number of iterations of the method: %" PetscInt_FMT "\n",its)); |
137 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(EPSGetST(eps,&st)); |
138 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(STGetKSP(st,&ksp)); |
139 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(KSPGetTotalIterations(ksp,&lits)); |
140 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(PetscPrintf(PETSC_COMM_WORLD," Number of linear iterations of the method: %" PetscInt_FMT "\n",lits)); |
141 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(EPSGetType(eps,&type)); |
142 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(PetscPrintf(PETSC_COMM_WORLD," Solution method: %s\n\n",type)); |
143 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(EPSGetDimensions(eps,&nev,NULL,NULL)); |
144 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(PetscPrintf(PETSC_COMM_WORLD," Number of requested eigenvalues: %" PetscInt_FMT "\n",nev)); |
145 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(EPSGetTolerances(eps,&tol,&maxit)); |
146 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(PetscPrintf(PETSC_COMM_WORLD," Stopping condition: tol=%.4g, maxit=%" PetscInt_FMT "\n",(double)tol,maxit)); |
147 | |||
148 | /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - | ||
149 | Display solution and clean up | ||
150 | - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */ | ||
151 | |||
152 | /* | ||
153 | Show detailed info unless -terse option is given by user | ||
154 | */ | ||
155 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(PetscOptionsHasName(NULL,NULL,"-terse",&terse)); |
156 |
5/8✓ Branch 0 taken 8 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 2 times.
✓ Branch 3 taken 6 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 2 times.
|
13 | if (terse) PetscCall(EPSErrorView(eps,EPS_ERROR_RELATIVE,NULL)); |
157 | else { | ||
158 | ✗ | PetscCall(PetscViewerPushFormat(PETSC_VIEWER_STDOUT_WORLD,PETSC_VIEWER_ASCII_INFO_DETAIL)); | |
159 | ✗ | PetscCall(EPSConvergedReasonView(eps,PETSC_VIEWER_STDOUT_WORLD)); | |
160 | ✗ | PetscCall(EPSErrorView(eps,EPS_ERROR_RELATIVE,PETSC_VIEWER_STDOUT_WORLD)); | |
161 | ✗ | PetscCall(PetscViewerPopFormat(PETSC_VIEWER_STDOUT_WORLD)); | |
162 | } | ||
163 | |||
164 | /* | ||
165 | Save eigenvectors, if requested | ||
166 | */ | ||
167 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(PetscOptionsGetString(NULL,NULL,"-evecs",filename,sizeof(filename),&evecs)); |
168 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(EPSGetConverged(eps,&nconv)); |
169 |
2/4✓ Branch 0 taken 8 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 8 times.
|
13 | if (nconv>0 && evecs) { |
170 | ✗ | PetscCall(PetscViewerBinaryOpen(PETSC_COMM_WORLD,filename,FILE_MODE_WRITE,&viewer)); | |
171 | ✗ | PetscCall(EPSIsHermitian(eps,&ishermitian)); | |
172 | ✗ | for (i=0;i<nconv;i++) { | |
173 | ✗ | PetscCall(EPSGetEigenvector(eps,i,xr,xi)); | |
174 | ✗ | PetscCall(VecView(xr,viewer)); | |
175 | #if !defined(PETSC_USE_COMPLEX) | ||
176 | ✗ | if (!ishermitian) PetscCall(VecView(xi,viewer)); | |
177 | #endif | ||
178 | } | ||
179 | ✗ | PetscCall(PetscViewerDestroy(&viewer)); | |
180 | } | ||
181 | |||
182 | /* | ||
183 | Free work space | ||
184 | */ | ||
185 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(EPSDestroy(&eps)); |
186 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(MatDestroy(&A)); |
187 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(MatDestroy(&B)); |
188 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(VecDestroy(&xr)); |
189 |
4/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 6 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✓ Branch 5 taken 2 times.
|
13 | PetscCall(VecDestroy(&xi)); |
190 |
1/8✗ Branch 0 not taken.
✓ Branch 1 taken 8 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.
|
13 | if (nini > 0) PetscCall(VecDestroyVecs(nini,&Iv)); |
191 |
1/8✗ Branch 0 not taken.
✓ Branch 1 taken 8 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.
|
13 | if (ncon > 0) PetscCall(VecDestroyVecs(ncon,&Cv)); |
192 |
3/6✓ Branch 0 taken 2 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 2 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
|
13 | PetscCall(SlepcFinalize()); |
193 | return 0; | ||
194 | } | ||
195 | |||
196 | /*TEST | ||
197 | |||
198 | test: | ||
199 | suffix: 1 | ||
200 | args: -f1 ${SLEPC_DIR}/share/slepc/datafiles/matrices/bfw62a.petsc -f2 ${SLEPC_DIR}/share/slepc/datafiles/matrices/bfw62b.petsc -eps_nev 4 -terse | ||
201 | requires: double !complex !defined(PETSC_USE_64BIT_INDICES) | ||
202 | |||
203 | test: | ||
204 | suffix: ciss_1 | ||
205 | args: -f1 ${DATAFILESPATH}/matrices/complex/mhd1280a.petsc -f2 ${DATAFILESPATH}/matrices/complex/mhd1280b.petsc -eps_type ciss -eps_ciss_usest 0 -eps_ciss_quadrule chebyshev -rg_type ring -rg_ring_center 0 -rg_ring_radius .49 -rg_ring_width 0.2 -rg_ring_startangle .25 -rg_ring_endangle .49 -terse -eps_max_it 1 | ||
206 | requires: double complex datafilespath !defined(PETSC_USE_64BIT_INDICES) | ||
207 | timeoutfactor: 2 | ||
208 | |||
209 | test: | ||
210 | suffix: 3 # test problem (A,A) | ||
211 | args: -f1 ${SLEPC_DIR}/share/slepc/datafiles/matrices/bfw62a.petsc -f2 ${SLEPC_DIR}/share/slepc/datafiles/matrices/bfw62a.petsc -eps_nev 4 -terse | ||
212 | requires: double !complex !defined(PETSC_USE_64BIT_INDICES) | ||
213 | |||
214 | TEST*/ | ||
215 |