| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex45.c | 78 PetscScalar Hx, Hy, Hz, HxHydHz, HyHzdHx, HxHzdHy; in ComputeRHS() local 85 Hx = 1.0 / (PetscReal)(mx - 1); in ComputeRHS() 88 HxHydHz = Hx * Hy / Hz; in ComputeRHS() 89 HxHzdHy = Hx * Hz / Hy; in ComputeRHS() 90 HyHzdHx = Hy * Hz / Hx; in ComputeRHS() 100 barray[k][j][i] = Hx * Hy * Hz; in ComputeRHS() 121 PetscScalar v[7], Hx, Hy, Hz, HxHydHz, HyHzdHx, HxHzdHy; in ComputeMatrix() local 128 Hx = 1.0 / (PetscReal)(mx - 1); in ComputeMatrix() 131 HxHydHz = Hx * Hy / Hz; in ComputeMatrix() 132 HxHzdHy = Hx * Hz / Hy; in ComputeMatrix() [all …]
|
| H A D | ex28.c | 48 PetscScalar Hx, val[2]; in ComputeInitialSolution() local 52 Hx = 2.0 * PETSC_PI / (PetscReal)mx; in ComputeInitialSolution() 58 val[0] = val[1] = PetscSinScalar((PetscScalar)i * Hx); in ComputeInitialSolution() 86 PetscScalar v[7], Hx; in ComputeMatrix() local 95 Hx = 2.0 * PETSC_PI / (PetscReal)mx; in ComputeMatrix() 97 lambda = 2.0 * Hx; in ComputeMatrix() 103 v[0] = Hx; in ComputeMatrix() 121 v[1] = Hx; in ComputeMatrix()
|
| H A D | ex66.c | 67 PetscScalar Hx, Hy, pi, uu, tt; in ComputeRHS() local 78 Hx = 1.0 / (PetscReal)M; in ComputeRHS() 84 …scCosScalar(uu * pi * ((PetscReal)i + 0.5) * Hx) * PetscCosScalar(tt * pi * ((PetscReal)j + 0.5) *… in ComputeRHS() 101 PetscScalar v[5], Hx, Hy, HydHx, HxdHy; in ComputeJacobian() local 109 Hx = 1.0 / (PetscReal)M; in ComputeJacobian() 111 HxdHy = Hx / Hy; in ComputeJacobian() 112 HydHx = Hy / Hx; in ComputeJacobian()
|
| H A D | ex29.c | 104 PetscScalar Hx, Hy, HydHx, HxdHy; in ComputeRHS() local 111 Hx = 1.0 / (PetscReal)(mx - 1); in ComputeRHS() 113 HxdHy = Hx / Hy; in ComputeRHS() 114 HydHx = Hy / Hx; in ComputeRHS() 120 …array[j][i] = PetscExpScalar(-((PetscReal)i * Hx) * ((PetscReal)i * Hx) / user->nu) * PetscExpScal… in ComputeRHS() 122 …r(-((PetscReal)i * Hx) * ((PetscReal)i * Hx) / user->nu) * PetscExpScalar(-((PetscReal)j * Hy) * (… in ComputeRHS() 159 PetscReal Hx, Hy, HydHx, HxdHy, rho; in ComputeMatrix() local 168 Hx = 1.0 / (PetscReal)(mx - 1); in ComputeMatrix() 170 HxdHy = Hx / Hy; in ComputeMatrix() 171 HydHx = Hy / Hx; in ComputeMatrix()
|
| H A D | ex32.c | 87 PetscScalar Hx, Hy; in ComputeRHS() local 94 Hx = 1.0 / (PetscReal)mx; in ComputeRHS() 99 …0.5) * Hx) * (((PetscReal)i + 0.5) * Hx) / user->nu) * PetscExpScalar(-(((PetscReal)j + 0.5) * Hy)… in ComputeRHS() 121 PetscScalar v[5], Hx, Hy, HydHx, HxdHy; in ComputeMatrix() local 128 Hx = 1.0 / (PetscReal)mx; in ComputeMatrix() 130 HxdHy = Hx / Hy; in ComputeMatrix() 131 HydHx = Hy / Hx; in ComputeMatrix()
|
| H A D | ex50.c | 71 PetscScalar Hx, Hy, pi, uu, tt; in ComputeRHS() local 82 Hx = 1.0 / (PetscReal)M; in ComputeRHS() 88 …scCosScalar(uu * pi * ((PetscReal)i + 0.5) * Hx) * PetscCosScalar(tt * pi * ((PetscReal)j + 0.5) *… in ComputeRHS() 105 PetscScalar v[5], Hx, Hy, HydHx, HxdHy; in ComputeJacobian() local 113 Hx = 1.0 / (PetscReal)M; in ComputeJacobian() 115 HxdHy = Hx / Hy; in ComputeJacobian() 116 HydHx = Hy / Hx; in ComputeJacobian()
|
| H A D | ex45f.F90 | 29 PetscScalar h, Hx, Hy 36 Hx = 1.0/real(mx - 1) 38 h = Hx*Hy 51 PetscScalar v(5), Hx, Hy 60 Hx = 1.0/real(mx - 1) 62 HxdHy = Hx/Hy 63 HydHx = Hy/Hx
|
| H A D | ex34.c | 32 PetscScalar Hx, Hy, Hz; in main() local 64 Hx = 1.0 / (PetscReal)mx; in main() 73 …array[k][j][i][d] -= PetscCosScalar(2 * PETSC_PI * (((PetscReal)i + 0.5) * Hx)) * PetscCosScalar(2… in main() 98 PetscScalar Hx, Hy, Hz; in ComputeRHS() local 106 Hx = 1.0 / (PetscReal)mx; in ComputeRHS() 115 …i + 0.5) * Hx)) * PetscCosScalar(2 * PETSC_PI * (((PetscReal)j + 0.5) * Hy)) * PetscCosScalar(2 * … in ComputeRHS() 136 PetscScalar v[7], Hx, Hy, Hz, HyHzdHx, HxHzdHy, HxHydHz; in ComputeMatrix() local 145 Hx = 1.0 / (PetscReal)mx; in ComputeMatrix() 148 HyHzdHx = Hy * Hz / Hx; in ComputeMatrix() 149 HxHzdHy = Hx * Hz / Hy; in ComputeMatrix() [all …]
|
| H A D | ex22f.F90 | 44 PetscScalar v(7), Hx, Hy, Hz 54 Hx = 1.0/real(mx - 1) 57 HxHydHz = Hx*Hy/Hz 58 HxHzdHy = Hx*Hz/Hy 59 HyHzdHx = Hy*Hz/Hx
|
| H A D | ex73.c | 943 PetscScalar Hx, Hy; in ComputeRHS_DMDA() local 951 Hx = 1.0 / (PetscReal)(mx - 1); in ComputeRHS_DMDA() 956 …r(-((PetscReal)i * Hx) * ((PetscReal)i * Hx) / user->nu) * PetscExpScalar(-((PetscReal)j * Hy) * (… in ComputeRHS_DMDA() 991 PetscReal Hx, Hy, HydHx, HxdHy, rho; in ComputeMatrix_DMDA() local 1001 Hx = 1.0 / (PetscReal)(mx - 1); in ComputeMatrix_DMDA() 1003 HxdHy = Hx / Hy; in ComputeMatrix_DMDA() 1004 HydHx = Hy / Hx; in ComputeMatrix_DMDA()
|
| /petsc/src/ksp/ksp/tests/ |
| H A D | ex32.c | 128 PetscScalar *v, *v_neighbor, Hx, Hy, Hz, HxHydHz, HyHzdHx, HxHzdHy; in ComputeMatrix() local 136 Hx = 1.0 / (PetscReal)(mx - 1); in ComputeMatrix() 139 HxHydHz = Hx * Hy / Hz; in ComputeMatrix() 140 HxHzdHy = Hx * Hz / Hy; in ComputeMatrix() 141 HyHzdHx = Hy * Hz / Hx; in ComputeMatrix()
|
| /petsc/src/tao/constrained/tutorials/ |
| H A D | tomographyADMM.c | 20 Mat A, ATA, H, Hx, D, Hz, DTD, HF; member 267 PetscCall(MatTransposeMatMult(user->A, user->A, MAT_INITIAL_MATRIX, PETSC_DETERMINE, &user->Hx)); in InitializeUserData() 272 PetscCall(MatAssemblyBegin(user->Hx, MAT_FINAL_ASSEMBLY)); in InitializeUserData() 273 PetscCall(MatAssemblyEnd(user->Hx, MAT_FINAL_ASSEMBLY)); in InitializeUserData() 298 PetscCall(MatDestroy(&user->Hx)); in DestroyContext() 342 PetscCall(TaoADMMSetMisfitHessianRoutine(tao, user->Hx, user->Hx, HessianMisfit, (void *)user)); in main()
|
| /petsc/src/mat/tests/ |
| H A D | ex129.c | 186 PetscScalar *v, *v_neighbor, Hx, Hy, Hz, HxHydHz, HyHzdHx, HxHzdHy, r1, r2; in ComputeMatrix() local 200 Hx = 1.0 / (PetscReal)(mx - 1); in ComputeMatrix() 203 HxHydHz = Hx * Hy / Hz; in ComputeMatrix() 204 HxHzdHy = Hx * Hz / Hy; in ComputeMatrix() 205 HyHzdHx = Hy * Hz / Hx; in ComputeMatrix()
|
| /petsc/src/tao/constrained/impls/admm/ |
| H A D | admm.h | 21 Mat Hx, Hxpre, Hz, Hzpre, ATA, BTB, JA, JApre, JB, JBpre; member
|
| H A D | admm.c | 353 …if (am->Hx && am->ops->misfithess) PetscCall(TaoSetHessian(am->subsolverX, am->Hx, am->Hx, SubHess… in TaoSolve_ADMM() 642 if (am->Hx) { in TaoDestroy_ADMM() 643 PetscCall(MatDestroy(&am->Hx)); in TaoDestroy_ADMM() 1142 PetscCall(MatDestroy(&am->Hx)); in TaoADMMSetMisfitHessianRoutine() 1143 am->Hx = H; in TaoADMMSetMisfitHessianRoutine()
|