From c9aec64fb9d9c20e4bf96b88c479383cf8f4ec5f Mon Sep 17 00:00:00 2001
From: Vaggelis Papoutsis <vaggelisp@gmail.com>
Date: Fri, 13 Dec 2019 20:22:23 +0200
Subject: [PATCH] BUG: NURBS3DVolume was getting a wrong set of points

findPointsInBox was receiving points in the Cartesian system,
not local coordinate system
---
 .../NURBS/NURBS3DVolume/NURBS3DVolume/NURBS3DVolume.C        | 5 ++++-
 1 file changed, 4 insertions(+), 1 deletion(-)

diff --git a/src/optimisation/adjointOptimisation/adjoint/parameterization/NURBS/NURBS3DVolume/NURBS3DVolume/NURBS3DVolume.C b/src/optimisation/adjointOptimisation/adjoint/parameterization/NURBS/NURBS3DVolume/NURBS3DVolume/NURBS3DVolume.C
index 5b03aa37fc..49d770cfef 100644
--- a/src/optimisation/adjointOptimisation/adjoint/parameterization/NURBS/NURBS3DVolume/NURBS3DVolume/NURBS3DVolume.C
+++ b/src/optimisation/adjointOptimisation/adjoint/parameterization/NURBS/NURBS3DVolume/NURBS3DVolume/NURBS3DVolume.C
@@ -1625,7 +1625,10 @@ Foam::tmp<Foam::vectorField> Foam::NURBS3DVolume::getPointsInBox()
     {
         findPointsInBox(localSystemCoordinates_);
     }
-    tmp<vectorField> pointsInBox(new vectorField(mesh_.points(), mapPtr_()));
+    tmp<vectorField> pointsInBox
+    (
+        new vectorField(localSystemCoordinates_, mapPtr_())
+    );
 
     return pointsInBox;
 }