fudaa-svn-commit Mailing List for Fudaa (Page 167)

SourceForge Headquarters 225 Broadway Suite 1600 San Diego, CA 92101 +1 (858) 422-6466

Revision: 4430
          http://fudaa.svn.sourceforge.net/fudaa/?rev=4430&view=rev
Author:   emmanuel_martin
Date:     2009-01-30 17:09:25 +0000 (Fri, 30 Jan 2009)

Log Message:
-----------
Tache #164 : "Edition 1D : Ne pas limiter ?\195?\160 3 le nombre de segments sur les profils"

Modified Paths:
--------------
    branches/FudaaModeleur_TC1Bis/fudaa_devel/fudaa/src/org/fudaa/fudaa/modeleur/modeleur1d/model/Bief.java
    branches/FudaaModeleur_TC1Bis/fudaa_devel/fudaa/src/org/fudaa/fudaa/modeleur/modeleur1d/model/ProfilContainer.java
    branches/FudaaModeleur_TC1Bis/fudaa_devel/fudaa/src/org/fudaa/fudaa/modeleur/modeleur1d/model/ProfilContainerAdapter.java
    branches/FudaaModeleur_TC1Bis/fudaa_devel/fudaa/src/org/fudaa/fudaa/modeleur/modeleur1d/view/VueTableau.java

Modified: branches/FudaaModeleur_TC1Bis/fudaa_devel/fudaa/src/org/fudaa/fudaa/modeleur/modeleur1d/model/Bief.java
===================================================================

--- branches/FudaaModeleur_TC1Bis/fudaa_devel/fudaa/src/org/fudaa/fudaa/modeleur/modeleur1d/model/Bief.java	2009-01-30 15:11:50 UTC (rev 4429)
+++ branches/FudaaModeleur_TC1Bis/fudaa_devel/fudaa/src/org/fudaa/fudaa/modeleur/modeleur1d/model/Bief.java	2009-01-30 17:09:25 UTC (rev 4430)
@@ -861,6 +861,16 @@
       else if(_att==GISAttributeConstants.INTERSECTION_LIMITE_STOCKAGE_DROITE)
         updateGeom(_idxGeom, limitesStockagsesModified_, 1, idxAttlsDroite_);
       else if(_att==GISAttributeConstants.INTERSECTIONS_LIGNES_DIRECTRICES)
+        /*
+         * Attention subtilit\xE9 : Quand l'attribut est
+         * GISAttributeConstants.INTERSECTIONS_LIGNES_DIRECTRICES c'est que
+         * l'\xE9v\xE9nement vient de la liste des index d'intersection avec une ligne
+         * directrice m\xE9moris\xE9e dans chaque profil. Or pour cette liste
+         * l'_idxGeom est en r\xE9alit\xE9 son index modifi\xE9. Vu que les lignes
+         * directrices sont ordonn\xE9es dans cette attribut, l'information qu'on
+         * re\xE7oit sous le nom _idxGeom est en faite l'index de la
+         * ligneDirectrice modifi\xE9e.
+         */
         updateLd(_idxGeom);
     }
 
@@ -870,7 +880,8 @@
         updateGeom(_idxGeom, rivesModified_, 1, idxAttRiveDroite_);
         updateGeom(_idxGeom, limitesStockagsesModified_, 0, idxAttlsGauche_);
         updateGeom(_idxGeom, limitesStockagsesModified_, 1, idxAttlsDroite_);
-        updateLd(_idxGeom);
+        for(int i=0;i<lignesDirectrices_.getNombre();i++)
+          updateLd(i);
       }
     }
     

Modified: branches/FudaaModeleur_TC1Bis/fudaa_devel/fudaa/src/org/fudaa/fudaa/modeleur/modeleur1d/model/ProfilContainer.java
===================================================================
--- branches/FudaaModeleur_TC1Bis/fudaa_devel/fudaa/src/org/fudaa/fudaa/modeleur/modeleur1d/model/ProfilContainer.java	2009-01-30 15:11:50 UTC (rev 4429)
+++ branches/FudaaModeleur_TC1Bis/fudaa_devel/fudaa/src/org/fudaa/fudaa/modeleur/modeleur1d/model/ProfilContainer.java	2009-01-30 17:09:25 UTC (rev 4430)
@@ -105,8 +105,6 @@
   /** Place l'intersection du profil avec la rive droite \xE0 l'index _idx. */
   public void setRiveDroite(int _idx, CtuluCommandContainer _cmd);
   
-  
-  
   /** Retourne le nombre de lignes directrices croisant le profil. */
   public int getNbLignesDirectrices();
   
@@ -127,4 +125,7 @@
   
   /** Retourne les index des lignes directrices pouvant \xEAtre d\xE9plac\xE9s au point indiqu\xE9.*/
   public int[] getAllowedMoveOfLignesDirectricesTo(int _idxPoint);
+  
+  /** Retourne vrai si le point indiqu\xE9 est un point de rupture. */
+  public boolean isRupturePoint(int _idx);
 }

Modified: branches/FudaaModeleur_TC1Bis/fudaa_devel/fudaa/src/org/fudaa/fudaa/modeleur/modeleur1d/model/ProfilContainerAdapter.java
===================================================================
--- branches/FudaaModeleur_TC1Bis/fudaa_devel/fudaa/src/org/fudaa/fudaa/modeleur/modeleur1d/model/ProfilContainerAdapter.java	2009-01-30 15:11:50 UTC (rev 4429)
+++ branches/FudaaModeleur_TC1Bis/fudaa_devel/fudaa/src/org/fudaa/fudaa/modeleur/modeleur1d/model/ProfilContainerAdapter.java	2009-01-30 17:09:25 UTC (rev 4430)
@@ -9,7 +9,10 @@
 
 import java.util.ArrayList;
 import java.util.BitSet;
