Lines Matching refs:adctx
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()
136 adctx->m = dofs * gxm * gym; in main()
137 adctx->n = dofs * gxm * gym; /* Number of dependent and independent variables */ in main()
155 if (adctx->sparse) { in main()
163 PetscCall(PetscMalloc1(adctx->n, &u_vec)); in main()
164 JP = (unsigned int **)malloc(adctx->m * sizeof(unsigned int *)); in main()
165 jac_pat(1, adctx->m, adctx->n, u_vec, JP, ctrl); in main()
166 if (adctx->sparse_view) PetscCall(PrintSparsity(MPI_COMM_WORLD, adctx->m, JP)); in main()
176 PetscCall(ISColoringGetIS(iscoloring, PETSC_USE_POINTER, &adctx->p, NULL)); in main()
179 PetscCall(AdolcMalloc2(adctx->n, adctx->p, &Seed)); in main()
186 PetscCall(AdolcMalloc2(adctx->m, adctx->p, &Rec)); in main()
187 PetscCall(GetRecoveryMatrix(Seed, JP, adctx->m, adctx->p, Rec)); in main()
190 adctx->Rec = Rec; in main()
191 for (i = 0; i < adctx->m; i++) free(JP[i]); in main()
200 adctx->p = adctx->n; in main()
201 PetscCall(AdolcMalloc2(adctx->n, adctx->p, &Seed)); in main()
202 PetscCall(Identity(adctx->n, Seed)); in main()
204 adctx->Seed = Seed; in main()
284 if (!adctx->no_an) { in main()
285 if (adctx->sparse) PetscCall(AdolcFree2(Rec)); in main()
289 PetscCall(PetscFree(adctx)); in main()
776 PetscCall(PetscAdolcComputeIJacobianLocalIDMass(1, A, u_vec, a, appctx->adctx)); in IJacobianAdolc()
1073 PetscCall(PetscAdolcComputeRHSJacobianLocal(1, A, u_vec, appctx->adctx)); in RHSJacobianAdolc()