67 #define MEHess(out, cf, etaB, etaD, dim, nc, cast) \  argument
75           PetscScalar c = (cf)[_m] * (etaB)[_ind] * (etaD)[_ind]; \
83 …aluate(PetscInt dim, PetscReal (*coordRange)[2], PetscInt nc, PetscScalar *cf, PetscScalar *cfWork…  in MultilinearEvaluate()  argument
120         MEdot(out, cf, etaB, (1 << dim), nc, (PetscScalar));  in MultilinearEvaluate()
124         MEdot(out, cf, etaB, (1 << dim), nc, PetscRealPart);  in MultilinearEvaluate()
140         for (k = 0; k < whol * nc; k++) cfWork[k] = cf[k];  in MultilinearEvaluate()
149             cf[k * nc + l]          = cfWork[2 * k * nc + l];  in MultilinearEvaluate()
150             cf[(k + half) * nc + l] = cfWork[(2 * k + 1) * nc + l];  in MultilinearEvaluate()
155             for (l = 0; l < nc; l++) out[l * dim + j] += d * e * cf[k * nc + l];  in MultilinearEvaluate()
159             for (l = 0; l < nc; l++) out[l * dim + j] += d * e * PetscRealPart(cf[k * nc + l]);  in MultilinearEvaluate()
168         MEHess(out, cf, etaB, etaD, dim, nc, (PetscScalar));  in MultilinearEvaluate()
172         MEHess(out, cf, etaB, etaD, dim, nc, PetscRealPart);  in MultilinearEvaluate()
471   PetscScalar *cv, *cf, *work;  in DMFieldCreateDA()  local
479   PetscCall(PetscMalloc3(nv, &cv, nv, &cf, nv, &work));  in DMFieldCreateDA()
481   for (i = 0; i < nv; i++) cf[i] = cv[i];  in DMFieldCreateDA()
483   dafield->cornerCoeffs = cf;  in DMFieldCreateDA()
491 …  for (k = 0; k < nc; k++) w[j * nc + k] = 0.5 * (cf[(2 * j + 1) * nc + k] - cf[(2 * j) * nc + k]);  in DMFieldCreateDA()
495 …  for (k = 0; k < nc; k++) w[j * nc + k] = 0.5 * (cf[(2 * j) * nc + k] + cf[(2 * j + 1) * nc + k]);  in DMFieldCreateDA()
497     for (j = 0; j < nv; j++) cf[j] = work[j];  in DMFieldCreateDA()