Home
last modified time | relevance | path

Searched refs:J2 (Results 1 – 4 of 4) sorted by relevance

/petsc/src/ksp/ksp/tutorials/ !
H A Dex29.c244 Mat J2; in ComputeMatrix() local
249 PetscCall(MatConvert(jac, MATIS, MAT_INITIAL_MATRIX, &J2)); in ComputeMatrix()
250 PetscCall(MatViewFromOptions(J2, NULL, "-view_conv")); in ComputeMatrix()
251 PetscCall(MatConvert(J2, jtype, MAT_INPLACE_MATRIX, &J2)); in ComputeMatrix()
253 PetscCall(MatSetOperation(J2, MATOP_VIEW, f)); in ComputeMatrix()
254 PetscCall(MatSetDM(J2, da)); in ComputeMatrix()
255 PetscCall(MatViewFromOptions(J2, NULL, "-view_conv_assembled")); in ComputeMatrix()
256 PetscCall(MatAXPY(J2, -1., jac, DIFFERENT_NONZERO_PATTERN)); in ComputeMatrix()
257 PetscCall(MatNorm(J2, NORM_FROBENIUS, &nrm)); in ComputeMatrix()
259 PetscCall(MatViewFromOptions(J2, NULL, "-view_conv_err")); in ComputeMatrix()
[all …]
H A Dex34.c283 Mat J2; in ComputeMatrix() local
288 PetscCall(MatConvert(jac, MATIS, MAT_INITIAL_MATRIX, &J2)); in ComputeMatrix()
289 PetscCall(MatViewFromOptions(J2, NULL, "-view_conv")); in ComputeMatrix()
290 PetscCall(MatConvert(J2, jtype, MAT_INPLACE_MATRIX, &J2)); in ComputeMatrix()
292 PetscCall(MatSetOperation(J2, MATOP_VIEW, f)); in ComputeMatrix()
293 PetscCall(MatSetDM(J2, da)); in ComputeMatrix()
294 PetscCall(MatViewFromOptions(J2, NULL, "-view_conv_assembled")); in ComputeMatrix()
295 PetscCall(MatAXPY(J2, -1., jac, DIFFERENT_NONZERO_PATTERN)); in ComputeMatrix()
296 PetscCall(MatNorm(J2, NORM_FROBENIUS, &nrm)); in ComputeMatrix()
298 PetscCall(MatViewFromOptions(J2, NULL, "-view_conv_err")); in ComputeMatrix()
[all …]
/petsc/src/binding/petsc4py/test/ !
H A Dtest_snes.py174 J2, P2, jac2 = self.snes.getJacobian()
175 self.assertEqual(J, J2)
176 self.assertEqual(J2, P2)
/petsc/src/ts/tutorials/ !
H A Dex30.c75 …const PetscReal J2 = (PetscSqr(a11 - a22) + PetscSqr(a22 - a33) + PetscSqr(a33 - a11)) / 6 + Pets… in Eigenvalues_Sym3x3() local
76 const PetscReal s = PetscSqrtReal(J2 / 3); in Eigenvalues_Sym3x3()
93 T[0] = T[0] - 2.0 * J2 / 3.0; in Eigenvalues_Sym3x3()
94 T[3] = T[3] - 2.0 * J2 / 3.0; in Eigenvalues_Sym3x3()
95 T[5] = T[5] - 2.0 * J2 / 3.0; in Eigenvalues_Sym3x3()
103 x[0] = -PetscSqrtReal(J2); in Eigenvalues_Sym3x3()
105 x[2] = PetscSqrtReal(J2); in Eigenvalues_Sym3x3()
113 const PetscReal sd = PetscSqrtReal(J2) * PetscSinReal(alpha); in Eigenvalues_Sym3x3()