From b1fe3cc51a56011a15442fd365ac4cba7727ddfa Mon Sep 17 00:00:00 2001 From: mattijs Date: Thu, 19 Feb 2009 11:15:32 +0000 Subject: [PATCH] processorFvPatchField now holds opposite value; not interpolate --- .../decomposePar/fvFieldDecomposer.C | 24 ++++++++++++------- 1 file changed, 15 insertions(+), 9 deletions(-) diff --git a/applications/utilities/parallelProcessing/decomposePar/fvFieldDecomposer.C b/applications/utilities/parallelProcessing/decomposePar/fvFieldDecomposer.C index 75b3ce779a..15eb039bfd 100644 --- a/applications/utilities/parallelProcessing/decomposePar/fvFieldDecomposer.C +++ b/applications/utilities/parallelProcessing/decomposePar/fvFieldDecomposer.C @@ -59,7 +59,6 @@ processorVolPatchFieldDecomposer addressing_(addressingSlice.size()), weights_(addressingSlice.size()) { - const scalarField& weights = mesh.weights().internalField(); const labelList& own = mesh.faceOwner(); const labelList& neighb = mesh.faceNeighbour(); @@ -72,15 +71,22 @@ processorVolPatchFieldDecomposer { // This is a regular face. it has been an internal face // of the original mesh and now it has become a face - // on the parallel boundary - addressing_[i].setSize(2); - weights_[i].setSize(2); + // on the parallel boundary. + // Give face the value of the neighbour. - addressing_[i][0] = own[ai]; - addressing_[i][1] = neighb[ai]; + addressing_[i].setSize(1); + weights_[i].setSize(1); + weights_[i][0] = 1.0; - weights_[i][0] = weights[ai]; - weights_[i][1] = 1.0 - weights[ai]; + if (addressingSlice[i] >= 0) + { + // I have the owner so use the neighbour value + addressing_[i][0] = neighb[ai]; + } + else + { + addressing_[i][0] = own[ai]; + } } else { @@ -89,7 +95,7 @@ processorVolPatchFieldDecomposer // do the interpolation properly (I would need to look // up the different (face) list of data), so I will // just grab the value from the owner cell - // + addressing_[i].setSize(1); weights_[i].setSize(1);