diff --git a/applications/test/syncTools/Test-syncTools.cxx b/applications/test/syncTools/Test-syncTools.cxx
index 05bcbf46e2..72e1eca1e7 100644
--- a/applications/test/syncTools/Test-syncTools.cxx
+++ b/applications/test/syncTools/Test-syncTools.cxx
@@ -525,6 +525,7 @@ public:
         }
     }
 };
+
 class edgePointTransformOp
 {
 public:
@@ -543,25 +544,24 @@ public:
             points1[i] = fld[i].second();
         }
 
-        pointField points0New;
-        pointField points1New;
         if (forward)
         {
-            points0New = vt.transformPosition(points0);
-            points1New = vt.transformPosition(points1);
+            vt.transformPositionList(points0);
+            vt.transformPositionList(points1);
         }
         else
         {
-            points0New = vt.invTransformPosition(points0);
-            points1New = vt.invTransformPosition(points1);
+            vt.invTransformPositionList(points0);
+            vt.invTransformPositionList(points1);
         }
 
         forAll(fld, i)
         {
-            fld[i] = PointPair(points0New[i], points1New[i]);
+            fld[i] = PointPair(points0[i], points1[i]);
         }
     }
 };
+
 class edgePointFlipOp
 {
 public:
@@ -572,6 +572,7 @@ public:
         return newVal;
     }
 };
