Return to minkowski_polynomials.adb CVS log

Up to [local] / OpenXM_contrib / PHC / Ada / Root_Counts / Dynlift

Annotation of OpenXM_contrib/PHC/Ada/Root_Counts/Dynlift/minkowski_polynomials.adb, Revision 1.1

1.1     ! maekawa     1: with Standard_Floating_Numbers;          use Standard_Floating_Numbers;
        !             2: with Standard_Complex_Numbers;           use Standard_Complex_Numbers;
        !             3: with Standard_Natural_Vectors;
        !             4: with Cayley_Embedding;                   use Cayley_Embedding;
        !             5: with Mixed_Volume_Computation;           use Mixed_Volume_Computation;
        !             6:
        !             7: package body Minkowski_Polynomials is
        !             8:
        !             9:   function Minkowski_Polynomial ( n,r : natural ) return Poly is
        !            10:
        !            11:     res : Poly := Null_Poly;
        !            12:     acc : Degrees := new Standard_Natural_Vectors.Vector'(1..r => 0);
        !            13:
        !            14:     procedure Generate_Monomials
        !            15:                   ( k,sum : in natural; deg : in out Degrees ) is
        !            16:
        !            17:     -- DESCRIPTION :
        !            18:     --   Generates all exponent vectors whose sum equals n.
        !            19:
        !            20:       t : Term;
        !            21:
        !            22:     begin
        !            23:       if k = r
        !            24:        then t.cf := Create(1.0);
        !            25:             deg(r) := n-sum;
        !            26:             t.dg := deg;
        !            27:             Add(res,t);
        !            28:        else for i in 0..(n-sum) loop
        !            29:               deg(k) := i;
        !            30:               Generate_Monomials(k+1,sum+i,deg);
        !            31:             end loop;
        !            32:       end if;
        !            33:     end Generate_Monomials;
        !            34:
        !            35:   begin
        !            36:     Generate_Monomials(1,0,acc);
        !            37:     Standard_Natural_Vectors.Clear
        !            38:       (Standard_Natural_Vectors.Link_to_Vector(acc));
        !            39:     return res;
        !            40:   end Minkowski_Polynomial;
        !            41:
        !            42:   function Convert ( dg : Degrees ) return Standard_Integer_Vectors.Vector is
        !            43:
        !            44:   -- DESCRIPTION :
        !            45:   --   Converts the degrees vector to a vector with integer numbers.
        !            46:
        !            47:     res : Standard_Integer_Vectors.Vector(dg'range);
        !            48:
        !            49:   begin
        !            50:     for i in res'range loop
        !            51:       res(i) := dg(i);
        !            52:     end loop;
        !            53:     return res;
        !            54:   end Convert;
        !            55:
        !            56:   procedure Minkowski_Polynomial
        !            57:                   ( p : in out Poly; t : in Triangulation; n : in natural;
        !            58:                     mix : in Vector; mixsub : out Mixed_Subdivision ) is
        !            59:
        !            60:     procedure Coefficient_Volume
        !            61:                   ( submix : in Vector; sub : in Mixed_Subdivision;
        !            62:                     vol : out natural ) is
        !            63:     begin
        !            64:       vol := Mixed_Volume(n,submix,sub);
        !            65:     end Coefficient_Volume;
        !            66:     procedure Coefficient_Volumes is
        !            67:       new Minkowski_Polynomial_Subdivisions(Coefficient_Volume);
        !            68:
        !            69:   begin
        !            70:     Coefficient_Volumes(p,t,n,mix,mixsub);
        !            71:   end Minkowski_Polynomial;
        !            72:
        !            73:   procedure Minkowski_Polynomial_Subdivisions
        !            74:                  ( p : in out Poly; t : in Triangulation; n : in natural;
        !            75:                    mix : in Vector; mixsub : out Mixed_Subdivision ) is
        !            76:
        !            77:     procedure Coefficient_Volume ( tt : in out Term; cont : out boolean ) is
        !            78:
        !            79:       wrkmix : Vector(mix'range) := Convert(tt.dg);
        !            80:       wrksub : Mixed_Subdivision := Extract_Mixed_Cells(n,wrkmix,t);
        !            81:       vol : natural;
        !            82:
        !            83:     begin
        !            84:       Deflate(n,wrksub);
        !            85:       Process(wrkmix,wrksub,vol);
        !            86:       tt.cf := Create(double_float(vol));
        !            87:       if wrkmix = mix
        !            88:        then mixsub := wrksub;
        !            89:        else Deep_Clear(wrksub);
        !            90:       end if;
        !            91:       cont := true;
        !            92:     end Coefficient_Volume;
        !            93:     procedure Coefficient_Volumes is new Changing_Iterator(Coefficient_Volume);
        !            94:
        !            95:   begin
        !            96:     Coefficient_Volumes(p);
        !            97:   end Minkowski_Polynomial_Subdivisions;
        !            98:
        !            99: end Minkowski_Polynomials;