Return to curves_into_grassmannian_io.adb CVS log

Up to [local] / OpenXM_contrib / PHC / Ada / Schubert

Annotation of OpenXM_contrib/PHC/Ada/Schubert/curves_into_grassmannian_io.adb, Revision 1.1

1.1     ! maekawa     1: with Characters_and_Numbers;             use Characters_and_Numbers;
        !             2: with Symbol_Table;                       use Symbol_Table;
        !             3: with Brackets;                           use Brackets;
        !             4: with Curves_into_Grassmannian;           use Curves_into_Grassmannian;
        !             5:
        !             6: package body Curves_into_Grassmannian_io is
        !             7:
        !             8:   procedure Set_up_Symbol_Table
        !             9:               ( m,p,q : in natural; top,bottom : in Bracket;
        !            10:                 addmix : in natural ) is
        !            11:
        !            12:     n : constant natural := Number_of_Variables(top,bottom) + 2 + addmix;
        !            13:     rws : constant natural := (m+p)*(q+1);
        !            14:     row,deg : natural;
        !            15:     sb : Symbol;
        !            16:
        !            17:   begin
        !            18:     if not Symbol_Table.Empty
        !            19:      then Symbol_Table.Clear;
        !            20:     end if;
        !            21:     Symbol_Table.Init(n);             -- initialization with #variables
        !            22:     sb := (sb'range => ' ');
        !            23:     sb(1) := 'x';
        !            24:     sb(4) := 's';
        !            25:     for j in 1..p loop                -- adding the rest columnwise
        !            26:       row := 0; deg := 0;
        !            27:       for i in 1..rws loop
        !            28:         row := row + 1;
        !            29:         if i >= top(j) and i <= bottom(j)
        !            30:          then sb(2) := Convert_Hexadecimal(row);
        !            31:               sb(3) := Convert_Hexadecimal(j);
        !            32:               sb(5) := Convert_Hexadecimal(deg);
        !            33:               Symbol_Table.Add(sb);
        !            34:         end if;
        !            35:         if i mod (m+p) = 0
        !            36:          then row := 0; deg := deg+1;
        !            37:         end if;
        !            38:       end loop;
        !            39:     end loop;
        !            40:     sb := (sb'range => ' ');
        !            41:     sb(1) := 's';
        !            42:     if addmix = 0
        !            43:      then Symbol_Table.Add(sb);       -- adding "s"
        !            44:      else sb(2) := '1';
        !            45:           Symbol_Table.Add(sb);       -- adding "s1"
        !            46:           sb(2) := '2';
        !            47:           Symbol_Table.Add(sb);       -- adding "s2"
        !            48:           sb(2) := ' ';
        !            49:     end if;
        !            50:     sb(1) := 't';
        !            51:     Symbol_Table.Add(sb);             -- adding "t"
        !            52:   end Set_up_Symbol_Table;
        !            53:
        !            54:   procedure One_Set_up_Symbol_Table
        !            55:               ( m,p,q : in natural; top,bottom : in Bracket ) is
        !            56:   begin
        !            57:     Set_up_Symbol_Table(m,p,q,top,bottom,0);
        !            58:   end One_Set_up_Symbol_Table;
        !            59:
        !            60:   procedure Two_Set_up_Symbol_Table
        !            61:               ( m,p,q : in natural; top,bottom : in Bracket ) is
        !            62:   begin
        !            63:     Set_up_Symbol_Table(m,p,q,top,bottom,1);
        !            64:   end Two_Set_up_Symbol_Table;
        !            65:
        !            66:   procedure Reduce_Symbols ( top,bottom : in Bracket; locmap : in Matrix ) is
        !            67:
        !            68:     ind : natural := 0;
        !            69:
        !            70:   begin
        !            71:     for j in locmap'range(2) loop
        !            72:       for i in top(j)..bottom(j) loop
        !            73:         ind := ind+1;
        !            74:         if locmap(i,j) = 1
        !            75:          then Symbol_Table.Remove(ind-j+1);
        !            76:         end if;
        !            77:       end loop;
        !            78:     end loop;
        !            79:   end Reduce_Symbols;
        !            80:
        !            81: end Curves_into_Grassmannian_io;