Changeset 17

Timestamp:

05/08/08 15:46:59 (8 months ago)

Author:

jpv

Message:

Ajout des équations de plans et des sections/separation
de solides

Files:

• trunk/solides.pro (modified) (15 diffs)
• trunk/tests/index.txt (modified) (1 diff)

trunk/solides.pro

@@ -1 +1 @@
-%!
-% PostScript prologue for pst-solides3d.tex.
-0, 2008/04/13
+1, 2008/05/08
-%
-%% COPYRIGHT 2008 by Jean-Paul Vignault, Manuel Luque
@@ -559 +559 @@
-
-%%%%% ### projtext ###
+%% syntaxe : str x0 y0 z0 [normal_vect] ultextp3d --> -
+%% syntaxe : str x0 y0 z0 [normal_vect] bool ultextp3d --> -
+%% syntaxe : str x0 y0 plantype ultextp3d --> -
+%% syntaxe : str x0 y0 plantype bool ultextp3d --> -
+%% syntaxe : str1 solid i str2 ultextp3d --> -
+%% syntaxe : str1 solid i str2 bool ultextp3d --> -
+%% syntaxe : str1 solid i alpha str2 bool ultextp3d --> -
- /initpr@jtext {
-5 dict begin
@@ -566 +573 @@
-      /mybool true def
-   } ifelse
-
-
+
-      /type_plan_proj true def
-
+
+
+
+
+
-      /z0 exch def
-      /y0 exch def
-      /x0 exch def
-
-
-
-
-
-
-
+
+
+
+
+
+
+
+
+
+
-      } {
-
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
-      } ifelse
-      %% y a-t-il un alpha
-      2 copy pop issolid {
-         /alpha 0 def
-      } {
-         /alpha exch def
-      } ifelse
-      /i exch def
-      /solid exch def
-      0 0
-   } ifelse
-} def
@@ -613 +633 @@
-/cltextp3d {initpr@jtext cltext_ closepr@jtext} def
-/bltextp3d {initpr@jtext bltext_ closepr@jtext} def
-b
+d
-/ubtextp3d {initpr@jtext ubtext_ closepr@jtext} def
-/cbtextp3d {initpr@jtext cbtext_ closepr@jtext} def
-/bbtextp3d {initpr@jtext bbtext_ closepr@jtext} def
-b
+d
-/uctextp3d {initpr@jtext uctext_ closepr@jtext} def
-/cctextp3d {initpr@jtext cctext_ closepr@jtext} def
-/bctextp3d {initpr@jtext bctext_ closepr@jtext} def
-b
+d
-/urtextp3d {initpr@jtext urtext_ closepr@jtext} def
-/crtextp3d {initpr@jtext crtext_ closepr@jtext} def
-/brtextp3d {initpr@jtext brtext_ closepr@jtext} def
-b
+d
-
-%%%%% ### currentppathtransform ###
@@ -797 +817 @@
-      /mybool true def
-   } ifelse
-
-
-
-
-
-
-
-
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
-%% syntaxe : solid i str bool solidprojpath --> -
-%% ou
@@ -2775 +2828 @@
-} def
-
-%% 
-%% 
-%% 
-%% 
-
-%% 
+
+
+
+
+
+
-/aretescachees true def
-%% 
+
-/solidgridOn {
-   /solidgrid true def
@@ -2832 +2885 @@
-
-%% syntaxe : solid table solidaddface -> -
+%% syntaxe : solid table (couleur) solidaddface -> -
-%% on ne se preoccupe pas des faces internes
-/solidaddface {
-
+
+
+
+
+
+
-   /table exch def
-   /solid exch def
@@ -2845 +2904 @@
-%   solid solidnombrefaces
-    solid [
-() 
+lac@uleur
-    ] solidputoutcolors 
-end
@@ -3592 +3651 @@
-} def
-
+/solidputcolors {
+3 dict begin
+   2 copy pop
+   isstring {
+      inoutputcolors
+   } {
+      outputcolors
+   } ifelse
+end
+} def
+
-%%%%% ### solidputhuecolors ###
-%% syntaxe : solid table solidputhuecolors --> -
@@ -3926 +3996 @@
-
-%%%%% ### solidlight ###
-%/lightintensity 1 def
-%/lightsrc {10 15 20} def
-
-/setlightintensity {
-   /lightintensity exch def
@@ -4275 +4342 @@
-
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%%%%                  plans affines                     %%%%
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+%%%%% ### planaffine ###
+%% plan : origine, base, range, ngrid
+%% [0 0 0 [1 0 0 0 1 0] [-3 3 -2 2] [1. 1.] ]
+
+/explan [0 0 0 [1 0 0 0 1 0 0 0 1] [-3 3 -2 2] [1. 1.] ] def
+
+%% syntaxe : any isplan --> bool
+/isplan {
+1 dict begin
+   /candidat exch def
+   candidat isarray {
+      candidat length 6 eq {
+         candidat 3 get isarray {
+            candidat 4 get isarray {
+               candidat 5 get isarray              
+            } {
+               false
+            } ifelse
+         } {
+            false
+         } ifelse
+      } {
+         false
+      } ifelse
+   } {
+      false
+   } ifelse
+end
+} def
+
+/newplanaffine {
+   [0 0 0 [1 0 0 0 1 0] [-3 3 -2 2] [1 1]] 
+} def
+
+%% syntaxe : plantype getorigine --> x y z
+/plangetorigine {
+1 dict begin
+   /plan exch def
+   plan isplan not {
+      (Erreur : mauvais type d argument dans plangetorigine) ==
+      Error
+   } if
+   plan 0 get 
+   plan 1 get 
+   plan 2 get 
+end
+} def
+
+%% syntaxe : plantype getbase --> [u v w]
+%% ou u, v et w vecteurs de R^3
+/plangetbase {
+1 dict begin
+   /plan exch def
+   plan isplan not {
+      (Erreur : mauvais type d argument dans plangetbase) ==
+      Error
+   } if
+   plan 3 get 
+end
+} def
+
+%% syntaxe : plantype getrange --> array
+%% ou array = [xmin xmax ymin ymax]
+/plangetrange {
+1 dict begin
+   /plan exch def
+   plan isplan not {
+      (Erreur : mauvais type d argument dans plangetrange) ==
+      Error
+   } if
+   plan 4 get 
+end
+} def
+
+%% syntaxe : plantype getngrid --> array
+%% ou array = [n1 n2]
+/plangetngrid {
+1 dict begin
+   /plan exch def
+   plan isplan not {
+      (Erreur : mauvais type d argument dans plangetngrid) ==
+      Error
+   } if
+   plan 5 get 
+end
+} def
+
+%% ===================
+
+%% syntaxe : plantype x y z putorigine --> -
+/planputorigine {
+4 dict begin
+   /z exch def
+   /y exch def
+   /x exch def
+   /plan exch def
+   plan isplan not {
+      (Erreur : mauvais type d argument dans planputorigine) ==
+      Error
+   } if
+   plan 0 x put 
+   plan 1 y put 
+   plan 2 z put 
+end
+} def
+
+%% syntaxe : plantype [u v w] putbase --> -
+%% ou u, v et w vecteurs de R^3
+/planputbase {
+2 dict begin
+   /base exch def
+   /plan exch def
+   plan isplan not {
+      (Erreur : mauvais type d argument dans planputbase) ==
+      Error
+   } if
+   plan 3 base put 
+end
+} def
+
+%% syntaxe : plantype array putrange --> -
+%% ou array = [xmin xmax ymin ymax]
+/planputrange {
+2 dict begin
+   /table exch def
+   /plan exch def
+   plan isplan not {
+      (Erreur : mauvais type d argument dans planputrange) ==
+      Error
+   } if
+   plan 4 table put 
+end
+} def
+
+%% syntaxe : plantype array putngrid --> -
+%% ou array = [n1 n2]
+/planputngrid {
+2 dict begin
+   /table exch def
+   /plan exch def
+   plan isplan not {
+      (Erreur : mauvais type d argument dans planputngrid) ==
+      quit
+   } if
+   plan 5 table put 
+end
+} def
+
+%% -3 3 -2 2 1. 1. newgrille
+%% drawsolid
+
+%orange
+
+%% plan : origine, base, range, ngrid
+
+%% syntaxe : plantype drawplanaffine --> -
+/drawplanaffine {
+5 dict begin
+   /plan exch def
+   plan plangetbase 
+   aload pop
+   /imK defpoint3d
+   /imJ defpoint3d
+   /imI defpoint3d
+   newpath
+      plan plangetrange plan plangetngrid aload pop  quadrillagexOy_
+      plan plangetorigine [imI imK] false planprojpath
+   stroke
+end
+} def
+
+
+%% %% syntaxe : [a b c d] (x0 y0 z0) alpha defeqplanaffine --> plantype
+%% %% plan defini par l equation ax+by+cz+d=0, 
+%% %% rotation de alpha autour de la normale (alpha est optionnel)
+%% %% origine (x0, y0, z0). l origine est optionnelle
+%% /defeqplanaffine {
+%% 5 dict begin
+%%    dup isarray {
+%%       /alpha 0 def
+%%    } {
+%%       dup isstring {
+%%          /alpha 0 def
+%%       } {
+%%          /alpha exch def
+%%       } ifelse
+%%    } ifelse
+%%    dup isstring {
+%%       cvx /origine exch def
+%%    } if
+%%    /table exch def
+%%    table length 4 ne {
+%%       (Erreur : mauvais type d argument dans defeqplanaffine) ==
+%%       Error
+%%    } if
+%%    table 0 get /a exch def
+%%    table 1 get /b exch def
+%%    table 2 get /c exch def
+%%    table 3 get /d exch def
+%%    /resultat newplanaffine def
+%%    [a b c alpha] normalvect_to_orthobase
+%%    /imK defpoint3d
+%%    /imJ defpoint3d
+%%    /imI defpoint3d
+%%    resultat [imI imJ imK] planputbase
+%%    currentdict /origine known {
+%%       origine /z exch def /y exch def /x exch def
+%%       a x mul b y mul add c z mul add d add 0 ne {
+%%          (Erreur : mauvaise origine dans defeqplanaffine) ==
+%%          Error
+%%       } if
+%%       resultat origine planputorigine
+%%    } {
+%%       c 0 ne {
+%%          resultat 0 0 d neg c div planputorigine
+%%       } {
+%%          a 0 ne {
+%%             resultat d neg a div 0 0 planputorigine
+%%          } {
+%%             resultat 0 d neg b div 0 planputorigine
+%%          } ifelse
+%%       } ifelse
+%%    } ifelse
+%%    resultat
+%% end
+%% } def
+
+%% /explan [0 0 0 [1 0 0 0 1 0 0 0 1] [-3 3 -2 2] [1 1] ] def
+%% explan drawplanaffine
+%% noir
+%% /explan [0 0 2 [1 0 0 0 1 0 0 0 1] [-3 3 -2 2] [1 .5] ] def
+%% explan drawplanaffine
+
+%% orange
+%% [0 0 1 -2] defeqplanaffine
+%% drawplanaffine
+%% noir
+%% [0 0 1 0] defeqplanaffine
+%% drawplanaffine
+%% bleu
+%% [1 1 1 0] (1 -1 0) defeqplanaffine
+%% drawplanaffine
+%% 
+
+/dessinebase {
+4 dict begin
+gsave
+   /V3 defpoint3d
+   /V2 defpoint3d
+   /V1 defpoint3d
+   /M0 defpoint3d
+   rouge
+   V3 newvecteur 
+   {M0 translatepoint3d} solidtransform
+   drawsolid**
+   bleu
+   V2 newvecteur 
+   {M0 translatepoint3d} solidtransform
+   drawsolid**
+   orange
+   V1 newvecteur 
+   {M0 translatepoint3d} solidtransform
+   drawsolid**
+grestore
+end
+} def
+
+%% syntaxe : solid i solidface2eqplan --> [a b c d]
+%% equation cartesienne de la face d'indice i du solide solid
+/solidface2eqplan {
+8 dict begin
+   /i exch def
+   /solid exch def
+   solid i solidnormaleface
+   /c exch def
+   /b exch def
+   /a exch def
+   solid 0 i solidgetsommetface
+   /z exch def
+   /y exch def
+   /x exch def
+   [a b c a x mul b y mul add c z mul add neg]
+end
+} def
+
+
+%% syntaxe : plantype newplan --> solid
+/newplan {
+5 dict begin
+   /lepl@n exch def
+   lepl@n plangetbase /@base exch def
+   @base 0 getp3d /@U defpoint3d
+   @base 1 getp3d /@V defpoint3d
+   lepl@n plangetorigine /@M defpoint3d
+   lepl@n plangetrange /@range exch def
+   lepl@n plangetngrid /@ngrid exch def
+   /@F {
+   2 dict begin
+      /@y exch def
+      /@x exch def
+      @U @x mulv3d
+      @V @y mulv3d
+      addv3d
+      @M addv3d
+   end
+   } def
+   @range aload pop @ngrid {@F} newsurfaceparametree
+end
+} def
+
+%% syntaxe : M eqplan --> real
+%% image de M par la fonction definie par l equation eqplan
+/pointeqplan {
+8 dict begin
+   /eqplan exch def
+   /@z exch def
+   /@y exch def
+   /@x exch def
+   /@a eqplan 0 get def
+   /@b eqplan 1 get def
+   /@c eqplan 2 get def
+   /@d eqplan 3 get def
+   @a @x mul @b @y mul add @c @z mul add @d add
+end
+} def
+
+%% syntaxe : [a b c d] (x0 y0 z0) alpha defeqplanaffine --> plantype
+%% plan defini par l equation ax+by+cz+d=0, 
+%% rotation de alpha autour de la normale (alpha est optionnel)
+%% origine (x0, y0, z0). l origine est optionnelle
+/eq2plan {
+5 dict begin
+   dup isarray {
+      /alpha 0 def
+   } {
+      dup isstring {
+         /alpha 0 def
+      } {
+         /alpha exch def
+      } ifelse
+   } ifelse
+   dup isstring {
+      cvx /origine exch def
+   } if
+   /table exch def
+   table length 4 ne {
+      (Erreur : mauvais type d argument dans eq2plan) ==
+      quit
+   } if
+   table 0 get /a exch def
+   table 1 get /b exch def
+   table 2 get /c exch def
+   table 3 get /d exch def
+   /resultat newplanaffine def
+   [a b c alpha] normalvect_to_orthobase
+   /imK defpoint3d
+   /imJ defpoint3d
+   /imI defpoint3d
+   resultat [imI imJ imK] planputbase
+   currentdict /origine known {
+      origine /z exch def /y exch def /x exch def
+      a x mul b y mul add c z mul add d add 0 ne {
+         (Erreur : mauvaise origine dans eq2plan) ==
+         quit
+      } if
+      resultat origine planputorigine
+   } {
+      c 0 ne {
+         resultat 0 0 d neg c div planputorigine
+      } {
+         a 0 ne {
+            resultat d neg a div 0 0 planputorigine
+         } {
+            b 0 ne {
+               resultat 0 d neg b div 0 planputorigine
+            } {
+               (Error dans eq2plan : (a,b,c) = (0,0,0)) ==
+            } ifelse
+         } ifelse
+      } ifelse
+   } ifelse
+   resultat
+end
+} def
+
+/points2eqplan {
+10 dict begin
+   /C defpoint3d
+   /B defpoint3d
+   /A defpoint3d
+   A B vecteur3d
+   A C vecteur3d
+   vectprod3d
+   normalize3d
+   /c exch def
+   /b exch def
+   /a exch def
+   A
+   /zA exch def
+   /yA exch def
+   /xA exch def
+   [a b c a xA mul b yA mul add c zA mul add neg]
+end
+} def
+
+%% /monplan 
+%% %[0 0 -2 [1 0 0 0 1 0 0 0 1] [-3 3 -2 2] [1. 1.]]
+%% [0 0 1 1] 30 eq2plan
+%% def
+%% 
+%% [0 0 1 -2] eq2plan newplan
+%% dup (blanc) outputcolors
+%% monplan newplan
+%% dup (blanc) outputcolors
+%% solidfuz
+%% drawsolid**
+%% monplan plangetorigine
+%% monplan plangetbase aload pop dessinebase
+
+%% syntaxe : x0 y0 z0 [normalvect] norm2plan
+/norm2plan {
+9 dict begin
+   normalvect_to_orthobase
+   /imK defpoint3d
+   /imJ defpoint3d
+   /imI defpoint3d
+   imK
+   /c exch def
+   /b exch def
+   /a exch def
+   /z0 exch def
+   /y0 exch def
+   /x0 exch def
+   [a b c a x0 mul b y0 mul add c z0 mul add neg] eq2plan
+   dup x0 y0 z0 planputorigine
+   dup [imI imJ] planputbase
+end
+} def
+
+%% syntaxe : plantype planxmarks
+/planxmarks {
+5 dict begin
+   /leplan exch def
+   leplan plangetrange aload pop
+   /ymax exch def
+   /ymin exch def
+   /xmax exch def
+   /xmin exch def
+   xmin truncate cvi xmkstep xmax truncate cvi {
+      dup 0 ne {
+         dup chaine cvs exch 0 leplan dctextp3d
+      } {
+         pop
+      } ifelse
+   } for
+end
+} def
+
+%% syntaxe : plantype planymarks
+/planymarks {
+5 dict begin
+   /leplan exch def
+   leplan plangetrange aload pop
+   /ymax exch def
+   /ymin exch def
+   /xmax exch def
+   /xmin exch def
+   ymin truncate cvi ymkstep ymax truncate cvi {
+      dup 0 ne {
+         dup chaine cvs exch 0 exch leplan cltextp3d
+      } {
+         pop
+      } ifelse
+   } for
+end
+} def
+
+%% syntaxe : plantype planmarks
+/planmarks {
+   dup planxmarks planymarks
+} def
+
+%% bleu
+%% [-3 3 -2 2] quadrillagexOy_
+%% stroke
+%% noir
+
+%% syntaxe : [xmin xmax ymin ymax] dx dy quadrillagexOy_
+/quadrillagexOy_ {
+4 dict begin
+   dup isarray {
+      /dx 1 def
+      /dy 1 def
+   } {
+      /dy exch def
+      dup isarray {
+         /dx dy def
+      } {
+         /dx exch def
+      } ifelse
+   } ifelse
+   /table exch def
+   table 0 get /xmin exch def
+   table 1 get /xmax exch def
+   table 2 get /ymin exch def
+   table 3 get /ymax exch def
+   ymin dy ymax {
+      /y exch def
+      xmin y smoveto
+      xmax y slineto
+   } for
+   xmin dx xmax {
+      /x exch def
+      x ymin smoveto
+      x ymax slineto
+   } for
+end
+} def
+
+%%%%% ### fin insertion ###
+
+
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-%%%%     operations sur des solides particuliers        %%%%
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+%%%%% ### solidplansection ###
+%% syntaxe : M eqplan --> real
+%% image de M par la fonction definie par l equation eqplan
+/pointeqplan {
+8 dict begin
+   /@qplan exch def
+   /@z exch def
+   /@y exch def
+   /@x exch def
+   /@a @qplan 0 get def
+   /@b @qplan 1 get def
+   /@c @qplan 2 get def
+   /@d @qplan 3 get def
+   @a @x mul @b @y mul add @c @z mul add @d add 
+end
+} def
+
+%% syntaxe : A B eqplan segment_inter_plan --> array true ou false
+%% array contient 1 point M si [AB] inter plan = {M}
+%% array contient les 2 points A et B si [AB] inter plan = [AB]
+/segment_inter_plan {
+4 dict begin
+   /plan exch def
+   plan aload pop
+   /d exch def
+   /c exch def
+   /b exch def
+   /a exch def
+   /B defpoint3d
+   /A defpoint3d
+   A 
+   /zA exch def
+   /yA exch def
+   /xA exch def
+   B 
+   /zB exch def
+   /yB exch def
+   /xB exch def
+   /imA a xA mul b yA mul add c zA mul add d add def
+   /imB a xB mul b yB mul add c zB mul add d add def
+   imA imB mul dup 0 gt {
+      %% pas d intersection
+      pop
+      false
+   } {
+      0 eq {
+         %% intersection en A ou en B
+         [ 
+            imA 0 eq {A} if 
+            imB 0 eq {B} if 
+         ] true
+      } {
+         %% intersection entre A et B
+         /k 
+            imA neg
+            xB xA sub a mul
+            yB yA sub b mul add
+            zB zA sub c mul add
+            dup 0 eq {
+               (Error dans segment_inter_plan) ==
+               quit
+            } if
+            div
+         def
+         [
+            A B vecteur3d
+            k mulv3d
+            A addv3d
+         ] true
+      } ifelse
+   } ifelse
+end
+} def
+
+%% syntaxe : solid i solidface2eqplan --> [a b c d]
+%% equation cartesienne de la face d'indice i du solide solid
+/solidface2eqplan {
+8 dict begin
+   /i exch def
+   /solid exch def
+   solid i solidnormaleface
+   /c exch def
+   /b exch def
+   /a exch def
+   solid 0 i solidgetsommetface
+   /z exch def
+   /y exch def
+   /x exch def
+   [a b c a x mul b y mul add c z mul add neg]
+end
+} def
+
+%% syntaxe : solid eqplan solidplansection --> solid2
+/solidplansection {
+10 dict begin
+   /eqplan exch def
+   dupsolid /result exch def
+   /solid exch def
+   /aenlever [] def
+%(nombre faces) ==
+%solid solidnombrefaces ==
+   0 1 solid solidnombrefaces 1 sub {
+      /i exch def
+      /lacouleur solid i solidgetfcolor def
+%(NOUVELLE FACE) == i ==
+      /F solid i solidgetface def %% table des indices
+      /n F length def %% nb d'aretes
+      /k1 -1 def
+      /k2 -1 def
+      /k3 -1 def
+      /k4 -1 def
+      /k3a -3 def
+      /k4a -3 def
+      %% pour chaque arete [AB]
+      0 1 n 1 sub {
+         /j exch def
+%(nvell arete) == j ==
+         solid j i solidgetsommetface /A defpoint3d
+         solid j 1 add n mod i solidgetsommetface /B defpoint3d
+         %% y a-t-il intersection
+         A B eqplan segment_inter_plan {
+%(inter) ==
+             dup length 6 eq {
+               /k1 -1 def
+               /k2 -1 def
+               /k3 -1 def
+               /k4 -1 def
+               /k3a -1 def
+               /k4a -1 def
+               exit %% c est deja scinde
+            } if
+            k1 0 lt {
+            %% 1ere intersection de la face
+               /k1 j def
+               result exch aload pop solidaddsommet
+               /k1a exch def
+            } {
+               k2 0 lt {
+               %% 2eme intersection de la face
+                  /k2 j def
+                  result exch aload pop solidaddsommet
+                  /k2a exch def
+               } {
+                  k3 0 lt {
+                  %% 3eme intersection de la face
+                     /k3 j def
+                     result exch aload pop solidaddsommet
+                     /k3a exch def
+                  } {
+                  %% 4eme intersection de la face
+                     /k4 j def
+                     result exch aload pop solidaddsommet
+                     /k4a exch def
+                  } ifelse
+               } ifelse
+            } ifelse
+         } if
+      } for
+      %% y a-t-il une coupe ?
+      k1 0 ge {
+%% (coupe) ==
+%% (n) == n ==
+%% k1 == k2 == k3 == k4 ==
+%% (a) ==
+%% k1a == k2a == k3a == k4a ==
+         k1a k2a eq k3 0 lt and {
+         } {
+         k1a k2a eq k3 0 ge and {
+            /k2a k3a def
+            /k2 k3 def
+         } if
+         k1a k3a eq k4 0 ge and {
+            /k2a k4a def
+            /k2 k4 def
+         } if
+            [
+               k1a F k1 1 add n mod get ne {
+                  k1a
+               } if
+               k1 1 add n mod 1 k2 {F exch get} for
+               k2a F k2 get ne {
+                  k2a
+               } if
+            ]
+            result exch lacouleur solidaddface
+            [
+               k2a F k2 1 add n mod get ne {
+                  k2a
+               } if
+               k2 1 add n ne {
+                  k2 1 add n mod 1 n 1 sub {F exch get} for
+               } if
+               0 1 k1 {F exch get} for
+               k1a F k1 get ne {
+                  k1a
+               } if
+            ]
+            result exch lacouleur solidaddface
+            /aenlever [aenlever aload pop i] store
+         } ifelse
+      } if
+   } for
+   result aenlever solidrmfaces
+   result 
+end    
+} def
+
+%% syntaxe : solid eqplan solideqplansepare --> solid1 solid2
+/solideqplansepare {
+10 dict begin
+   /eqplan exch def
+   eqplan solidplansection
+   /solid exch def
+   /n solid solidnombrefaces def
+   /lesneg [] def
+   /lespos [] def
+   0 1 n 1 sub {
+      /i exch def
+      solid i solidcentreface /G defpoint3d
+      G eqplan pointeqplan dup 0 gt {
+         pop
+         /lespos [lespos aload pop i] store
+      } {
+         0 lt {
+            /lesneg [lesneg aload pop i] store
+         } {
+            /lesneg [lesneg aload pop i] store
+            /lespos [lespos aload pop i] store
+         } ifelse
+      } ifelse
+   } for
+   solid
+   dupsolid dup lesneg solidrmfaces 
+   /result1 exch def
+   dupsolid dup lespos solidrmfaces
+   /result2 exch def
+   pop
+   /sommetspos [] def
+   /sommetsneg [] def
+   %% pour chaque face du cote negatif
+   0 1 lesneg length 1 sub {
+      lesneg exch get /i exch def
+      /F solid i solidgetface def
+      %% pour chaque sommet de cette face
+      0 1 F length 1 sub {
+         /j exch def
+         /sommet F j get def
+         %% si le sommet n'est pas encore note
+         sommet sommetsneg in not {
+            %% et s'il est isole, on peut l'enlever
+            result1 sommet solidsommetsadjsommet length 0 eq {
+               /sommetsneg [sommetsneg aload pop sommet] store
+            } if
+         } {
+            pop
+         } ifelse
+      } for
+   } for
+   sommetsneg bubblesort reverse {result1 exch solidrmsommet} apply
+
+   %% pour chaque face du cote positif
+   0 1 lespos length 1 sub {
+      lespos exch get /i exch def
+      /F solid i solidgetface def
+      %% pour chaque sommet de cette face
+      0 1 F length 1 sub {
+         /j exch def
+         /sommet F j get def
+         %% si le sommet n'est pas encore note
+         sommet sommetspos in not {
+            %% et s'il est isole, on peut l'enlever
+            result2 sommet solidsommetsadjsommet length 0 eq {
+               /sommetspos [sommetspos aload pop sommet] store
+            } if
+         } {
+            pop
+         } ifelse
+      } for
+   } for
+   sommetspos bubblesort reverse {result2 exch solidrmsommet} apply
+
+   result1 result2
+end
+} def
-
-%%%%% ### dualpolyedre ###
@@ -4881 +5759 @@
-   [[/nx /ny] [6 6] [6 8] [10 8] [16 12] [16 36]] gestionsolidmode
-
+   nx isinteger not {
+       %% alors nx est un dx
+       /nx xmax xmin sub nx div cvi store
+   } if
+   ny isinteger not {
+       %% alors ny est un dy
+       /ny ymax ymin sub ny div cvi store
+   } if
+   /dy ymax ymin sub ny div def %% le pas sur y
+   /dx xmax xmin sub nx div def %% le pas sur x
+
-   %% ny = nb de meridiens
-   %% nx = nb d horizontales
-
-   /r exch def
-   nx isinteger not {
-       %% alors nx est un dx
-       /nx r nx div cvi store
-   } if
-   ny isinteger not {
-       %% alors ny est un dy
-       /ny 360 ny div cvi store
-   } if
-   /dy r ny div def %% le pas sur y
-   /dx r nx div def %% le pas sur x
-
-   /F [
-       %% 1er etage
@@ -6200 +7078 @@
-   /lafonction str cvx def
-   /laderivee str (') append cvx def
-
+%%
-   /tmax exch def
-   /tmin exch def
@@ -6213 +7091 @@
-      %% definition du repere de Frenet (k1, k2, k3) au point f(a)
-      a0 lafonction /M defpoint3d
-
-
+
+
+
+
+
+
+
+
+
+
-      k1 baseplannormal /K3 defpoint3d /K2 defpoint3d
-%      M k1 k2 K3 dessinebase
-%      a0 laderivee2nd normalize3d /k2 defpoint3d
+
+      %% projete orthogonal du dernier rayon sur le plan actuel
+      %% (normal a la vitesse)
-      K2 tub@dernierk2 aload pop K2 scalprod3d mulv3d 
-      K3 tub@dernierk2 aload pop K3 scalprod3d mulv3d addv3d /k2 defpoint3d
+%      M k1 K2 K3 dessinebase
-      k1 norme3d 0 eq {
-         tub@dernierk1 aload pop /k1 defpoint3d
@@ -6237 +7126 @@
-      } ifelse
-      k3 k1 vectprod3d normalize3d /k2 defpoint3d
+%%      M k1 k2 k3 dessinebase
-      /tub@dernierk2 [k2] store
-      /@n 360 N div def %% le pas angulaire

trunk/tests/index.txt

@@ -8 +8 @@
-test_08 : ecriture de fichiers dat
-test_09 : lecture de fichier dat
+test_10 : section et separation de solide
+test_11 : equations de plan