From: eap <eap@opencascade.com>
Date: Tue, 16 Nov 2010 07:05:07 +0000 (+0000)
Subject: 0021074: EDF 1682 MESH: Bug with start and end length hypothesis
X-Git-Tag: tag_Before_ViscLayers_PrototypeVersion~20
X-Git-Url: http://git.salome-platform.org/gitweb/?p=modules%2Fsmesh.git;a=commitdiff_plain;h=89028932c04f2bad1141bf8c2d0a7073973221c5;ds=sidebyside

0021074: EDF 1682 MESH: Bug with start and end length hypothesis
  fix compensateError()
---

diff --git a/src/StdMeshers/StdMeshers_Regular_1D.cxx b/src/StdMeshers/StdMeshers_Regular_1D.cxx
index 6771e53ed..ce6742579 100644
--- a/src/StdMeshers/StdMeshers_Regular_1D.cxx
+++ b/src/StdMeshers/StdMeshers_Regular_1D.cxx
@@ -377,10 +377,14 @@ static void compensateError(double a1, double an,
     }
 
     double q  = dUn / ( nPar - 1 );
-    if ( !adjustNeighbors2an ) {
-      for ( itU = theParams.rbegin(), i = 1; i < nPar; itU++, i++ ) {
+    if ( !adjustNeighbors2an )
+    {
+      q = Abs( dUn / ( Utgt - Un )); // factor of segment length change
+      for ( itU = theParams.rbegin(), i = 1; i < nPar; i++ ) {
+        double prevU = *itU;
         (*itU) += dUn;
-        dUn -= q;
+        ++itU;
+        dUn = q * (*itU - prevU) * (prevU-U1)/(Un-U1);
       }
     }
     else {