Return to reduction_of_nonsquare_systems.adb CVS log

Up to [local] / OpenXM_contrib / PHC / Ada / Homotopy

Annotation of OpenXM_contrib/PHC/Ada/Homotopy/reduction_of_nonsquare_systems.adb, Revision 1.1.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_Polynomials;       use Standard_Complex_Polynomials;
                      4: with Reduction_of_Polynomials;           use Reduction_of_Polynomials;
                      5:
                      6: package body Reduction_of_Nonsquare_Systems is
                      7:
                      8:   function Random_Square ( p : in Poly_Sys ) return Poly_Sys is
                      9:
                     10:     m : constant natural := Number_of_Unknowns(p(p'first));
                     11:     res : Poly_Sys(1..m);
                     12:
                     13:   begin
                     14:     for i in res'range loop
                     15:       Copy(p(i),res(i));
                     16:       for j in m+1..p'last loop
                     17:         declare
                     18:           a : Complex_Number := Random1;
                     19:           tmp : Poly := a*p(j);
                     20:         begin
                     21:           Add(res(i),tmp);
                     22:           Clear(tmp);
                     23:         end;
                     24:       end loop;
                     25:     end loop;
                     26:     return res;
                     27:   end Random_Square;
                     28:
                     29:   function Reduced_Square ( p : in Poly_Sys ) return Poly_Sys is
                     30:
                     31:     m : constant natural := Number_of_Unknowns(p(p'first));
                     32:     res : Poly_Sys(1..m);
                     33:
                     34:   begin
                     35:     for i in res'range loop
                     36:       Copy(p(i),res(i));
                     37:       for j in m+1..p'last loop
                     38:         declare
                     39:           tmp : Poly := Rpoly(res(i),p(j));
                     40:         begin
                     41:           Copy(tmp,res(i)); Clear(tmp);
                     42:         end;
                     43:       end loop;
                     44:     end loop;
                     45:     return res;
                     46:   end Reduced_Square;
                     47:
                     48: end Reduction_of_Nonsquare_Systems;