void CPrimalGrid::SetCoord_CG(su2double **val_coord) { unsigned short iDim, iNode, NodeFace, iFace; AD::StartPreacc(); AD::SetPreaccIn(val_coord, GetnNodes(), nDim); for (iDim = 0; iDim < nDim; iDim++) { Coord_CG[iDim] = 0.0; for (iNode = 0; iNode < GetnNodes(); iNode++) Coord_CG[iDim] += val_coord[iNode][iDim]/su2double(GetnNodes()); } for (iFace = 0; iFace < GetnFaces(); iFace++) for (iDim = 0; iDim < nDim; iDim++) { Coord_FaceElems_CG[iFace][iDim] = 0.0; for (iNode = 0; iNode < GetnNodesFace(iFace); iNode++) { NodeFace = GetFaces(iFace, iNode); Coord_FaceElems_CG[iFace][iDim] += val_coord[NodeFace][iDim]/su2double(GetnNodesFace(iFace)); } } AD::SetPreaccOut(Coord_CG, nDim); AD::SetPreaccOut(Coord_FaceElems_CG, GetnFaces(), nDim); AD::EndPreacc(); }
void CPrimalGrid::SetCG(double **val_coord) { unsigned short iDim, iNode, NodeFace, iFace; for (iDim = 0; iDim < nDim; iDim++) { Coord_CG[iDim] = 0.0; for (iNode = 0; iNode < GetnNodes(); iNode++) Coord_CG[iDim] += val_coord[iNode][iDim]/double(GetnNodes()); } for (iFace = 0; iFace < GetnFaces(); iFace++) for (iDim = 0; iDim < nDim; iDim++) { Coord_FaceElems_CG[iFace][iDim] = 0.0; for (iNode = 0; iNode < GetnNodesFace(iFace); iNode++) { NodeFace = GetFaces(iFace, iNode); Coord_FaceElems_CG[iFace][iDim] += val_coord[NodeFace][iDim]/double(GetnNodesFace(iFace)); } } }