Changeset 29

Timestamp:

05/25/08 10:09:46 (8 months ago)

Author:

jpv

Message:

Modif de la section par un plan pour ajouter la face de coupe

Files:

• trunk/solides.pro (modified) (15 diffs)
• trunk/tests/test_10.pdf (modified) (previous)
• trunk/tests/test_10.tex (modified) (1 diff)

trunk/solides.pro

@@ -1 +1 @@
-%!
-% PostScript prologue for pst-solides3d.tex.
-5, 2008/05/18
+6, 2008/05/25
-%
-%% COPYRIGHT 2008 by Jean-Paul Vignault, Manuel Luque
@@ -2852 +2852 @@
-} def
-
+%% syntaxe : solid solidfacesreverse -> -
+/solidfacesreverse {
+5 dict begin
+   /solid exch def
+   /n solid solidnombrefaces def
+   0 1 n 1 sub {
+      /i exch def
+      /F solid i solidgetface reverse def
+      /m F length def
+      solid i [F aload pop m 0 roll ] solidputface
+   } for
+end
+} def
+
-/solidputfcolors {
-   2 exch put
@@ -2915 +2929 @@
-} def
-
-%% %% syntaxe : solid array solidputoutfaces --> -
-%% /solidputoutfaces {
-%% 4 dict begin
-%%    /facesexternes exch def
-%%    /solid exch def
-%%    /n1 facesexternes length def
-%%    /IO solid solidgetinouttable def
-%%    /facesinternes solid solidgetinfaces def
-%%    /n2 facesinternes length def
-%%    solid
-%%       [facesexternes aload pop facesinternes aload pop]
-%%       solidputfaces
-%%    IO 0 0 put
-%%    IO 1 n1 1 sub put
-%%    IO 2 n1 put
-%%    IO 3 n1 n2 add 1 sub put
-%% end
-%% } def
-
-%% syntaxe : solid array solidputoutfaces --> -
-/solidputoutfaces {
@@ -3011 +3006 @@
-} def
-
-%/solidintersectiontype -1 def
-%/solidintersectionlinewidth [1 2] def
-%/solidintersectioncolor [(bleu) (rouge)] def
-%/solidintersectionplan [[1 0 .5 0] [1 0 .5 -1]] def
-
-%% /tracelignedeniveau? false def
-%% /hauteurlignedeniveau 1 def
@@ -3021 +3011 @@
-%% /linewidthlignedeniveau 4 def
-
-/solidgrid true def
-/aretescachees true def
-/defaultsolidmode 2 def
-/solidgridOn {
-   /solidgrid true def
@@ -3336 +3323 @@
-      vecteurnormal angle3doriente
-   end} papply3d
-
-
+
+
+
-end
-} def
@@ -5486 +5474 @@
-} def
-
+%% syntaxe : array1 arrayrmdouble --> array2
+%% remplace 2 elts identiques consecutifs par 1 elt
+/arrayrmdouble {
+5 dict begin
+   /table exch def
+   /result [table 0 get] def
+   /j 0 def
+   1 1 table length 1 sub {
+      /i exch def
+      table i get
+      result j get
+      eq not {
+         /result [result aload pop table i get] store
+         /j j 1 add store
+      } if
+   } for
+   result
+end
+} def
+
-%% syntaxe : solid eqplan solidplansection --> solid2
-/solidplansection {
@@ -5493 +5501 @@
-   /solid exch def
-   /aenlever [] def
+   /indnouveauxsommets [] def
-%(nombre faces) ==
-%solid solidnombrefaces ==
@@ -5523 +5532 @@
-               /k3a -1 def
-               /k4a -1 def
+               dup 0 getp3d /A defpoint3d
+               1 getp3d /B defpoint3d
+               result A solidaddsommet /a1 exch def
+               result B solidaddsommet /a2 exch def
+               /indnouveauxsommets [
+                  indnouveauxsommets aload pop a1 a2
+               ] store
-               exit %% c est deja scinde
-            } if
@@ -5560 +5576 @@
-%% k1a == k2a == k3a == k4a ==
-         k1a k2a eq k3 0 lt and {
+            %% 1 pt d'intersection
-         } {
-         k1a k2a eq k3 0 ge and {
+            %% 2 pts d'intersection
-            /k2a k3a def
-            /k2 k3 def
-         } if
-         k1a k3a eq k4 0 ge and {
+            %% 2 pts d'intersection
-            /k2a k4a def
-            /k2 k4 def
@@ -5579 +5598 @@
-            ]
-            result exch lacouleur solidaddface
+            /indnouveauxsommets [indnouveauxsommets aload pop k1a k2a] store
-            [
-               k2a F k2 1 add n mod get ne {
@@ -5597 +5617 @@
-   } for
-   result aenlever solidrmfaces
+
+   indnouveauxsommets bubblesort arrayrmdouble
+   /indnouveauxsommets exch def
+   %% maintenant, on ajoute la face de plan de coupe
+   /nouveauxsommets [
+      0 1 indnouveauxsommets length 1 sub {
+         /i exch def
+         result indnouveauxsommets i get solidgetsommet
+      } for
+   ] def
+
+   0 0 0 eqplan pointeqplan 0 eq {
+      /ptref {0 1 1} def
+   } {
+      /ptref {0 0 0} def
+   } ifelse
+
+
+   %% restera a traiter le cas limite ou la nouvelle face existe deja
+   %% tester si max(indicestries) < nb sommets avant section
+   nouveauxsommets ptref ordonnepoints3d
+   /indicestries exch def
+   /nvelleface [
+      0 1 indicestries length 1 sub {
+         /m exch def
+         indnouveauxsommets indicestries m get get
+      } for
+   ] def
+   /F result solidgetfaces def
+   /FC result solidgetfcolors def
+   /IO result solidgetinouttable def
+   IO dup 1 get 1 add 1 exch put
+   result [nvelleface F aload pop] solidputfaces
+   result [lacouleur FC aload pop] solidputfcolors
-   result 
+   result
-end    
-} def
@@ -5608 +5663 @@
-   /solid exch def
-   /n solid solidnombrefaces def
-
-
+
+
+
+
+
+
-   0 1 n 1 sub {
-      /i exch def
@@ -5620 +5679 @@
-            /lesneg [lesneg aload pop i] store
-         } {
- 
- 
+%
+%
-         } ifelse
-      } ifelse
-   } for
-   solid
- 
+
-   /result1 exch def
-   dupsolid dup lespos solidrmfaces
-   /result2 exch def
-   pop
+   %% on rajoute maintenant la face du plan de coupe
+   result1 facecoupe couleurfacecoupe solidaddface
+   %% et on verifie l'orientation
+   result1 dup solidnombrefaces 1 sub solidnormaleface 
+   result1 dup solidnombrefaces 1 sub solidcentreface addv3d
+   eqplan pointeqplan 0 gt {
+      %% l'orientation est mauvaise
+      result1 dup solidnombrefaces 1 sub solidrmface 
+      result2 facecoupe couleurfacecoupe solidaddface
+      result1 facecoupe reverse couleurfacecoupe solidaddface
+   } {
+      %% l'orientation est ok
+      result2 facecoupe reverse solidaddface
+   } ifelse
+   
-   /sommetspos [] def
-   /sommetsneg [] def
@@ -5674 +5748 @@
-   } for
-   sommetspos bubblesort reverse {result2 exch solidrmsommet} apply
-
+
-   result1 result2
-end

trunk/tests/test_10.tex

@@ -26 +26 @@
-\psframe(-5,-4)(6,9)
-\codejps{
-{
--2 1 2 [5 12] newcylindre
+dup (jaune) (vert) inoutputcolors
-[1 1 1 0] solideqplansepare
-
-
-
+
-exch
-{0 0 60 rotateOpoint3d} solidtransform
-{0 2 0 translatepoint3d} solidtransform
-
-
-
+
-drawsolid**
-}