Return to generic_polynomial_systems.adb CVS log

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

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

1.1     ! maekawa     1: with unchecked_deallocation;
        !             2:
        !             3: package body Generic_Polynomial_Systems is
        !             4:
        !             5: -- COPYING :
        !             6:
        !             7:   procedure Copy ( p : in Poly_Sys; q : in out Poly_Sys ) is
        !             8:   begin
        !             9:     for i in p'range loop
        !            10:       Copy(p(i),q(i));
        !            11:     end loop;
        !            12:   end Copy;
        !            13:
        !            14: -- ARITHMETIC OPERATIONS :
        !            15:
        !            16:   function "+" ( p,q : Poly_Sys ) return Poly_Sys is
        !            17:
        !            18:     res : Poly_Sys(p'range);
        !            19:
        !            20:   begin
        !            21:     for i in p'range loop
        !            22:       res(i) := p(i)+q(i);
        !            23:     end loop;
        !            24:     return res;
        !            25:   end "+";
        !            26:
        !            27:   function "-" ( p,q : Poly_Sys ) return Poly_Sys is
        !            28:
        !            29:     res : Poly_Sys(p'range);
        !            30:
        !            31:   begin
        !            32:     for i in p'range loop
        !            33:       res(i) := p(i)-q(i);
        !            34:     end loop;
        !            35:     return res;
        !            36:   end "-";
        !            37:
        !            38:   function "-" ( p : Poly_Sys ) return Poly_Sys is
        !            39:
        !            40:     res : Poly_Sys(p'range);
        !            41:
        !            42:   begin
        !            43:     for i in p'range loop
        !            44:       res(i) := -p(i);
        !            45:     end loop;
        !            46:     return res;
        !            47:   end "-";
        !            48:
        !            49:   function "*" ( a : number; p : Poly_Sys ) return Poly_Sys is
        !            50:
        !            51:     res : Poly_Sys(p'range);
        !            52:
        !            53:   begin
        !            54:     for k in p'range loop
        !            55:       res(k) := a*p(k);
        !            56:     end loop;
        !            57:     return res;
        !            58:   end "*";
        !            59:
        !            60:   function "*" ( p : Poly_Sys; a : number ) return Poly_Sys is
        !            61:   begin
        !            62:     return a*p;
        !            63:   end "*";
        !            64:
        !            65:   procedure Add ( p : in out Poly_Sys; q : in Poly_Sys ) is
        !            66:   begin
        !            67:     for i in p'range loop
        !            68:       Add(p(i),q(i));
        !            69:     end loop;
        !            70:   end Add;
        !            71:
        !            72:   procedure Sub ( p : in out Poly_Sys; q : in Poly_Sys ) is
        !            73:   begin
        !            74:     for i in p'range loop
        !            75:       Sub(p(i),q(i));
        !            76:     end loop;
        !            77:   end Sub;
        !            78:
        !            79:   procedure Min ( p : in out Poly_Sys ) is
        !            80:   begin
        !            81:     for i in p'range loop
        !            82:       Min(p(i));
        !            83:     end loop;
        !            84:   end Min;
        !            85:
        !            86:   procedure Mul ( p : in out Poly_Sys; a : in number ) is
        !            87:   begin
        !            88:     for k in p'range loop
        !            89:       Mul(p(k),a);
        !            90:     end loop;
        !            91:   end Mul;
        !            92:
        !            93:   function Diff ( p : Poly_Sys; i : natural ) return Poly_Sys is
        !            94:
        !            95:     res : Poly_Sys(p'range);
        !            96:
        !            97:   begin
        !            98:     for j in p'range loop
        !            99:       res(j) := Diff(p(j),i);
        !           100:     end loop;
        !           101:     return res;
        !           102:   end Diff;
        !           103:
        !           104:   procedure Diff ( p : in out Poly_Sys; i : in natural ) is
        !           105:   begin
        !           106:     for j in p'range loop
        !           107:       Diff(p(j),i);
        !           108:     end loop;
        !           109:   end Diff;
        !           110:
        !           111: -- DESTRUCTORS :
        !           112:
        !           113:   procedure Clear ( p : in out Poly_Sys ) is
        !           114:   begin
        !           115:     for i in p'range loop
        !           116:       Clear(p(i));
        !           117:     end loop;
        !           118:   end Clear;
        !           119:
        !           120:   procedure Clear ( p : in out Link_to_Poly_Sys ) is
        !           121:
        !           122:     procedure free is new unchecked_deallocation(Poly_Sys,Link_to_Poly_Sys);
        !           123:
        !           124:   begin
        !           125:     if p /= null
        !           126:      then Clear(p.all);
        !           127:     end if;
        !           128:     free(p);
        !           129:   end Clear;
        !           130:
        !           131: end Generic_Polynomial_Systems;