+import java.util.HashSet;
+import java.util.Iterator;
 import java.util.List;
+import java.util.Set;
 
 import org.fudaa.ctulu.CtuluCommand;
 import org.fudaa.ctulu.CtuluCommandComposite;
@@ -22,10 +25,10 @@
 import org.fudaa.ctulu.gis.GISAttributeModelIntegerList;
 import org.fudaa.ctulu.gis.GISCoordinateSequenceContainerInterface;
 import org.fudaa.ctulu.gis.GISGeometryFactory;
-import org.fudaa.ctulu.gis.GISPolyligne;
 import org.fudaa.ctulu.gis.GISZoneCollection;
 import org.fudaa.ctulu.gis.GISZoneCollectionLigneBrisee;
 import org.fudaa.fudaa.commun.FudaaLib;
+import org.fudaa.fudaa.modeleur.LibUtils;
 
 import com.vividsolutions.jts.geom.Coordinate;
 import com.vividsolutions.jts.geom.CoordinateSequence;
@@ -97,6 +100,21 @@
     }
   }
   
+  /**
+   * Marge d'acceptation avant de consid\xE9r\xE9 un point comme \xE9tant un point de
+   * rupture. La valeur est la diff\xE9rence des coefficients directeurs des
+   * segments avant et apr\xE8s le point.
+   * En donnant la valeur en radian de l'angle minimal : tolerance_=Math.tan(#ValeurDegr\xE9#)
+   *    Avec #ValeurDegr\xE9# C [0;PI/2[
+   * En donnant la valeur en radian de l'angle maximal : tolerance_=Math.tan(Math.PI-#ValeurDegr\xE9#)
+   *    Avec #ValeurDegr\xE9# C ]PI/2;PI]
+   * En donnant la valeur en degr\xE9 de l'angle minimal : tolerance_=Math.tan(#ValeurDegr\xE9#*Math.PI/2/90)
+   *    Avec #ValeurDegr\xE9# C [0;90[
+   * En donnant la valeur en degr\xE9 de l'angle maximal : tolerance_=Math.tan((180-#ValeurDegr\xE9#)*Math.PI/2/90)
+   *    Avec #ValeurDegr\xE9# C ]90;180]
+   */
+  static final double tolerance_=Math.tan(89*Math.PI/2/90);
+  
   /** Le biefContainer contenant le bief contenant les profils utilis\xE9s. */
   private BiefContainer biefContainer_;
   /** Le profil selectionn\xE9. */
@@ -106,17 +124,22 @@
   // Caches \\
   private int idxZMax_;
   private int idxZMin_;
-  private double curviligneDecalage_; // d\xE9calage de l'axe hydraulique.
-  // Le profil est potentiellement coup\xE9 en trois, les deux indices qui suivent
-  // indique l'index de chacune de ces ruptures.
-  private int idxRupture1_;
-  private int idxRupture2_;
+  /** D\xE9calage de l'axe hydraulique. */
+  private double curviligneDecalage_;
+  /** Les points de ruptures du profils sont contenus dans cet ensemble. */
+  private Set<Integer> idxRuptures_=new HashSet<Integer>();
 
   public ProfilContainerAdapter(BiefContainer _biefContainer) {
     biefContainer_=_biefContainer;
     biefContainer_.addBiefContainerListener(this);
   }
 
+  public boolean isRupturePoint(int _idx) {
+    if(_idx<0||_idx>=getNbPoint())
+      throw new IllegalArgumentException(FudaaLib.getS("_idx d'appartient pas au profil."));
+    return idxRuptures_.contains(_idx);
+  }
+  
   /**
    * Change le profil selectionn\xE9.
    * @throws DataGeometryException 
@@ -202,7 +225,7 @@
    * Permet de choisir la g\xE9om\xE9trie sur lequel l'instance va travailler.
    * Hypoth\xE8se importante : les points de la g\xE9om\xE9tries sont correctement
    * ordonn\xE9s (les points sont ordonn\xE9s en ordre croissant de leur abscisse
-   * curviligne). La veracit\xE9 de cette hypoth\xE8se n'est PAS v\xE9rifi\xE9 dans setData.
+   * curviligne). La veracit\xE9 de cette hypoth\xE8se n'est PAS v\xE9rifi\xE9 dans updateData.
    * 
    * @param _zone la zone contenant la g\xE9om\xE9trie
    * @param _idxGeom l'indice de la g\xE9om\xE9trie dans _zone
@@ -213,6 +236,7 @@
       z_=null;
       curv_=null;
       curviligneDecalage_=0;
+      idxRuptures_.clear();
     }
     else if (biefContainer_.getZoneProfils()!=null&&idxProfilSelected_>=0&&idxProfilSelected_<biefContainer_.getZoneProfils().getNbGeometries()&&biefContainer_.getZoneProfils().getAttributeIsZ()!=null) {
       try {
@@ -249,11 +273,8 @@
             throw new ProfilContainerException(FudaaLib.getS("Le profil se coupe lui m\xEAme. Il n'est donc pas valide."));
         }
         
-        // / D\xE9termination des points de ruptures et verification qu'il y en a
-        // au maximun deux. \\\
-        idxRupture1_=-1;
-        idxRupture2_=-1;
-        double upsilon=0.0001; // Marge d'erreur des calculs sur double
+        /// D\xE9termination des points de ruptures. \\\
+        idxRuptures_.clear();
         double coefDirecteur=0;
         // coefDirX : Indique si le coefficient directeur est calcul\xE9 sur les x
         // ou sur les y (au d\xE9nominateur).
@@ -288,18 +309,22 @@
             newCoefDirecteur=0;
           }
           // D\xE9tection des ruptures
-          if (coefDirX!=newCoefDirX||Math.abs(coefDirecteur-newCoefDirecteur)>upsilon) {
-            if (idxRupture1_==-1)
-              idxRupture1_=idx-1;
-            else if (idxRupture2_==-1)
-              idxRupture2_=idx-1;
-            else
-              throw new ProfilContainerException(FudaaLib.getS("Il y a plus que deux ruptures dans le profil."));
-          }
+          if (coefDirX!=newCoefDirX||Math.abs(coefDirecteur-newCoefDirecteur)>tolerance_)
+            idxRuptures_.add(idx-1);
           // Pr\xE9paration de l'it\xE9ration suivante
           coefDirecteur=newCoefDirecteur;
           coefDirX=newCoefDirX;
         }
+        // Lin\xE9arisation du profil
+        if(idxRuptures_.size()>0) {
+          int[] idxRuptures=new int[idxRuptures_.size()+2];
+          idxRuptures[0]=0;
+          Iterator<Integer> it=idxRuptures_.iterator();
+          for(int i=0;i<idxRuptures_.size();i++)
+            idxRuptures[i+1]=it.next();
+          idxRuptures[idxRuptures.length-1]=seq.size()-1;
+          LibUtils.linearisePolyligne(biefContainer_.getZoneProfils(), idxProfilSelected_, idxRuptures, null);
+        }
       }
       catch(ProfilContainerException _exp) {
         // Remise dans un \xE9tat coh\xE9rent
@@ -307,6 +332,7 @@
         z_=null;
         curv_=null;
         curviligneDecalage_=0;
+        idxRuptures_.clear();
         fireProfilContainerDataModified();
         // Propagation de l'exception
         throw _exp;
@@ -701,18 +727,18 @@
      * m\xEAme raison. Pour ne pas trop perdre en performance, seuls les
      * coordonn\xE9es modifi\xE9es sont clon\xE9es.
      */
