「セルcellやパッチpatchを構成する面faceそれぞれの構成点pointを取得する方法」の編集履歴(バックアップ)一覧はこちら
追加された行は緑色になります。
削除された行は赤色になります。
*セルcellやパッチpatchを構成する面faceそれぞれの構成点pointを取得する方法
**セルcellを構成する面faceそれぞれの構成点pointを取得する方法
forAll(U, celli)
{
labelList poLabel = mesh.cellPoints(celli);
forAll(poLabel, ind) {
point pos = thePoints[poLabel[ind]];
Info << "cell[" << celli << "]=(" << pos.x() << "," << pos.y() << "," << pos.z() << ")" << endl;
}//forAll
}//forAll
**パッチpatchを構成する面faceそれぞれの構成点pointを取得する方法
const faceList& theFaces = mesh.faces();
const pointField& thePoints = mesh.points();
forAll(U.boundaryField(), patchi)
{
fvPatchVectorField& pU = U.boundaryField()[patchi];
const fvPatch& bm = mesh.boundary()[patchi];
label indexS = bm.patch().start();
forAll(pU, facei)
{
const pointField facepo = (theFaces[facei+indexS]).points(thePoints);
forAll(facepo, fapindex) {
point pos = facepo[fapindex];
Info << "patch[" << patchi << "]-"face[" << facei << "]=(" << pos.x() << "," << pos.y() << "," << pos.z() << ")" << endl;
}//forAll
}//forAll
}//forAll
<2014/09/20sat>
[EOF]&counter()