| 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[] = "Use the matrix exponential to compute rightmost eigenvalues.\n\n" | ||
| 12 | "Same problem as ex9.c but with explicitly created matrix. The command line options are:\n" | ||
| 13 | " -n <n>, where <n> = block dimension of the 2x2 block matrix.\n" | ||
| 14 | " -L <L>, where <L> = bifurcation parameter.\n" | ||
| 15 | " -alpha <alpha>, -beta <beta>, -delta1 <delta1>, -delta2 <delta2>,\n" | ||
| 16 | " where <alpha> <beta> <delta1> <delta2> = model parameters.\n\n"; | ||
| 17 | |||
| 18 | #include <slepceps.h> | ||
| 19 | #include <slepcmfn.h> | ||
| 20 | |||
| 21 | /* | ||
| 22 | This example computes the eigenvalues with largest real part of the | ||
| 23 | following matrix | ||
| 24 | |||
| 25 | A = [ tau1*T+(beta-1)*I alpha^2*I | ||
| 26 | -beta*I tau2*T-alpha^2*I ], | ||
| 27 | |||
| 28 | where | ||
| 29 | |||
| 30 | T = tridiag{1,-2,1} | ||
| 31 | h = 1/(n+1) | ||
| 32 | tau1 = delta1/(h*L)^2 | ||
| 33 | tau2 = delta2/(h*L)^2 | ||
| 34 | |||
| 35 | but it builds A explicitly, as opposed to ex9.c | ||
| 36 | */ | ||
| 37 | |||
| 38 | /* Routines for shell spectral transformation */ | ||
| 39 | PetscErrorCode STApply_Exp(ST,Vec,Vec); | ||
| 40 | PetscErrorCode STBackTransform_Exp(ST,PetscInt,PetscScalar*,PetscScalar*); | ||
| 41 | PetscErrorCode STDestroy_Exp(void**); | ||
| 42 | |||
| 43 | 60 | int main(int argc,char **argv) | |
| 44 | { | ||
| 45 | 60 | Mat A; /* operator matrix */ | |
| 46 | 60 | EPS eps; /* eigenproblem solver context */ | |
| 47 | 60 | ST st; /* spectral transformation context */ | |
| 48 | 60 | MFN mfn; /* matrix function solver object to compute exp(A)*v */ | |
| 49 | 60 | FN f; | |
| 50 | 60 | EPSType type; | |
| 51 | 60 | PetscScalar alpha,beta,tau1,tau2,delta1,delta2,L,h; | |
| 52 | 60 | PetscInt n=30,i,Istart,Iend,nev; | |
| 53 | 60 | PetscBool isShell,terse; | |
| 54 | |||
| 55 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
60 | PetscFunctionBeginUser; |
| 56 |
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.
|
60 | PetscCall(SlepcInitialize(&argc,&argv,NULL,help)); |
| 57 | #if defined(PETSC_HAVE_COMPLEX) | ||
| 58 |
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.
|
60 | PetscCall(PetscOptionsGetInt(NULL,NULL,"-n",&n,NULL)); |
| 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.
|
60 | PetscCall(PetscPrintf(PETSC_COMM_WORLD,"\nBrusselator wave model with matrix exponential, n=%" PetscInt_FMT "\n\n",n)); |
| 60 | |||
| 61 | /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - | ||
| 62 | Generate the matrix | ||
| 63 | - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */ | ||
| 64 | |||
| 65 | 60 | alpha = 2.0; | |
| 66 | 60 | beta = 5.45; | |
| 67 | 60 | delta1 = 0.008; | |
| 68 | 60 | delta2 = 0.004; | |
| 69 | 60 | L = 0.51302; | |
| 70 | |||
| 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.
|
60 | PetscCall(PetscOptionsGetScalar(NULL,NULL,"-L",&L,NULL)); |
| 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.
|
60 | PetscCall(PetscOptionsGetScalar(NULL,NULL,"-alpha",&alpha,NULL)); |
| 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.
|
60 | PetscCall(PetscOptionsGetScalar(NULL,NULL,"-beta",&beta,NULL)); |
| 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.
|
60 | PetscCall(PetscOptionsGetScalar(NULL,NULL,"-delta1",&delta1,NULL)); |
| 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.
|
60 | PetscCall(PetscOptionsGetScalar(NULL,NULL,"-delta2",&delta2,NULL)); |
| 76 | |||
| 77 | 60 | h = 1.0 / (PetscReal)(n+1); | |
| 78 | 60 | tau1 = delta1 / ((h*L)*(h*L)); | |
| 79 | 60 | tau2 = delta2 / ((h*L)*(h*L)); | |
| 80 | |||
| 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.
|
60 | PetscCall(MatCreate(PETSC_COMM_WORLD,&A)); |
| 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.
|
60 | PetscCall(MatSetSizes(A,PETSC_DECIDE,PETSC_DECIDE,2*n,2*n)); |
| 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.
|
60 | PetscCall(MatSetFromOptions(A)); |
| 84 | |||
| 85 |
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.
|
60 | PetscCall(MatGetOwnershipRange(A,&Istart,&Iend)); |
| 86 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
7140 | for (i=Istart;i<Iend;i++) { |
| 87 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
7080 | if (i<n) { /* upper blocks */ |
| 88 |
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.
|
3540 | if (i>0) PetscCall(MatSetValue(A,i,i-1,tau1,INSERT_VALUES)); |
| 89 |
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.
|
3540 | if (i<n-1) PetscCall(MatSetValue(A,i,i+1,tau1,INSERT_VALUES)); |
| 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.
|
3540 | PetscCall(MatSetValue(A,i,i,-2.0*tau1+beta-1.0,INSERT_VALUES)); |
| 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.
|
3540 | PetscCall(MatSetValue(A,i,i+n,alpha*alpha,INSERT_VALUES)); |
| 92 | } else { /* lower blocks */ | ||
| 93 |
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.
|
3540 | if (i>n) PetscCall(MatSetValue(A,i,i-1,tau2,INSERT_VALUES)); |
| 94 |
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.
|
3540 | if (i<2*n-1) PetscCall(MatSetValue(A,i,i+1,tau2,INSERT_VALUES)); |
| 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.
|
3540 | PetscCall(MatSetValue(A,i,i,-2.0*tau2-alpha*alpha,INSERT_VALUES)); |
| 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.
|
7080 | PetscCall(MatSetValue(A,i,i-n,-beta,INSERT_VALUES)); |
| 97 | } | ||
| 98 | } | ||
| 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.
|
60 | PetscCall(MatAssemblyBegin(A,MAT_FINAL_ASSEMBLY)); |
| 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.
|
60 | PetscCall(MatAssemblyEnd(A,MAT_FINAL_ASSEMBLY)); |
| 101 | |||
| 102 | /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - | ||
| 103 | Create the eigensolver and set various options | ||
| 104 | - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */ | ||
| 105 | |||
| 106 |
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.
|
60 | PetscCall(EPSCreate(PETSC_COMM_WORLD,&eps)); |
| 107 |
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.
|
60 | PetscCall(EPSSetOperators(eps,A,NULL)); |
| 108 |
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.
|
60 | PetscCall(EPSSetProblemType(eps,EPS_NHEP)); |
| 109 |
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.
|
60 | PetscCall(EPSSetWhichEigenpairs(eps,EPS_LARGEST_REAL)); |
| 110 |
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.
|
60 | PetscCall(EPSGetST(eps,&st)); |
| 111 |
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.
|
60 | PetscCall(STSetType(st,STSHELL)); |
| 112 |
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.
|
60 | PetscCall(EPSSetFromOptions(eps)); |
| 113 | |||
| 114 | /* | ||
| 115 | Initialize shell spectral transformation | ||
| 116 | */ | ||
| 117 |
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.
|
60 | PetscCall(PetscObjectTypeCompare((PetscObject)st,STSHELL,&isShell)); |
| 118 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
60 | if (isShell) { |
| 119 | |||
| 120 | /* Create the MFN object to be used by the spectral transform */ | ||
| 121 |
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(MFNCreate(PETSC_COMM_WORLD,&mfn)); |
| 122 |
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(MFNSetOperator(mfn,A)); |
| 123 |
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(MFNGetFN(mfn,&f)); |
| 124 |
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(FNSetType(f,FNEXP)); |
| 125 |
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(FNSetScale(f,0.03,1.0)); /* this can be set with -fn_scale */ |
| 126 |
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(MFNSetFromOptions(mfn)); |
| 127 | |||
| 128 | /* Set callback functions */ | ||
| 129 |
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(STShellSetApply(st,STApply_Exp)); |
| 130 |
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(STShellSetBackTransform(st,STBackTransform_Exp)); |
| 131 |
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(STShellSetContext(st,mfn)); |
| 132 |
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(STShellSetContextDestroy(st,STDestroy_Exp)); |
| 133 |
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(PetscObjectSetName((PetscObject)st,"STEXP")); |
| 134 | } | ||
| 135 | |||
| 136 | /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - | ||
| 137 | Solve the eigensystem | ||
| 138 | - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */ | ||
| 139 | |||
| 140 |
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.
|
60 | PetscCall(EPSSolve(eps)); |
| 141 |
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.
|
60 | PetscCall(EPSGetType(eps,&type)); |
| 142 |
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.
|
60 | PetscCall(PetscPrintf(PETSC_COMM_WORLD," Solution method: %s\n\n",type)); |
| 143 |
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.
|
60 | PetscCall(EPSGetDimensions(eps,&nev,NULL,NULL)); |
| 144 |
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.
|
60 | PetscCall(PetscPrintf(PETSC_COMM_WORLD," Number of requested eigenvalues: %" PetscInt_FMT "\n",nev)); |
| 145 | |||
| 146 | /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - | ||
| 147 | Display solution and clean up | ||
| 148 | - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */ | ||
| 149 | |||
| 150 | /* show detailed info unless -terse option is given by user */ | ||
| 151 |
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.
|
60 | PetscCall(PetscOptionsHasName(NULL,NULL,"-terse",&terse)); |
| 152 |
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.
|
60 | if (terse) PetscCall(EPSErrorView(eps,EPS_ERROR_RELATIVE,NULL)); |
| 153 | else { | ||
| 154 | ✗ | PetscCall(PetscViewerPushFormat(PETSC_VIEWER_STDOUT_WORLD,PETSC_VIEWER_ASCII_INFO_DETAIL)); | |
| 155 | ✗ | PetscCall(EPSConvergedReasonView(eps,PETSC_VIEWER_STDOUT_WORLD)); | |
| 156 | ✗ | PetscCall(EPSErrorView(eps,EPS_ERROR_RELATIVE,PETSC_VIEWER_STDOUT_WORLD)); | |
| 157 | ✗ | PetscCall(PetscViewerPopFormat(PETSC_VIEWER_STDOUT_WORLD)); | |
| 158 | } | ||
| 159 |
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.
|
60 | PetscCall(EPSDestroy(&eps)); |
| 160 |
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.
|
60 | PetscCall(MatDestroy(&A)); |
| 161 | #else | ||
| 162 | SETERRQ(PETSC_COMM_WORLD,PETSC_ERR_SUP,"This example requires C99 complex numbers"); | ||
| 163 | #endif | ||
| 164 |
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.
|
60 | PetscCall(SlepcFinalize()); |
| 165 | return 0; | ||
| 166 | } | ||
| 167 | |||
| 168 | /* ------------------------------------------------------------------- */ | ||
| 169 | /* | ||
| 170 | STBackTransform_Exp - Undoes the exp(A) transformation by taking logarithms. | ||
| 171 | |||
| 172 | Input Parameters: | ||
| 173 | + st - spectral transformation context | ||
| 174 | - n - number of eigenvalues to transform | ||
| 175 | |||
| 176 | Input/Output Parameters: | ||
| 177 | + eigr - pointer to real part of eigenvalues | ||
| 178 | - eigi - pointer to imaginary part of eigenvalues | ||
| 179 | */ | ||
| 180 | 5890 | PetscErrorCode STBackTransform_Exp(ST st,PetscInt n,PetscScalar *eigr,PetscScalar *eigi) | |
| 181 | { | ||
| 182 | #if defined(PETSC_HAVE_COMPLEX) | ||
| 183 | 5890 | PetscInt j; | |
| 184 | 5890 | MFN mfn; | |
| 185 | 5890 | FN fn; | |
| 186 | 5890 | PetscScalar tau,eta; | |
| 187 | #if !defined(PETSC_USE_COMPLEX) | ||
| 188 | 1785 | PetscComplex theta,lambda; | |
| 189 | #endif | ||
| 190 | |||
| 191 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
5890 | PetscFunctionBeginUser; |
| 192 |
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.
|
5890 | PetscCall(STShellGetContext(st,&mfn)); |
| 193 |
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.
|
5890 | PetscCall(MFNGetFN(mfn,&fn)); |
| 194 |
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.
|
5890 | PetscCall(FNGetScale(fn,&tau,&eta)); |
| 195 |
2/2✓ Branch 0 taken 10 times.
✓ Branch 1 taken 10 times.
|
17710 | for (j=0;j<n;j++) { |
| 196 | #if defined(PETSC_USE_COMPLEX) | ||
| 197 | 8230 | eigr[j] = PetscLogComplex(eigr[j]/eta)/tau; | |
| 198 | #else | ||
| 199 | 3590 | theta = PetscCMPLX(eigr[j],eigi[j])/eta; | |
| 200 | 3590 | lambda = PetscLogComplex(theta)/tau; | |
| 201 | 3590 | eigr[j] = PetscRealPartComplex(lambda); | |
| 202 | 3590 | eigi[j] = PetscImaginaryPartComplex(lambda); | |
| 203 | #endif | ||
| 204 | } | ||
| 205 |
5/12✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 2 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
✗ Branch 10 not taken.
✗ Branch 11 not taken.
|
1178 | PetscFunctionReturn(PETSC_SUCCESS); |
| 206 | #else | ||
| 207 | return 0; | ||
| 208 | #endif | ||
| 209 | } | ||
| 210 | |||
| 211 | /* | ||
| 212 | STApply_Exp - Applies the operator exp(tau*A) to a given vector using an MFN object. | ||
| 213 | |||
| 214 | Input Parameters: | ||
| 215 | + st - spectral transformation context | ||
| 216 | - x - input vector | ||
| 217 | |||
| 218 | Output Parameter: | ||
| 219 | . y - output vector | ||
| 220 | */ | ||
| 221 | 570 | PetscErrorCode STApply_Exp(ST st,Vec x,Vec y) | |
| 222 | { | ||
| 223 | 570 | MFN mfn; | |
| 224 | |||
| 225 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
570 | PetscFunctionBeginUser; |
| 226 |
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.
|
570 | PetscCall(STShellGetContext(st,&mfn)); |
| 227 |
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.
|
570 | PetscCall(MFNSolve(mfn,x,y)); |
| 228 |
5/12✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 2 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
✗ Branch 10 not taken.
✗ Branch 11 not taken.
|
114 | PetscFunctionReturn(PETSC_SUCCESS); |
| 229 | } | ||
| 230 | |||
| 231 | /* | ||
| 232 | STDestroy_Exp - This routine destroys the shell ST context. | ||
| 233 | |||
| 234 | Input Parameter: | ||
| 235 | . ctx - user-defined spectral transformation context | ||
| 236 | */ | ||
| 237 | 10 | PetscErrorCode STDestroy_Exp(void **ctx) | |
| 238 | { | ||
| 239 | 10 | MFN mfn = (MFN)*ctx; | |
| 240 | |||
| 241 |
1/2✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
|
10 | PetscFunctionBeginUser; |
| 242 |
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(MFNDestroy(&mfn)); |
| 243 |
5/12✓ Branch 0 taken 2 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 2 times.
✓ Branch 4 taken 2 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 2 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 2 times.
✗ Branch 10 not taken.
✗ Branch 11 not taken.
|
2 | PetscFunctionReturn(PETSC_SUCCESS); |
| 244 | } | ||
| 245 | |||
| 246 | /*TEST | ||
| 247 | |||
| 248 | testset: | ||
| 249 | args: -eps_nev 4 -mfn_ncv 16 -terse | ||
| 250 | requires: c99_complex !single | ||
| 251 | filter: sed -e "s/-2/+2/g" | ||
| 252 | output_file: output/ex36_1.out | ||
| 253 | test: | ||
| 254 | suffix: 1 | ||
| 255 | requires: !__float128 | ||
| 256 | test: | ||
| 257 | suffix: 1_quad | ||
| 258 | args: -eps_tol 1e-14 | ||
| 259 | requires: __float128 | ||
| 260 | |||
| 261 | test: | ||
| 262 | suffix: 2 | ||
| 263 | args: -n 56 -eps_nev 2 -st_type sinvert -eps_target -390 -eps_target_magnitude -eps_type power | ||
| 264 | args: -eps_power_shift_type {{constant rayleigh}} -eps_two_sided {{0 1}} -eps_tol 1e-14 -terse | ||
| 265 | requires: c99_complex !single | ||
| 266 | filter: sed -e "s/[+-]0\.0*i//g" | ||
| 267 | |||
| 268 | test: | ||
| 269 | suffix: 3 | ||
| 270 | args: -n 100 -st_type sinvert -eps_type ciss -rg_type ellipse -rg_ellipse_center 0 -rg_ellipse_radius 6 -eps_all -eps_tol 1e-6 -terse | ||
| 271 | requires: c99_complex !single | ||
| 272 | filter: sed -e "s/-3.37036-3.55528i, -3.37036+3.55528i/-3.37036+3.55528i, -3.37036-3.55528i/" -e "s/-1.79853-3.03216i, -1.79853+3.03216i/-1.79853+3.03216i, -1.79853-3.03216i/" -e "s/-0.67471-2.52856i, -0.67471+2.52856i/-0.67471+2.52856i, -0.67471-2.52856i/" -e "s/0.00002-2.13950i, 0.00002+2.13950i/0.00002+2.13950i, 0.00002-2.13950i/" | ||
| 273 | |||
| 274 | TEST*/ | ||
| 275 |