Lines Matching refs:J
35 PetscScalar **J; in PetscAdolcComputeRHSJacobian() local
38 PetscCall(AdolcMalloc2(m, p, &J)); in PetscAdolcComputeRHSJacobian()
39 if (adctx->Seed) fov_forward(tag, m, n, p, u_vec, adctx->Seed, NULL, J); in PetscAdolcComputeRHSJacobian()
40 else jacobian(tag, m, n, u_vec, J); in PetscAdolcComputeRHSJacobian()
42 PetscCall(RecoverJacobian(A, INSERT_VALUES, m, p, adctx->Rec, J, NULL)); in PetscAdolcComputeRHSJacobian()
46 … if (fabs(J[i][j]) > 1.e-16) PetscCall(MatSetValues(A, 1, &i, 1, &j, &J[i][j], INSERT_VALUES)); in PetscAdolcComputeRHSJacobian()
50 PetscCall(AdolcFree2(J)); in PetscAdolcComputeRHSJacobian()
73 PetscScalar **J; in PetscAdolcComputeRHSJacobianLocal() local
76 PetscCall(AdolcMalloc2(m, p, &J)); in PetscAdolcComputeRHSJacobianLocal()
77 if (adctx->Seed) fov_forward(tag, m, n, p, u_vec, adctx->Seed, NULL, J); in PetscAdolcComputeRHSJacobianLocal()
78 else jacobian(tag, m, n, u_vec, J); in PetscAdolcComputeRHSJacobianLocal()
80 PetscCall(RecoverJacobianLocal(A, INSERT_VALUES, m, p, adctx->Rec, J, NULL)); in PetscAdolcComputeRHSJacobianLocal()
84 …if (fabs(J[i][j]) > 1.e-16) PetscCall(MatSetValuesLocal(A, 1, &i, 1, &j, &J[i][j], INSERT_VALUES)); in PetscAdolcComputeRHSJacobianLocal()
88 PetscCall(AdolcFree2(J)); in PetscAdolcComputeRHSJacobianLocal()
112 PetscScalar **J; in PetscAdolcComputeIJacobian() local
115 PetscCall(AdolcMalloc2(m, p, &J)); in PetscAdolcComputeIJacobian()
118 if (adctx->Seed) fov_forward(tag1, m, n, p, u_vec, adctx->Seed, NULL, J); in PetscAdolcComputeIJacobian()
119 else jacobian(tag1, m, n, u_vec, J); in PetscAdolcComputeIJacobian()
122 PetscCall(RecoverJacobian(A, INSERT_VALUES, m, p, adctx->Rec, J, NULL)); in PetscAdolcComputeIJacobian()
126 … if (fabs(J[i][j]) > 1.e-16) PetscCall(MatSetValues(A, 1, &i, 1, &j, &J[i][j], INSERT_VALUES)); in PetscAdolcComputeIJacobian()
134 if (adctx->Seed) fov_forward(tag2, m, n, p, u_vec, adctx->Seed, NULL, J); in PetscAdolcComputeIJacobian()
135 else jacobian(tag2, m, n, u_vec, J); in PetscAdolcComputeIJacobian()
137 PetscCall(RecoverJacobian(A, ADD_VALUES, m, p, adctx->Rec, J, &a)); in PetscAdolcComputeIJacobian()
141 if (fabs(J[i][j]) > 1.e-16) { in PetscAdolcComputeIJacobian()
142 J[i][j] *= a; in PetscAdolcComputeIJacobian()
143 PetscCall(MatSetValues(A, 1, &i, 1, &j, &J[i][j], ADD_VALUES)); in PetscAdolcComputeIJacobian()
150 PetscCall(AdolcFree2(J)); in PetscAdolcComputeIJacobian()
172 PetscScalar **J; in PetscAdolcComputeIJacobianIDMass() local
175 PetscCall(AdolcMalloc2(m, p, &J)); in PetscAdolcComputeIJacobianIDMass()
178 if (adctx->Seed) fov_forward(tag, m, n, p, u_vec, adctx->Seed, NULL, J); in PetscAdolcComputeIJacobianIDMass()
179 else jacobian(tag, m, n, u_vec, J); in PetscAdolcComputeIJacobianIDMass()
182 PetscCall(RecoverJacobian(A, INSERT_VALUES, m, p, adctx->Rec, J, NULL)); in PetscAdolcComputeIJacobianIDMass()
186 … if (fabs(J[i][j]) > 1.e-16) PetscCall(MatSetValues(A, 1, &i, 1, &j, &J[i][j], INSERT_VALUES)); in PetscAdolcComputeIJacobianIDMass()
192 PetscCall(AdolcFree2(J)); in PetscAdolcComputeIJacobianIDMass()
217 PetscScalar **J; in PetscAdolcComputeIJacobianLocal() local
220 PetscCall(AdolcMalloc2(m, p, &J)); in PetscAdolcComputeIJacobianLocal()
223 if (adctx->Seed) fov_forward(tag1, m, n, p, u_vec, adctx->Seed, NULL, J); in PetscAdolcComputeIJacobianLocal()
224 else jacobian(tag1, m, n, u_vec, J); in PetscAdolcComputeIJacobianLocal()
226 PetscCall(RecoverJacobianLocal(A, INSERT_VALUES, m, p, adctx->Rec, J, NULL)); in PetscAdolcComputeIJacobianLocal()
230 …if (fabs(J[i][j]) > 1.e-16) PetscCall(MatSetValuesLocal(A, 1, &i, 1, &j, &J[i][j], INSERT_VALUES)); in PetscAdolcComputeIJacobianLocal()
238 if (adctx->Seed) fov_forward(tag2, m, n, p, u_vec, adctx->Seed, NULL, J); in PetscAdolcComputeIJacobianLocal()
239 else jacobian(tag2, m, n, u_vec, J); in PetscAdolcComputeIJacobianLocal()
241 PetscCall(RecoverJacobianLocal(A, ADD_VALUES, m, p, adctx->Rec, J, &a)); in PetscAdolcComputeIJacobianLocal()
245 if (fabs(J[i][j]) > 1.e-16) { in PetscAdolcComputeIJacobianLocal()
246 J[i][j] *= a; in PetscAdolcComputeIJacobianLocal()
247 PetscCall(MatSetValuesLocal(A, 1, &i, 1, &j, &J[i][j], ADD_VALUES)); in PetscAdolcComputeIJacobianLocal()
254 PetscCall(AdolcFree2(J)); in PetscAdolcComputeIJacobianLocal()
276 PetscScalar **J; in PetscAdolcComputeIJacobianLocalIDMass() local
279 PetscCall(AdolcMalloc2(m, p, &J)); in PetscAdolcComputeIJacobianLocalIDMass()
282 if (adctx->Seed) fov_forward(tag, m, n, p, u_vec, adctx->Seed, NULL, J); in PetscAdolcComputeIJacobianLocalIDMass()
283 else jacobian(tag, m, n, u_vec, J); in PetscAdolcComputeIJacobianLocalIDMass()
285 PetscCall(RecoverJacobianLocal(A, INSERT_VALUES, m, p, adctx->Rec, J, NULL)); in PetscAdolcComputeIJacobianLocalIDMass()
289 …if (fabs(J[i][j]) > 1.e-16) PetscCall(MatSetValuesLocal(A, 1, &i, 1, &j, &J[i][j], INSERT_VALUES)); in PetscAdolcComputeIJacobianLocalIDMass()
295 PetscCall(AdolcFree2(J)); in PetscAdolcComputeIJacobianLocalIDMass()
322 PetscScalar **J, *concat, **S; in PetscAdolcComputeRHSJacobianP() local
326 PetscCall(AdolcMalloc2(m, p, &J)); in PetscAdolcComputeRHSJacobianP()
334 fov_forward(tag, m, n + p, p, concat, S, NULL, J); in PetscAdolcComputeRHSJacobianP()
341 if (fabs(J[i][j]) > 1.e-16) PetscCall(MatSetValues(A, 1, &i, 1, &j, &J[i][j], INSERT_VALUES)); in PetscAdolcComputeRHSJacobianP()
344 PetscCall(AdolcFree2(J)); in PetscAdolcComputeRHSJacobianP()
366 PetscScalar **J, *concat, **S; in PetscAdolcComputeRHSJacobianPLocal() local
370 PetscCall(AdolcMalloc2(m, p, &J)); in PetscAdolcComputeRHSJacobianPLocal()
378 fov_forward(tag, m, n + p, p, concat, S, NULL, J); in PetscAdolcComputeRHSJacobianPLocal()
385 …if (fabs(J[i][j]) > 1.e-16) PetscCall(MatSetValuesLocal(A, 1, &i, 1, &j, &J[i][j], INSERT_VALUES)); in PetscAdolcComputeRHSJacobianPLocal()
388 PetscCall(AdolcFree2(J)); in PetscAdolcComputeRHSJacobianPLocal()
415 PetscScalar **J; in PetscAdolcComputeIJacobianAndDiagonalLocal() local
418 PetscCall(AdolcMalloc2(m, p, &J)); in PetscAdolcComputeIJacobianAndDiagonalLocal()
421 if (adctx->Seed) fov_forward(tag1, m, n, p, u_vec, adctx->Seed, NULL, J); in PetscAdolcComputeIJacobianAndDiagonalLocal()
422 else jacobian(tag1, m, n, u_vec, J); in PetscAdolcComputeIJacobianAndDiagonalLocal()
424 PetscCall(RecoverDiagonalLocal(diag, INSERT_VALUES, m, adctx->rec, J, NULL)); in PetscAdolcComputeIJacobianAndDiagonalLocal()
427 … if (fabs(J[i][i]) > 1.e-16) PetscCall(VecSetValuesLocal(diag, 1, &i, &J[i][i], INSERT_VALUES)); in PetscAdolcComputeIJacobianAndDiagonalLocal()
434 if (adctx->Seed) fov_forward(tag2, m, n, p, u_vec, adctx->Seed, NULL, J); in PetscAdolcComputeIJacobianAndDiagonalLocal()
435 else jacobian(tag2, m, n, u_vec, J); in PetscAdolcComputeIJacobianAndDiagonalLocal()
437 PetscCall(RecoverDiagonalLocal(diag, ADD_VALUES, m, adctx->rec, J, NULL)); in PetscAdolcComputeIJacobianAndDiagonalLocal()
440 if (fabs(J[i][i]) > 1.e-16) { in PetscAdolcComputeIJacobianAndDiagonalLocal()
441 J[i][i] *= a; in PetscAdolcComputeIJacobianAndDiagonalLocal()
442 PetscCall(VecSetValuesLocal(diag, 1, &i, &J[i][i], ADD_VALUES)); in PetscAdolcComputeIJacobianAndDiagonalLocal()
448 PetscCall(AdolcFree2(J)); in PetscAdolcComputeIJacobianAndDiagonalLocal()