Return to random_product_system_io.adb CVS log

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

Annotation of OpenXM_contrib/PHC/Ada/Root_Counts/Product/random_product_system_io.adb, Revision 1.1.1.1

1.1       maekawa     1: with integer_io;
                      2: with Standard_Floating_Numbers;          use Standard_Floating_Numbers;
                      3: with Standard_Floating_Numbers_io;       use Standard_Floating_Numbers_io;
                      4: with Standard_Complex_Numbers;           use Standard_Complex_Numbers;
                      5: with Standard_Complex_Vectors;           use Standard_Complex_Vectors;
                      6: with Standard_Natural_Vectors;
                      7: with Symbol_Table,Symbol_Table_io;
                      8: with Standard_Complex_Polynomials;       use Standard_Complex_Polynomials;
                      9: with Standard_Complex_Polynomials_io;    use Standard_Complex_Polynomials_io;
                     10: with Standard_Complex_Poly_Systems;      use Standard_Complex_Poly_Systems;
                     11: with Standard_Complex_Poly_Systems_io;   use Standard_Complex_Poly_Systems_io;
                     12: with Random_Product_System;
                     13:
                     14: package body Random_Product_System_io is
                     15:
                     16:   procedure get ( n : out natural ) is
                     17:   begin
                     18:     Random_Product_System_io.get(Standard_Input,n);
                     19:   end get;
                     20:
                     21:   procedure get ( file : in file_type; n : out natural ) is
                     22:
                     23:     nn : natural;
                     24:
                     25:   begin
                     26:     integer_io.get(file,nn); n := nn;
                     27:     Random_Product_System.Init(nn);
                     28:     declare
                     29:       h : Vector(0..nn);
                     30:       pp : Poly;
                     31:       d : degrees := new Standard_Natural_Vectors.Vector'(1..nn => 0);
                     32:       stop : boolean;
                     33:     begin
                     34:       get(file,nn,pp); Clear(pp);
                     35:       for i in 1..nn loop
                     36:         stop := false;
                     37:         while not stop loop
                     38:           get(file,nn,pp);
                     39:           stop := (pp = Null_Poly);
                     40:           exit when stop;
                     41:           h(0) := Coeff(pp,d);
                     42:           for j in 1..nn loop
                     43:             d(j) := 1;
                     44:             h(j) := Coeff(pp,d);
                     45:             d(j) := 0;
                     46:           end loop;
                     47:           Random_Product_System.Add_Hyperplane(i,h);
                     48:         end loop;
                     49:       end loop;
                     50:     end;
                     51:   end get;
                     52:
                     53:   procedure put ( n,fore,after,exp : in natural ) is
                     54:   begin
                     55:     Random_Product_System_io.put(Standard_Output,n,fore,after,exp);
                     56:   end put;
                     57:
                     58:   procedure put ( file : in file_type; n,fore,after,exp : in natural ) is
                     59:
                     60:     h : Vector(0..n);
                     61:
                     62:     procedure Write_Number ( file : in file_type; x : in Complex_Number ) is
                     63:     begin
                     64:       if IMAG_PART(x) + 1.0 = 1.0
                     65:        then put(file,REAL_PART(x),fore,after,exp);
                     66:        else put(file,'(');
                     67:             put(file,REAL_PART(x),fore,after,exp);
                     68:             put(file,'+');
                     69:             put(file,IMAG_PART(x),fore,after,exp);
                     70:             put(file,')');
                     71:       end if;
                     72:     end Write_Number;
                     73:
                     74:   begin
                     75:     for i in 1..n loop
                     76:       put(file,"The hyperplanes for the "); integer_io.put(file,i,1);
                     77:       put_line(file,"th equation :");
                     78:       for j in 1..Random_Product_System.Number_Of_Hyperplanes(i) loop
                     79:         h := Random_Product_System.Get_Hyperplane(i,j);
                     80:         put(file,' ');
                     81:         for k in 1..n loop
                     82:          Write_Number(file,h(k));
                     83:          put(file,'*');
                     84:          Symbol_Table_io.put(file,Symbol_Table.Get(k));
                     85:          put(file," + ");
                     86:         end loop;
                     87:         Write_Number(file,h(0));
                     88:         new_line(file);
                     89:       end loop;
                     90:     end loop;
                     91:   end put;
                     92:
                     93: end Random_Product_System_io;