[EDF24091] : Use of area per cell field was missing in computation of Force in Torseu...

[tools/paravisaddons_common.git] / src / TorseurCIH / plugin / TorseurCIHModule / vtkTorseurCIH.cxx

diff --git a/src/TorseurCIH/plugin/TorseurCIHModule/vtkTorseurCIH.cxx b/src/TorseurCIH/plugin/TorseurCIHModule/vtkTorseurCIH.cxx
index 9ac9f26ba214562506005fff7e7882c9ed21e8ad..f9f43b93dd028e0ab5d9273c7338775f4ecbeb45 100644 (file)
--- a/src/TorseurCIH/plugin/TorseurCIHModule/vtkTorseurCIH.cxx
+++ b/src/TorseurCIH/plugin/TorseurCIHModule/vtkTorseurCIH.cxx
@@ -702,7 +702,10 @@ vtkSmartPointer<vtkTable> ComputeTorseurCIH(vtkUnstructuredGrid* usgIn)
     F_y = ForceBuilder({ 3, 1, 5 }, matrix, eqn);
     F_z = ForceBuilder({ 4, 5, 2 }, matrix, eqn);
   }
+  //
+  //
   MCAuto<DataArrayDouble> F(DataArrayDouble::Meld({ F_x, F_y, F_z }));
+  F->multiplyEqual(area_vector);
   double ZeForce[3], normalFace[3];
   F->accumulate(ZeForce);
   eqn->accumulate(normalFace);