59 …ReadStgInflow(const MPI_Comm comm, const char path[PETSC_MAX_PATH_LEN], StgShur14Context stg_ctx) {  in ReadStgInflow()  argument
70   CeedScalar  rij[6][stg_ctx->nprofs];  in ReadStgInflow()
71   CeedScalar *wall_dist              = &stg_ctx->data[stg_ctx->offsets.wall_dist];  in ReadStgInflow()
72   CeedScalar *eps                    = &stg_ctx->data[stg_ctx->offsets.eps];  in ReadStgInflow()
73   CeedScalar *lt                     = &stg_ctx->data[stg_ctx->offsets.lt];  in ReadStgInflow()
74 …CeedScalar(*ubar)[stg_ctx->nprofs] = (CeedScalar(*)[stg_ctx->nprofs]) & stg_ctx->data[stg_ctx->off…  in ReadStgInflow()
76   for (PetscInt i = 0; i < stg_ctx->nprofs; i++) {  in ReadStgInflow()
100 …CeedScalar(*cij)[stg_ctx->nprofs] = (CeedScalar(*)[stg_ctx->nprofs]) & stg_ctx->data[stg_ctx->offs…  in ReadStgInflow()
101   PetscCall(CalcCholeskyDecomp(comm, stg_ctx->nprofs, rij, cij));  in ReadStgInflow()
116 …e ReadStgRand(const MPI_Comm comm, const char path[PETSC_MAX_PATH_LEN], StgShur14Context stg_ctx) {  in ReadStgRand()  argument
127   CeedScalar *phi                     = &stg_ctx->data[stg_ctx->offsets.phi];  in ReadStgRand()
128 …CeedScalar(*d)[stg_ctx->nmodes]     = (CeedScalar(*)[stg_ctx->nmodes]) & stg_ctx->data[stg_ctx->of…  in ReadStgRand()
129 …CeedScalar(*sigma)[stg_ctx->nmodes] = (CeedScalar(*)[stg_ctx->nmodes]) & stg_ctx->data[stg_ctx->of…  in ReadStgRand()
131   for (PetscInt i = 0; i < stg_ctx->nmodes; i++) {  in ReadStgRand()
163                                  StgShur14Context *stg_ctx) {  in GetStgContextData()  argument
168   const PetscBool need_rand = (!(*stg_ctx)->mean_only) || (*stg_ctx)->use_fluctuating_IC;  in GetStgContextData()
179     *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()
202     PetscCall(ReadStgRand(comm, stg_rand_path, *stg_ctx));  in GetStgContextData()
203     CeedScalar *kappa     = &(*stg_ctx)->data[(*stg_ctx)->offsets.kappa];  in GetStgContextData()
204     CeedScalar *wall_dist = &(*stg_ctx)->data[(*stg_ctx)->offsets.wall_dist];  in GetStgContextData()
205     CeedScalar *lt        = &(*stg_ctx)->data[(*stg_ctx)->offsets.lt];  in GetStgContextData()
208     CeedPragmaSIMD for (PetscInt i = 0; i < (*stg_ctx)->nprofs; i++) {  in GetStgContextData()
214 …CeedPragmaSIMD for (PetscInt i = 0; i < (*stg_ctx)->nmodes; i++) kappa[i] = kmin * pow((*stg_ctx)-…  in GetStgContextData()