void Foam::sampledPatch::print(Ostream& os) const { os << "sampledPatch: " << name() << " :" << " patch:" << patchName() << " faces:" << faces().size() << " points:" << points().size(); }
Ostream& processorBoundaryPatch::operator<<(Ostream& os) const { os << patchName() << nl << token::BEGIN_BLOCK << nl << " type " << patchType() << token::END_STATEMENT << nl << " nFaces " << patchSize() << token::END_STATEMENT << nl << " startFace " << patchStart() << token::END_STATEMENT << nl << " myProcNo " << myProcNo_ << token::END_STATEMENT << nl << " neighbProcNo " << neighbProcNo_ << token::END_STATEMENT << nl << token::END_BLOCK << endl; return os; }
int main(int argc, char *argv[]) { # include <OpenFOAM/addRegionOption.H> timeSelector::addOptions(); argList::validArgs.append("fieldName"); argList::validArgs.append("patchName"); # include <OpenFOAM/setRootCase.H> # include <OpenFOAM/createTime.H> instantList timeDirs = timeSelector::select0(runTime, args); # include <OpenFOAM/createNamedMesh.H> word fieldName(args.additionalArgs()[0]); word patchName(args.additionalArgs()[1]); forAll(timeDirs, timeI) { runTime.setTime(timeDirs[timeI], timeI); Info<< "Time = " << runTime.timeName() << endl; IOobject fieldHeader ( fieldName, runTime.timeName(), mesh, IOobject::MUST_READ ); // Check field exists if (fieldHeader.headerOk()) { mesh.readUpdate(); label patchi = mesh.boundaryMesh().findPatchID(patchName); if (patchi < 0) { FatalError << "Unable to find patch " << patchName << nl << exit(FatalError); } // Give patch area if (isA<cyclicPolyPatch>(mesh.boundaryMesh()[patchi])) { Info<< " Cyclic patch vector area: " << nl; label nFaces = mesh.boundaryMesh()[patchi].size(); vector sum1 = vector::zero; vector sum2 = vector::zero; for (label i=0; i<nFaces/2; i++) { sum1 += mesh.Sf().boundaryField()[patchi][i]; sum2 += mesh.Sf().boundaryField()[patchi][i+nFaces/2]; } reduce(sum1, sumOp<vector>()); reduce(sum2, sumOp<vector>()); Info<< " - half 1 = " << sum1 << ", " << mag(sum1) << nl << " - half 2 = " << sum2 << ", " << mag(sum2) << nl << " - total = " << (sum1 + sum2) << ", " << mag(sum1 + sum2) << endl; Info<< " Cyclic patch area magnitude = " << gSum(mesh.magSf().boundaryField()[patchi])/2.0 << endl; } else { Info<< " Area vector of patch " << patchName << '[' << patchi << ']' << " = " << gSum(mesh.Sf().boundaryField()[patchi]) << endl; Info<< " Area magnitude of patch " << patchName << '[' << patchi << ']' << " = " << gSum(mesh.magSf().boundaryField()[patchi]) << endl; } // Read field and calc integral if (fieldHeader.headerClassName() == volScalarField::typeName) { Info<< " Reading " << volScalarField::typeName << " " << fieldName << endl; volScalarField field(fieldHeader, mesh); Info<< " Integral of " << fieldName << " over vector area of patch " << patchName << '[' << patchi << ']' << " = " << gSum ( mesh.Sf().boundaryField()[patchi] *field.boundaryField()[patchi] ) << nl; Info<< " Integral of " << fieldName << " over area magnitude of patch " << patchName << '[' << patchi << ']' << " = " << gSum ( mesh.magSf().boundaryField()[patchi] *field.boundaryField()[patchi] ) << nl; } else if ( fieldHeader.headerClassName() == surfaceScalarField::typeName ) { Info<< " Reading " << surfaceScalarField::typeName << " " << fieldName << endl; surfaceScalarField field(fieldHeader, mesh); scalar sumField = gSum(field.boundaryField()[patchi]); Info<< " Integral of " << fieldName << " over patch " << patchName << '[' << patchi << ']' << " = " << sumField << nl; } else { FatalError << "Only possible to integrate " << volScalarField::typeName << "s " << "and " << surfaceScalarField::typeName << "s" << nl << exit(FatalError); } } else { Info<< " No field " << fieldName << endl; } Info<< endl; }