AMI improvements

• Added new faceAreaWeightAMI2D AMIMethod:

  • performs intersection using a new 2D triangle class;

  • candidate face matches set using an AABBTree method (vs advancing front for faceAreaWeightAMI).

  • Use by setting the AMIMethod entry when specifying the AMI in the constant/polyMesh/boundary file, e.g.

    AMI
    {
        type            cyclicACMI;
        AMIMethod       faceAreaWeightAMI2D; // new method
        Cbb             0.1; // optional coefficient
        nFaces          1000;
        startFace       100000;
        matchTolerance  0.0001;
        transform       noOrdering;
        neighbourPatch  AMI1;
        nonOverlapPatch AMI1_non_overlap;
    }

  • The optional Cbb coeffcient controls the size of the bounding box used when looking for candidate pairs; the value of 0.1 is the default and worked well for a large range of test cases. For badly matched AMI patches this may need to be increased.

• Deprecated the partialFaceAreaWeightAMI class - primarily used by ACMI:

  • functionality now offered by the AMI variants.

• NOTE: both AABBTree and octree variants included for review; ONLY 1 METHOD TO BE ADDED

EDIT: octree method remoed

src/OpenFOAM/primitives/Vector2D/Vector2D.H

@@ -6,7 +6,7 @@
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
     Copyright (C) 2011-2016 OpenFOAM Foundation
-    Copyright (C) 2018-2020 OpenCFD Ltd.
+    Copyright (C) 2018-2021 OpenCFD Ltd.
 -------------------------------------------------------------------------------
 License
     This file is part of OpenFOAM.
@@ -116,13 +116,18 @@ public:
         //- Access to the vector y component
         inline Cmpt& y();
 
-
         //- Normalise the vector by its magnitude
         inline Vector2D<Cmpt>& normalise();
 
-
         //- Perp dot product (dot product with perpendicular vector)
         inline scalar perp(const Vector2D<Cmpt>& b) const;
+
+        //- Return true if vector is within tol
+        inline bool isClose
+        (
+            const Vector2D<Cmpt>& b,
+            const scalar tol = 1e-10
+        ) const;
 };