Return to generic_jacobian_matrices.adb CVS log

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

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

1.1     ! maekawa     1: package body Generic_Jacobian_Matrices is
        !             2:
        !             3: -- CREATORS :
        !             4:
        !             5:   function Create ( p : Poly_Sys ) return Jaco_Mat is
        !             6:
        !             7:     res : Jaco_Mat(p'range,1..Number_of_Unknowns(p(p'first)));
        !             8:
        !             9:   begin
        !            10:     for i in res'range(1) loop
        !            11:       for j in res'range(2) loop
        !            12:         res(i,j) := Diff(p(i),j);
        !            13:       end loop;
        !            14:     end loop;
        !            15:     return res;
        !            16:   end Create;
        !            17:
        !            18:   function Create ( j : Jaco_Mat ) return Eval_Jaco_Mat is
        !            19:
        !            20:     res : Eval_Jaco_Mat(j'range(1),j'range(2));
        !            21:
        !            22:   begin
        !            23:     for k in j'range(1) loop
        !            24:       for l in j'range(2) loop
        !            25:         res(k,l) := Create(j(k,l));
        !            26:       end loop;
        !            27:     end loop;
        !            28:     return res;
        !            29:   end Create;
        !            30:
        !            31:   procedure Create ( p : Poly_Sys;
        !            32:                      j : out Eval_Coeff_Jaco_Mat; m : out Mult_Factors ) is
        !            33:
        !            34:     nb : constant natural := Number_of_Unknowns(p(p'first));
        !            35:     nbk : natural;
        !            36:
        !            37:   begin
        !            38:     for k in p'range loop
        !            39:       nbk := Number_of_Terms(p(k));
        !            40:       for l in 1..nb loop
        !            41:         declare
        !            42:           mkl : Vector(1..nbk);
        !            43:         begin
        !            44:           Diff(p(k),l,j(k,l),mkl);
        !            45:           m(k,l) := new Vectors.Vector'(mkl);
        !            46:         end;
        !            47:       end loop;
        !            48:     end loop;
        !            49:   end Create;
        !            50:
        !            51: -- EVALUATORS :
        !            52:
        !            53:   function Eval ( j : Jaco_Mat; x : Vector ) return Matrix is
        !            54:
        !            55:     m : Matrix(j'range(1),j'range(2));
        !            56:
        !            57:   begin
        !            58:     for k in j'range(1) loop
        !            59:       for l in j'range(2) loop
        !            60:         m(k,l) := Eval(Poly(j(k,l)),x);
        !            61:       end loop;
        !            62:     end loop;
        !            63:     return m;
        !            64:   end Eval;
        !            65:
        !            66:   function Eval ( j : Eval_Jaco_Mat; x : Vector ) return Matrix is
        !            67:
        !            68:     m : Matrix(j'range(1),j'range(2));
        !            69:
        !            70:   begin
        !            71:     for k in j'range(1) loop
        !            72:       for l in j'range(2) loop
        !            73:         m(k,l) := Eval(Eval_Poly(j(k,l)),x);
        !            74:       end loop;
        !            75:     end loop;
        !            76:     return m;
        !            77:   end Eval;
        !            78:
        !            79:   function Eval ( j : Eval_Coeff_Jaco_Mat; m : Mult_Factors;
        !            80:                   c : VecVec; x : Vector ) return Matrix is
        !            81:
        !            82:     res : Matrix(j'range(1),j'range(2));
        !            83:
        !            84:   begin
        !            85:     for k in j'range(1) loop
        !            86:       declare
        !            87:         cm : Vector(c(k)'range);
        !            88:       begin
        !            89:         for l in j'range(2) loop
        !            90:           for i in cm'range loop
        !            91:             cm(i) := m(k,l)(i)*c(k)(i);
        !            92:           end loop;
        !            93:           res(k,l) := Eval(Eval_Coeff_Poly(j(k,l)),cm,x);
        !            94:         end loop;
        !            95:       end;
        !            96:     end loop;
        !            97:     return res;
        !            98:   end Eval;
        !            99:
        !           100: -- DESTRUCTORS :
        !           101:
        !           102:   procedure Clear ( j : in out Jaco_Mat ) is
        !           103:   begin
        !           104:     for k in j'range(1) loop
        !           105:       for l in j'range(2) loop
        !           106:         Clear(j(k,l));
        !           107:       end loop;
        !           108:     end loop;
        !           109:   end Clear;
        !           110:
        !           111:   procedure Clear ( j : in out Eval_Jaco_Mat ) is
        !           112:   begin
        !           113:     for k in j'range(1) loop
        !           114:       for l in j'range(2) loop
        !           115:         Clear(j(k,l));
        !           116:       end loop;
        !           117:     end loop;
        !           118:   end Clear;
        !           119:
        !           120:   procedure Clear ( j : in out Eval_Coeff_Jaco_Mat ) is
        !           121:   begin
        !           122:     for k in j'range(1) loop
        !           123:       for l in j'range(2) loop
        !           124:         Clear(j(k,l));
        !           125:       end loop;
        !           126:     end loop;
        !           127:   end Clear;
        !           128:
        !           129:   procedure Clear ( m : in out Mult_Factors ) is
        !           130:   begin
        !           131:     for k in m'range(1) loop
        !           132:       for l in m'range(2) loop
        !           133:         Clear(m(k,l));
        !           134:       end loop;
        !           135:     end loop;
        !           136:   end Clear;
        !           137:
        !           138: end Generic_Jacobian_Matrices;