WebSVN – wimsdev – Diff – /trunk/wims/public_html/scripts/slib/geo2D/hyptiling

     (nombre de sommets) pour chacun des sommets d'une face (pavage de Catalan);\
     </li><li>\
     cas de type 2: transitifs sur l'ensemble des sommets et le premier paramètre est\
     la configuration de sommets [d_1, ... , d_n],\
     c'est-à-dire la suite des types de faces (nombre de sommets) ayant un sommet commun.\
-  </li></ul>\
+  </li><li>Le cas de type 3 montre le pavage et son dual sur le même dessin.</li></ul>\
   Dans le premier cas, la tuile de base est un polygone tangentiel\
     (polygone convexe admettant un cercle inscrit) dont les angles\
     sont la suite des 2*pi/d_1, ..., 2*pi/d_n.\
   Dans le deuxième cas, le pavage est dual du pavage de type 1 de configuration de faces\
   [d_1, ... , d_n].<br>\
   le polygone est un polygone à d_1 côtés égaux et d_1 angles égaux à 2*pi/d_2.<br>\
    Ces pavages existent si sum_i=1^n 1/d_i<(n-2)/2 et si les d_i sont pairs.\
   Il semble qu'on puisse affaiblir cette deuxième condition en supposant que\
   pour tout i et j, d_i est pair ou d_{i-j}=d_{i+j}.<br>\
   Le pavage (qui est théoriquement infini) est limité au cercle de rayon 1-eps (eps>0.04).\
-  Cependant, le nombre de polygones dessinés est limité à 500.\
+  Cependant, le nombre de polygones dessinés est limité à 400.\
   Si eps est un entier, seules les esp premières tuiles sont dessinées.\
   <br>\
   L'option depl change le centrage initial du polygone de base: entre autre,\
   depl=1: le centre du polygone de départ est 0,\
-  depl=0 le cercle est tangent au point 0 au polygone.
+  depl=0 :le cercle est tangent au point 0 au polygone.
   type=1: catalan
   type=2: dual de catalan
   type=3: les deux
   !!!avec |a|>|b|; elle envoie 0 sur b/\bar{a} et le point 1 de l'horizon sur\
   !!!a/\bar{a}.
 slib_example=[4,8,10],0.1,html type=1\
 [4,8,10],0.01,html color=blue depl=1\
 [4,8,10],0.01,html color=blue type=2\
-[8,8,4,8],0.01,html color=blue type=1\
+[4,8,10],0.01,html color=[blue,green] type=3\
-[8,8,4,8],0.01,html color=blue type=2\
+[8,8,4,8],0.01,html type=2\
-[4,8,6,4,6,6],0.01,html\
+[8,8,4,8],0.01,html \
+[4,8,6,4,6,6],0.01,html color=blue type=1\
+[4,8,6,4,6,6],0.01,html color=green type=2\
+[4,8,6,4,6,6],0.01,html color=[blue,green] type=3\
-[34,5,34],0.1,html fill=[1]\
+[34,5,34],0.1,html color=green fill=[1]\
 [6,8,10],0.1,html color=[blue,grey] fill\
 [6,8,10],0.1,html fill=[1,5,7]\
 [4,8,9,8,6],0.01,html
 !exit
 :proc
-!reset slib_Test slib_type slib_bound slib_size slib_color slib_dessin slib_dessin1 slib_dessin2 slib_dessin3 slib_depl slib_fill slib_dual
+!reset slib_Test slib_type slib_type0 slib_bound slib_size slib_color slib_dessin slib_dessin1 slib_dessin2 slib_dessin3 slib_depl slib_fill slib_dual
 !distribute items $wims_read_parm into slib_data,slib_bound,slib_option
 !set slib_gpprog=hyptiling
 !if $slib_bound<0.01
   !!reset slib_bound
 !endif
 !set slib_size=!getopt size in $slib_option
 !default slib_size=600
 !set slib_color=!getopt color in $slib_option
 !set slib_color=!declosing $slib_color
 !set slib_type=!getopt type in $slib_option
-!default slib_type=1
-!if $slib_type=1
-  !if fill iswordof $slib_option
-    !set slib_pol=fhyppolygon
-    !default slib_color=black,white
-  !else
-    !set slib_pol=hyppolygon
-    !default slib_color=black
+!default slib_color=black
-  !endif
-!endif
-!if $slib_type=2
-  !set slib_pol=hyprays
-  !default slib_color=black,black
-!endif
 !set slib_color=$slib_color,$(slib_color[1])
 !set slib_color=$(slib_color[1,2])
 !set slib_fill=!getopt fill in $slib_option
 !set slib_fill=!declosing $slib_fill
