From 6897a85cc34cd0727761aaed3eb5f2cddae3c186 Mon Sep 17 00:00:00 2001
From: Henry Weller <http://cfd.direct>
Date: Wed, 21 Apr 2021 21:21:14 +0100
Subject: [PATCH] volPointInterpolation::pointConstraints: Resize
 patchPatchPointConstraints_ if necessary

---
 .../volPointInterpolation/pointConstraints.C     | 16 ++++++++++++----
 1 file changed, 12 insertions(+), 4 deletions(-)

diff --git a/src/finiteVolume/interpolation/volPointInterpolation/pointConstraints.C b/src/finiteVolume/interpolation/volPointInterpolation/pointConstraints.C
index 4e3e1e5204..508120c938 100644
--- a/src/finiteVolume/interpolation/volPointInterpolation/pointConstraints.C
+++ b/src/finiteVolume/interpolation/volPointInterpolation/pointConstraints.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     | Website:  https://openfoam.org
-    \\  /    A nd           | Copyright (C) 2013-2018 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2013-2021 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -90,8 +90,7 @@ void Foam::pointConstraints::makePatchPatchAddressing()
     Map<label> patchPatchPointSet(2*nPatchPatchPoints);
 
     // Constraints (initialised to unconstrained)
-    patchPatchPointConstraints_.setSize(nPatchPatchPoints);
-    patchPatchPointConstraints_ = pointConstraint();
+    patchPatchPointConstraints_.setSize(nPatchPatchPoints, pointConstraint());
 
     // From constraint index to mesh point
     labelList patchPatchPoints(nPatchPatchPoints);
@@ -251,6 +250,16 @@ void Foam::pointConstraints::makePatchPatchAddressing()
                     constraintI = iter();
                 }
 
+                // Extend the patchPatchPointConstraints_ array if necessary
+                if (patchPatchPointConstraints_.size() <= constraintI)
+                {
+                    patchPatchPointConstraints_.setSize
+                    (
+                        constraintI + 1,
+                        pointConstraint()
+                    );
+                }
+
                 // Combine (new or existing) constraint with one
                 // on coupled.
                 patchPatchPointConstraints_[constraintI].combine
@@ -262,7 +271,6 @@ void Foam::pointConstraints::makePatchPatchAddressing()
     }
 
 
-
     nPatchPatchPoints = pppi;
     patchPatchPoints.setSize(nPatchPatchPoints);
     patchPatchPointConstraints_.setSize(nPatchPatchPoints);