Lines Matching refs:phys
146 Physics phys = (Physics)ctx; in PhysicsBoundary_Advect_Inflow() local
147 Physics_Advect *advect = (Physics_Advect *)phys->data; in PhysicsBoundary_Advect_Inflow()
161 …cScalar *xR, PetscInt numConstants, const PetscScalar constants[], PetscScalar *flux, Physics phys) in PhysicsRiemann_Advect() argument
163 Physics_Advect *advect = (Physics_Advect *)phys->data; in PhysicsRiemann_Advect()
213 Physics phys = (Physics)ctx; in PhysicsSolution_Advect() local
214 Physics_Advect *advect = (Physics_Advect *)phys->data; in PhysicsSolution_Advect()
252 Physics phys = (Physics)ctx; in PhysicsFunctional_Advect() local
253 Physics_Advect *advect = (Physics_Advect *)phys->data; in PhysicsFunctional_Advect()
257 PetscCall(PhysicsSolution_Advect(mod, time, x, yexact, phys)); in PhysicsFunctional_Advect()
263 static PetscErrorCode SetUpBC_Advect(DM dm, PetscDS prob, Physics phys) in SetUpBC_Advect() argument
271 …nflowids), inflowids, 0, 0, NULL, (PetscVoidFn *)PhysicsBoundary_Advect_Inflow, NULL, phys, NULL)); in SetUpBC_Advect()
272 …lowids), outflowids, 0, 0, NULL, (PetscVoidFn *)PhysicsBoundary_Advect_Outflow, NULL, phys, NULL)); in SetUpBC_Advect()
276 static PetscErrorCode PhysicsCreate_Advect(Model mod, Physics phys, PetscOptionItems PetscOptionsOb… in PhysicsCreate_Advect() argument
281 phys->field_desc = PhysicsFields_Advect; in PhysicsCreate_Advect()
282 phys->riemann = (PetscRiemannFn *)(PetscVoidFn *)PhysicsRiemann_Advect; in PhysicsCreate_Advect()
284 phys->data = advect; in PhysicsCreate_Advect()
301 phys->maxspeed = Norm2Real(tilted->wind); in PhysicsCreate_Advect()
314 phys->maxspeed = 3.; /* radius of mesh, kludge */ in PhysicsCreate_Advect()
320 PetscCall(ModelSolutionSetDefault(mod, PhysicsSolution_Advect, phys)); in PhysicsCreate_Advect()
322 …FunctionalRegister(mod, "Solution", &advect->functional.Solution, PhysicsFunctional_Advect, phys)); in PhysicsCreate_Advect()
323 …(ModelFunctionalRegister(mod, "Error", &advect->functional.Error, PhysicsFunctional_Advect, phys)); in PhysicsCreate_Advect()
361 Physics phys = (Physics)ctx; in PhysicsFunctional_SW() local
362 Physics_SW *sw = (Physics_SW *)phys->data; in PhysicsFunctional_SW()
376 static PetscErrorCode SetUpBC_SW(DM dm, PetscDS prob, Physics phys) in SetUpBC_SW() argument
383 …Y_LENGTH(wallids), wallids, 0, 0, NULL, (PetscVoidFn *)PhysicsBoundary_SW_Wall, NULL, phys, NULL)); in SetUpBC_SW()
388 static PetscErrorCode CreateQFunctionContext_SW(Physics phys, Ceed ceed, CeedQFunctionContext *qfCt… in CreateQFunctionContext_SW() argument
390 Physics_SW *in = (Physics_SW *)phys->data; in CreateQFunctionContext_SW()
423 static PetscErrorCode PhysicsCreate_SW(Model mod, Physics phys, PetscOptionItems PetscOptionsObject) in PhysicsCreate_SW() argument
429 phys->field_desc = PhysicsFields_SW; in PhysicsCreate_SW()
431 phys->data = sw; in PhysicsCreate_SW()
448 phys->riemann = (PetscRiemannFn *)(PetscVoidFn *)PhysicsRiemann_SW; in PhysicsCreate_SW()
451 phys->maxspeed = PetscSqrtReal(2.0 * sw->gravity); /* Mach 1 for depth of 2 */ in PhysicsCreate_SW()
453 PetscCall(ModelSolutionSetDefault(mod, PhysicsSolution_SW, phys)); in PhysicsCreate_SW()
454 …scCall(ModelFunctionalRegister(mod, "Height", &sw->functional.Height, PhysicsFunctional_SW, phys)); in PhysicsCreate_SW()
455 …PetscCall(ModelFunctionalRegister(mod, "Speed", &sw->functional.Speed, PhysicsFunctional_SW, phys)… in PhysicsCreate_SW()
456 …scCall(ModelFunctionalRegister(mod, "Energy", &sw->functional.Energy, PhysicsFunctional_SW, phys)); in PhysicsCreate_SW()
476 Physics phys = (Physics)ctx; in PhysicsSolution_Euler() local
477 Physics_Euler *eu = (Physics_Euler *)phys->data; in PhysicsSolution_Euler()
530 if (c > phys->maxspeed) phys->maxspeed = c; in PhysicsSolution_Euler()
540 Physics phys = (Physics)ctx; in PhysicsBoundary_Euler_Wall() local
541 Physics_Euler *eu = (Physics_Euler *)phys->data; in PhysicsBoundary_Euler_Wall()
562 Physics phys = (Physics)ctx; in PhysicsFunctional_Euler() local
563 Physics_Euler *eu = (Physics_Euler *)phys->data; in PhysicsFunctional_Euler()
577 static PetscErrorCode SetUpBC_Euler(DM dm, PetscDS prob, Physics phys) in SetUpBC_Euler() argument
579 Physics_Euler *eu = (Physics_Euler *)phys->data; in SetUpBC_Euler()
586 …ENGTH(wallids), wallids, 0, 0, NULL, (PetscVoidFn *)PhysicsBoundary_Euler_Wall, NULL, phys, NULL)); in SetUpBC_Euler()
589 …ENGTH(wallids), wallids, 0, 0, NULL, (PetscVoidFn *)PhysicsBoundary_Euler_Wall, NULL, phys, NULL)); in SetUpBC_Euler()
595 static PetscErrorCode CreateQFunctionContext_Euler(Physics phys, Ceed ceed, CeedQFunctionContext *q… in CreateQFunctionContext_Euler() argument
597 Physics_Euler *in = (Physics_Euler *)phys->data; in CreateQFunctionContext_Euler()
630 static PetscErrorCode PhysicsCreate_Euler(Model mod, Physics phys, PetscOptionItems PetscOptionsObj… in PhysicsCreate_Euler() argument
635 phys->field_desc = PhysicsFields_Euler; in PhysicsCreate_Euler()
636 phys->riemann = (PetscRiemannFn *)(PetscVoidFn *)PhysicsRiemann_Euler_Godunov; in PhysicsCreate_Euler()
638 phys->data = eu; in PhysicsCreate_Euler()
660 phys->riemann = (PetscRiemannFn *)(PetscVoidFn *)PhysicsRiemann_Euler; in PhysicsCreate_Euler()
694 phys->maxspeed = 0.; /* will get set in solution */ in PhysicsCreate_Euler()
695 PetscCall(ModelSolutionSetDefault(mod, PhysicsSolution_Euler, phys)); in PhysicsCreate_Euler()
696 …PetscCall(ModelFunctionalRegister(mod, "Speed", &eu->monitor.Speed, PhysicsFunctional_Euler, phys)… in PhysicsCreate_Euler()
697 …scCall(ModelFunctionalRegister(mod, "Energy", &eu->monitor.Energy, PhysicsFunctional_Euler, phys)); in PhysicsCreate_Euler()
698 …Call(ModelFunctionalRegister(mod, "Density", &eu->monitor.Density, PhysicsFunctional_Euler, phys)); in PhysicsCreate_Euler()
699 …ll(ModelFunctionalRegister(mod, "Momentum", &eu->monitor.Momentum, PhysicsFunctional_Euler, phys)); in PhysicsCreate_Euler()
700 …ll(ModelFunctionalRegister(mod, "Pressure", &eu->monitor.Pressure, PhysicsFunctional_Euler, phys)); in PhysicsCreate_Euler()
1220 Physics phys = mod->physics; in Transfer() local
1224 …PetscCallMPI(MPIU_Allreduce(&phys->maxspeed, &mod->maxspeed, 1, MPIU_REAL, MPIU_MAX, PetscObjectCo… in Transfer()
1243 Physics phys; in main() local
1265 phys = mod->physics; in main()
1320 PetscCall(PetscMemzero(phys, sizeof(struct _n_Physics))); in main()
1321 PetscCall((*physcreate)(mod, phys, PetscOptionsObject)); in main()
1323 …for (phys->nfields = 0, phys->dof = 0; phys->field_desc[phys->nfields].name; phys->nfields++) phys… in main()
1324 …PetscCheck(phys->dof > 0, comm, PETSC_ERR_ARG_WRONGSTATE, "Physics '%s' did not set dof", physname… in main()
1402 PetscCall(PetscFVSetNumComponents(fvm, phys->dof)); in main()
1407 for (f = 0, dof = 0; f < phys->nfields; f++) { in main()
1408 PetscInt newDof = phys->field_desc[f].dof; in main()
1411 PetscCall(PetscFVSetComponentName(fvm, dof, phys->field_desc[f].name)); in main()
1418 …PetscCall(PetscSNPrintf(compName, sizeof(compName), "%s_%" PetscInt_FMT, phys->field_desc[f].name,… in main()
1431 PetscCall((*mod->setupbc)(dm, prob, phys)); in main()
1526 …PetscCallMPI(MPIU_Allreduce(&phys->maxspeed, &mod->maxspeed, 1, MPIU_REAL, MPIU_MAX, PetscObjectCo… in main()