Return to ts_durker.adb CVS log

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

Annotation of OpenXM_contrib/PHC/Ada/Root_Counts/Implift/ts_durker.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;          use Standard_Floating_Numbers;
                      3: with Standard_Floating_Numbers_io;       use Standard_Floating_Numbers_io;
                      4: with Standard_Complex_Numbers;           use Standard_Complex_Numbers;
                      5: with Standard_Complex_Numbers_io;        use Standard_Complex_Numbers_io;
                      6: with Standard_Complex_Vectors;           use Standard_Complex_Vectors;
                      7: with Standard_Random_Vectors;            use Standard_Random_Vectors;
                      8: with Durand_Kerner;
                      9:
                     10: procedure ts_durker is
                     11:
                     12: -- DESCRIPTION :
                     13: --   Test on the solver for polynomial equations in one variable.
                     14:
                     15:   procedure Read ( cv : in out Vector ) is
                     16:   begin
                     17:     for i in cv'range loop
                     18:       put(' '); put(i,1); put(" : ");
                     19:       get(cv(i));
                     20:     end loop;
                     21:   end Read;
                     22:
                     23:   procedure Write ( step : in natural; z,res : in Vector ) is
                     24:   begin
                     25:     put("Output after step  "); put(step,1); put_line(" :");
                     26:     put_line
                     27:    ("------------------------------------------------------------------------");
                     28:     put_line
                     29:    ("|    APPROXIMATED ROOTS                        |     RESIDUALS         |");
                     30:     put_line
                     31:    ("------------------------------------------------------------------------");
                     32:     for i in z'range loop
                     33:       put("| "); put(z(i)); put(" | "); put(AbsVal(res(i))); put(" |");
                     34:       new_line;
                     35:     end loop;
                     36:     put_line
                     37:    ("------------------------------------------------------------------------");
                     38:   end Write;
                     39:
                     40:   procedure dk is new Durand_Kerner(Write);
                     41:
                     42:   procedure Main is
                     43:
                     44:     n,nb : natural;
                     45:
                     46:   begin
                     47:     put("Give the degree of p(x) : "); get(n);
                     48:     declare
                     49:       p : Vector(0..n);
                     50:       z,res : Vector(1..n);
                     51:       max : natural;
                     52:       eps : double_float;
                     53:     begin
                     54:       put_line("p(x) = a_0 + a_1*x + a_2*x^2 + ... + a_n*x^n");
                     55:       put("Give "); put(n+1,1); put_line(" complex coefficients of p(x) :");
                     56:       Read(p);
                     57:       z := Random_Vector(n);    -- first approximation is random vector
                     58:       res := z;
                     59:       put("Give the maximum number of steps : "); get(max);
                     60:       put("Give the required accuracy       : "); get(eps);
                     61:       dk(p,z,res,max,eps,nb);
                     62:     end;
                     63:   end Main;
                     64:
                     65: begin
                     66:   Main;
                     67: end ts_durker;
                     68: