Lines Matching refs:dimR
237 PetscInt dim, dimR, offset; in DMFieldEvaluate_DS() local
258 PetscCall(DMGetDimension(field->dm, &dimR)); in DMFieldEvaluate_DS()
268 …alloc3(gatherSize * dim, &cellPoints, gatherMax * dim, &coordsReal, gatherMax * dimR, &coordsRef)); in DMFieldEvaluate_DS()
269 …PetscCall(PetscMalloc4(gatherMax * dimR, &v, gatherMax * dimR * dimR, &J, gatherMax * dimR * dimR,… in DMFieldEvaluate_DS()
297 PetscCall(PetscMalloc1(dimR * nq, &quadPoints)); in DMFieldEvaluate_DS()
298 for (p = 0; p < dimR * nq; p++) quadPoints[p] = coordsRef[p]; in DMFieldEvaluate_DS()
299 PetscCall(PetscQuadratureSetData(quad, dimR, 0, nq, quadPoints, NULL)); in DMFieldEvaluate_DS()
323 for (d = 0; d < dimR; d++) { in DMFieldEvaluate_DS()
325 …for (e = 0; e < dimR; e++) vs[d] += invJ[dimR * dimR * p + e * dimR + d] * cD[(nc * p + g) * dimR … in DMFieldEvaluate_DS()
327 for (d = 0; d < dimR; d++) cD[(nc * p + g) * dimR + d] = vs[d]; in DMFieldEvaluate_DS()
336 for (d = 0; d < dimR; d++) { in DMFieldEvaluate_DS()
338 …for (e = 0; e < dimR; e++) v[d] += invJ[dimR * dimR * p + e * dimR + d] * cD[(nc * p + g) * dimR +… in DMFieldEvaluate_DS()
340 for (d = 0; d < dimR; d++) cD[(nc * p + g) * dimR + d] = v[d]; in DMFieldEvaluate_DS()
351 for (g = 0; g < nc * dimR; g++) { in DMFieldEvaluate_DS()
354 for (d = 0; d < dimR; d++) { in DMFieldEvaluate_DS()
356 …for (e = 0; e < dimR; e++) vs[d] += invJ[dimR * dimR * p + e * dimR + d] * cH[(nc * dimR * p + g) … in DMFieldEvaluate_DS()
358 for (d = 0; d < dimR; d++) cH[(nc * dimR * p + g) * dimR + d] = vs[d]; in DMFieldEvaluate_DS()
361 for (f = 0; f < dimR; f++) { in DMFieldEvaluate_DS()
364 for (d = 0; d < dimR; d++) { in DMFieldEvaluate_DS()
366 …for (e = 0; e < dimR; e++) vs[d] += invJ[dimR * dimR * p + e * dimR + d] * cH[((nc * p + g) * dimR… in DMFieldEvaluate_DS()
368 for (d = 0; d < dimR; d++) cH[((nc * p + g) * dimR + d) * dimR + f] = vs[d]; in DMFieldEvaluate_DS()
377 for (g = 0; g < nc * dimR; g++) { in DMFieldEvaluate_DS()
378 for (d = 0; d < dimR; d++) { in DMFieldEvaluate_DS()
380 …for (e = 0; e < dimR; e++) v[d] += invJ[dimR * dimR * p + e * dimR + d] * cH[(nc * dimR * p + g) *… in DMFieldEvaluate_DS()
382 for (d = 0; d < dimR; d++) cH[(nc * dimR * p + g) * dimR + d] = v[d]; in DMFieldEvaluate_DS()
385 for (f = 0; f < dimR; f++) { in DMFieldEvaluate_DS()
386 for (d = 0; d < dimR; d++) { in DMFieldEvaluate_DS()
388 …for (e = 0; e < dimR; e++) v[d] += invJ[dimR * dimR * p + e * dimR + d] * cH[((nc * p + g) * dimR … in DMFieldEvaluate_DS()
390 for (d = 0; d < dimR; d++) cH[((nc * p + g) * dimR + d) * dimR + f] = v[d]; in DMFieldEvaluate_DS()