| 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 the solution of a HEP without calling EPSSetFromOptions (based on ex1.c).\n\n" | ||
| 12 | "The command line options are:\n" | ||
| 13 | " -n <n>, where <n> = number of grid subdivisions = matrix dimension.\n" | ||
| 14 | " -type <eps_type> = eps type to test.\n\n"; | ||
| 15 | |||
| 16 | #include <slepceps.h> | ||
| 17 | |||
| 18 | 111 | int main(int argc,char **argv) | |
| 19 | { | ||
| 20 | 111 | Mat A; /* problem matrix */ | |
| 21 | 111 | EPS eps; /* eigenproblem solver context */ | |
| 22 | 111 | PetscReal tol=1000*PETSC_MACHINE_EPSILON; | |
| 23 | 111 | PetscInt n=30,i,Istart,Iend,nev; | |
| 24 | 111 | PetscBool flg,gd2; | |
| 25 | 111 | char epstype[30] = "krylovschur"; | |
| 26 | |||
| 27 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
111 | PetscFunctionBeginUser; |
| 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.
|
111 | PetscCall(SlepcInitialize(&argc,&argv,NULL,help)); |
| 29 | |||
| 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.
|
111 | PetscCall(PetscOptionsGetInt(NULL,NULL,"-n",&n,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.
|
111 | PetscCall(PetscOptionsGetString(NULL,NULL,"-type",epstype,sizeof(epstype),NULL)); |
| 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.
|
111 | PetscCall(PetscPrintf(PETSC_COMM_WORLD,"\n1-D Laplacian Eigenproblem, n=%" PetscInt_FMT "\n\n",n)); |
| 33 | |||
| 34 | /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - | ||
| 35 | Compute the operator matrix that defines the eigensystem, Ax=kx | ||
| 36 | - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */ | ||
| 37 | |||
| 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.
|
111 | PetscCall(MatCreate(PETSC_COMM_WORLD,&A)); |
| 39 |
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.
|
111 | PetscCall(MatSetSizes(A,PETSC_DECIDE,PETSC_DECIDE,n,n)); |
| 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.
|
111 | PetscCall(MatSetFromOptions(A)); |
| 41 | |||
| 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.
|
111 | PetscCall(MatGetOwnershipRange(A,&Istart,&Iend)); |
| 43 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
3441 | for (i=Istart;i<Iend;i++) { |
| 44 |
6/8✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ 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.
|
3330 | if (i>0) PetscCall(MatSetValue(A,i,i-1,-1.0,INSERT_VALUES)); |
| 45 |
6/8✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ 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.
|
3330 | if (i<n-1) PetscCall(MatSetValue(A,i,i+1,-1.0,INSERT_VALUES)); |
| 46 |
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.
|
3330 | PetscCall(MatSetValue(A,i,i,2.0,INSERT_VALUES)); |
| 47 | } | ||
| 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.
|
111 | PetscCall(MatAssemblyBegin(A,MAT_FINAL_ASSEMBLY)); |
| 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.
|
111 | PetscCall(MatAssemblyEnd(A,MAT_FINAL_ASSEMBLY)); |
| 50 | |||
| 51 | /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - | ||
| 52 | Create the eigensolver and set various options | ||
| 53 | - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */ | ||
| 54 | /* | ||
| 55 | Create eigensolver context | ||
| 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.
|
111 | PetscCall(EPSCreate(PETSC_COMM_WORLD,&eps)); |
| 58 | |||
| 59 | /* | ||
| 60 | Set operators. In this case, it is a standard eigenvalue problem | ||
| 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.
|
111 | PetscCall(EPSSetOperators(eps,A,NULL)); |
| 63 |
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.
|
111 | PetscCall(EPSSetProblemType(eps,EPS_HEP)); |
| 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.
|
111 | PetscCall(EPSSetDimensions(eps,4,PETSC_DETERMINE,PETSC_DETERMINE)); |
| 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.
|
111 | PetscCall(EPSSetTolerances(eps,tol,PETSC_CURRENT)); |
| 66 | |||
| 67 | /* | ||
| 68 | Set solver parameters at runtime | ||
| 69 | */ | ||
| 70 |
3/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 taken 2 times.
|
111 | PetscCall(PetscStrcmp(epstype,"gd2",&flg)); |
| 71 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
111 | if (flg) { |
| 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.
|
10 | PetscCall(EPSSetType(eps,EPSGD)); |
| 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.
|
10 | PetscCall(EPSGDSetDoubleExpansion(eps,PETSC_TRUE)); |
| 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.
|
10 | PetscCall(EPSGDGetDoubleExpansion(eps,&gd2)); /* not used */ |
| 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.
|
101 | } else PetscCall(EPSSetType(eps,epstype)); |
| 76 |
3/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 taken 2 times.
|
111 | PetscCall(PetscStrcmp(epstype,"jd",&flg)); |
| 77 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
111 | if (flg) { |
| 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.
|
10 | PetscCall(EPSSetWhichEigenpairs(eps,EPS_TARGET_MAGNITUDE)); |
| 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.
|
10 | PetscCall(EPSSetTarget(eps,4.0)); |
| 80 | } | ||
| 81 |
3/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 taken 2 times.
|
111 | PetscCall(PetscStrcmp(epstype,"lanczos",&flg)); |
| 82 |
6/8✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ 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.
|
111 | if (flg) PetscCall(EPSLanczosSetReorthog(eps,EPS_LANCZOS_REORTHOG_LOCAL)); |
| 83 |
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.
|
111 | PetscCall(PetscObjectTypeCompareAny((PetscObject)eps,&flg,EPSRQCG,EPSLOBPCG,"")); |
| 84 |
6/8✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
✓ 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.
|
111 | if (flg) PetscCall(EPSSetWhichEigenpairs(eps,EPS_SMALLEST_REAL)); |
| 85 | |||
| 86 | /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - | ||
| 87 | Solve the eigensystem | ||
| 88 | - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */ | ||
| 89 | |||
| 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.
|
111 | PetscCall(EPSSolve(eps)); |
| 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.
|
111 | PetscCall(EPSGetDimensions(eps,&nev,NULL,NULL)); |
| 92 |
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.
|
111 | PetscCall(PetscPrintf(PETSC_COMM_WORLD," Number of requested eigenvalues: %" PetscInt_FMT "\n",nev)); |
| 93 | |||
| 94 | /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - | ||
| 95 | Display solution and clean up | ||
| 96 | - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */ | ||
| 97 | |||
| 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.
|
111 | PetscCall(EPSErrorView(eps,EPS_ERROR_RELATIVE,NULL)); |
| 99 |
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.
|
111 | PetscCall(EPSDestroy(&eps)); |
| 100 |
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.
|
111 | PetscCall(MatDestroy(&A)); |
| 101 |
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.
|
111 | PetscCall(SlepcFinalize()); |
| 102 | return 0; | ||
| 103 | } | ||
| 104 | |||
| 105 | /*TEST | ||
| 106 | |||
| 107 | testset: | ||
| 108 | filter: sed -e "s/3.95905/3.95906/" | ||
| 109 | output_file: output/test4_1.out | ||
| 110 | test: | ||
| 111 | suffix: 1 | ||
| 112 | args: -type {{krylovschur subspace arnoldi lanczos gd jd gd2 lapack}} | ||
| 113 | test: | ||
| 114 | suffix: 1_arpack | ||
| 115 | args: -type arpack | ||
| 116 | requires: arpack | ||
| 117 | test: | ||
| 118 | suffix: 1_primme | ||
| 119 | args: -type primme | ||
| 120 | requires: primme | ||
| 121 | |||
| 122 | test: | ||
| 123 | suffix: 2 | ||
| 124 | args: -type {{rqcg lobpcg}} | ||
| 125 | |||
| 126 | TEST*/ | ||
| 127 |