Return to matrix_indeterminates.adb CVS log

Up to [local] / OpenXM_contrib / PHC / Ada / Math_Lib / Polynomials

Annotation of OpenXM_contrib/PHC/Ada/Math_Lib/Polynomials/matrix_indeterminates.adb, Revision 1.1

1.1     ! maekawa     1: with Characters_and_Numbers;             use Characters_and_Numbers;
        !             2: with Standard_Complex_Numbers;           use Standard_Complex_Numbers;
        !             3: with Standard_Natural_Vectors;
        !             4:
        !             5: package body Matrix_Indeterminates is
        !             6:
        !             7:   procedure Initialize_Symbols ( n,d : in natural ) is
        !             8:   begin
        !             9:     Symbol_Table.Init(n*d);
        !            10:     for i in 1..n loop
        !            11:       for j in 1..d loop
        !            12:         declare
        !            13:           s : Symbol := X_ij(i,j);
        !            14:         begin
        !            15:           Symbol_Table.Add(s);
        !            16:         end;
        !            17:       end loop;
        !            18:     end loop;
        !            19:   end Initialize_Symbols;
        !            20:
        !            21:   function X_ij ( i,j : natural ) return Symbol is
        !            22:
        !            23:     res : Symbol;
        !            24:
        !            25:   begin
        !            26:     res(1) := 'x';
        !            27:     res(2) := Convert_Hexadecimal(i);
        !            28:     res(3) := Convert_Hexadecimal(j);
        !            29:     for i in 4..res'last loop
        !            30:       res(i) := ' ';
        !            31:     end loop;
        !            32:     return res;
        !            33:   end X_ij;
        !            34:
        !            35:   function Monomial ( n,d,i,j : natural ) return Poly is
        !            36:
        !            37:     res : Poly;
        !            38:     t : Term;
        !            39:
        !            40:   begin
        !            41:     t.cf := Create(1.0);
        !            42:     t.dg := new Standard_Natural_Vectors.Vector'(1..n*d => 0);
        !            43:     t.dg((i-1)*d+j) := 1;
        !            44:     res := Create(t);
        !            45:     Clear(t.dg);
        !            46:     return res;
        !            47:   end Monomial;
        !            48:
        !            49:   procedure Reduce_Symbols ( locmap : in Matrix ) is
        !            50:
        !            51:     d : constant natural := locmap'length(2);
        !            52:
        !            53:   begin
        !            54:     for i in reverse locmap'range(1) loop
        !            55:       for j in reverse locmap'range(2) loop
        !            56:         if locmap(i,j) /= 2
        !            57:          then Symbol_Table.Remove(X_ij(i,j));
        !            58:         end if;
        !            59:       end loop;
        !            60:     end loop;
        !            61:   end Reduce_Symbols;
        !            62:
        !            63:   procedure Clear_Symbols is
        !            64:   begin
        !            65:     Symbol_Table.Clear;
        !            66:   end Clear_Symbols;
        !            67:
        !            68: end Matrix_Indeterminates;