Return to exponent_vectors.adb CVS log

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

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

1.1     ! maekawa     1: package body Exponent_Vectors is
        !             2:
        !             3: -- CREATORS :
        !             4:
        !             5:   function Create ( p : Standard_Complex_Laur_Polys.Poly ) return VecVec is
        !             6:
        !             7:     use Standard_Complex_Laur_Polys;
        !             8:     res : VecVec(1..Number_of_Terms(p));
        !             9:     ind : natural := 0;
        !            10:
        !            11:     procedure Add_Exponent ( t : in Term; continue : out boolean ) is
        !            12:     begin
        !            13:       ind := ind + 1;
        !            14:       res(ind) := new Vector(t.dg'range);
        !            15:       for i in t.dg'range loop
        !            16:         res(ind)(i) := t.dg(i);
        !            17:       end loop;
        !            18:       continue := true;
        !            19:     end Add_Exponent;
        !            20:     procedure Add_Exponents is new Visiting_Iterator(Add_Exponent);
        !            21:
        !            22:   begin
        !            23:     Add_Exponents(p);
        !            24:     return res;
        !            25:   end Create;
        !            26:
        !            27:   function Create ( p : Standard_Complex_Polynomials.Poly ) return VecVec is
        !            28:
        !            29:     use Standard_Complex_Polynomials;
        !            30:     res : VecVec(1..Number_of_Terms(p));
        !            31:     ind : natural := 0;
        !            32:
        !            33:     procedure Add_Exponent ( t : in Term; continue : out boolean ) is
        !            34:     begin
        !            35:       ind := ind + 1;
        !            36:       res(ind) := new Vector(t.dg'range);
        !            37:       for i in t.dg'range loop
        !            38:         res(ind)(i) := t.dg(i);
        !            39:       end loop;
        !            40:       continue := true;
        !            41:     end Add_Exponent;
        !            42:     procedure Add_Exponents is new Visiting_Iterator(Add_Exponent);
        !            43:
        !            44:   begin
        !            45:     Add_Exponents(p);
        !            46:     return res;
        !            47:   end Create;
        !            48:
        !            49:   function Create ( p : Poly_Sys ) return Exponent_Vectors_Array is
        !            50:
        !            51:     res : Exponent_Vectors_Array(p'range);
        !            52:
        !            53:   begin
        !            54:     for i in p'range loop
        !            55:       declare
        !            56:         cpi : constant VecVec := Create(p(i));
        !            57:       begin
        !            58:         res(i) := new VecVec(cpi'range);
        !            59:         for j in cpi'range loop
        !            60:           res(i)(j) := cpi(j);
        !            61:         end loop;
        !            62:       end;  -- a detour for GNAT 3.07
        !            63:      -- res(i) := new VecVec'(Create(p(i)));
        !            64:     end loop;
        !            65:     return res;
        !            66:   end Create;
        !            67:
        !            68:   function Create ( p : Laur_Sys ) return Exponent_Vectors_Array is
        !            69:
        !            70:     res : Exponent_Vectors_Array(p'range);
        !            71:
        !            72:   begin
        !            73:     for i in p'range loop
        !            74:       declare
        !            75:         cpi : constant VecVec := Create(p(i));
        !            76:       begin
        !            77:         res(i) := new VecVec(cpi'range);
        !            78:         for j in cpi'range loop
        !            79:           res(i)(j) := cpi(j);
        !            80:         end loop;
        !            81:       end;  -- a detour for GNAT 3.07
        !            82:      -- res(i) := new VecVec'(Create(p(i)));
        !            83:     end loop;
        !            84:     return res;
        !            85:   end Create;
        !            86:
        !            87: -- SELECTOR :
        !            88:
        !            89:   function Position ( ev : VecVec; v : Vector ) return integer is
        !            90:   begin
        !            91:     for i in ev'range loop
        !            92:       if Equal(ev(i).all,v)
        !            93:        then return i;
        !            94:       end if;
        !            95:     end loop;
        !            96:     return ev'last+1;
        !            97:   end Position;
        !            98:
        !            99: -- EVALUATORS :
        !           100:
        !           101:   function Eval ( e : Vector; c : Complex_Number;
        !           102:                   x : Standard_Complex_Vectors.Vector )
        !           103:                 return Complex_Number is
        !           104:
        !           105:     res : Complex_Number := c;
        !           106:
        !           107:   begin
        !           108:     for i in e'range loop
        !           109:       for j in 1..e(i) loop
        !           110:         res := res*x(i);
        !           111:       end loop;
        !           112:       for j in 1..-e(i) loop
        !           113:         res := res/x(i);
        !           114:       end loop;
        !           115:     end loop;
        !           116:     return res;
        !           117:   end Eval;
        !           118:
        !           119:   function Eval ( ev : VecVec; c,x : Standard_Complex_Vectors.Vector )
        !           120:                 return Complex_Number is
        !           121:
        !           122:     res : Complex_Number := Eval(ev(ev'first).all,c(c'first),x);
        !           123:
        !           124:   begin
        !           125:     for i in c'first+1..c'last loop
        !           126:       res := res + Eval(ev(i).all,c(i),x);
        !           127:     end loop;
        !           128:     return res;
        !           129:   end Eval;
        !           130:
        !           131:   function Eval ( ev : Exponent_Vectors_Array;
        !           132:                   c : Standard_Complex_VecVecs.VecVec;
        !           133:                   x : Standard_Complex_Vectors.Vector )
        !           134:                 return Standard_Complex_Vectors.Vector is
        !           135:
        !           136:     res : Standard_Complex_Vectors.Vector(x'range);
        !           137:
        !           138:   begin
        !           139:     for i in ev'range loop
        !           140:       res(i) := Eval(ev(i).all,c(i).all,x);
        !           141:     end loop;
        !           142:     return res;
        !           143:   end Eval;
        !           144:
        !           145: -- DESTRUCTORS :
        !           146:
        !           147:   procedure Clear ( v : in out Exponent_Vectors_Array ) is
        !           148:   begin
        !           149:     for i in v'range loop
        !           150:       Deep_Clear(v(i));
        !           151:     end loop;
        !           152:   end Clear;
        !           153:
        !           154: end Exponent_Vectors;