-!set slib_type=!getopt type in $slib_option
 !set slib_depl=!getopt depl in $slib_option
 !default slib_depl=0
 !default slib_type=1
 !set slib_pqr=!declosing $slib_data
 !set slib_pqr=!exec pari [$slib_pqr]
 !if regular notin $slib_option
   !set slib_n=!itemcnt $slib_pqr
   !set slib_testexist=!exec pari a=[$slib_pqr];test=vecsum(vector($slib_n,i,1/a[i]))>=($slib_n-2)/2;a=concat(a,a[1]);for(i=2,$slib_n,if(a[i]%2==1, test+=a[i-1]!=a[i+1])); test
   !if $slib_testexist>0
     !set slib_out=error $slib_testexist>0
     !exit
   !endif
 !endif
 !if regular isin $slib_option
   !if $[1/$(slib_pqr[1])+1/$(slib_pqr[2])]>=1/2
     slib_out=Error regular
     !exit
   !readproc gp/$slib_gpprog.gp
   !set slib_hyptiling=!exec pari tikz=0;$(slib_header_$slib_gpprog);
 !!endif
 slib_depl:0 pas au centre
 au centre
+!default slib_type=1
 -> 1,2
 -> 3,4
-!set slib_hyptiling=!exec pari wwww=wims_catalan([$slib_pqr],$slib_bound,$slib_depl,$slib_type);default(realprecision,3);wwww
+- > le faire en deux fois.
-!set slib_data=!trim $(slib_hyptiling[2*$slib_type-1,2*$slib_type])
+!set slib_type0=$slib_type
+!if $slib_type=3
-!! depends only on $slib_data -- could be an "offhypdraw"
+  !set slib_type0=1,2
+!else
-!set slib_data=!declosing $slib_data
+  !set slib_type0=$slib_type
+!endif
+!set slib_hyptiling=!exec pari wwww=wims_catalan([$slib_pqr],$slib_bound,$slib_depl,$(slib_type0[1]));default(realprecision,3);wwww
-!distribute items $slib_data into slib_som,slib_face
+!set slib_data=!trim $(slib_hyptiling[2*$(slib_type0[1])-1,2*$(slib_type0[1])])
-slib_som=!declosing $slib_som
-slib_face=!declosing $slib_face
-slib_facecnt=!itemcnt $slib_face
-!for slib_k=1 to $slib_facecnt
-  !reset slib_tmp slib_test slib_col
-  !set slib_ff=!declosing $(slib_face[$slib_k])
-  !if $slib_type=2
+!if $(slib_type0[1])=1
-    !set slib_ffcnt=!itemcnt $slib_ff
+  !if fill iswordof $slib_option
-    !set slib_ff=$slib_ffcnt,$slib_ff,1
+    !set slib_pol=fhyppolygon
-  !endif
-  !set slib_ffcnt=!itemcnt $slib_ff
-  !set slib_col=$(slib_color[1+$(slib_ff[$[$(slib_ff[1])+2]])])
-  !if $slib_k isin $slib_fill
-    !set slib_pol1=f$slib_pol
+    !default slib_color=black,white
   !else
-    !set slib_pol1=$slib_pol
+    !set slib_pol=hyppolygon
-  !endif
-  !set slib_tmp=$slib_pol1 $(slib_col)
+    !set slib_col=$(slib_color[1])
-  !set slib_tmp=!append item $(slib_som[$(slib_ff[2..$(slib_ff[1])+1]);]) to $slib_tmp
-  !set slib_tmp=!replace internal ; by , in $slib_tmp
-  !if $slib_k<500
-    !set slib_dessin=!append line $slib_tmp to $slib_dessin
-  !else
-    !goto out
-    !if $slib_k<800
-      !set slib_dessin1=!append line $slib_tmp to $slib_dessin1
-    !else
-      !if $slib_k<2400
-        !set slib_dessin2=!append line $slib_tmp to $slib_dessin2
-      !else
-        !set slib_dessin3=!append line $slib_tmp to $slib_dessin3
-      !endif
-    !endif
   !endif
+!else
+  !set slib_pol=hyprays
+  !set slib_col=$(slib_color[1])
-:end
+!endif
+!readproc slib/geo2D/offhyptiling
+!if $slib_type=3
+  !set slib_type0=2
+  !set slib_hyptiling=!exec pari wwww=wims_catalan([$slib_pqr],$slib_bound,$slib_depl,$(slib_type0));default(realprecision,3);wwww
+  !set slib_data=!trim $(slib_hyptiling[2*$(slib_type0)-1,2*$(slib_type0)])
+  !set slib_pol=hyprays
+  !set slib_col=$(slib_color[2])
+  !readproc slib/geo2D/offhyptiling
-!next
+!endif
-:out
+!! depends only on $slib_data -- could be an "offhypdraw"
 !set slib_range=-1,1
 !set slib_out=xrange $slib_range\
 yrange $slib_range\
 circles black,0,0,1\
 $slib_dessin\
-$slib_dessin1
+$slib_dessin1\
 $slib_dessin2\
 $slib_dessin3
 !if url iswordof $slib_option or html iswordof $slib_option
   !set insdraw_size=$slib_size,$slib_size
   !insdraw $slib_out
   !set slib_out=$ins_url
   !if html iswordof $slib_option

Subversion Repositories wimsdev

(root)/trunk/wims/public_html/scripts/slib/geo2D/hyptiling – Rev 18529 → 18530