ENH: Integration updates

This commit is contained in:
Andrew Heather
2016-04-26 14:34:27 +01:00
parent 9a0bd5831a
commit 90ee8d433f
4 changed files with 19 additions and 68 deletions

View File

@ -2,11 +2,7 @@
========= | ========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | \\ / O peration |
<<<<<<< HEAD:applications/utilities/postProcessing/graphics/PV4Readers/PV4FoamReader/vtkPV4Foam/vtkPV4FoamVolFields.H
\\ / A nd | Copyright (C) 2011-2014 OpenFOAM Foundation
=======
\\ / A nd | Copyright (C) 2011-2016 OpenFOAM Foundation \\ / A nd | Copyright (C) 2011-2016 OpenFOAM Foundation
>>>>>>> foundation:applications/utilities/postProcessing/graphics/PVReaders/PVFoamReader/vtkPVFoam/vtkPVFoamVolFields.H
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
@ -48,18 +44,10 @@ InClass
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
template<class Type> template<class Type>
<<<<<<< HEAD:applications/utilities/postProcessing/graphics/PV4Readers/PV4FoamReader/vtkPV4Foam/vtkPV4FoamVolFields.H void Foam::vtkPVFoam::convertVolField
void Foam::vtkPV4Foam::convertVolField
( (
const PtrList<PrimitivePatchInterpolation<primitivePatch>>& ppInterpList, const PtrList<PrimitivePatchInterpolation<primitivePatch>>& ppInterpList,
const GeometricField<Type, fvPatchField, volMesh>& tf, const GeometricField<Type, fvPatchField, volMesh>& tf,
=======
void Foam::vtkPVFoam::convertVolFields
(
const fvMesh& mesh,
const PtrList<PrimitivePatchInterpolation<primitivePatch>>& ppInterpList,
const IOobjectList& objects,
>>>>>>> foundation:applications/utilities/postProcessing/graphics/PVReaders/PVFoamReader/vtkPVFoam/vtkPVFoamVolFields.H
const bool interpFields, const bool interpFields,
vtkMultiBlockDataSet* output vtkMultiBlockDataSet* output
) )
@ -81,27 +69,7 @@ void Foam::vtkPVFoam::convertVolFields
( (
volPointInterpolation::New(mesh).interpolate(tf).ptr() volPointInterpolation::New(mesh).interpolate(tf).ptr()
); );
<<<<<<< HEAD:applications/utilities/postProcessing/graphics/PV4Readers/PV4FoamReader/vtkPV4Foam/vtkPV4FoamVolFields.H
} }
=======
// Interpolated field (demand driven)
autoPtr<GeometricField<Type, pointPatchField, pointMesh>> ptfPtr;
if (interpFields)
{
if (debug)
{
Info<< "convertVolFieldBlock interpolating:" << tf.name()
<< endl;
}
ptfPtr.reset
(
volPointInterpolation::New(tf.mesh()).interpolate(tf).ptr()
);
}
>>>>>>> foundation:applications/utilities/postProcessing/graphics/PVReaders/PVFoamReader/vtkPVFoam/vtkPVFoamVolFields.H
// Convert activated internalMesh regions // Convert activated internalMesh regions
convertVolFieldBlock convertVolFieldBlock
@ -169,7 +137,6 @@ void Foam::vtkPVFoam::convertVolFields
tmp<Field<Type>> tpptf tmp<Field<Type>> tpptf
( (
<<<<<<< HEAD:applications/utilities/postProcessing/graphics/PV4Readers/PV4FoamReader/vtkPV4Foam/vtkPV4FoamVolFields.H
fvPatchField<Type>(p, tf).patchInternalField() fvPatchField<Type>(p, tf).patchInternalField()
); );
@ -181,22 +148,6 @@ void Foam::vtkPVFoam::convertVolFields
arrayRangePatches_, arrayRangePatches_,
datasetNo datasetNo
); );
=======
isType<emptyFvPatchField<Type>>(ptf)
||
(
reader_->GetExtrapolatePatches()
&& !polyPatch::constraintType(patches[patchId].type())
)
)
{
fvPatch p(ptf.patch().patch(), tf.mesh().boundary());
tmp<Field<Type>> tpptf
(
fvPatchField<Type>(p, tf).patchInternalField()
);
>>>>>>> foundation:applications/utilities/postProcessing/graphics/PVReaders/PVFoamReader/vtkPVFoam/vtkPVFoamVolFields.H
if (interpFields) if (interpFields)
{ {
@ -310,7 +261,7 @@ void Foam::vtkPVFoam::convertVolFields
template<class Type> template<class Type>
void Foam::vtkPV4Foam::convertVolFields void Foam::vtkPVFoam::convertVolFields
( (
const fvMesh& mesh, const fvMesh& mesh,
const PtrList<PrimitivePatchInterpolation<primitivePatch>>& ppInterpList, const PtrList<PrimitivePatchInterpolation<primitivePatch>>& ppInterpList,
@ -345,7 +296,7 @@ void Foam::vtkPV4Foam::convertVolFields
template<class Type> template<class Type>
void Foam::vtkPV4Foam::convertDimFields void Foam::vtkPVFoam::convertDimFields
( (
const fvMesh& mesh, const fvMesh& mesh,
const PtrList<PrimitivePatchInterpolation<primitivePatch>>& ppInterpList, const PtrList<PrimitivePatchInterpolation<primitivePatch>>& ppInterpList,

View File

@ -305,7 +305,7 @@ void Foam::meshRefinement::weightedSum
} }
sum.setSize(meshPoints.size()); sum.setSize(meshPoints.size());
sum = Type(zero); sum = Type(Zero);
forAll(edges, edgeI) forAll(edges, edgeI)
{ {

View File

@ -284,7 +284,8 @@ void Foam::cyclicPeriodicAMIPolyPatch::resetAMI
ownStr.reset(new OBJstream(dir/name() + postfix)); ownStr.reset(new OBJstream(dir/name() + postfix));
neiStr.reset(new OBJstream(dir/neighbPatch().name() + postfix)); neiStr.reset(new OBJstream(dir/neighbPatch().name() + postfix));
Info<< "cyclicPeriodicAMIPolyPatch::resetAMI : patch:" << name() InfoInFunction
<< "patch:" << name()
<< " writing accumulated AMI to " << ownStr().name() << " writing accumulated AMI to " << ownStr().name()
<< " and " << neiStr().name() << endl; << " and " << neiStr().name() << endl;
} }
@ -368,7 +369,8 @@ void Foam::cyclicPeriodicAMIPolyPatch::resetAMI
if (debug) if (debug)
{ {
Info<< "cyclicPeriodicAMIPolyPatch::resetAMI : patch:" << name() InfoInFunction
<< "patch:" << name()
<< " srcSum:" << srcSum << " srcSum:" << srcSum
<< " tgtSum:" << tgtSum << " tgtSum:" << tgtSum
<< " direction:" << direction << " direction:" << direction
@ -391,8 +393,8 @@ void Foam::cyclicPeriodicAMIPolyPatch::resetAMI
if (debug) if (debug)
{ {
Info<< "cyclicPeriodicAMIPolyPatch::resetAMI : patch:" InfoInFunction
<< name() << "patch:" << name()
<< " moving this side from:" << " moving this side from:"
<< gAverage(thisPatch.points()) << gAverage(thisPatch.points())
<< " to:" << gAverage(thisPoints) << endl; << " to:" << gAverage(thisPoints) << endl;
@ -402,8 +404,8 @@ void Foam::cyclicPeriodicAMIPolyPatch::resetAMI
if (debug) if (debug)
{ {
Info<< "cyclicPeriodicAMIPolyPatch::resetAMI : patch:" InfoInFunction
<< name() << "patch:" << name()
<< " appending weights with untransformed slave side" << " appending weights with untransformed slave side"
<< endl; << endl;
} }
@ -421,8 +423,8 @@ void Foam::cyclicPeriodicAMIPolyPatch::resetAMI
if (debug) if (debug)
{ {
Info<< "cyclicPeriodicAMIPolyPatch::resetAMI : patch:" InfoInFunction
<< name() << "patch:" << name()
<< " moving neighbour side from:" << " moving neighbour side from:"
<< gAverage(nbrPatch.points()) << gAverage(nbrPatch.points())
<< " to:" << gAverage(nbrPoints) << endl; << " to:" << gAverage(nbrPoints) << endl;
@ -457,7 +459,8 @@ void Foam::cyclicPeriodicAMIPolyPatch::resetAMI
if (debug) if (debug)
{ {
Info<< "cyclicPeriodicAMIPolyPatch::resetAMI : patch:" << name() InfoInFunction
<< "patch:" << name()
<< " iteration:" << iter << " iteration:" << iter
<< " srcSum:" << srcSum << " srcSum:" << srcSum
<< " tgtSum:" << tgtSum << " tgtSum:" << tgtSum

View File

@ -95,10 +95,7 @@ void kinematicSingleLayer::resetPrimaryRegionSourceTerms()
void kinematicSingleLayer::transferPrimaryRegionThermoFields() void kinematicSingleLayer::transferPrimaryRegionThermoFields()
{ {
if (debug) DebugInFunction << endl;
{
InfoInFunction << endl;
}
// Update fields from primary region via direct mapped // Update fields from primary region via direct mapped
// (coupled) boundary conditions // (coupled) boundary conditions
@ -928,7 +925,7 @@ void kinematicSingleLayer::postEvolveRegion()
{ {
if (debug) if (debug)
{ {
Info<< "kinematicSingleLayer::postEvolveRegion()" << endl; InfoInFunction << endl;
} }
// Reset source terms for next time integration // Reset source terms for next time integration