Return to m_homogeneous_start_systems.adb CVS log

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

Annotation of OpenXM_contrib/PHC/Ada/Root_Counts/Product/m_homogeneous_start_systems.adb, Revision 1.1

1.1     ! maekawa     1: with Standard_Complex_Numbers;           use Standard_Complex_Numbers;
        !             2: with Standard_Random_Numbers;            use Standard_Random_Numbers;
        !             3: with Standard_Complex_Vectors;           use Standard_Complex_Vectors;
        !             4: with Standard_Integer_Matrices;
        !             5: with Sets_of_Unknowns;                   use Sets_of_Unknowns;
        !             6: with Degrees_in_Sets_of_Unknowns;        use Degrees_in_Sets_of_Unknowns;
        !             7: with Random_Product_System;
        !             8:
        !             9: package body m_Homogeneous_Start_Systems is
        !            10:
        !            11:   procedure Create_Random_Hyperplanes ( index,n,d : in natural; s : in Set ) is
        !            12:   begin
        !            13:     for i in 1..d loop
        !            14:       declare
        !            15:         h : Standard_Complex_Vectors.Vector(0..n);
        !            16:       begin
        !            17:         h(0) := Random1;
        !            18:         for j in 1..Dimension(s) loop
        !            19:           if Is_In(s,j)
        !            20:            then h(j) := Random1;
        !            21:            else h(j) := Create(0.0);
        !            22:           end if;
        !            23:         end loop;
        !            24:         Random_Product_System.Add_Hyperplane(index,h);
        !            25:       end;
        !            26:     end loop;
        !            27:   end Create_Random_Hyperplanes;
        !            28:
        !            29:   procedure Create_Random_System
        !            30:               ( n,m : natural; z : partition;
        !            31:                 d : Standard_Integer_Matrices.Matrix ) is
        !            32:
        !            33:   begin
        !            34:     for j in 1..m loop
        !            35:       for i in 1..n loop
        !            36:         Create_Random_Hyperplanes(i,n,d(i,j),z(j));
        !            37:       end loop;
        !            38:     end loop;
        !            39:   end Create_Random_System;
        !            40:
        !            41:   procedure m_Homogeneous_Start_System
        !            42:                  ( p : in Poly_Sys; z : in partition;
        !            43:                    q : out Poly_Sys; qsols : in out Solution_List ) is
        !            44:
        !            45:     n : constant natural := p'length;
        !            46:     m : constant natural := z'last;
        !            47:     d : constant Standard_Integer_Matrices.Matrix := Degree_Table(p,z);
        !            48:     nl : natural := 0;
        !            49:
        !            50:   begin
        !            51:     Random_Product_System.Init(n);
        !            52:     Create_Random_System(n,m,z,d);
        !            53:     Random_Product_System.Solve(qsols,nl);
        !            54:     q := Random_Product_System.Polynomial_System;
        !            55:   end m_Homogeneous_Start_System;
        !            56:
        !            57: end m_Homogeneous_Start_Systems;