+
 void testEdgeFlip2(const polyMesh& mesh, Random& rndGen)
 {
     Info<< nl << "Testing edge-wise (oriented) data synchronisation." << endl;
@@ -759,13 +760,11 @@ class pointListOps
 {
 public:
 
-    void operator()(pointList& lhs, const pointList& rhs) const
+    void operator()(UList<point>& lhs, const UList<point>& rhs) const
     {
         forAll(lhs, i)
         {
-            point& l = lhs[i];
-            const point& r = rhs[i];
-            maxMagSqrEqOp<vector>()(l, r);
+            maxMagSqrEqOp<vector>()(lhs[i], rhs[i]);
         }
     }
 
@@ -773,38 +772,31 @@ public:
     (
         const vectorTensorTransform& vt,
         const bool forward,
-        List<pointList>& fld
+        UList<pointList>& fld
     ) const
     {
         if (forward)
         {
-            for (auto& elems : fld)
+            for (auto& pts : fld)
             {
-                for (auto& elem : elems)
-                {
-                    elem = vt.transformPosition(elem);
-                }
+                vt.transformPositionList(pts);
             }
         }
         else
         {
-            for (auto& elems : fld)
+            for (auto& pts : fld)
             {
-                for (auto& elem : elems)
-                {
-                    elem = vt.invTransformPosition(elem);
-                }
+                vt.invTransformPositionList(pts);
             }
         }
     }
 
     //- Transform patch-based field
-    void operator()(const coupledPolyPatch& cpp, List<pointList>& fld) const
+    void operator()(const coupledPolyPatch& cpp, UList<pointList>& fld) const
     {
         forAll(fld, facei)
         {
-            pointList& pts = fld[facei];
-            for (auto& pt : pts)
+            for (auto& pt : fld[facei])
             {
                 cpp.transformPosition(pt, facei);
             }
diff --git a/applications/utilities/mesh/manipulation/checkMesh/checkGeometry.C b/applications/utilities/mesh/manipulation/checkMesh/checkGeometry.C
index b8008a71b4..4619f3b446 100644
--- a/applications/utilities/mesh/manipulation/checkMesh/checkGeometry.C
+++ b/applications/utilities/mesh/manipulation/checkMesh/checkGeometry.C
@@ -39,8 +39,7 @@ Foam::label Foam::findOppositeWedge
          && isA<wedgePolyPatch>(patches[patchi])
         )
         {
-            const wedgePolyPatch& pp =
-                refCast<const wedgePolyPatch>(patches[patchi]);
+            const auto& pp = refCast<const wedgePolyPatch>(patches[patchi]);
 
             // Calculate (cos of) angle to wpp (not pp!) centre normal
             scalar ppCosAngle = wpp.centreNormal() & pp.n();
@@ -80,8 +79,7 @@ bool Foam::checkWedges
     {
         if (patches[patchi].size() && isA<wedgePolyPatch>(patches[patchi]))
         {
-            const wedgePolyPatch& pp =
-                refCast<const wedgePolyPatch>(patches[patchi]);
+            const auto& pp = refCast<const wedgePolyPatch>(patches[patchi]);
 
             scalar wedgeAngle = acos(pp.cosAngle());
 
@@ -105,7 +103,7 @@ bool Foam::checkWedges
                 return true;
             }
 
-            const wedgePolyPatch& opp =
+            const auto& opp =
                 refCast<const wedgePolyPatch>(patches[oppositePatchi]);
 
 
@@ -275,25 +273,27 @@ namespace Foam
         void operator()
         (
             const coupledPolyPatch& cpp,
-            List<pointField>& pts
+            UList<pointField>& pts
         ) const
         {
-            // Each element of pts is all the points in the face. Convert into
-            // lists of size cpp to transform.
+            // Each element of pts is all the points in the face.
+            // Convert into lists of size cpp to transform.
 
             List<pointField> newPts(pts.size());
             forAll(pts, facei)
             {
-                newPts[facei].setSize(pts[facei].size());
+                newPts[facei].resize(pts[facei].size());
             }
 
+            pointField ptsAtIndex(pts.size());
+
             label index = 0;
             while (true)
             {
                 label n = 0;
 
                 // Extract for every face the i'th position
-                pointField ptsAtIndex(pts.size(), Zero);
+                ptsAtIndex = Zero;
                 forAll(cpp, facei)
                 {
                     const pointField& facePts = pts[facei];
@@ -326,7 +326,8 @@ namespace Foam
                 index++;
             }
 
-            pts.transfer(newPts);
+            // Transfer newPts -> pts
+            std::move(newPts.begin(), newPts.end(), pts.begin());
         }
     };
 }
@@ -352,10 +353,7 @@ bool Foam::checkCoupledPoints
     {
         if (patches[patchi].coupled())
         {
-            const coupledPolyPatch& cpp = refCast<const coupledPolyPatch>
-            (
-                patches[patchi]
-            );
+            const auto& cpp = refCast<const coupledPolyPatch>(patches[patchi]);
 
             forAll(cpp, i)
             {
@@ -387,8 +385,7 @@ bool Foam::checkCoupledPoints
     {
         if (patches[patchi].coupled())
         {
-            const coupledPolyPatch& cpp =
-                refCast<const coupledPolyPatch>(patches[patchi]);
+            const auto& cpp = refCast<const coupledPolyPatch>(patches[patchi]);
 
             if (cpp.owner())
             {
@@ -1023,7 +1020,7 @@ Foam::label Foam::checkGeometry
         {
             if (isA<cyclicAMIPolyPatch>(pbm[patchi]))
             {
-                const cyclicAMIPolyPatch& cpp =
+                const auto& cpp =
                     refCast<const cyclicAMIPolyPatch>(pbm[patchi]);
 
                 if (cpp.owner())
@@ -1058,7 +1055,7 @@ Foam::label Foam::checkGeometry
 
                         if (isA<cyclicACMIPolyPatch>(pbm[patchi]))
                         {
-                            const cyclicACMIPolyPatch& pp =
+                            const auto& pp =
                                 refCast<const cyclicACMIPolyPatch>(pbm[patchi]);
 
                             scalarField mergedMask;
@@ -1108,8 +1105,9 @@ Foam::label Foam::checkGeometry
 
                         if (isA<cyclicACMIPolyPatch>(pbm[patchi]))
                         {
-                            const cyclicACMIPolyPatch& pp =
+                            const auto& pp =
                                 refCast<const cyclicACMIPolyPatch>(pbm[patchi]);
+
                             scalarField mergedMask;
                             globalFaces().gather
                             (
diff --git a/src/OpenFOAM/meshes/polyMesh/mapPolyMesh/mapDistribute/mapDistribute.C b/src/OpenFOAM/meshes/polyMesh/mapPolyMesh/mapDistribute/mapDistribute.C
index 1a8cc8ea39..f0dd047945 100644
--- a/src/OpenFOAM/meshes/polyMesh/mapPolyMesh/mapDistribute/mapDistribute.C
+++ b/src/OpenFOAM/meshes/polyMesh/mapPolyMesh/mapDistribute/mapDistribute.C
@@ -40,96 +40,6 @@ namespace Foam
 
 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
 
-template<>
-void Foam::mapDistribute::transform::operator()
-(
-    const vectorTensorTransform&,
-    const bool,
-    List<label>&
-) const {}
-
-template<>
-void Foam::mapDistribute::transform::operator()
-(
-    const coupledPolyPatch&,
-    UList<label>&
-) const {}
-
-template<>
-void Foam::mapDistribute::transform::operator()
-(
-    const coupledPolyPatch&,
-    Map<label>&
-) const {}
-
-template<>
-void Foam::mapDistribute::transform::operator()
-(
-    const coupledPolyPatch&,
-    EdgeMap<label>&
-) const {}
-
-
-template<>
-void Foam::mapDistribute::transform::operator()
-(
-    const vectorTensorTransform&,
-    const bool,
-    List<scalar>&
-) const {}
-
-template<>
-void Foam::mapDistribute::transform::operator()
-(
-    const coupledPolyPatch&,
-    UList<scalar>&
-) const {}
-
-template<>
-void Foam::mapDistribute::transform::operator()
-(
-    const coupledPolyPatch&,
-    Map<scalar>&
-) const {}
-
-template<>
-void Foam::mapDistribute::transform::operator()
-(
-    const coupledPolyPatch&,
-    EdgeMap<scalar>&
-) const {}
-
-
-template<>
-void Foam::mapDistribute::transform::operator()
-(
-    const vectorTensorTransform&,
-    const bool,
-    List<bool>&
-) const {}
-
-template<>
-void Foam::mapDistribute::transform::operator()
-(
-    const coupledPolyPatch&,
-    UList<bool>&
-) const {}
-
-template<>
-void Foam::mapDistribute::transform::operator()
-(
-    const coupledPolyPatch&,
-    Map<bool>&
-) const {}
-
-template<>
-void Foam::mapDistribute::transform::operator()
-(
-    const coupledPolyPatch&,
-    EdgeMap<bool>&
-) const {}
-
-
 void Foam::mapDistribute::printLayout(Ostream& os) const
 {
     mapDistributeBase::printLayout(os);
diff --git a/src/OpenFOAM/meshes/polyMesh/mapPolyMesh/mapDistribute/mapDistribute.H b/src/OpenFOAM/meshes/polyMesh/mapPolyMesh/mapDistribute/mapDistribute.H
index 92b98611c9..95c92a999a 100644
--- a/src/OpenFOAM/meshes/polyMesh/mapPolyMesh/mapDistribute/mapDistribute.H
+++ b/src/OpenFOAM/meshes/polyMesh/mapPolyMesh/mapDistribute/mapDistribute.H
@@ -219,11 +219,14 @@ public:
             (
                 const vectorTensorTransform& vt,
                 const bool forward,
-                List<Type>& fld
+                UList<Type>& fld
             ) const
             {
-                const tensor T(forward ? vt.R() : vt.R().T());
-                transformList(T, fld);
+                if constexpr (is_rotational_vectorspace_v<Type>)
+                {
+                    const tensor rot(forward ? vt.R() : vt.R().T());
+                    transformList(rot, fld);
+                }
             }
 
             template<class Type>
@@ -234,9 +237,13 @@ public:
                 List<List<Type>>& flds
             ) const
             {
-                for (List<Type>& fld : flds)
+                if constexpr (is_rotational_vectorspace_v<Type>)
                 {
-                    operator()(vt, forward, fld);
+                    const tensor rot(forward ? vt.R() : vt.R().T());
+                    for (auto& fld : flds)
+                    {
+                        transformList(rot, fld);
+                    }
                 }
             }
 
@@ -244,9 +251,12 @@ public:
             template<class Type>
             void operator()(const coupledPolyPatch& cpp, UList<Type>& fld) const
             {
-                if (!cpp.parallel())
+                if constexpr (is_rotational_vectorspace_v<Type>)
                 {
-                    transformList(cpp.forwardT(), fld);
+                    if (!cpp.parallel())
+                    {
+                        transformList(cpp.forwardT(), fld);
+                    }
                 }
             }
 
@@ -255,9 +265,12 @@ public:
             void operator()(const coupledPolyPatch& cpp, Container<Type>& map)
             const
             {
-                if (!cpp.parallel())
+                if constexpr (is_rotational_vectorspace_v<Type>)
                 {
-                    transformList(cpp.forwardT(), map);
+                    if (!cpp.parallel())
+                    {
+                        transformList(cpp.forwardT(), map);
+                    }
                 }
             }
         };
@@ -271,17 +284,16 @@ public:
             (
                 const vectorTensorTransform& vt,
                 const bool forward,
-                List<point>& fld
+                UList<point>& fld
             ) const
             {
-                pointField pfld(std::move(fld));
                 if (forward)
                 {
-                    fld = vt.transformPosition(pfld);
+                    vt.transformPositionList(fld);
                 }
                 else
                 {
-                    fld = vt.invTransformPosition(pfld);
+                    vt.invTransformPositionList(fld);
                 }
             }
 
@@ -292,9 +304,19 @@ public:
                 List<List<point>>& flds
             ) const
             {
-                for (List<point>& fld : flds)
+                if (forward)
                 {
-                    operator()(vt, forward, fld);
+                    for (auto& fld : flds)
+                    {
+                        vt.transformPositionList(fld);
+                    }
+                }
+                else
+                {
+                    for (auto& fld : flds)
+                    {
+                        vt.invTransformPositionList(fld);
+                    }
                 }
             }
 
@@ -682,86 +704,6 @@ public:
 };
 
 
-// Template specialisation for primitives that do not need transform
-template<>
-void mapDistribute::transform::operator()
-(
-    const vectorTensorTransform&,
-    const bool,
-    List<label>&
-) const;
-template<>
-void mapDistribute::transform::operator()
-(
-    const coupledPolyPatch&,
-    UList<label>&
-) const;
-template<>
-void mapDistribute::transform::operator()
-(
-    const coupledPolyPatch&,
-    Map<label>&
-) const;
-template<>
-void mapDistribute::transform::operator()
-(
-    const coupledPolyPatch&,
-    EdgeMap<label>&
-) const;
-
-template<>
-void mapDistribute::transform::operator()
-(
-    const coupledPolyPatch&,
-    UList<scalar>&
-) const;
-template<>
-void mapDistribute::transform::operator()
-(
-    const vectorTensorTransform&,
-    const bool,
-    List<scalar>&
-) const;
-template<>
-void mapDistribute::transform::operator()
-(
-    const coupledPolyPatch&,
-    Map<scalar>&
-) const;
-template<>
-void mapDistribute::transform::operator()
-(
-    const coupledPolyPatch&,
-    EdgeMap<scalar>&
-) const;
-
-template<>
-void mapDistribute::transform::operator()
-(
-    const coupledPolyPatch& cpp,
-    UList<bool>& fld
-) const;
-template<>
-void mapDistribute::transform::operator()
-(
-    const vectorTensorTransform&,
-    const bool,
-    List<bool>&
-) const;
-template<>
-void mapDistribute::transform::operator()
-(
-    const coupledPolyPatch&,
-    Map<bool>&
-) const;
-template<>
-void mapDistribute::transform::operator()
-(
-    const coupledPolyPatch&,
-    EdgeMap<bool>&
-) const;
-
-
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
 } // End namespace Foam
diff --git a/src/OpenFOAM/meshes/polyMesh/syncTools/dummyTransform.H b/src/OpenFOAM/meshes/polyMesh/syncTools/dummyTransform.H
index 355c5f9e9c..2dd2362744 100644
--- a/src/OpenFOAM/meshes/polyMesh/syncTools/dummyTransform.H
+++ b/src/OpenFOAM/meshes/polyMesh/syncTools/dummyTransform.H
@@ -6,7 +6,7 @@
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
     Copyright (C) 2011-2016 OpenFOAM Foundation
-    Copyright (C) 2020 OpenCFD Ltd.
+    Copyright (C) 2020-2025 OpenCFD Ltd.
 -------------------------------------------------------------------------------
 License
     This file is part of OpenFOAM.
@@ -47,17 +47,27 @@ namespace Foam
 class dummyTransform
 {
 public:
+
     template<class T>
     void operator()
     (
         const vectorTensorTransform& vt,
         const bool forward,
-        List<T>& fld
+        UList<T>& fld
+    ) const
+    {}
+
+    template<class Type>
+    void operator()
+    (
+        const vectorTensorTransform& vt,
+        const bool forward,
+        List<List<Type>>& flds
     ) const
     {}
 
     template<class T>
-    void operator()(const coupledPolyPatch& cpp, Field<T>& fld) const
+    void operator()(const coupledPolyPatch& cpp, UList<T>& fld) const
     {}
 
     template<class T, template<class> class Container>
diff --git a/src/OpenFOAM/primitives/globalIndexAndTransform/vectorTensorTransform/vectorTensorTransform.C b/src/OpenFOAM/primitives/globalIndexAndTransform/vectorTensorTransform/vectorTensorTransform.C
index 8aa2ffc20d..a79d9e0e01 100644
--- a/src/OpenFOAM/primitives/globalIndexAndTransform/vectorTensorTransform/vectorTensorTransform.C
+++ b/src/OpenFOAM/primitives/globalIndexAndTransform/vectorTensorTransform/vectorTensorTransform.C
@@ -6,6 +6,7 @@
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
     Copyright (C) 2011-2016 OpenFOAM Foundation
+    Copyright (C) 2025 OpenCFD Ltd.
 -------------------------------------------------------------------------------
 License
     This file is part of OpenFOAM.
@@ -60,6 +61,13 @@ Foam::vectorTensorTransform::vectorTensorTransform(Istream& is)
 
 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 
+void Foam::vectorTensorTransform::checkRotation(const scalar tol)
+{
+    // Detect identity and zero rotations
+    hasR_ = !(R_.is_identity(tol) || Foam::mag(R_) < tol);
+}
+
+
 Foam::word Foam::name(const vectorTensorTransform& s)
 {
     OStringStream buf;
@@ -70,32 +78,6 @@ Foam::word Foam::name(const vectorTensorTransform& s)
 }
 
 
-template<>
-Foam::tmp<Foam::Field<bool>> Foam::vectorTensorTransform::transform
-(
-    const Field<bool>& fld
-) const
-{
-    return fld;
-}
-template<>
-Foam::tmp<Foam::Field<Foam::label>> Foam::vectorTensorTransform::transform
-(
-    const Field<label>& fld
-) const
-{
-    return fld;
-}
-template<>
-Foam::tmp<Foam::Field<Foam::scalar>> Foam::vectorTensorTransform::transform
-(
-    const Field<scalar>& fld
-) const
-{
-    return fld;
-}
-
-
 // * * * * * * * * * * * * * * * IOstream Operators  * * * * * * * * * * * * //
 
 Foam::Istream& Foam::operator>>(Istream& is, vectorTensorTransform& tr)
diff --git a/src/OpenFOAM/primitives/globalIndexAndTransform/vectorTensorTransform/vectorTensorTransform.H b/src/OpenFOAM/primitives/globalIndexAndTransform/vectorTensorTransform/vectorTensorTransform.H
index 3a7aa26a72..258b29b4af 100644
--- a/src/OpenFOAM/primitives/globalIndexAndTransform/vectorTensorTransform/vectorTensorTransform.H
+++ b/src/OpenFOAM/primitives/globalIndexAndTransform/vectorTensorTransform/vectorTensorTransform.H
@@ -6,7 +6,7 @@
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
     Copyright (C) 2011-2016 OpenFOAM Foundation
-    Copyright (C) 2020 OpenCFD Ltd.
+    Copyright (C) 2020-2025 OpenCFD Ltd.
 -------------------------------------------------------------------------------
 License
     This file is part of OpenFOAM.
@@ -34,12 +34,12 @@ Description
 SourceFiles
     vectorTensorTransformI.H
     vectorTensorTransform.C
-    vectorTensorTransformTemplates.C
+    vectorTensorTransform.txx
 
 \*---------------------------------------------------------------------------*/
 
-#ifndef vectorTensorTransform_H
-#define vectorTensorTransform_H
+#ifndef Foam_vectorTensorTransform_H
+#define Foam_vectorTensorTransform_H
 
 #include "tensor.H"
 #include "word.H"
@@ -126,39 +126,57 @@ public:
 
     // Member Functions
 
-        // Access
+    // Access
 
-            inline const vector& t() const;
+        //- True if it has a non-identity rotation tensor
+        bool hasR() const noexcept { return hasR_; }
 
-            inline const tensor& R() const;
+        //- The translation vector
+        const vector& t() const noexcept { return t_; }
 
-            inline bool hasR() const;
+        //- The (forward) rotation tensor
+        const tensor& R() const noexcept { return R_; }
 
 
-        // Edit
+    // Edit
 
-            inline vector& t();
+        //- Non-const access to the translation vector
+        vector& t() noexcept { return t_; }
 
-            inline tensor& R();
+        //- Non-const access to the rotation tensor. Sets hasR = true
+        inline tensor& R() noexcept;
+
+        //- Test for identity rotation and set hasR accordingly
+        void checkRotation(const scalar tol = ROOTVSMALL);
 
 
-        // Transform
+    // Transform
 
-            //- Transform the given position
-            inline vector transformPosition(const vector& v) const;
+        //- Transform the given position
+        inline vector transformPosition(const vector& v) const;
 
-            //- Transform the given pointField
-            inline pointField transformPosition(const pointField& pts) const;
+        //- Transform the given field of points
+        inline pointField transformPosition(const pointField& pts) const;
 
-            //- Inverse transform the given position
-            inline vector invTransformPosition(const vector& v) const;
+        //- Inplace transform the given points
+        inline void transformPositionList(UList<point>& pts) const;
 
-            //- Inverse transform the given pointField
-            inline pointField invTransformPosition(const pointField& pts) const;
+        //- Inverse transform the given position
+        inline vector invTransformPosition(const vector& v) const;
 
-            //- Transform the given field
-            template<class Type>
-            tmp<Field<Type>> transform(const Field<Type>&) const;
+        //- Inverse transform the given field of points
+        inline pointField invTransformPosition(const pointField& pts) const;
+
+        //- Inplace inverse transform the given points
+        inline void invTransformPositionList(UList<point>& pts) const;
+
+        //- Transform the given field
+        template<class Type>
+        tmp<Field<Type>> transform(const Field<Type>& fld) const;
+
+        //- Inplace transform the given field
+        template<class Type>
+        void transformList(UList<Type>& fld) const;
 
 
     // Member Operators
@@ -196,18 +214,6 @@ inline vectorTensorTransform inv(const vectorTensorTransform& tr);
 //- Return a string representation of a vectorTensorTransform
 word name(const vectorTensorTransform&);
 
-//- Template specialisations
-template<>
-tmp<Field<bool>> vectorTensorTransform::transform(const Field<bool>&) const;
-
-template<>
-tmp<Field<label>> vectorTensorTransform::transform(const Field<label>&) const;
-
-template<>
-tmp<Field<scalar>> vectorTensorTransform::transform(const Field<scalar>&)
-const;
-
-
 // * * * * * * * * * * * * * * * Global Operators  * * * * * * * * * * * * * //
 
 inline bool operator==
@@ -262,7 +268,7 @@ inline vectorTensorTransform operator&
 #include "vectorTensorTransformI.H"
 
 #ifdef NoRepository
-    #include "vectorTensorTransformTemplates.C"
+    #include "vectorTensorTransform.txx"
 #endif
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
diff --git a/src/OpenFOAM/primitives/globalIndexAndTransform/vectorTensorTransform/vectorTensorTransformTemplates.C b/src/OpenFOAM/primitives/globalIndexAndTransform/vectorTensorTransform/vectorTensorTransform.txx
similarity index 75%
rename from src/OpenFOAM/primitives/globalIndexAndTransform/vectorTensorTransform/vectorTensorTransformTemplates.C
rename to src/OpenFOAM/primitives/globalIndexAndTransform/vectorTensorTransform/vectorTensorTransform.txx
index aeeab30c08..0f72f8c146 100644
--- a/src/OpenFOAM/primitives/globalIndexAndTransform/vectorTensorTransform/vectorTensorTransformTemplates.C
+++ b/src/OpenFOAM/primitives/globalIndexAndTransform/vectorTensorTransform/vectorTensorTransform.txx
@@ -6,6 +6,7 @@
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
     Copyright (C) 2011-2016 OpenFOAM Foundation
+    Copyright (C) 2025 OpenCFD Ltd.
 -------------------------------------------------------------------------------
 License
     This file is part of OpenFOAM.
@@ -33,9 +34,16 @@ Foam::tmp<Foam::Field<Type>> Foam::vectorTensorTransform::transform
     const Field<Type>& fld
 ) const
 {
-    if (hasR_)
+    if constexpr (is_rotational_vectorspace_v<Type>)
     {
-        return R() & fld;
+        if (hasR_)
+        {
+            return (R() & fld);
+        }
+        else
+        {
+            return fld;
+        }
     }
     else
     {
@@ -44,4 +52,23 @@ Foam::tmp<Foam::Field<Type>> Foam::vectorTensorTransform::transform
 }
 
 
+template<class Type>
+void Foam::vectorTensorTransform::transformList
+(
+    UList<Type>& fld
+) const
+{
+    if constexpr (is_rotational_vectorspace_v<Type>)
+    {
+        if (hasR_)
+        {
+            for (auto& val : fld)
+            {
+                val = (R() & val);
+            }
+        }
+    }
+}
+
+
 // ************************************************************************* //
diff --git a/src/OpenFOAM/primitives/globalIndexAndTransform/vectorTensorTransform/vectorTensorTransformI.H b/src/OpenFOAM/primitives/globalIndexAndTransform/vectorTensorTransform/vectorTensorTransformI.H
index c1c7cf58f7..a2747bc59c 100644
--- a/src/OpenFOAM/primitives/globalIndexAndTransform/vectorTensorTransform/vectorTensorTransformI.H
+++ b/src/OpenFOAM/primitives/globalIndexAndTransform/vectorTensorTransform/vectorTensorTransformI.H
@@ -6,6 +6,7 @@
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
     Copyright (C) 2011-2016 OpenFOAM Foundation
+    Copyright (C) 2025 OpenCFD Ltd.
 -------------------------------------------------------------------------------
 License
     This file is part of OpenFOAM.
@@ -66,35 +67,10 @@ inline Foam::vectorTensorTransform::vectorTensorTransform(const tensor& R)
 
 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 
-inline const Foam::vector& Foam::vectorTensorTransform::t() const
-{
-    return t_;
-}
-
-
-inline const Foam::tensor& Foam::vectorTensorTransform::R() const
-{
-    return R_;
-}
-
-
-inline bool Foam::vectorTensorTransform::hasR() const
-{
-    return hasR_;
-}
-
-
-inline Foam::vector& Foam::vectorTensorTransform::t()
-{
-    return t_;
-}
-
-
-inline Foam::tensor& Foam::vectorTensorTransform::R()
+inline Foam::tensor& Foam::vectorTensorTransform::R() noexcept
 {
     // Assume that non-const access to R changes it from I, so set
     // hasR to true
-
     hasR_ = true;
 
     return R_;
@@ -136,6 +112,28 @@ inline Foam::pointField Foam::vectorTensorTransform::transformPosition
 }
 
 
+inline void Foam::vectorTensorTransform::transformPositionList
+(
+    UList<point>& pts
+) const
+{
+    if (hasR_)
+    {
+        for (auto& p : pts)
+        {
+            p = t() + (R() & p);
+        }
+    }
+    else
+    {
+        for (auto& p : pts)
+        {
+            p += t();
+        }
+    }
+}
+
+
 inline Foam::vector Foam::vectorTensorTransform::invTransformPosition
 (
     const vector& v
@@ -171,6 +169,30 @@ inline Foam::pointField Foam::vectorTensorTransform::invTransformPosition
 }
 
 
+inline void Foam::vectorTensorTransform::invTransformPositionList
+(
+    UList<point>& pts
+) const
+{
+    if (hasR_)
+    {
+        const auto rot = R().T();
+
+        for (auto& p : pts)
+        {
+            p = (rot & (p - t()));
+        }
+    }
+    else
+    {
+        for (auto& p : pts)
+        {
+            p -= t();
+        }
+    }
+}
+
+
 // * * * * * * * * * * * * * * * Member Operators  * * * * * * * * * * * * * //
 
 inline void Foam::vectorTensorTransform::operator&=
diff --git a/src/mesh/snappyHexMesh/meshRefinement/weightedPosition.C b/src/mesh/snappyHexMesh/meshRefinement/weightedPosition.C
index 5db17e07f9..909365c7b3 100644
--- a/src/mesh/snappyHexMesh/meshRefinement/weightedPosition.C
+++ b/src/mesh/snappyHexMesh/meshRefinement/weightedPosition.C
@@ -62,7 +62,7 @@ void Foam::weightedPosition::getPoints
     List<point>& out
 )
 {
-    out.setSize(in.size());
+    out.resize_nocopy(in.size());
     forAll(in, i)
     {
         out[i] = in[i].second();
@@ -78,11 +78,24 @@ void Foam::weightedPosition::getPoints
 void Foam::weightedPosition::setPoints
 (
     const UList<point>& in,
-    List<weightedPosition>& out
+    UList<weightedPosition>& out
 )
 {
-    out.setSize(in.size());
-    forAll(in, i)
+    #ifdef FULLDEBUG
+    if (in.size() != out.size())
+    {
+        FatalErrorInFunction
+            << "Mismatch in size. Number of points:" << in.size()
+            << " != number of weighted positions:" << out.size()
+            << exit(FatalError);
+    }
+
+    const label len = out.size();
+    #else
+    const label len = Foam::min(in.size(), out.size());
+    #endif
+
+    for (label i = 0; i < len; ++i)
     {
         out[i].second() = out[i].first()*in[i];
     }
@@ -104,7 +117,7 @@ void Foam::weightedPosition::operator()
 (
     const vectorTensorTransform& vt,
     const bool forward,
-    List<weightedPosition>& fld
+    UList<weightedPosition>& fld
 ) const
 {
     pointField pfld;
@@ -112,11 +125,11 @@ void Foam::weightedPosition::operator()
 
     if (forward)
     {
-        pfld = vt.transformPosition(pfld);
+        vt.transformPositionList(pfld);
     }
     else
     {
-        pfld = vt.invTransformPosition(pfld);
+        vt.invTransformPositionList(pfld);
     }
 
     setPoints(pfld, fld);
@@ -130,7 +143,7 @@ void Foam::weightedPosition::operator()
     List<List<weightedPosition>>& flds
 ) const
 {
-    for (List<weightedPosition>& fld : flds)
+    for (auto& fld : flds)
     {
         operator()(vt, forward, fld);
     }
diff --git a/src/mesh/snappyHexMesh/meshRefinement/weightedPosition.H b/src/mesh/snappyHexMesh/meshRefinement/weightedPosition.H
index a56fe80d25..1ca5eaf06e 100644
--- a/src/mesh/snappyHexMesh/meshRefinement/weightedPosition.H
+++ b/src/mesh/snappyHexMesh/meshRefinement/weightedPosition.H
@@ -125,7 +125,7 @@ public:
             static void setPoints
             (
                 const UList<point>& in,
-                List<weightedPosition>& out
+                UList<weightedPosition>& out  // [in,out]
             );
 
 
@@ -133,13 +133,13 @@ public:
         static void plusEqOp(weightedPosition& x, const weightedPosition& y);
 
 
-        // Transformation operators
+        // Transformation Operators
 
             void operator()
             (
                 const vectorTensorTransform& vt,
                 const bool forward,
-                List<weightedPosition>& fld
+                UList<weightedPosition>& fld
             ) const;
 
             void operator()
diff --git a/src/meshTools/algorithms/PatchEdgeFaceWave/PatchEdgeFaceWave.H b/src/meshTools/algorithms/PatchEdgeFaceWave/PatchEdgeFaceWave.H
index 0a59fb7ac8..df5feab901 100644
--- a/src/meshTools/algorithms/PatchEdgeFaceWave/PatchEdgeFaceWave.H
+++ b/src/meshTools/algorithms/PatchEdgeFaceWave/PatchEdgeFaceWave.H
@@ -380,19 +380,11 @@ public:
         List<Type>& fld
     ) const
     {
-        if (forward)
+        const tensor rot(forward ? vt.R() : vt.R().T());
+
+        for (auto& val : fld)
         {
-            forAll(fld, i)
-            {
-                fld[i].transform(mesh_, patch_, vt.R(), tol_, td_);
-            }
-        }
-        else
-        {
-            forAll(fld, i)
-            {
-                fld[i].transform(mesh_, patch_, vt.R().T(), tol_, td_);
-            }
+            val.transform(mesh_, patch_, rot, tol_, td_);
         }
     }
 };