Return to ts_binsolve.adb CVS log

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

Annotation of OpenXM_contrib/PHC/Ada/Root_Counts/Implift/ts_binsolve.adb, Revision 1.1.1.1

1.1       maekawa     1: with text_io,integer_io;                 use text_io,integer_io;
                      2: with Standard_Floating_Numbers_io;       use Standard_Floating_Numbers_io;
                      3: with Standard_Complex_Numbers;           use Standard_Complex_Numbers;
                      4: with Standard_Complex_Numbers_io;        use Standard_Complex_Numbers_io;
                      5: with Standard_Complex_Vectors;           use Standard_Complex_Vectors;
                      6: with Standard_Complex_Vectors_io;        use Standard_Complex_Vectors_io;
                      7: with Standard_Integer_VecVecs;           use Standard_Integer_VecVecs;
                      8: with Standard_Integer_VecVecs_io;        use Standard_Integer_VecVecs_io;
                      9: with Standard_Complex_Solutions;         use Standard_Complex_Solutions;
                     10: with Standard_Complex_Solutions_io;      use Standard_Complex_Solutions_io;
                     11: with Binomial_System_Solvers;            use Binomial_System_Solvers;
                     12:
                     13: procedure ts_binsolve is
                     14:
                     15: -- DESCRIPTION :
                     16: --   Test on the solver of binomial systems.
                     17:
                     18:   ans : character;
                     19:   n : natural;
                     20:
                     21:   procedure put_Bar is
                     22:   begin
                     23:     for i in 1..50 loop
                     24:       put('*');
                     25:     end loop;
                     26:     new_line;
                     27:   end put_Bar;
                     28:
                     29:   procedure put ( n : in natural; v : in Vector; sols : in Solution_List ) is
                     30:
                     31:     tmp : Solution_List;
                     32:     s : Solution(n);
                     33:     nb : natural;
                     34:
                     35:   begin
                     36:     tmp := sols;
                     37:     nb := 1;
                     38:     while not Is_Null(tmp) loop
                     39:       s := Head_Of(tmp).all;
                     40:       put_Bar;
                     41:       put("Solution number "); put(nb,1); put_line(" :"); put(s);
                     42:       put("The residual : "); put(REAL_PART(v(nb))); new_line;
                     43:       nb := nb + 1;
                     44:       put_Bar;
                     45:       tmp := Tail_Of(tmp);
                     46:     end loop;
                     47:   end put;
                     48:
                     49: begin
                     50:   loop
                     51:     put("Give the dimension : "); get(n);
                     52:     declare
                     53:       vv : VecVec(1..n);
                     54:       cv : Standard_Complex_Vectors.Vector(1..n);
                     55:       sols : Solution_List;
                     56:     begin
                     57:       put("Give "); put(n,1); put(' '); put_line("integer vectors :");
                     58:       get(n,vv);
                     59:       put_line("Give a vector of nonzero constants : ");
                     60:       get(cv);
                     61:       Solve(vv,cv,n,sols);
                     62:       declare
                     63:        nb : natural := Length_Of(sols);
                     64:        res : Standard_Complex_Vectors.Vector(1..nb);
                     65:       begin
                     66:        Residuals(vv,cv,n,sols,res);
                     67:         put("There are "); put(nb,1); put_line(" solutions :");
                     68:         put(n,res,sols);
                     69:       end;
                     70:       Clear(sols);
                     71:     end;
                     72:     put("Do you want more tests ? (y/n) "); get(ans);
                     73:     exit when ans /= 'y';
                     74:   end loop;
                     75: end ts_binsolve;