Go to the documentation of this file.
44 fieldTableName_(iF.name()),
50 startSampledValues_(0),
69 fieldTableName_(iF.name()),
70 setAverage_(
dict.getOrDefault(
"setAverage",
false)),
71 perturb_(
dict.getOrDefault(
"perturb", 1
e-5)),
83 startSampledValues_(0),
90 if (
dict.found(
"offset"))
97 mapMethod_ !=
"planarInterpolation"
98 && mapMethod_ !=
"nearest"
102 <<
"mapMethod should be one of 'planarInterpolation'"
106 dict.readIfPresent(
"fieldTableName", fieldTableName_);
108 if (
dict.found(
"value"))
137 fieldTableName_(ptf.fieldTableName_),
138 setAverage_(ptf.setAverage_),
139 perturb_(ptf.perturb_),
140 mapMethod_(ptf.mapMethod_),
143 startSampleTime_(-1),
144 startSampledValues_(0),
147 endSampledValues_(0),
149 offset_(ptf.offset_.clone())
161 fieldTableName_(ptf.fieldTableName_),
162 setAverage_(ptf.setAverage_),
163 perturb_(ptf.perturb_),
164 mapMethod_(ptf.mapMethod_),
165 mapperPtr_(ptf.mapperPtr_),
166 sampleTimes_(ptf.sampleTimes_),
167 startSampleTime_(ptf.startSampleTime_),
168 startSampledValues_(ptf.startSampledValues_),
169 startAverage_(ptf.startAverage_),
170 endSampleTime_(ptf.endSampleTime_),
171 endSampledValues_(ptf.endSampledValues_),
172 endAverage_(ptf.endAverage_),
173 offset_(ptf.offset_.clone())
186 fieldTableName_(ptf.fieldTableName_),
187 setAverage_(ptf.setAverage_),
188 perturb_(ptf.perturb_),
189 mapMethod_(ptf.mapMethod_),
190 mapperPtr_(ptf.mapperPtr_),
191 sampleTimes_(ptf.sampleTimes_),
192 startSampleTime_(ptf.startSampleTime_),
193 startSampledValues_(ptf.startSampledValues_),
194 startAverage_(ptf.startAverage_),
195 endSampleTime_(ptf.endSampleTime_),
196 endSampledValues_(ptf.endSampledValues_),
197 endAverage_(ptf.endAverage_),
198 offset_(ptf.offset_.clone())
211 if (startSampledValues_.size())
213 startSampledValues_.autoMap(m);
214 endSampledValues_.autoMap(m);
218 startSampleTime_ = -1;
233 refCast<const timeVaryingMappedFixedValuePointPatchField<Type>>(ptf);
235 startSampledValues_.
rmap(tiptf.startSampledValues_, addr);
236 endSampledValues_.rmap(tiptf.endSampledValues_, addr);
240 startSampleTime_ = -1;
248 const Time& time = this->db().
time();
251 if (startSampleTime_ == -1 && endSampleTime_ == -1)
253 const polyMesh& pMesh = this->
patch().boundaryMesh().mesh()();
277 polyMesh::meshSubDir,
292 /this->
patch().name()
313 && mapMethod_ !=
"planarInterpolation"
330 const fileName samplePointsDir = samplePointsFile.
path();
331 sampleTimes_ = Time::findTimes(samplePointsDir);
335 Info<<
"timeVaryingMappedFixedValuePointPatchField : In directory "
336 << samplePointsDir <<
" found times "
337 << pointToPointPlanarInterpolation::timeNames(sampleTimes_)
346 bool foundTime = mapperPtr_().findTime
358 <<
"Cannot find starting sampling values for current time "
360 <<
"Have sampling values for times "
361 << pointToPointPlanarInterpolation::timeNames(sampleTimes_) <<
nl
364 <<
"\n on patch " << this->
patch().name()
365 <<
" of field " << fieldTableName_
372 if (lo != startSampleTime_)
374 startSampleTime_ = lo;
376 if (startSampleTime_ == endSampleTime_)
381 Pout<<
"checkTable : Setting startValues to (already read) "
383 /this->
patch().name()
384 /sampleTimes_[startSampleTime_].name()
387 startSampledValues_ = endSampledValues_;
388 startAverage_ = endAverage_;
394 Pout<<
"checkTable : Reading startValues from "
396 /this->
patch().name()
397 /sampleTimes_[lo].name()
406 /this->
patch().name()
407 /sampleTimes_[startSampleTime_].name()
424 startAverage_ = vals.
average();
427 if (vals.size() != mapperPtr_().sourceSize())
430 <<
"Number of values (" << vals.size()
431 <<
") differs from the number of points ("
432 << mapperPtr_().sourceSize()
436 startSampledValues_ = mapperPtr_().interpolate(vals);
440 if (hi != endSampleTime_)
444 if (endSampleTime_ == -1)
449 Pout<<
"checkTable : Clearing endValues" <<
endl;
451 endSampledValues_.clear();
457 Pout<<
"checkTable : Reading endValues from "
459 /this->
patch().name()
460 /sampleTimes_[endSampleTime_].name()
469 /this->
patch().name()
470 /sampleTimes_[endSampleTime_].name()
491 if (vals.size() != mapperPtr_().sourceSize())
494 <<
"Number of values (" << vals.size()
495 <<
") differs from the number of points ("
496 << mapperPtr_().sourceSize()
500 endSampledValues_ = mapperPtr_().interpolate(vals);
520 if (endSampleTime_ == -1)
525 Pout<<
"updateCoeffs : Sampled, non-interpolated values"
526 <<
" from start time:"
527 << sampleTimes_[startSampleTime_].
name() <<
nl;
531 wantedAverage = startAverage_;
535 scalar start = sampleTimes_[startSampleTime_].value();
536 scalar
end = sampleTimes_[endSampleTime_].value();
538 scalar
s = (this->db().time().value()-start)/(
end-start);
542 Pout<<
"updateCoeffs : Sampled, interpolated values"
543 <<
" between start time:"
544 << sampleTimes_[startSampleTime_].
name()
545 <<
" and end time:" << sampleTimes_[endSampleTime_].
name()
546 <<
" with weight:" <<
s <<
endl;
549 this->
operator==((1-s)*startSampledValues_ +
s*endSampledValues_);
550 wantedAverage = (1-
s)*startAverage_ +
s*endAverage_;
563 Pout<<
"updateCoeffs :"
564 <<
" actual average:" << averagePsi
565 <<
" wanted average:" << wantedAverage
569 if (
mag(averagePsi) < VSMALL)
572 const Type offset = wantedAverage - averagePsi;
575 Pout<<
"updateCoeffs :"
576 <<
" offsetting with:" << offset <<
endl;
582 const scalar scale =
mag(wantedAverage)/
mag(averagePsi);
586 Pout<<
"updateCoeffs :"
587 <<
" scaling with:" << scale <<
endl;
596 const scalar t = this->db().time().timeOutputValue();
602 Pout<<
"updateCoeffs : set fixedValue to min:" <<
gMin(*
this)
603 <<
" max:" <<
gMax(*
this)
625 this->internalField().
name(),
632 "planarInterpolation",
638 offset_->writeData(os);
int debug
Static debugging option.
vectorField pointField
pointField is a vectorField.
virtual const pointField & points() const
Return raw points.
Ostream & writeEntryIfDifferent(const word &key, const T &value1, const T &value2)
Write a keyword/value entry only when the two values differ.
Defines the attributes of an object for which implicit objectRegistry management is supported,...
A simple wrapper around bool so that it can be read as a word: true/false, on/off,...
Class to control time during OpenFOAM simulations that is also the top-level objectRegistry.
A class for handling words, derived from Foam::string.
A class for handling file names.
A primitive field of type <T> with automated input and output.
gmvFile<< "tracers "<< particles.size()<< nl;for(const passiveParticle &p :particles){ gmvFile<< p.position().x()<< " ";}gmvFile<< nl;for(const passiveParticle &p :particles){ gmvFile<< p.position().y()<< " ";}gmvFile<< nl;for(const passiveParticle &p :particles){ gmvFile<< p.position().z()<< " ";}gmvFile<< nl;forAll(lagrangianScalarNames, i){ word name=lagrangianScalarNames[i];IOField< scalar > s(IOobject(name, runTime.timeName(), cloud::prefix, mesh, IOobject::MUST_READ, IOobject::NO_WRITE))
static std::string path(const std::string &str)
Return directory path name (part before last /)
static constexpr const zero Zero
Global zero (0)
Type gAverage(const FieldField< Field, Type > &f)
virtual void updateCoeffs()
Update the coefficients associated with the patch field.
static std::string name(const std::string &str)
Return basename (part beyond last /), including its extension.
const Time & time() const
Return time.
const fileName & facesInstance() const
Return the current instance directory for faces.
Interpolates between two sets of unstructured points using 2D Delaunay triangulation....
Ostream & endl(Ostream &os)
Add newline and flush stream.
const Type & value() const
Return const reference to value.
Basic pointPatch represents a set of points from the mesh.
prefixOSstream Pout
An Ostream wrapper for parallel output to std::cout.
Abstract base class for point-mesh patch fields.
Mesh consisting of general polyhedral cells.
Foam::pointPatchFieldMapper.
const fileName & pointsInstance() const
Return the current instance directory for points.
A FixedValue boundary condition for pointField.
Generic templated field type.
tmp< faMatrix< Type > > operator==(const faMatrix< Type > &, const faMatrix< Type > &)
Like IOField but falls back to raw IFstream if no header found. Optionally reads average value....
messageStream Info
Information stream (uses stdout - output is on the master only)
A time-varying form of a mapped fixed value boundary condition.
word name(const complex &c)
Return string representation of complex.
virtual const fileName & name() const
Return the name of the stream.
virtual void rmap(const pointPatchField< Type > &, const labelList &)
Reverse map the given PointPatchField onto.
virtual void write(Ostream &) const
Write.
gmvFile<< "tracers "<< particles.size()<< nl;for(const passiveParticle &p :particles){ gmvFile<< p.position().x()<< ' ';}gmvFile<< nl;for(const passiveParticle &p :particles){ gmvFile<< p.position().y()<< ' ';}gmvFile<< nl;for(const passiveParticle &p :particles){ gmvFile<< p.position().z()<< ' ';}gmvFile<< nl;for(const word &name :lagrangianScalarNames){ IOField< scalar > fld(IOobject(name, runTime.timeName(), cloud::prefix, mesh, IOobject::MUST_READ, IOobject::NO_WRITE))
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
constexpr auto end(C &c) -> decltype(c.end())
Return iterator to the end of the container c.
timeVaryingMappedFixedValuePointPatchField(const pointPatch &, const DimensionedField< Type, pointMesh > &)
Construct from patch and internal field.
errorManipArg< error, int > exit(error &err, const int errNo=1)
tmp< DimensionedField< TypeR, GeoMesh > > New(const tmp< DimensionedField< TypeR, GeoMesh >> &tdf1, const word &name, const dimensionSet &dimensions)
Global function forwards to reuseTmpDimensionedField::New.
pointField points0(pointIOField(IOobject("points", mesh.time().constant(), polyMesh::meshSubDir, mesh, IOobject::MUST_READ, IOobject::NO_WRITE, false)))
const Type & average() const
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
const std::string patch
OpenFOAM patch number as a std::string.
fileName caseConstant() const
virtual void autoMap(const pointPatchFieldMapper &)
Map (and resize as needed) from self given a mapping object.
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
const dimensionedScalar e
Elementary charge.
void write(vtk::formatter &fmt, const Type &val, const label n=1)
Component-wise write of a value (N times)
#define FatalIOErrorInFunction(ios)
Report an error message using Foam::FatalIOError.
An Ostream is an abstract base class for all output systems (streams, files, token lists,...
Type gMin(const FieldField< Field, Type > &f)
void checkTable()
Find boundary data inbetween current time and interpolate.
const word & constant() const
Return constant name.
Type gMax(const FieldField< Field, Type > &f)
string evaluate(const std::string &s, size_t pos=0, size_t len=std::string::npos)
Field with dimensions and associated with geometry type GeoMesh which is used to size the field and a...