=================================================================== RCS file: /home/cvs/OpenXM/src/kan96xx/Doc/gfan.sm1,v retrieving revision 1.13 retrieving revision 1.18 diff -u -p -r1.13 -r1.18 --- OpenXM/src/kan96xx/Doc/gfan.sm1 2009/08/26 04:54:17 1.13 +++ OpenXM/src/kan96xx/Doc/gfan.sm1 2009/12/11 02:09:09 1.18 @@ -1,9 +1,25 @@ -% $OpenXM: OpenXM/src/kan96xx/Doc/gfan.sm1,v 1.12 2005/07/07 07:53:37 takayama Exp $ +% $OpenXM: OpenXM/src/kan96xx/Doc/gfan.sm1,v 1.17 2009/09/04 11:13:11 takayama Exp $ % cp cone.sm1 $OpenXM_HOME/src/kan96xx/Doc/gfan.sm1 -% $Id: gfan.sm1,v 1.13 2009/08/26 04:54:17 takayama Exp $ +% $Id: gfan.sm1,v 1.18 2009/12/11 02:09:09 takayama Exp $ % iso-2022-jp %%Ref: @s/2004/08/21-note.pdf +%% gfan.sm1 works only for polymake 2.0 Use webservice of 2.0. +[(gfan) +[ + (gfan.sm1 is a package to compute global and local Grobner fans.) + (See R.Bahloul and N.Takayama, arxiv, math.AG/0412044 and references as to algorithms.) + (At the beginning of the source code gfan.sm1, there are sample inputs cone.sample and cone.sample2.) + ( ) + (gfan.sm1 works only with polymake 2.0. We provide a web service of computing ) + (with polymake 2.0. /@@@polymake.web 1 def is set by default in gfan.sm1.) + (See changelog-ja.tex as to details on the difference between 2.0 and later versions.) + ( ) + ( cone.Wt cone.Lpt {vertices in the output} are weights on the rays of the Grobner cone.) + ( cone.L gives a basis of the linearity space.) +] +] putUsages + %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %% Two examples are given below to get a global Grobner fan and %% a local Grobner fan ; cone.sample and cone.sample2 @@ -250,9 +266,13 @@ dhcone.printGrobnerFan % If you use the cgi/polymake on the net, then uncomment out. %/doPolymake {doPolymake.OoHG} def (Using doPolymake.OoHG ) message %/polymake.start {polymake.start.OoHG} def (Using polymake.start.OoHG ) message +/@@@polymake.web 1 def %% Choose it automatically. -[(which) (polymake)] oxshell tag 0 eq { - (Polymake is not installed in this system.) message +[(which) (polymake)] oxshell tag 0 eq +@@@polymake.web 1 eq +or +{ + (Polymake is not installed in this system or @@@polymake.web is set.) message /doPolymake {doPolymake.OoHG} def (Using doPolymake.OoHG ) message /polymake.start {polymake.start.OoHG} def @@ -1020,7 +1040,7 @@ def %> /getConeInfo { /arg1 set - [/ww /g /ceq /ceq2 /cdim /mmc /mmL /rr /ineq /ppt] pushVariables + [/ww /g /ceq /ceq2 /cdim /mmc /mmL /rr /ineq /ppt /rr0 /mm0 /mm1] pushVariables [ /ceq arg1 def ceq pruneZeroVector /ceq set @@ -1029,13 +1049,22 @@ def (Monomial ideal is not accepted as an input.) cone_ir_input } { } ifelse + /mm1 + ( Use [(keep_tmp_files) 1] oxshell to check the input to polymake2tfb. See /tmp or $TMP ) + def + ceq genPo2 /ceq2 set % ceq2 は polymake.data(polymake.INEQUALITIES(...)) 形式 % polymake で ceq2 の次元の計算. /getConeInfo.ceq ceq def /getConeInfo.ceq2 ceq2 def cone.debug { (Calling polymake DIM.) message } { } ifelse - [(DIM) ceq2] doPolymake 1 get /rr set + [(DIM) ceq2] doPolymake /rr0 set + % rr0 2 get message + rr0 2 get 1 get 0 get /mm0 set + mm0 length 0 eq { } + { [mm0 mm1] cat error } ifelse + rr0 1 get /rr set cone.debug {(Done.) message } { } ifelse % test5 には次のコメントとりさる. 上の行をコメントアウト. % test5.data tfbToTree /rr set @@ -1050,13 +1079,24 @@ def % FACETS を持っていないなら再度計算する. % POINTED, NOT__POINTED も得られる cone.debug { (Calling polymake FACETS.) message } { } ifelse - [(FACETS) ceq2] doPolymake 1 get /rr set + [(FACETS) ceq2] doPolymake /rr0 set + + % rr0 2 get message + rr0 2 get 1 get 0 get /mm0 set + mm0 length 0 eq { } + { [mm0 mm1] cat error } ifelse + + rr0 1 get /rr set cone.debug { (Done.) message } { } ifelse } { } ifelse rr (VERTICES) getNode tag 0 eq { (internal error: VERTICES is not found.) error - } { } ifelse + } { + rr (VERTICES) getNode + (UNDEF) getNode tag 0 eq { } + { (internal error: VERTICES is UNDEF. See rr. Set /@@@polymake.web 1 def) error } ifelse + } ifelse /cone.getConeInfo.rr1 rr def @@ -1626,7 +1666,7 @@ def /vlist arg1 def /wlist arg2 def wlist length vlist length eq { - } { (cone_wtowv: length of the argument must be the same.) error} ifelse + } { (cone_wtowv: length of the argument must be the same. Please check the values of cone.vlist cone.vv cone.type parametrizeWeightSpace) error} ifelse wlist to_int32 /wlist set [