From 6e818675629627ce1fc3e78fe77c483084a9e21f Mon Sep 17 00:00:00 2001
From: mattijs <mattijs>
Date: Mon, 2 Jul 2012 17:18:52 +0100
Subject: [PATCH] ENH: meshToMesh: allow deeper patch octree
---
.../meshToMesh/calculateMeshToMeshAddressing.C | 8 +++++---
1 file changed, 5 insertions(+), 3 deletions(-)
diff --git a/src/sampling/meshToMeshInterpolation/meshToMesh/calculateMeshToMeshAddressing.C b/src/sampling/meshToMeshInterpolation/meshToMesh/calculateMeshToMeshAddressing.C
index fda5525f095..83716c95c4f 100644
--- a/src/sampling/meshToMeshInterpolation/meshToMesh/calculateMeshToMeshAddressing.C
+++ b/src/sampling/meshToMeshInterpolation/meshToMesh/calculateMeshToMeshAddressing.C
@@ -171,13 +171,15 @@ void Foam::meshToMesh::calcAddressing()
wallBb.max() + vector(typDim, typDim, typDim)
);
+ // Note: allow more levels than in meshSearch. Assume patch
+ // is not as big as all boundary faces
indexedOctree<treeDataFace> oc
(
treeDataFace(false, fromPatch),
shiftedBb, // overall search domain
- 8, // maxLevel
+ 12, // maxLevel
10, // leafsize
- 3.0 // duplicity
+ 6.0 // duplicity
);
const vectorField::subField centresToBoundary =
@@ -185,7 +187,7 @@ void Foam::meshToMesh::calcAddressing()
boundaryAddressing_[patchi].setSize(toPatch.size());
- scalar distSqr = sqr(GREAT);
+ scalar distSqr = sqr(wallBb.mag());
forAll(toPatch, toi)
{
--
GitLab