| /petsc/src/ts/tutorials/autodiff/adolc-utils/ |
| H A D | drivers.cxx | 33 AdolcCtx *adctx = (AdolcCtx *)ctx; in PetscAdolcComputeRHSJacobian() local
34 PetscInt i, j, m = adctx->m, n = adctx->n, p = adctx->p; in PetscAdolcComputeRHSJacobian()
39 if (adctx->Seed) fov_forward(tag, m, n, p, u_vec, adctx->Seed, NULL, J); in PetscAdolcComputeRHSJacobian()
41 if (adctx->sparse) { in PetscAdolcComputeRHSJacobian()
42 PetscCall(RecoverJacobian(A, INSERT_VALUES, m, p, adctx->Rec, J, NULL)); in PetscAdolcComputeRHSJacobian()
71 AdolcCtx *adctx = (AdolcCtx *)ctx; in PetscAdolcComputeRHSJacobianLocal() local
72 PetscInt i, j, m = adctx->m, n = adctx->n, p = adctx->p; in PetscAdolcComputeRHSJacobianLocal()
77 if (adctx->Seed) fov_forward(tag, m, n, p, u_vec, adctx->Seed, NULL, J); in PetscAdolcComputeRHSJacobianLocal()
79 if (adctx->sparse) { in PetscAdolcComputeRHSJacobianLocal()
80 PetscCall(RecoverJacobianLocal(A, INSERT_VALUES, m, p, adctx->Rec, J, NULL)); in PetscAdolcComputeRHSJacobianLocal()
[all …]
|
| /petsc/src/ts/tutorials/autodiff/ |
| H A D | ex16adj_tl.cxx | 31 AdolcCtx *adctx; member
101 PetscCall(PetscMalloc1(user->adctx->n, &J)); in RHSJacobian()
106 for (i = 0; i < user->adctx->m; i++) { in RHSJacobian()
107 …for (j = 0; j < user->adctx->n; j++) PetscCall(MatSetValues(A, 1, &i, 1, &j, &J[i][j], INSERT_VALU… in RHSJacobian()
159 …for (i = 0; i < user->adctx->m; i++) PetscCall(MatSetValues(A, 1, &i, 1, &j, &J[i][user->adctx->n]… in RHSJacobianP()
198 AdolcCtx *adctx; in main() local
212 PetscCall(PetscNew(&adctx)); in main()
215 adctx->m = 2; in main()
216 adctx->n = 2; in main()
217 adctx->p = 2; in main()
[all …]
|
| H A D | ex16adj.cxx | 27 AdolcCtx *adctx; member
126 PetscCall(PetscAdolcComputeRHSJacobian(1, A, x, user->adctx)); in RHSJacobian()
141 PetscCall(PetscAdolcComputeRHSJacobianP(3, A, x, &user->mu, user->adctx)); in RHSJacobianP()
178 AdolcCtx *adctx; in main() local
192 PetscCall(PetscNew(&adctx)); in main()
195 adctx->m = 2; in main()
196 adctx->n = 2; in main()
197 adctx->p = 2; in main()
198 adctx->num_params = 1; in main()
199 user.adctx = adctx; in main()
[all …]
|
| H A D | adr_ex5adj.cxx | 58 AdolcCtx *adctx; /* Automatic differentation support */ member
78 AdolcCtx *adctx; in main() local
88 PetscCall(PetscNew(&adctx)); in main()
91 …appctx.aijpc = PETSC_FALSE, adctx->no_an = PETSC_FALSE, adctx->sparse = PETSC_FALSE, adctx->sparse… in main()
92 adctx->sparse_view_done = PETSC_FALSE; in main()
94 PetscCall(PetscOptionsGetBool(NULL, NULL, "-no_annotation", &adctx->no_an, NULL)); in main()
96 PetscCall(PetscOptionsGetBool(NULL, NULL, "-adolc_sparse", &adctx->sparse, NULL)); in main()
97 PetscCall(PetscOptionsGetBool(NULL, NULL, "-adolc_sparse_view", &adctx->sparse_view, NULL)); in main()
102 appctx.adctx = adctx; in main()
134 if (!adctx->no_an) { in main()
[all …]
|
| H A D | adr_ex1.cxx | 19 AdolcCtx *adctx; /* Automatic differentiation support */ member
144 PetscCall(PetscAdolcComputeIJacobian(1, 2, A, u, a, appctx->adctx)); in IJacobian()
179 AdolcCtx *adctx; in main() local
190 PetscCall(PetscNew(&adctx)); in main()
191 adctx->m = n; in main()
192 adctx->n = n; in main()
193 adctx->p = n; in main()
194 ctx.adctx = adctx; in main()
283 PetscCall(PetscFree(adctx)); in main()
|
| H A D | ex16opt_ic.cxx | 32 AdolcCtx *adctx; member
98 PetscCall(PetscAdolcComputeRHSJacobian(1, A, x, user->adctx)); in RHSJacobian()
131 AdolcCtx *adctx; in main() local
147 PetscCall(PetscNew(&adctx)); in main()
152 adctx->m = 2; in main()
153 adctx->n = 2; in main()
154 adctx->p = 2; in main()
155 user.adctx = adctx; in main()
266 PetscCall(PetscFree(adctx)); in main()
|
| H A D | adr_ex5adj_mf.cxx | 33 AdolcCtx *adctx; /* Automatic differentation support */ member
118 PetscCall(PetscNew(&appctx.adctx)); in main()
120 PetscCall(PetscFree(appctx.adctx)); in main()
|