63 …ReadStgInflow(const MPI_Comm comm, const char path[PETSC_MAX_PATH_LEN], StgShur14Context stg_ctx) {  in ReadStgInflow()  argument
74   CeedScalar  rij[6][stg_ctx->nprofs];  in ReadStgInflow()
75   CeedScalar *wall_dist              = &stg_ctx->data[stg_ctx->offsets.wall_dist];  in ReadStgInflow()
76   CeedScalar *eps                    = &stg_ctx->data[stg_ctx->offsets.eps];  in ReadStgInflow()
77   CeedScalar *lt                     = &stg_ctx->data[stg_ctx->offsets.lt];  in ReadStgInflow()
78 …CeedScalar(*ubar)[stg_ctx->nprofs] = (CeedScalar(*)[stg_ctx->nprofs]) & stg_ctx->data[stg_ctx->off…  in ReadStgInflow()
80   for (PetscInt i = 0; i < stg_ctx->nprofs; i++) {  in ReadStgInflow()
104 …CeedScalar(*cij)[stg_ctx->nprofs] = (CeedScalar(*)[stg_ctx->nprofs]) & stg_ctx->data[stg_ctx->offs…  in ReadStgInflow()
105   PetscCall(CalcCholeskyDecomp(comm, stg_ctx->nprofs, rij, cij));  in ReadStgInflow()
120 …e ReadStgRand(const MPI_Comm comm, const char path[PETSC_MAX_PATH_LEN], StgShur14Context stg_ctx) {  in ReadStgRand()  argument
131   CeedScalar *phi                     = &stg_ctx->data[stg_ctx->offsets.phi];  in ReadStgRand()
132 …CeedScalar(*d)[stg_ctx->nmodes]     = (CeedScalar(*)[stg_ctx->nmodes]) & stg_ctx->data[stg_ctx->of…  in ReadStgRand()
133 …CeedScalar(*sigma)[stg_ctx->nmodes] = (CeedScalar(*)[stg_ctx->nmodes]) & stg_ctx->data[stg_ctx->of…  in ReadStgRand()
135   for (PetscInt i = 0; i < stg_ctx->nmodes; i++) {  in ReadStgRand()
167                                  StgShur14Context *stg_ctx) {  in GetStgContextData()  argument
180     *temp_ctx                   = **stg_ctx;  in GetStgContextData()
194     PetscCall(PetscFree(*stg_ctx));  in GetStgContextData()
195     PetscCall(PetscMalloc(temp_ctx->total_bytes, stg_ctx));  in GetStgContextData()
196     **stg_ctx = *temp_ctx;  in GetStgContextData()
200   PetscCall(ReadStgInflow(comm, stg_inflow_path, *stg_ctx));  in GetStgContextData()
201   PetscCall(ReadStgRand(comm, stg_rand_path, *stg_ctx));  in GetStgContextData()
204     CeedScalar *kappa     = &(*stg_ctx)->data[(*stg_ctx)->offsets.kappa];  in GetStgContextData()
205     CeedScalar *wall_dist = &(*stg_ctx)->data[(*stg_ctx)->offsets.wall_dist];  in GetStgContextData()
206     CeedScalar *lt        = &(*stg_ctx)->data[(*stg_ctx)->offsets.lt];  in GetStgContextData()
209     CeedPragmaSIMD for (PetscInt i = 0; i < (*stg_ctx)->nprofs; i++) {  in GetStgContextData()
215 …CeedPragmaSIMD for (PetscInt i = 0; i < (*stg_ctx)->nmodes; i++) { kappa[i] = kmin * pow((*stg_ctx…  in GetStgContextData()