512   RotCtx   *rc  = (RotCtx *)ctx;  in DMPlexBasisTransformSetUp_Rotation_Internal()  local
513   PetscInt  dim = rc->dim;  in DMPlexBasisTransformSetUp_Rotation_Internal()
517   PetscCall(PetscMalloc2(PetscSqr(dim), &rc->R, PetscSqr(dim), &rc->RT));  in DMPlexBasisTransformSetUp_Rotation_Internal()
520     c1       = PetscCosReal(rc->alpha);  in DMPlexBasisTransformSetUp_Rotation_Internal()
521     s1       = PetscSinReal(rc->alpha);  in DMPlexBasisTransformSetUp_Rotation_Internal()
522     rc->R[0] = c1;  in DMPlexBasisTransformSetUp_Rotation_Internal()
523     rc->R[1] = s1;  in DMPlexBasisTransformSetUp_Rotation_Internal()
524     rc->R[2] = -s1;  in DMPlexBasisTransformSetUp_Rotation_Internal()
525     rc->R[3] = c1;  in DMPlexBasisTransformSetUp_Rotation_Internal()
526     PetscCall(PetscArraycpy(rc->RT, rc->R, PetscSqr(dim)));  in DMPlexBasisTransformSetUp_Rotation_Internal()
527     DMPlex_Transpose2D_Internal(rc->RT);  in DMPlexBasisTransformSetUp_Rotation_Internal()
530     c1       = PetscCosReal(rc->alpha);  in DMPlexBasisTransformSetUp_Rotation_Internal()
531     s1       = PetscSinReal(rc->alpha);  in DMPlexBasisTransformSetUp_Rotation_Internal()
532     c2       = PetscCosReal(rc->beta);  in DMPlexBasisTransformSetUp_Rotation_Internal()
533     s2       = PetscSinReal(rc->beta);  in DMPlexBasisTransformSetUp_Rotation_Internal()
534     c3       = PetscCosReal(rc->gamma);  in DMPlexBasisTransformSetUp_Rotation_Internal()
535     s3       = PetscSinReal(rc->gamma);  in DMPlexBasisTransformSetUp_Rotation_Internal()
536     rc->R[0] = c1 * c3 - c2 * s1 * s3;  in DMPlexBasisTransformSetUp_Rotation_Internal()
537     rc->R[1] = c3 * s1 + c1 * c2 * s3;  in DMPlexBasisTransformSetUp_Rotation_Internal()
538     rc->R[2] = s2 * s3;  in DMPlexBasisTransformSetUp_Rotation_Internal()
539     rc->R[3] = -c1 * s3 - c2 * c3 * s1;  in DMPlexBasisTransformSetUp_Rotation_Internal()
540     rc->R[4] = c1 * c2 * c3 - s1 * s3;  in DMPlexBasisTransformSetUp_Rotation_Internal()
541     rc->R[5] = c3 * s2;  in DMPlexBasisTransformSetUp_Rotation_Internal()
542     rc->R[6] = s1 * s2;  in DMPlexBasisTransformSetUp_Rotation_Internal()
543     rc->R[7] = -c1 * s2;  in DMPlexBasisTransformSetUp_Rotation_Internal()
544     rc->R[8] = c2;  in DMPlexBasisTransformSetUp_Rotation_Internal()
545     PetscCall(PetscArraycpy(rc->RT, rc->R, PetscSqr(dim)));  in DMPlexBasisTransformSetUp_Rotation_Internal()
546     DMPlex_Transpose3D_Internal(rc->RT);  in DMPlexBasisTransformSetUp_Rotation_Internal()
556   RotCtx *rc = (RotCtx *)ctx;  in DMPlexBasisTransformDestroy_Rotation_Internal()  local
559   PetscCall(PetscFree2(rc->R, rc->RT));  in DMPlexBasisTransformDestroy_Rotation_Internal()
560   PetscCall(PetscFree(rc));  in DMPlexBasisTransformDestroy_Rotation_Internal()
566   RotCtx *rc = (RotCtx *)ctx;  in DMPlexBasisTransformGetMatrix_Rotation_Internal()  local
571     *A = rc->R;  in DMPlexBasisTransformGetMatrix_Rotation_Internal()
573     *A = rc->RT;  in DMPlexBasisTransformGetMatrix_Rotation_Internal()
874   RotCtx  *rc;  in DMPlexCreateBasisRotation()  local
879   PetscCall(PetscMalloc1(1, &rc));  in DMPlexCreateBasisRotation()
880   dm->transformCtx       = rc;  in DMPlexCreateBasisRotation()
884   rc->dim                = cdim;  in DMPlexCreateBasisRotation()
885   rc->alpha              = alpha;  in DMPlexCreateBasisRotation()
886   rc->beta               = beta;  in DMPlexCreateBasisRotation()
887   rc->gamma              = gamma;  in DMPlexCreateBasisRotation()