Lines Matching refs:d
111 PetscInt vStart, vEnd, v, d, coordSize; in ChangeCoordinates() local
134 for (d = 0; d < spaceDim; ++d) coords[off + d] = vertexCoords[(v - vStart) * spaceDim + d]; in ChangeCoordinates()
149 PetscInt d, i, j; in CheckFEMGeometry() local
153 for (d = 0; d < spaceDim; ++d) { in CheckFEMGeometry()
154 if (v0[d] != v0Ex[d]) { in CheckFEMGeometry()
179 PetscInt d; in CheckFVMGeometry() local
183 for (d = 0; d < spaceDim; ++d) { in CheckFVMGeometry()
185 …d], centroidEx[d]) < tol, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Cell %" PetscInt_FMT ", Invalid centro… in CheckFVMGeometry()
186 …d], normalEx[d]) < tol, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Cell %" PetscInt_FMT ", Invalid normal[%… in CheckFVMGeometry()
198 PetscInt coneSize, f, dim, cdim, d; in CheckGaussLaw() local
213 for (d = 0; d < cdim; ++d) integral[d] += sgn * area * normal[d]; in CheckGaussLaw()
215 for (d = 0; d < cdim; ++d) in CheckGaussLaw()
216 …d]) < tol, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Cell %" PetscInt_FMT " Surface integral for component… in CheckGaussLaw()
272 PetscInt d, e, f, p; in CheckCell() local
278 for (d = 0; d < cdim; ++d) PetscCall(PetscRandomGetValue(r, &trans[d])); in CheckCell()
306 for (d = 0; d < cdim; ++d) { in CheckCell()
307 v0ExT[d] = v0Ex[d]; in CheckCell()
309 JExT[d * cdim + e] = JEx[d * cdim + e]; in CheckCell()
310 invJExT[d * cdim + e] = invJEx[d * cdim + e]; in CheckCell()
313 for (d = 0; d < cdim; ++d) { in CheckCell()
314 v0ExT[d] *= scale; in CheckCell()
315 v0ExT[d] += PetscRealPart(trans[d]); in CheckCell()
318 JExT[d * cdim + e] *= scale; in CheckCell()
319 invJExT[d * cdim + e] /= scale; in CheckCell()
321 if (d < dim) detJExT *= scale; in CheckCell()
324 for (d = 0; d < cdim; ++d) { in CheckCell()
325 for (e = 0, rot[d] = 0.0; e < cdim; ++e) rot[d] += R[d * cdim + e] * v0ExT[e]; in CheckCell()
327 for (d = 0; d < cdim; ++d) v0ExT[d] = rot[d]; in CheckCell()
328 for (d = 0; d < cdim; ++d) { in CheckCell()
330 …for (f = 0, rotM[d * cdim + e] = 0.0; f < cdim; ++f) rotM[d * cdim + e] += R[d * cdim + f] * JExT[… in CheckCell()
333 for (d = 0; d < cdim; ++d) { in CheckCell()
334 for (e = 0; e < cdim; ++e) JExT[d * cdim + e] = rotM[d * cdim + e]; in CheckCell()
336 for (d = 0; d < cdim; ++d) { in CheckCell()
338 …for (f = 0, rotM[d * cdim + e] = 0.0; f < cdim; ++f) rotM[d * cdim + e] += invJExT[d * cdim + f] *… in CheckCell()
341 for (d = 0; d < cdim; ++d) { in CheckCell()
342 for (e = 0; e < cdim; ++e) invJExT[d * cdim + e] = rotM[d * cdim + e]; in CheckCell()
347 for (d = 0; d < cdim; ++d) { in CheckCell()
348 centroidExT[d] = centroidEx[d]; in CheckCell()
349 normalExT[d] = normalEx[d]; in CheckCell()
351 for (d = 0; d < cdim; ++d) { in CheckCell()
352 centroidExT[d] *= scale; in CheckCell()
353 centroidExT[d] += PetscRealPart(trans[d]); in CheckCell()
354 normalExT[d] /= scale; in CheckCell()
356 if (d < dim) volExT *= scale; in CheckCell()
359 for (d = 0; d < cdim; ++d) { in CheckCell()
360 for (e = 0, rot[d] = 0.0; e < cdim; ++e) rot[d] += R[d * cdim + e] * centroidExT[e]; in CheckCell()
362 for (d = 0; d < cdim; ++d) centroidExT[d] = rot[d]; in CheckCell()
363 for (d = 0; d < cdim; ++d) { in CheckCell()
364 for (e = 0, rot[d] = 0.0; e < cdim; ++e) rot[d] += R[d * cdim + e] * normalExT[e]; in CheckCell()
366 for (d = 0; d < cdim; ++d) normalExT[d] = rot[d]; in CheckCell()
367 for (d = 0, norm = 0.0; d < cdim; ++d) norm += PetscSqr(normalExT[d]); in CheckCell()
370 for (d = 0; d < cdim; ++d) normalExT[d] /= norm; in CheckCell()
372 for (d = 0; d < cdim; ++d) { in CheckCell()
374 newCoords[p * cdim + d] *= scale; in CheckCell()
375 newCoords[p * cdim + d] += trans[d]; in CheckCell()
379 for (d = 0; d < cdim; ++d) { in CheckCell()
380 …for (e = 0, rot[d] = 0.0; e < cdim; ++e) rot[d] += R[d * cdim + e] * PetscRealPart(newCoords[p * c… in CheckCell()
382 for (d = 0; d < cdim; ++d) newCoords[p * cdim + d] = rot[d]; in CheckCell()
397 for (d = 0; d < cdim; ++d) { in CheckCell()
398 faceCentroidExT[d] = faceCentroidEx[off + d]; in CheckCell()
399 faceNormalExT[d] = faceNormalEx[off + d]; in CheckCell()
401 for (d = 0; d < cdim; ++d) { in CheckCell()
402 faceCentroidExT[d] *= scale; in CheckCell()
403 faceCentroidExT[d] += PetscRealPart(trans[d]); in CheckCell()
404 faceNormalExT[d] /= scale; in CheckCell()
406 if (d < dim - 1) faceVolExT *= scale; in CheckCell()
408 for (d = 0; d < cdim; ++d) { in CheckCell()
409 … for (e = 0, rot[d] = 0.0; e < cdim; ++e) rot[d] += R[d * cdim + e] * faceCentroidExT[e]; in CheckCell()
411 for (d = 0; d < cdim; ++d) faceCentroidExT[d] = rot[d]; in CheckCell()
412 for (d = 0; d < cdim; ++d) { in CheckCell()
413 for (e = 0, rot[d] = 0.0; e < cdim; ++e) rot[d] += R[d * cdim + e] * faceNormalExT[e]; in CheckCell()
415 for (d = 0; d < cdim; ++d) faceNormalExT[d] = rot[d]; in CheckCell()
416 for (d = 0, norm = 0.0; d < cdim; ++d) norm += PetscSqr(faceNormalExT[d]); in CheckCell()
418 for (d = 0; d < cdim; ++d) faceNormalExT[d] /= norm; in CheckCell()
653 PetscInt dim, d, cStart, cEnd, cEndInterior, c; in main() local
672 for (d = 0; d < dim; ++d) { in main()
673 if (d > 0) PetscCall(PetscPrintf(PETSC_COMM_SELF, ", ")); in main()
674 PetscCall(PetscPrintf(PETSC_COMM_SELF, "%12.2g", (double)cg->centroid[d])); in main()