Lines Matching refs:J
454 PetscCall(DMFieldEvaluateFE(field, pointIS, quad, PETSC_REAL, g->v, g->J, NULL)); in DMFieldCreateFEGeom()
461 PetscReal J[16] = {0}; in DMFieldCreateFEGeom() local
464 for (k = 0; k < dim; k++) J[j * dE + k] = g->J[i * dE * dim + j * dim + k]; in DMFieldCreateFEGeom()
469 for (k = 0; k < dE; k++) J[j * dE + k] = (j == k) ? 1. : 0.; in DMFieldCreateFEGeom()
474 PetscReal norm = PetscSqrtReal(J[0] * J[0] + J[2] * J[2]); in DMFieldCreateFEGeom()
476 J[1] = -J[2] / norm; in DMFieldCreateFEGeom()
477 J[3] = J[0] / norm; in DMFieldCreateFEGeom()
479 PetscReal inorm = 1. / PetscSqrtReal(J[0] * J[0] + J[3] * J[3] + J[6] * J[6]); in DMFieldCreateFEGeom()
480 PetscReal x = J[0] * inorm; in DMFieldCreateFEGeom()
481 PetscReal y = J[3] * inorm; in DMFieldCreateFEGeom()
482 PetscReal z = J[6] * inorm; in DMFieldCreateFEGeom()
487 J[1] = -y; in DMFieldCreateFEGeom()
488 J[2] = -z; in DMFieldCreateFEGeom()
489 J[4] = 1. - y * y * inv1pX; in DMFieldCreateFEGeom()
490 J[5] = -y * z * inv1pX; in DMFieldCreateFEGeom()
491 J[7] = -y * z * inv1pX; in DMFieldCreateFEGeom()
492 J[8] = 1. - z * z * inv1pX; in DMFieldCreateFEGeom()
496 J[1] = z; in DMFieldCreateFEGeom()
497 J[2] = y; in DMFieldCreateFEGeom()
498 J[4] = -y * z * inv1mX; in DMFieldCreateFEGeom()
499 J[5] = 1. - y * y * inv1mX; in DMFieldCreateFEGeom()
500 J[7] = 1. - z * z * inv1mX; in DMFieldCreateFEGeom()
501 J[8] = -y * z * inv1mX; in DMFieldCreateFEGeom()
508 J[2] = J[3] * J[7] - J[6] * J[4]; in DMFieldCreateFEGeom()
509 J[5] = J[6] * J[1] - J[0] * J[7]; in DMFieldCreateFEGeom()
510 J[8] = J[0] * J[4] - J[3] * J[1]; in DMFieldCreateFEGeom()
512 inorm = 1. / PetscSqrtReal(J[2] * J[2] + J[5] * J[5] + J[8] * J[8]); in DMFieldCreateFEGeom()
514 J[2] *= inorm; in DMFieldCreateFEGeom()
515 J[5] *= inorm; in DMFieldCreateFEGeom()
516 J[8] *= inorm; in DMFieldCreateFEGeom()
519 for (j = 0; j < dE * dE; j++) g->J[i * dE * dE + j] = J[j]; in DMFieldCreateFEGeom()