Return to bablphc.adb CVS log

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

Annotation of OpenXM_contrib/PHC/Ada/Main/bablphc.adb, Revision 1.1.1.1

1.1       maekawa     1: with text_io;                            use text_io;
                      2: with Communications_with_User;           use Communications_with_User;
                      3: with Timing_Package;                     use Timing_Package;
                      4: with Standard_Complex_Poly_Systems;      use Standard_Complex_Poly_Systems;
                      5: with Standard_Complex_Poly_Systems_io;   use Standard_Complex_Poly_Systems_io;
                      6: with Standard_Complex_Solutions;         use Standard_Complex_Solutions;
                      7: with Scaling;                            use Scaling;
                      8: with Black_Box_Root_Counting;
                      9: with Black_Polynomial_Continuations;     use Black_Polynomial_Continuations;
                     10:
                     11: procedure bablphc ( infilename,outfilename : in string ) is
                     12:
                     13:   procedure Read_System ( file : in out file_type; filename : in string;
                     14:                           p : in out Link_to_Poly_Sys ) is
                     15:   begin
                     16:     if filename /= ""
                     17:      then Open_Input_File(file,filename);
                     18:           get(file,p);
                     19:     end if;
                     20:   exception
                     21:     when others => put_line("Something is wrong with argument file...");
                     22:                    p := null; return;
                     23:   end Read_System;
                     24:
                     25:   procedure Timing_Summary ( file : in file_type;
                     26:                              roco,hoco,poco,total : in duration ) is
                     27:
                     28:     b0 : constant string :=
                     29:      "  ---------------------------------------------------------------------";
                     30:     b1 : constant string :=
                     31:      "  |                    TIMING INFORMATION SUMMARY                     |";
                     32:     b2 : constant string :=
                     33:      "  |   root counts  |  start system  |  continuation  |   total time   |";
                     34:
                     35:   begin
                     36:     put_line(file,b0);
                     37:     put_line(file,b1);
                     38:     put_line(file,b0);
                     39:     put_line(file,b2);
                     40:     put_line(file,b0);
                     41:     put(file,"  | ");
                     42:     print_hms(file,roco); put(file," | ");
                     43:     print_hms(file,hoco); put(file," | ");
                     44:     print_hms(file,poco); put(file," | ");
                     45:     print_hms(file,total); put_line(file," |");
                     46:     put_line(file,b0);
                     47:   end Timing_Summary;
                     48:
                     49:   procedure Main is
                     50:
                     51:     timer : timing_widget;
                     52:     infile,outfile : file_type;
                     53:     p,q : Link_to_Poly_Sys;
                     54:     rc : natural;
                     55:     sols : Solution_List;
                     56:     roco,hoco,poco,total : duration;
                     57:
                     58:   begin
                     59:     Read_System(infile,infilename,p);
                     60:     if p = null
                     61:      then new_line; get(p);
                     62:     end if;
                     63:     Create_Output_File(outfile,outfilename);
                     64:     put(outfile,p'last,p.all);
                     65:     q := new Poly_Sys(p'range);
                     66:     tstart(timer);
                     67:     declare
                     68:       pp : Poly_Sys(p'range) := p.all;
                     69:     begin
                     70:       Black_Box_Root_Counting(outfile,pp,rc,q.all,sols,roco,hoco);
                     71:       if rc /= 0
                     72:        then
                     73:          Scale(pp);
                     74:          Black_Box_Polynomial_Continuation(outfile,pp,q.all,sols,poco);
                     75:        end if;
                     76:     end;
                     77:     tstop(timer);
                     78:     total := Elapsed_User_Time(timer);
                     79:     Close(infile);
                     80:     new_line(outfile);
                     81:     print_times(outfile,timer,"Solving the polynomial system");
                     82:     new_line(outfile);
                     83:     Timing_Summary(outfile,roco,hoco,poco,total);
                     84:     Close(outfile);
                     85:   end Main;
                     86:
                     87: begin
                     88:   Main;
                     89: end bablphc;