45 void Foam::projectEdge::findNearest
49 pointConstraint& constraint
57 List<pointConstraint> boundaryConstraint(1);
68 near = boundaryNear[0];
69 constraint = boundaryConstraint[0];
74 constraint = pointConstraint();
99 if (surfaces_[i] == -1)
102 <<
"Cannot find surface " << names[i] <<
" in geometry"
116 points_[start_] +
lambda*(points_[end_] - points_[start_])
124 findNearest(start, near, constraint);
135 static label eIter = 0;
144 Info<<
"Writing lines from straight-line start points"
145 <<
" to projected points to " << debugStr().name() <<
endl;
152 const point& startPt = points_[start_];
153 const point& endPt = points_[end_];
154 const vector d = endPt-startPt;
159 points[i] = startPt+lambdas[i]*d;
164 const label maxIter = 10;
167 const scalar relTol = 0.1;
168 const scalar absTol = 1
e-4;
170 scalar initialResidual = 0;
172 for (
label iter = 0; iter < maxIter; iter++)
189 if (lambdas[0] < small)
193 if (lambdas.
last() > 1 - small)
198 if (debugStr.
valid())
216 projLambdas /= projLambdas.
last();
227 interpolator.
valueWeights(lambdas[i], indices, weights);
232 predicted += weights[indexi]*
points[indices[indexi]];
234 residual[i] = predicted -
points[i];
237 const scalar scalarResidual =
sum(
mag(residual));
241 Pout<<
"Iter:" << iter <<
" initialResidual:" << initialResidual
242 <<
" residual:" << scalarResidual <<
endl;
245 if (scalarResidual < absTol*0.5*lambdas.
size())
251 initialResidual = scalarResidual;
253 else if (scalarResidual/initialResidual < relTol)
258 if (debugStr.
valid())
#define forAll(list, i)
Loop across all elements in list.
Macros for easy insertion into run-time selection tables.
An Istream is an abstract base class for all input systems (streams, files, token lists etc)....
void size(const label)
Override size to be inconsistent with allocated storage.
void setSize(const label)
Reset size of List.
OFstream which keeps track of vertices.
T & last()
Return the last element of the list.
An auto-pointer similar to the STL auto_ptr but with automatic casting to a reference to the type and...
bool valid() const
Return true if the autoPtr valid (ie, the pointer is set)
void reset(T *=nullptr)
If object pointer already set, delete object and set to given.
Define a curved edge that is parameterised for 0<lambda<1 between the start and end point.
const pointField & points_
A list of keyword definitions, which are a keyword followed by any number of values (e....
virtual bool valueWeights(const scalar t, labelList &indices, scalarField &weights) const
Calculate weights and indices to calculate t from samples.
Accumulates point constraints through successive applications of the applyConstraint function.
Defines the edge from the projection onto a surface (single surface) or intersection of two surfaces.
projectEdge(const dictionary &dict, const label index, const searchableSurfaces &geometry, const pointField &points, Istream &)
Construct from Istream setting pointsList.
virtual point position(const scalar) const
Return the point positions corresponding to the curve parameters.
static void findNearest(const PtrList< searchableSurface > &, const labelList &surfacesToTest, const pointField &, const scalarField &nearestDistSqr, labelList &surfaces, List< pointIndexHit > &)
Find nearest. Return -1 (and a miss()) or surface and nearest.
Container for searchableSurfaces.
label findSurfaceID(const word &name) const
Find index of surface. Return -1 if not found.
A class for managing temporary objects.
T & ref() const
Return non-const reference or generate a fatal error.
#define FatalIOErrorInFunction(ios)
Report an error message using Foam::FatalIOError.
dimensionedScalar lambda(viscosity->lookup("lambda"))
errorManipArg< error, int > exit(error &err, const int errNo=1)
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
Ostream & endl(Ostream &os)
Add newline and flush stream.
line< point, const point & > linePointRef
Line using referred points.
dimensioned< Type > sum(const DimensionedField< Type, GeoMesh > &df)
vectorField pointField
pointField is a vectorField.
vector point
Point is a vector.
Field< scalar > scalarField
Specialisation of Field<T> for scalar.
dimensioned< scalar > mag(const dimensioned< Type > &)
defineTypeNameAndDebug(combustionModel, 0)
prefixOSstream Pout(cout, "Pout")
word name(const complex &)
Return a string representation of a complex.
dimensioned< scalar > magSqr(const dimensioned< Type > &)
addToRunTimeSelectionTable(ensightPart, ensightPartCells, istream)
Unit conversion functions.