===================================================================
RCS file: /home/cvs/OpenXM/src/kan96xx/Doc/gfan.sm1,v
retrieving revision 1.11
retrieving revision 1.13
diff -u -p -r1.11 -r1.13
--- OpenXM/src/kan96xx/Doc/gfan.sm1	2005/07/07 06:07:46	1.11
+++ OpenXM/src/kan96xx/Doc/gfan.sm1	2009/08/26 04:54:17	1.13
@@ -1,6 +1,6 @@
-%  $OpenXM: OpenXM/src/kan96xx/Doc/gfan.sm1,v 1.10 2005/07/07 01:31:21 takayama Exp $
+%  $OpenXM: OpenXM/src/kan96xx/Doc/gfan.sm1,v 1.12 2005/07/07 07:53:37 takayama Exp $
 % cp cone.sm1 $OpenXM_HOME/src/kan96xx/Doc/gfan.sm1
-% $Id: gfan.sm1,v 1.11 2005/07/07 06:07:46 takayama Exp $
+% $Id: gfan.sm1,v 1.13 2009/08/26 04:54:17 takayama Exp $
 % iso-2022-jp
 %%Ref:  @s/2004/08/21-note.pdf
 
@@ -415,6 +415,13 @@ cone.comment message
 %>
 /cone.DhH  0 def
 
+%<
+% Global
+% gbCheck をするか? しないと結果はあやふや. しかしメモリ exhaust は防げる.
+% 使うときは /cone.epsilon,  /cone.epsilon.limit を十分小さくしておく.
+%>
+/cone.do_gbCheck 1 def
+
 % Default の cone.gb の定義. 各プログラムで再度定義してもよい.
 /cone.gb {
   cone.DhH { 
@@ -1017,6 +1024,11 @@ def
   [
      /ceq arg1 def 
      ceq pruneZeroVector /ceq set
+
+     ceq length 0 eq {
+       (Monomial ideal is not accepted as an input.) cone_ir_input
+     } { } ifelse
+
      ceq genPo2 /ceq2 set
      % ceq2 は polymake.data(polymake.INEQUALITIES(...)) 形式
      % polymake で ceq2 の次元の計算.
@@ -1871,6 +1883,11 @@ def
 % note: 2004.9.2
     cone (facetsv) getNode 2 get facet_i get /v set
     cone (facets)  getNode 2 get facet_i get /f set
+
+    v length 0 eq {
+       (The codimension of the linarity space of the Grobner cone seems to be 1 or 0.) cone_ir_input
+     } { } ifelse
+
     /vp v 0 get def
     1 1 v length 1 sub {
       /ii set
@@ -2048,15 +2065,20 @@ def
 %>
 /cone.gb_Dh {
   /arg2 set /arg1 set
-  [/ff /ww /gg] pushVariables
+  [/ff /ww /gg /gbopt] pushVariables
   [
     /ff arg1 def
     /ww arg2 def
     [(AutoReduce) 1] system_variable
     [cone.vv ring_of_differential_operators
      [ww] weight_vector 0] define_ring
-    [ff {toString .} map] ff getAttributeList setAttributeList 
-    groebner 0 get /gg set
+    %(---) messagen ff getAttributeList message
+    ff getAttributeList tag 0 eq {/gbopt [ ] def }
+    {
+       /gbopt ff getAttributeList def
+    } ifelse
+   [ff {toString .} map gbopt] 
+    groebner 0 get /gg set   %% groenber は attribute を受け付けない.
     /cone.gb_Dh.g gg def
     /arg1 gg def
   ] pop
@@ -3541,8 +3563,15 @@ def   
     %gNew が newWeight での GB か check. Yes なら reduced basis へ.
     %No なら null を戻す. 
 %%Ref: note @s/2005/06/30-note-gfan.pdf
-    gNew [(gbCheck) 1] setAttributeList newWeight 
-       cone.gb (gb) getAttribute 
+    cone.do_gbCheck not { 
+       (Warning! gbCheck is skipped.) message
+    } { 
+       (Doing gbCheck.) message
+    } ifelse
+    cone.do_gbCheck {
+     gNew [(gbCheck) 1] setAttributeList newWeight 
+        cone.gb (gb) getAttribute 
+    } { 1 } ifelse
     1 eq {
      gNew [(reduceOnly) 1] setAttributeList newWeight cone.gb /arg1 set
     }{ /arg1 null def } ifelse
@@ -3824,4 +3853,21 @@ def   
  ] pop
  popVariables
  arg1
+} def
+
+%%change
+/cone_ir_input {
+  /arg1 set 
+  [/msg ] pushVariables
+  [
+    /msg arg1 def
+    (---------------) message
+    msg message
+    (  ) message 
+    (Please also refer to the value of the variables cone.getConeInfo.rr0) message
+    ( cone.getConeInfo.rr1 cone.Lp cone.cinit) message 
+    $ cone.cinit (FACETS) getNode ::  $  message
+    (We are sorry that we cannot accept this input.) error
+  ] pop
+  popVariables
 } def