Return to ts_fewsolve.adb CVS log

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

Annotation of OpenXM_contrib/PHC/Ada/Root_Counts/Implift/ts_fewsolve.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 Standard_Floating_Numbers;          use Standard_Floating_Numbers;
                      4: with Numbers_io;                         use Numbers_io;
                      5: with Standard_Complex_Poly_Systems;      use Standard_Complex_Poly_Systems;
                      6: with Standard_Complex_Poly_Systems_io;   use Standard_Complex_Poly_Systems_io;
                      7: with Standard_Complex_Laur_Systems;      use Standard_Complex_Laur_Systems;
                      8: with Standard_Poly_Laur_Convertors;      use Standard_Poly_Laur_Convertors;
                      9: with Standard_Complex_Solutions;         use Standard_Complex_Solutions;
                     10: with Standard_Complex_Solutions_io;      use Standard_Complex_Solutions_io;
                     11: with Standard_Root_Refiners;             use Standard_Root_Refiners;
                     12: with Transforming_Laurent_Systems;       use Transforming_Laurent_Systems;
                     13: with Fewnomial_System_Solvers;           use Fewnomial_System_Solvers;
                     14:
                     15: procedure ts_fewsolve is
                     16:
                     17: -- DESCRIPTION :
                     18: --   Test on the resolution of fewnomial systems.
                     19:
                     20:   ans : character;
                     21:   lp : Link_to_Poly_Sys;
                     22:
                     23:   procedure Call_Root_Refiner ( file : in file_type; p : in Poly_Sys;
                     24:                                 sols : in out Solution_List ) is
                     25:
                     26:     epsxa,epsfa,tolsing : double_float;
                     27:     numit : natural;
                     28:     max : constant natural := 5;
                     29:
                     30:   begin
                     31:     epsxa := 1.0E-14;
                     32:     epsfa := 1.0E-14;
                     33:     tolsing := 1.0E-08;
                     34:     Reporting_Root_Refiner(file,p,sols,epsxa,epsfa,tolsing,numit,max,false);
                     35:   end Call_Root_Refiner;
                     36:
                     37: begin
                     38:   loop
                     39:     get(lp);
                     40:     declare
                     41:       l : Laur_Sys(lp'range) := Polynomial_to_Laurent_System(lp.all);
                     42:       sols : Solution_List;
                     43:       fail : boolean;
                     44:     begin
                     45:       Shift(l);
                     46:       if not Is_Fewnomial_System(l)
                     47:        then put_line("The system has too many terms ...");
                     48:        else Solve(l,sols,fail);
                     49:             if fail
                     50:              then put_line("A failure occured ...");
                     51:              else put_line("These are the solutions : ");
                     52:                   Call_Root_Refiner(Standard_Output,lp.all,sols);
                     53:             end if;
                     54:       end if;
                     55:       Clear(sols); Clear(lp); Clear(l);
                     56:     end;
                     57:     put("Do you want more tests ? (y/n) "); Ask_Yes_or_No(ans);
                     58:     exit when ans /= 'y';
                     59:   end loop;
                     60: end ts_fewsolve;