-    // Les nouvelles valeurs du points \xE0 d\xE9placer
-    double newX;
-    double newY;
-    // Delta pour savoir si des doubles sont \xE9gaux
-    double tolerance=0.000001;
+//    // Les nouvelles valeurs du points \xE0 d\xE9placer
+//    double newX;
+//    double newY;
+//    // Delta pour savoir si des doubles sont \xE9gaux
+//    double tolerance=0.000001;
     
     // Verifie que le point reste encadr\xE9 par le(s) m\xEAme(s) point(s)
     if((_idxPoint==0&&_value>=curv_.get(1))||(_idxPoint==seq.size()-1&&_value<=curv_.get(seq.size()-2))||(_idxPoint!=0&&_idxPoint!=seq.size()-1&&(_value<=curv_.get(_idxPoint-1)||_value>=curv_.get(_idxPoint+1))))
       throw new ProfilContainerException(FudaaLib.getS("Il est interdit de d\xE9placer un point au dela des points l'encadrant."));
     
     // Cas o\xF9 l'index n'est pas sur un point de rupture \\
-    if (_idxPoint!=idxRupture1_&&_idxPoint!=idxRupture2_) {
+    if (!idxRuptures_.contains(_idxPoint)) {
       // Le point reste encadr\xE9 par les deux m\xEAme points => le signe de xa-xb et
       // de ya-yb ne change pas
       Coordinate newCoord=UtilsProfil1d.getCoordinateXY(seq, _value);
@@ -726,202 +752,202 @@
         curv_.set(_idxPoint, _value);
     }
     
-    // Cas complexe : l'index est sur un point de rupture \\
-    else {
-      /* 
-       * Abstract :
-       * Pour trouver le point restectant la concervation de labcisse curviligne, est 
-       * utilis\xE9 par la suite une m\xE9thode d'intersection de deux cercles.
-       * - Soient A, B, C, trois points non confondus. Soit la polyligne A, C, B.
-       * On veut construire C' tq cruv(C')=c et que curv(B) reste constant.
-       * Donc :
-       *  dist(AC)+dist(CB)=dist(AC')+dist(C'B)
-       * Ce probl\xE8me peut se r\xE9duire \xE0 la recherche de l'intersection de deux cercles :
-       * C1=cercle(centre=A, rayon=c-curv(A)) avec C2=cercle(centre=B, rayon=curv(B)-c)
-       * 
-       * Cas 1 : ya-yb != 0:
-       *  La r\xE9solution de ce probl\xE8me donne l'\xE9quation :
-       *    x\xB2[M\xB2+1] + x[2yaM-2NM-2xa] + xa\xB2+ya\xB2+N\xB2-rayonC1\xB2-2yaN=0
-       *    avec N=(rayonC2\xB2-rayonC1\xB2-xb\xB2+xa\xB2-yb\xB2+ya\xB2)/(2(ya-yb))
-       *    et M=(xa-xb)/(ya-yb)
-       *  C'est une \xE9quation du second degr\xE9.
-       *  Pour la r\xE9solution de l'\xE9quation, nous nommerons :
-       *    equA=[M\xB2+1]
-       *    equB=[2yaM-2NM-2xa]
-       *    equC=xa\xB2+ya\xB2+N\xB2-rayonC1\xB2-2yaN
-       *  La r\xE9solution de l'\xE9quation r\xE9sultante donne deux solutions (c1' et c2').
-       *  Pour trouver celle qui nous interesse (si elle existe) la contrainte suivante doit \xEAtre v\xE9rifi\xE9e :
-       *    vecteurNorm\xE9(C, Proj(C, AB)) = vecteurNorm\xE9(CX', Proj(CX', AB))
-       *    A noter que Proj(C1', AB)=Proj(C2', AB) de plus pour simplifier cette r\xE9solution sera faite sur z=0.
-       * 
-       * Cas 2 : ya-yb = 0 et xa-xb != 0:
-       *  La r\xE9solution de ce prol\xE8me donne l'\xE9quation :
-       *    y\xB2-2yay+xa+K\xB2-2xaK-rayonC1\xB2=0
-       *    avec K=(rayonC2\xB2-rayonC1\xB2-xb\xB2+xa\xB2)/(2*(xa-xb))
-       *    # k est \xE9galement la valeur de x des deux solutions potentielles
-       *  C'est un \xE9quation du second degr\xE9.
-       *  Pour la r\xE9solution de l'aquation, nous nommerons :
-       *    equA=1
-       *    equB=2ya
-       *    equC=K\xB2-2xaK-rayonC1\xB2
-       *  La r\xE9solution de l'\xE9quation r\xE9sultante donne deux solutions (c1' et c2').
-       *  Pour trouver celle qui nous interesse (si elle existe) la contrainte suivante doit \xEAtre v\xE9rifi\xE9e :
-       *    vecteurNorm\xE9(C, Proj(C, AB)) = vecteurNorm\xE9(CX', Proj(CX', AB))
-       *    A noter que Proj(C1', AB)=Proj(C2', AB) de plus pour simplifier cette r\xE9solution sera faite sur z=0.
-       * 
-       * Cas 3 : ya-yb = 0 et xa-xb=0:
-       *  quel qu'en soit le r\xE9sultat, c'est incoh\xE9rent pour notre contexte.
-       *  
-       * De plus on verifie que le r\xE9sultat ne provoque pas de croisement avec les autres axes
-       * du profil.
-       * 
-       * Une fois ce point trouv\xE9, on va \xE9galement d\xE9placer tous les points des deux axes adjacents,
-       * pour que le profil concerve sa forme g\xE9n\xE9rale. Ce r\xE9aligement se faire en concervant les 
-       * abscisses curvilignes de chacun des points.
-       * Pour y parvenir le nouveau coefficient directeur des axes est calcul\xE9 et appliqu\xE9 au rapport de
-       * l'abcisse curviligne avec l'abscisse du point extr\xE9mit\xE9 de l'axe.
-       * 
-       * Derni\xE8re petite note pour la fin, le cas dit 'simple' pr\xE9c\xE9dent
-       * (d\xE9placement d'un point qui n'est pas un point de rupture) n'est qu'un
-       * cas particulier de la r\xE9solution ci dessus (la solution de l'\xE9quation
-       * du second degr\xE9 serait unique). On pourrait donc l'enlever, mais je pense que
-       * s\xE9parer ces deux cas ne fait pas de mal pour une \xE9ventuelle relecture
-       * et/ou modification de ce code.
-       */
-      
-      // Application de la m\xE9thode d\xE9crite ci dessus \\
-      
-      // Valuation des variables d'entr\xE9es n\xE9c\xE9ssaires
-      final int idxPtdeb=_idxPoint==idxRupture1_?0:idxRupture1_;
-      final int idxPtFin=_idxPoint==idxRupture2_||idxRupture2_==-1?seq.size()-1:idxRupture2_;
-      final Coordinate a=seq.getCoordinate(idxPtdeb);
-      final Coordinate b=seq.getCoordinate(idxPtFin);
-      final Coordinate c=seq.getCoordinate(_idxPoint);
-      final double rayonC1=_value-curv_.get(idxPtdeb);
-      final double rayonC2=curv_.get(idxPtFin)-_value;
-      // Solutions
-      double x1;
-      double y1;
-      double x2;
-      double y2;
-      // Cas 1
-      if(Math.abs(a.y-b.y)>tolerance){
-        // Variables interm\xE9diaires
-        double n=(rayonC2*rayonC2-rayonC1*rayonC1-b.x*b.x+a.x*a.x-b.y*b.y+a.y*a.y)/(2*(a.y-b.y));
-        double m=(a.x-b.x)/(a.y-b.y);
-        double equA=m*m+1;
-        double equB=2*(a.y*m-n*m-a.x);
-        double equC=a.x*a.x+a.y*a.y+n*n-rayonC1*rayonC1-2*a.y*n;
-        // Variables de r\xE9solution
-        double deltaTmp=equB*equB-4*equA*equC;
-        if(deltaTmp<0)
-          throw new ProfilContainerException(FudaaLib.getS("Ce d\xE9placement n'est pas possible sans violer la contrainte de constance de la longueur du profil."));
-        double delta=Math.sqrt(equB*equB-4*equA*equC);
-        x1=(-equB+delta)/(2*equA);
-        y1=n-x1*m;
-        x2=(-equB-delta)/(2*equA);
-        y2=n-x2*m;
-      }
-      // Cas 2
-      else if(Math.abs(a.x-b.x)>tolerance){
-        // Variables interm\xE9diaires
-        double k=(rayonC2*rayonC2-rayonC1*rayonC1-b.x*b.x+a.x*a.x)/(2*(a.x-b.x));
-        double equA=1;
-        double equB=2*a.y;
-        double equC=k*k-2*a.x*k-rayonC1*rayonC1;
-        // Variables de r\xE9solution
-        double deltaTmp=equB*equB-4*equA*equC;
-        if(deltaTmp<0)
-          throw new ProfilContainerException(FudaaLib.getS("Ce d\xE9placement n'est pas possible sans violer la contrainte de constance de la longueur du profil."));
-        double delta=Math.sqrt(equB*equB-4*equA*equC);
-        y1=(-equB+delta)/(2*equA);
-        x1=k;
-        y2=(-equB-delta)/(2*equA);
-        x2=k;
-      }
-      // Cas 3
-      else
-        throw new ProfilContainerException(FudaaLib.getS("Ce d\xE9placement n'est pas valide. Ce cas ne devrait jamais arriver. CSI_1"));
-      // Recherche du r\xE9sultat nous interessant dans ce cas ci
-      Coordinate aZ0=new Coordinate(a.x, a.y, 0);
-      Coordinate bZ0=new Coordinate(b.x, b.y, 0);
-      Coordinate cZ0=new Coordinate(c.x, c.y, 0);
-      Coordinate vecReference=UtilsProfil1d.vec(cZ0, UtilsProfil1d.proj(cZ0, aZ0, bZ0));
-      Coordinate cProj=UtilsProfil1d.proj(new Coordinate(x1, y1, 0), aZ0, bZ0);
-      assert(UtilsProfil1d.egal(cProj, UtilsProfil1d.proj(new Coordinate(x2, y2, 0), aZ0, bZ0), 0.0001, false));
-      if(UtilsProfil1d.egal(vecReference, UtilsProfil1d.vec(new Coordinate(x1, y1, 0), cProj), 0.0001, true)){
-        newX=x1;
-        newY=y1;
-      }
-      else if(UtilsProfil1d.egal(vecReference, UtilsProfil1d.vec(new Coordinate(x2, y2, 0), cProj), 0.0001, true)) {
-        newX=x2;
-        newY=y2;
-      }
-      else
-        throw new ProfilContainerException(FudaaLib.getS("Ce d\xE9placement n'est pas valide. Ce cas ne devrait jamais arriver. CSI_2"));
-
-      // Verifie que le point trouv\xE9 ne coupe pas un autre axe du profil \\
-      boolean noCut=true;
-      Coordinate newP=new Coordinate(newX, newY);
-      GISPolyligne axe1=GISGeometryFactory.INSTANCE.createLineString(new Coordinate[]{a, newP});
-      GISPolyligne axe2=GISGeometryFactory.INSTANCE.createLineString(new Coordinate[]{newP, b});
-      if(idxPtdeb>0) {
-        /*
-         * Test de croisement entre les deux axes modifi\xE9s et l'axe avant. l'axe
-         * a, newP a normalement un point d'intersectino (a).
-         */
-        GISPolyligne axeTest=GISGeometryFactory.INSTANCE.createLineString(new Coordinate[]{new Coordinate(coords[0]), new Coordinate(coords[idxPtdeb])});
-        noCut=axe1.intersection(axeTest).getNumPoints()==1;
-        noCut=noCut&&!axe2.intersects(axeTest);
-      }
-      if(idxPtFin<seq.size()-1) {
-        /*
-         * Test de croisement entre les deux axes modifi\xE9s et l'axe avant. l'axe
-         * newP, b a normalement un point d'intersectino (b).
-         */
-        GISPolyligne axeTest=GISGeometryFactory.INSTANCE.createLineString(new Coordinate[]{new Coordinate(coords[idxPtFin]), new Coordinate(coords[coords.length-1])});
-        noCut=noCut&&axe2.intersection(axeTest).getNumPoints()==1;
-        noCut=noCut&&!axe1.intersects(axeTest);
-      }
-      // Traitement de l'\xE9ventuelle d\xE9tection d'un croisement
-      if(!noCut)
-        throw new ProfilContainerException(FudaaLib.getS("Ce d\xE9placement est impossible sans obtenir un profil dont des axes se coupent."));
-      // Enregistrement du nouveau point C' et de la nouvelle valeur curviligne
-      coords[_idxPoint]=new Coordinate(newX, newY, coords[_idxPoint].z);
-      // Mise a jour de la table des valeurs curvilignes
-      curv_.set(_idxPoint, _value);
-      
-      // R\xE9alignement des points des axes adjacents \\
-      for(int inc=0;inc<2;inc++){
-        // D\xE9termination des points de d\xE9but et de fin d'axe
-        int idxDeb;
-        int idxFin;
-        if(inc==0) { // premi\xE8re it\xE9ration, axe : debut, C'
-          idxDeb=_idxPoint==idxRupture1_?0:idxRupture1_;
-          idxFin=_idxPoint;
-        }
-        else { // seconde it\xE9ration, axe : C', fin
-          idxDeb=_idxPoint;
-          idxFin=_idxPoint==idxRupture2_||idxRupture2_==-1 ? seq.size()-1:idxRupture2_;
-        }
-        // R\xE9alignement sur l'axe
-        if (Math.abs(coords[idxFin].x-coords[idxDeb].x)>tolerance) {
-          double sign=Math.signum(coords[idxFin].x-coords[idxDeb].x);
-          double newCoefDirecteur=(coords[idxFin].y-coords[idxDeb].y)/(coords[idxFin].x-coords[idxDeb].x);
-          for (int i=idxDeb+1; i<idxFin; i++) {
-            coords[i]=new Coordinate(0, 0, coords[i].z);
-            coords[i].x=coords[idxDeb].x+sign*(curv_.get(i)-curv_.get(idxDeb))/Math.sqrt(1+newCoefDirecteur*newCoefDirecteur);
-            coords[i].y=coords[idxDeb].y+newCoefDirecteur*(coords[i].x-coords[idxDeb].x);
-          }
-        }
-        else { // Cas du bout de profil vertical
-          double sign=Math.signum(coords[idxFin].y-coords[idxDeb].y);
-          for (int i=idxDeb+1; i<idxFin; i++)
-            coords[i]=new Coordinate(coords[idxDeb].x, coords[idxDeb].y+sign*(curv_.get(i)-curv_.get(idxDeb)),
-                coords[i].z);
-        }
-      }
-    }
+//    // Cas complexe : l'index est sur un point de rupture \\
+//    else {
+//      /* 
+//       * Abstract :
+//       * Pour trouver le point restectant la concervation de labcisse curviligne, est 
+//       * utilis\xE9 par la suite une m\xE9thode d'intersection de deux cercles.
+//       * - Soient A, B, C, trois points non confondus. Soit la polyligne A, C, B.
+//       * On veut construire C' tq cruv(C')=c et que curv(B) reste constant.
+//       * Donc :
+//       *  dist(AC)+dist(CB)=dist(AC')+dist(C'B)
+//       * Ce probl\xE8me peut se r\xE9duire \xE0 la recherche de l'intersection de deux cercles :
+//       * C1=cercle(centre=A, rayon=c-curv(A)) avec C2=cercle(centre=B, rayon=curv(B)-c)
+//       * 
+//       * Cas 1 : ya-yb != 0:
+//       *  La r\xE9solution de ce probl\xE8me donne l'\xE9quation :
+//       *    x\xB2[M\xB2+1] + x[2yaM-2NM-2xa] + xa\xB2+ya\xB2+N\xB2-rayonC1\xB2-2yaN=0
+//       *    avec N=(rayonC2\xB2-rayonC1\xB2-xb\xB2+xa\xB2-yb\xB2+ya\xB2)/(2(ya-yb))
+//       *    et M=(xa-xb)/(ya-yb)
+//       *  C'est une \xE9quation du second degr\xE9.
+//       *  Pour la r\xE9solution de l'\xE9quation, nous nommerons :
+//       *    equA=[M\xB2+1]
+//       *    equB=[2yaM-2NM-2xa]
+//       *    equC=xa\xB2+ya\xB2+N\xB2-rayonC1\xB2-2yaN
+//       *  La r\xE9solution de l'\xE9quation r\xE9sultante donne deux solutions (c1' et c2').
+//       *  Pour trouver celle qui nous interesse (si elle existe) la contrainte suivante doit \xEAtre v\xE9rifi\xE9e :
+//       *    vecteurNorm\xE9(C, Proj(C, AB)) = vecteurNorm\xE9(CX', Proj(CX', AB))
+//       *    A noter que Proj(C1', AB)=Proj(C2', AB) de plus pour simplifier cette r\xE9solution sera faite sur z=0.
+//       * 
+//       * Cas 2 : ya-yb = 0 et xa-xb != 0:
+//       *  La r\xE9solution de ce prol\xE8me donne l'\xE9quation :
+//       *    y\xB2-2yay+xa+K\xB2-2xaK-rayonC1\xB2=0
+//       *    avec K=(rayonC2\xB2-rayonC1\xB2-xb\xB2+xa\xB2)/(2*(xa-xb))
+//       *    # k est \xE9galement la valeur de x des deux solutions potentielles
+//       *  C'est un \xE9quation du second degr\xE9.
+//       *  Pour la r\xE9solution de l'aquation, nous nommerons :
+//       *    equA=1
+//       *    equB=2ya
+//       *    equC=K\xB2-2xaK-rayonC1\xB2
+//       *  La r\xE9solution de l'\xE9quation r\xE9sultante donne deux solutions (c1' et c2').
+//       *  Pour trouver celle qui nous interesse (si elle existe) la contrainte suivante doit \xEAtre v\xE9rifi\xE9e :
+//       *    vecteurNorm\xE9(C, Proj(C, AB)) = vecteurNorm\xE9(CX', Proj(CX', AB))
+//       *    A noter que Proj(C1', AB)=Proj(C2', AB) de plus pour simplifier cette r\xE9solution sera faite sur z=0.
+//       * 
+//       * Cas 3 : ya-yb = 0 et xa-xb=0:
+//       *  quel qu'en soit le r\xE9sultat, c'est incoh\xE9rent pour notre contexte.
+//       *  
+//       * De plus on verifie que le r\xE9sultat ne provoque pas de croisement avec les autres axes
+//       * du profil.
+//       * 
+//       * Une fois ce point trouv\xE9, on va \xE9galement d\xE9placer tous les points des deux axes adjacents,
+//       * pour que le profil concerve sa forme g\xE9n\xE9rale. Ce r\xE9aligement se faire en concervant les 
+//       * abscisses curvilignes de chacun des points.
+//       * Pour y parvenir le nouveau coefficient directeur des axes est calcul\xE9 et appliqu\xE9 au rapport de
+//       * l'abcisse curviligne avec l'abscisse du point extr\xE9mit\xE9 de l'axe.
+//       * 
+//       * Derni\xE8re petite note pour la fin, le cas dit 'simple' pr\xE9c\xE9dent
+//       * (d\xE9placement d'un point qui n'est pas un point de rupture) n'est qu'un
+//       * cas particulier de la r\xE9solution ci dessus (la solution de l'\xE9quation
+//       * du second degr\xE9 serait unique). On pourrait donc l'enlever, mais je pense que
+//       * s\xE9parer ces deux cas ne fait pas de mal pour une \xE9ventuelle relecture
+//       * et/ou modification de ce code.
+//       */
+//      
+//      // Application de la m\xE9thode d\xE9crite ci dessus \\
+//      
+//      // Valuation des variables d'entr\xE9es n\xE9c\xE9ssaires
+//      final int idxPtdeb=_idxPoint==idxRupture1_?0:idxRupture1_;
+//      final int idxPtFin=_idxPoint==idxRupture2_||idxRupture2_==-1?seq.size()-1:idxRupture2_;
+//      final Coordinate a=seq.getCoordinate(idxPtdeb);
+//      final Coordinate b=seq.getCoordinate(idxPtFin);
+//      final Coordinate c=seq.getCoordinate(_idxPoint);
+//      final double rayonC1=_value-curv_.get(idxPtdeb);
+//      final double rayonC2=curv_.get(idxPtFin)-_value;
+//      // Solutions
+//      double x1;
+//      double y1;
+//      double x2;
+//      double y2;
+//      // Cas 1
+//      if(Math.abs(a.y-b.y)>tolerance){
+//        // Variables interm\xE9diaires
+//        double n=(rayonC2*rayonC2-rayonC1*rayonC1-b.x*b.x+a.x*a.x-b.y*b.y+a.y*a.y)/(2*(a.y-b.y));
+//        double m=(a.x-b.x)/(a.y-b.y);
+//        double equA=m*m+1;
+//        double equB=2*(a.y*m-n*m-a.x);
+//        double equC=a.x*a.x+a.y*a.y+n*n-rayonC1*rayonC1-2*a.y*n;
+//        // Variables de r\xE9solution
+//        double deltaTmp=equB*equB-4*equA*equC;
+//        if(deltaTmp<0)
+//          throw new ProfilContainerException(FudaaLib.getS("Ce d\xE9placement n'est pas possible sans violer la contrainte de constance de la longueur du profil."));
+//        double delta=Math.sqrt(equB*equB-4*equA*equC);
+//        x1=(-equB+delta)/(2*equA);
+//        y1=n-x1*m;
+//        x2=(-equB-delta)/(2*equA);
+//        y2=n-x2*m;
+//      }
+//      // Cas 2
+//      else if(Math.abs(a.x-b.x)>tolerance){
+//        // Variables interm\xE9diaires
+//        double k=(rayonC2*rayonC2-rayonC1*rayonC1-b.x*b.x+a.x*a.x)/(2*(a.x-b.x));
+//        double equA=1;
+//        double equB=2*a.y;
+//        double equC=k*k-2*a.x*k-rayonC1*rayonC1;
+//        // Variables de r\xE9solution
+//        double deltaTmp=equB*equB-4*equA*equC;
+//        if(deltaTmp<0)
+//          throw new ProfilContainerException(FudaaLib.getS("Ce d\xE9placement n'est pas possible sans violer la contrainte de constance de la longueur du profil."));
+//        double delta=Math.sqrt(equB*equB-4*equA*equC);
+//        y1=(-equB+delta)/(2*equA);
+//        x1=k;
+//        y2=(-equB-delta)/(2*equA);
+//        x2=k;
+//      }
+//      // Cas 3
+//      else
+//        throw new ProfilContainerException(FudaaLib.getS("Ce d\xE9placement n'est pas valide. Ce cas ne devrait jamais arriver. CSI_1"));
+//      // Recherche du r\xE9sultat nous interessant dans ce cas ci
+//      Coordinate aZ0=new Coordinate(a.x, a.y, 0);
+//      Coordinate bZ0=new Coordinate(b.x, b.y, 0);
+//      Coordinate cZ0=new Coordinate(c.x, c.y, 0);
+//      Coordinate vecReference=UtilsProfil1d.vec(cZ0, UtilsProfil1d.proj(cZ0, aZ0, bZ0));
+//      Coordinate cProj=UtilsProfil1d.proj(new Coordinate(x1, y1, 0), aZ0, bZ0);
+//      assert(UtilsProfil1d.egal(cProj, UtilsProfil1d.proj(new Coordinate(x2, y2, 0), aZ0, bZ0), 0.0001, false));
+//      if(UtilsProfil1d.egal(vecReference, UtilsProfil1d.vec(new Coordinate(x1, y1, 0), cProj), 0.0001, true)){
+//        newX=x1;
+//        newY=y1;
+//      }
+//      else if(UtilsProfil1d.egal(vecReference, UtilsProfil1d.vec(new Coordinate(x2, y2, 0), cProj), 0.0001, true)) {
+//        newX=x2;
+//        newY=y2;
+//      }
+//      else
+//        throw new ProfilContainerException(FudaaLib.getS("Ce d\xE9placement n'est pas valide. Ce cas ne devrait jamais arriver. CSI_2"));
+//
+//      // Verifie que le point trouv\xE9 ne coupe pas un autre axe du profil \\
+//      boolean noCut=true;
+//      Coordinate newP=new Coordinate(newX, newY);
+//      GISPolyligne axe1=GISGeometryFactory.INSTANCE.createLineString(new Coordinate[]{a, newP});
+//      GISPolyligne axe2=GISGeometryFactory.INSTANCE.createLineString(new Coordinate[]{newP, b});
+//      if(idxPtdeb>0) {
+//        /*
+//         * Test de croisement entre les deux axes modifi\xE9s et l'axe avant. l'axe
+//         * a, newP a normalement un point d'intersectino (a).
+//         */
+//        GISPolyligne axeTest=GISGeometryFactory.INSTANCE.createLineString(new Coordinate[]{new Coordinate(coords[0]), new Coordinate(coords[idxPtdeb])});
+//        noCut=axe1.intersection(axeTest).getNumPoints()==1;
+//        noCut=noCut&&!axe2.intersects(axeTest);
+//      }
+//      if(idxPtFin<seq.size()-1) {
+//        /*
+//         * Test de croisement entre les deux axes modifi\xE9s et l'axe avant. l'axe
+//         * newP, b a normalement un point d'intersectino (b).
+//         */
+//        GISPolyligne axeTest=GISGeometryFactory.INSTANCE.createLineString(new Coordinate[]{new Coordinate(coords[idxPtFin]), new Coordinate(coords[coords.length-1])});
+//        noCut=noCut&&axe2.intersection(axeTest).getNumPoints()==1;
+//        noCut=noCut&&!axe1.intersects(axeTest);
+//      }
+//      // Traitement de l'\xE9ventuelle d\xE9tection d'un croisement
+//      if(!noCut)
+//        throw new ProfilContainerException(FudaaLib.getS("Ce d\xE9placement est impossible sans obtenir un profil dont des axes se coupent."));
+//      // Enregistrement du nouveau point C' et de la nouvelle valeur curviligne
+//      coords[_idxPoint]=new Coordinate(newX, newY, coords[_idxPoint].z);
+//      // Mise a jour de la table des valeurs curvilignes
+//      curv_.set(_idxPoint, _value);
+//      
+//      // R\xE9alignement des points des axes adjacents \\
+//      for(int inc=0;inc<2;inc++){
+//        // D\xE9termination des points de d\xE9but et de fin d'axe
+//        int idxDeb;
+//        int idxFin;
+//        if(inc==0) { // premi\xE8re it\xE9ration, axe : debut, C'
+//          idxDeb=_idxPoint==idxRupture1_?0:idxRupture1_;
+//          idxFin=_idxPoint;
+//        }
+//        else { // seconde it\xE9ration, axe : C', fin
+//          idxDeb=_idxPoint;
+//          idxFin=_idxPoint==idxRupture2_||idxRupture2_==-1 ? seq.size()-1:idxRupture2_;
+//        }
+//        // R\xE9alignement sur l'axe
+//        if (Math.abs(coords[idxFin].x-coords[idxDeb].x)>tolerance) {
+//          double sign=Math.signum(coords[idxFin].x-coords[idxDeb].x);
+//          double newCoefDirecteur=(coords[idxFin].y-coords[idxDeb].y)/(coords[idxFin].x-coords[idxDeb].x);
+//          for (int i=idxDeb+1; i<idxFin; i++) {
+//            coords[i]=new Coordinate(0, 0, coords[i].z);
+//            coords[i].x=coords[idxDeb].x+sign*(curv_.get(i)-curv_.get(idxDeb))/Math.sqrt(1+newCoefDirecteur*newCoefDirecteur);
+//            coords[i].y=coords[idxDeb].y+newCoefDirecteur*(coords[i].x-coords[idxDeb].x);
+//          }
+//        }
+//        else { // Cas du bout de profil vertical
+//          double sign=Math.signum(coords[idxFin].y-coords[idxDeb].y);
+//          for (int i=idxDeb+1; i<idxFin; i++)
+//            coords[i]=new Coordinate(coords[idxDeb].x, coords[idxDeb].y+sign*(curv_.get(i)-curv_.get(idxDeb)),
+//                coords[i].z);
+//        }
+//      }
+//    }
     CtuluCommandComposite cmd=new CtuluCommandComposite("Changement abs. curviligne");
     // Enregistrement des nouvelles coordonn\xE9es \\
     biefContainer_.getZoneProfils().setCoordinateSequence(idxProfilSelected_, GISGeometryFactory.INSTANCE.getCoordinateSequenceFactory().create(coords), cmd);

Modified: branches/FudaaModeleur_TC1Bis/fudaa_devel/fudaa/src/org/fudaa/fudaa/modeleur/modeleur1d/view/VueTableau.java
===================================================================
--- branches/FudaaModeleur_TC1Bis/fudaa_devel/fudaa/src/org/fudaa/fudaa/modeleur/modeleur1d/view/VueTableau.java	2009-01-30 15:11:50 UTC (rev 4429)
+++ branches/FudaaModeleur_TC1Bis/fudaa_devel/fudaa/src/org/fudaa/fudaa/modeleur/modeleur1d/view/VueTableau.java	2009-01-30 17:09:25 UTC (rev 4430)
@@ -112,7 +112,7 @@
       if(column==0)
         return false;
       else
-        return true;
+        return !data_.isRupturePoint(row);
     }
     
     public Class<?> getColumnClass(int columnIndex) {


This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site.




2008	Jan	Feb	Mar	Apr	May	Jun (46)	Jul (37)	Aug (154)	Sep (140)	Oct (132)	Nov (104)	Dec (67)
2009	Jan (113)	Feb (73)	Mar (102)	Apr (106)	May (114)	Jun (67)	Jul (116)	Aug (48)	Sep (108)	Oct (296)	Nov (56)	Dec (53)
2010	Jan (95)	Feb (31)	Mar (40)	Apr (12)	May (10)	Jun (27)	Jul (19)	Aug (81)	Sep (48)	Oct (45)	Nov (40)	Dec (7)
2011	Jan (16)	Feb (32)	Mar (55)	Apr (38)	May (24)	Jun (37)	Jul (12)	Aug (16)	Sep (97)	Oct (90)	Nov (101)	Dec (89)
2012	Jan (18)	Feb (2)	Mar (54)	Apr (69)	May (97)	Jun (131)	Jul (70)	Aug (48)	Sep (48)	Oct (119)	Nov (194)	Dec (60)
2013	Jan (73)	Feb (35)	Mar (42)	Apr (28)	May (7)	Jun (17)	Jul (27)	Aug (10)	Sep (38)	Oct (12)	Nov (4)	Dec (16)
2014	Jan (33)	Feb (37)	Mar (19)	Apr (3)	May (47)	Jun (26)	Jul (22)	Aug	Sep (10)	Oct (30)	Nov (24)	Dec (19)
2015	Jan (13)	Feb (16)	Mar (36)	Apr (19)	May (5)	Jun (5)	Jul (3)	Aug (11)	Sep (22)	Oct (7)	Nov (14)	Dec
2016	Jan	Feb (26)	Mar (13)	Apr (61)	May	Jun	Jul (4)	Aug	Sep	Oct (27)	Nov (14)	Dec (21)
2017	Jan (30)	Feb (4)	Mar (3)	Apr (5)	May (69)	Jun (29)	Jul (1)	Aug	Sep	Oct	Nov	Dec

fudaa-svn-commit Mailing List for Fudaa (Page 167)

https://gitlab.com/fudaa

fudaa-svn-commit — svn commit with diffs