33   for (PetscInt i = 0; i < nprofs; i++) {  in CalcCholeskyDecomp()  local
34     Cij[0][i] = sqrt(Rij[0][i]);  in CalcCholeskyDecomp()
35     Cij[3][i] = Rij[3][i] / Cij[0][i];  in CalcCholeskyDecomp()
36     Cij[1][i] = sqrt(Rij[1][i] - Square(Cij[3][i]));  in CalcCholeskyDecomp()
37     Cij[4][i] = Rij[4][i] / Cij[0][i];  in CalcCholeskyDecomp()
38     Cij[5][i] = (Rij[5][i] - Cij[3][i] * Cij[4][i]) / Cij[1][i];  in CalcCholeskyDecomp()
39     Cij[2][i] = sqrt(Rij[2][i] - Square(Cij[4][i]) - Square(Cij[5][i]));  in CalcCholeskyDecomp()
41     PetscCheck(!isnan(Cij[0][i]) && !isnan(Cij[1][i]) && !isnan(Cij[2][i]), comm, PETSC_ERR_FP,  in CalcCholeskyDecomp()
42 … at profile point %" PetscInt_FMT ". Either STGInflow has non-SPD matrix or contains nan.", i + 1);  in CalcCholeskyDecomp()
76   for (PetscInt i = 0; i < stg_ctx->nprofs; i++) {  in ReadStgInflow()  local
80 …does not have correct number of columns (%d instead of %" PetscInt_FMT ")", i, path, ndims, dims[1…  in ReadStgInflow()
82     wall_dist[i] = (CeedScalar)atof(array[0]);  in ReadStgInflow()
83     ubar[0][i]   = (CeedScalar)atof(array[1]);  in ReadStgInflow()
84     ubar[1][i]   = (CeedScalar)atof(array[2]);  in ReadStgInflow()
85     ubar[2][i]   = (CeedScalar)atof(array[3]);  in ReadStgInflow()
86     rij[0][i]    = (CeedScalar)atof(array[4]);  in ReadStgInflow()
87     rij[1][i]    = (CeedScalar)atof(array[5]);  in ReadStgInflow()
88     rij[2][i]    = (CeedScalar)atof(array[6]);  in ReadStgInflow()
89     rij[3][i]    = (CeedScalar)atof(array[7]);  in ReadStgInflow()
90     rij[4][i]    = (CeedScalar)atof(array[8]);  in ReadStgInflow()
91     rij[5][i]    = (CeedScalar)atof(array[9]);  in ReadStgInflow()
92     lt[i]        = (CeedScalar)atof(array[12]);  in ReadStgInflow()
93     eps[i]       = (CeedScalar)atof(array[13]);  in ReadStgInflow()
95 …PetscCheck(wall_dist[i] >= 0, comm, PETSC_ERR_FILE_UNEXPECTED, "Distance to wall in %s cannot be n…  in ReadStgInflow()
96 …PetscCheck(lt[i] >= 0, comm, PETSC_ERR_FILE_UNEXPECTED, "Turbulent length scale in %s cannot be ne…  in ReadStgInflow()
97 …PetscCheck(eps[i] >= 0, comm, PETSC_ERR_FILE_UNEXPECTED, "Turbulent dissipation in %s cannot be ne…  in ReadStgInflow()
131   for (PetscInt i = 0; i < stg_ctx->nmodes; i++) {  in ReadStgRand()  local
135 …does not have correct number of columns (%d instead of %" PetscInt_FMT ")", i, path, ndims, dims[1…  in ReadStgRand()
137     d[0][i]     = (CeedScalar)atof(array[0]);  in ReadStgRand()
138     d[1][i]     = (CeedScalar)atof(array[1]);  in ReadStgRand()
139     d[2][i]     = (CeedScalar)atof(array[2]);  in ReadStgRand()
140     phi[i]      = (CeedScalar)atof(array[3]);  in ReadStgRand()
141     sigma[0][i] = (CeedScalar)atof(array[4]);  in ReadStgRand()
142     sigma[1][i] = (CeedScalar)atof(array[5]);  in ReadStgRand()
143     sigma[2][i] = (CeedScalar)atof(array[6]);  in ReadStgRand()
208     CeedPragmaSIMD for (PetscInt i = 0; i < (*stg_ctx)->nprofs; i++) {  in GetStgContextData()  local
209       le = PetscMin(2 * wall_dist[i], 3 * lt[i]);  in GetStgContextData()
214 …CeedPragmaSIMD for (PetscInt i = 0; i < (*stg_ctx)->nmodes; i++) kappa[i] = kmin * pow((*stg_ctx)-…  in GetStgContextData()  local
281     for (PetscInt i = 0; i < 3; i++) domain_size[i] = domain_max[i] - domain_min[i];  in SetupStg()  local
342       for (int i = 0; i < 4; i++) comps[i] = i;  in SetupStrongStg()  local
347       for (int i = 0; i < 4; i++) comps[i] = i + 1;  in SetupStrongStg()  local
352       for (int i = 0; i < 4; i++) comps[i] = i + 1;  in SetupStrongStg()  local