15   Physics_MR phys;  in PhysicsContext_MR()  local
20   PetscCall(PetscMalloc1(1, &phys));  in PhysicsContext_MR()
21   PetscCall(ProcessPhysics_MR(comm, phys, *units));  in PhysicsContext_MR()
23   CeedQFunctionContextSetData(*ctx, CEED_MEM_HOST, CEED_COPY_VALUES, sizeof(*phys), phys);  in PhysicsContext_MR()
24   PetscCall(PetscFree(phys));  in PhysicsContext_MR()
33   Physics_MR  phys, phys_smoother;  in PhysicsSmootherContext_MR()  local
48     CeedQFunctionContextGetData(ctx, CEED_MEM_HOST, &phys);  in PhysicsSmootherContext_MR()
50     PetscCall(PetscMemcpy(phys_smoother, phys, sizeof(*phys)));  in PhysicsSmootherContext_MR()
51     CeedQFunctionContextRestoreData(ctx, &phys);  in PhysicsSmootherContext_MR()
65 PetscErrorCode ProcessPhysics_MR(MPI_Comm comm, Physics_MR phys, Units units) {  in ProcessPhysics_MR()  argument
67   phys->mu_1      = -1;  in ProcessPhysics_MR()
68   phys->mu_2      = -1;  in ProcessPhysics_MR()
69   phys->lambda    = -1;  in ProcessPhysics_MR()
78 …PetscCall(PetscOptionsScalar("-mu_1", "Material Property mu_1", NULL, phys->mu_1, &phys->mu_1, NUL…  in ProcessPhysics_MR()
79 …if (phys->mu_1 < 0) SETERRQ(PETSC_COMM_SELF, PETSC_ERR_SUP, "Mooney-Rivlin model requires non-nega…  in ProcessPhysics_MR()
81 …PetscCall(PetscOptionsScalar("-mu_2", "Material Property mu_2", NULL, phys->mu_2, &phys->mu_2, NUL…  in ProcessPhysics_MR()
82 …if (phys->mu_2 < 0) SETERRQ(PETSC_COMM_SELF, PETSC_ERR_SUP, "Mooney-Rivlin model requires non-nega…  in ProcessPhysics_MR()
86   phys->lambda = 2 * (phys->mu_1 + phys->mu_2) * nu / (1 - 2 * nu);  in ProcessPhysics_MR()
103   phys->mu_1 *= units->Pascal;  in ProcessPhysics_MR()
104   phys->mu_2 *= units->Pascal;  in ProcessPhysics_MR()
105   phys->lambda *= units->Pascal;  in ProcessPhysics_MR()