Return to bablvali.adb CVS log

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

Annotation of OpenXM_contrib/PHC/Ada/Main/bablvali.adb, Revision 1.1

1.1     ! maekawa     1: with text_io,integer_io;                 use text_io,integer_io;
        !             2: with Timing_Package;                     use Timing_Package;
        !             3: with Communications_with_User;           use Communications_with_User;
        !             4: with File_Scanning;                      use File_Scanning;
        !             5: with Standard_Floating_Numbers;          use Standard_Floating_Numbers;
        !             6: with Standard_Floating_Numbers_io;       use Standard_Floating_Numbers_io;
        !             7: with Numbers_io;                         use Numbers_io;
        !             8: with Standard_Complex_Poly_Systems;      use Standard_Complex_Poly_Systems;
        !             9: with Standard_Complex_Poly_Systems_io;   use Standard_Complex_Poly_Systems_io;
        !            10: with Standard_Complex_Solutions;         use Standard_Complex_Solutions;
        !            11: with Standard_Complex_Solutions_io;      use Standard_Complex_Solutions_io;
        !            12: with Standard_Root_Refiners;             use Standard_Root_Refiners;
        !            13:
        !            14: procedure bablvali ( infilename,outfilename : in string ) is
        !            15:
        !            16:   timer : timing_widget;
        !            17:   lp : Link_to_Poly_Sys;
        !            18:   infile,outfile : file_type;
        !            19:   ans : character;
        !            20:   n : natural;
        !            21:   sysonfile,found : boolean;
        !            22:   sols : Solution_List;
        !            23:
        !            24:   procedure Read_System ( file : in out file_type; filename : in string ) is
        !            25:   begin
        !            26:     if filename /= ""
        !            27:      then Open(file,in_file,filename);
        !            28:           get(file,lp);
        !            29:           sysonfile := true;
        !            30:      else sysonfile := false;
        !            31:     end if;
        !            32:   exception
        !            33:     when others => put_line("Something is wrong with argument file...");
        !            34:                    sysonfile := false;
        !            35:                    lp := null; return;
        !            36:   end Read_System;
        !            37:
        !            38:   procedure Refine_Roots is
        !            39:
        !            40:     epsxa,epsfa : constant double_float := 10.0**(-8);
        !            41:     tolsing : constant double_float := 10.0**(-8);
        !            42:     maxit : constant natural := 3;
        !            43:     numb : natural := 0;
        !            44:     refsols : Solution_List;
        !            45:
        !            46:   begin
        !            47:     new_line(outfile);
        !            48:     put_line(outfile,"ROOT REFINING PARAMETERS");
        !            49:     put(outfile,"  tolerance for error on the root : ");
        !            50:     put(outfile,epsxa,2,3,3); new_line(outfile);
        !            51:     put(outfile,"  tolerance for residual          : ");
        !            52:     put(outfile,epsfa,2,3,3); new_line(outfile);
        !            53:     put(outfile,"  tolerance for singular roots    : ");
        !            54:     put(outfile,tolsing,2,3,3); new_line(outfile);
        !            55:     put(outfile,"  maximum number of iterations    : ");
        !            56:     put(outfile,maxit,2); new_line(outfile);
        !            57:     tstart(timer);
        !            58:     Reporting_Root_Refiner
        !            59:       (outfile,lp.all,sols,refsols,epsxa,epsfa,tolsing,numb,maxit,false);
        !            60:     tstop(timer);
        !            61:     new_line(outfile);
        !            62:     put_line(outfile,"THE SOLUTIONS :");
        !            63:     put(outfile,Length_Of(refsols),Head_Of(refsols).n,refsols);
        !            64:     new_line(outfile);
        !            65:     print_times(outfile,timer,"Root refining");
        !            66:   end Refine_Roots;
        !            67:
        !            68: begin
        !            69:   Read_System(infile,infilename);
        !            70:   if lp = null
        !            71:    then
        !            72:      new_line;
        !            73:      put("Is the system on file ? (y/n) ");
        !            74:      Ask_Yes_or_No(ans);
        !            75:      if ans = 'y'
        !            76:       then put_line("Reading the name of the input file.");
        !            77:            Read_Name_and_Open_File(infile);
        !            78:            get(infile,lp);
        !            79:            sysonfile := true;
        !            80:       else put("Give the dimension : "); get(n);
        !            81:            lp := new Poly_Sys(1..n);
        !            82:            put("Give "); put(n,1); put(" "); put(n,1);
        !            83:            put_line("-variate polynomials :");
        !            84:            get(n,lp.all);
        !            85:            skip_line;
        !            86:            sysonfile := false;
        !            87:      end if;
        !            88:   end if;
        !            89:
        !            90:   Create_Output_File(outfile,outfilename);
        !            91:   put(outfile,lp.all);
        !            92:
        !            93:   if sysonfile
        !            94:    then Scan_and_Skip(infile,"SOLUTIONS",found);
        !            95:         if found
        !            96:          then get(infile,sols);
        !            97:         end if;
        !            98:         Close(infile);
        !            99:    else found := false;
        !           100:   end if;
        !           101:   if not found
        !           102:    then new_line; Read(sols);
        !           103:   end if;
        !           104:
        !           105:   Refine_Roots;
        !           106: end bablvali;