Return to mainphc.adb CVS log

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

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

1.1       maekawa     1: with text_io,integer_io;                 use text_io,integer_io;
                      2: with Communications_with_User;           use Communications_with_User;
                      3: with Timing_Package;                     use Timing_Package;
                      4: with Standard_Complex_Numbers;           use Standard_Complex_Numbers;
                      5: with Standard_Complex_Vectors;           use Standard_Complex_Vectors;
                      6: with Standard_Complex_Polynomials;       use Standard_Complex_Polynomials;
                      7: with Standard_Complex_Poly_Systems;      use Standard_Complex_Poly_Systems;
                      8: with Standard_Complex_Poly_Systems_io;   use Standard_Complex_Poly_Systems_io;
                      9: with Standard_Complex_Solutions;         use Standard_Complex_Solutions;
                     10: with Homotopy;
                     11: with Drivers_for_Scaling;                use Drivers_for_Scaling;
                     12: with Drivers_for_Reduction;              use Drivers_for_Reduction;
                     13: with Drivers_for_Poly_Continuation;      use Drivers_for_Poly_Continuation;
                     14: with Driver_for_Root_Counts;
                     15: with Drivers_for_Homotopy_Creation;      use Drivers_for_Homotopy_Creation;
                     16: with Driver_for_Root_Refining;
                     17: with Bye_Bye_Message;
                     18:
                     19: procedure mainphc ( infilename,outfilename : in string ) is
                     20:
                     21:   outpt : file_type;
                     22:   lp : Link_to_Poly_Sys := null;
                     23:   timer : timing_widget;
                     24:
                     25:   procedure Display_Options is
                     26:
                     27:   -- DESCRIPTION :
                     28:   --   Displays an overview of all options on screen.
                     29:
                     30:     o : array(1..8) of string(1..65);
                     31:
                     32:   begin
                     33:     put_line("Running in full mode.  Note also the following options:");
                     34:     o(1) := "  phc -s : Equation and variable Scaling on system and solutions ";
                     35:     o(2) := "  phc -d : Linear and nonlinear Reduction w.r.t. the total degree";
                     36:     o(3) := "  phc -r : Root counting and Construction of start systems       ";
                     37:     o(4) := "  phc -m : Mixed-Volume Computation by four lifting strategies   ";
                     38:     o(5) := "  phc -p : Polynomial Continuation by a homotopy in one parameter";
                     39:     o(6) := "  phc -v : Validation, refinement and purification of solutions  ";
                     40:     o(7) := "  phc -e : SAGBI/Pieri homotopies to intersect linear subspaces  ";
                     41:     o(8) := "  phc -b : Batch or black-box processing                         ";
                     42:     for i in o'range loop
                     43:       put_line(o(i));
                     44:     end loop;
                     45:   end Display_Options;
                     46:
                     47:   procedure Read_System ( filename : in string ) is
                     48:
                     49:     file : file_type;
                     50:     n : natural;
                     51:
                     52:   begin
                     53:     if filename /= ""
                     54:      then Open_Input_File(file,filename);
                     55:           get(file,n);
                     56:           lp := new Poly_Sys(1..n);
                     57:           get(file,n,lp.all);
                     58:           Close(file);
                     59:     end if;
                     60:   exception
                     61:     when others =>
                     62:       new_line;
                     63:       put("Could not open file with name "); put_line(filename);
                     64:       lp := null; return;
                     65:   end Read_System;
                     66:
                     67: begin
                     68:   new_line;
                     69:   Display_Options;
                     70:
                     71: -- READ THE INPUT SYSTEM AND OUTPUT FILE
                     72:
                     73:   Read_System(infilename);
                     74:   if lp = null
                     75:    then new_line; get(lp);
                     76:   end if;
                     77:
                     78:   Create_Output_File(outpt,outfilename);
                     79:   put(outpt,lp'last,lp.all);
                     80:
                     81:   tstart(timer);
                     82:
                     83:   declare
                     84:
                     85:     p,q,scalp,projp : Poly_Sys(lp'range);
                     86:     target : Complex_Number;
                     87:     basis,roco : natural;
                     88:     scalvec : Link_to_Vector;
                     89:     sols : Solution_List;
                     90:     proj : boolean;
                     91:
                     92:   begin
                     93:
                     94:     p := lp.all;  Copy(p,scalp);
                     95:
                     96:    -- PREPROCESSING : SCALING AND REDUCTION
                     97:
                     98:     Driver_for_Scaling(outpt,scalp,basis,scalvec);
                     99:     Driver_for_Reduction(outpt,scalp,roco,true);
                    100:
                    101:    -- APPLY ROOT COUNTING METHODS TO CONSTRUCT A START SYSTEM
                    102:
                    103:     Copy(scalp,projp);
                    104:     Driver_for_Root_Counts(outpt,projp,q,true,sols,roco);
                    105:
                    106:     if Length_Of(sols) > 0
                    107:      then
                    108:
                    109:       -- CONSTRUCTION OF THE HOMOTOPY
                    110:
                    111:        Driver_for_Homotopy_Construction(outpt,projp,q,sols,target);
                    112:
                    113:       -- CONTINUATION
                    114:
                    115:        proj := (Number_of_Unknowns(p(p'first)) > p'last);
                    116:        if Head_Of(sols).t /= Create(0.0)
                    117:         then Set_Continuation_Parameter(sols,Create(0.0));
                    118:        end if;
                    119:        Driver_for_Polynomial_Continuation(outpt,sols,proj,target);
                    120:
                    121:       -- ROOT REFINING
                    122:
                    123:        Driver_for_Root_Refining(outpt,scalp,p,basis,scalvec,sols);
                    124:      end if;
                    125:
                    126:   end;
                    127:
                    128:   tstop(timer);
                    129:   new_line(outpt);
                    130:   print_times(outpt,timer,"solving the polynomial system");
                    131:   new_line(outpt);
                    132:   put(outpt,Bye_Bye_Message);
                    133:   Close(outpt);
                    134:
                    135: end mainphc;