Return to ts_expvec.adb CVS log

Up to [local] / OpenXM_contrib / PHC / Ada / Math_Lib / Polynomials

Annotation of OpenXM_contrib/PHC/Ada/Math_Lib/Polynomials/ts_expvec.adb, Revision 1.1.1.1

1.1       maekawa     1: with text_io,integer_io;                 use text_io,integer_io;
                      2: with Symbol_Table;
                      3: with Standard_Complex_Numbers;           use Standard_Complex_Numbers;
                      4: with Standard_Complex_Numbers_io;        use Standard_Complex_Numbers_io;
                      5: with Standard_Natural_Vectors;
                      6: with Standard_Complex_Vectors;           use Standard_Complex_Vectors;
                      7: with Standard_Complex_Vectors_io;        use Standard_Complex_Vectors_io;
                      8: with Standard_Integer_VecVecs;
                      9: with Standard_Integer_VecVecs_io;        use Standard_Integer_VecVecs_io;
                     10: with Standard_Complex_Polynomials;       use Standard_Complex_Polynomials;
                     11: with Standard_Complex_Polynomials_io;    use Standard_Complex_Polynomials_io;
                     12: with Standard_Complex_Poly_Functions;    use Standard_Complex_Poly_Functions;
                     13: with Exponent_Vectors;                   use Exponent_Vectors;
                     14:
                     15: procedure ts_expvec is
                     16:
                     17: -- DESCRIPTION :
                     18: --   This routine provides basic testing routines for complex polynomials.
                     19:
                     20:   procedure Read ( m : out natural;
                     21:                    q : out Standard_Complex_Polynomials.Poly ) is
                     22:
                     23:   -- DESCRIPTION :
                     24:   --   Tests the input/output of a polynomial in several variables
                     25:   --   and with complex coefficients.
                     26:
                     27:     n : natural;
                     28:     p : Standard_Complex_Polynomials.Poly;
                     29:
                     30:   begin
                     31:     put("Give the number of variables : "); get(n);
                     32:     Symbol_Table.Init(n);
                     33:     put_line("Give a polynomial (terminate with ;) : "); get(p);
                     34:     put_line("Your polynomial : "); put(p); new_line;
                     35:     Symbol_Table.Clear;
                     36:     q := p;
                     37:     m := n;
                     38:   end Read;
                     39:
                     40:   function Coeff ( p : Standard_Complex_Polynomials.Poly;
                     41:                    e : Standard_Integer_VecVecs.VecVec )
                     42:                  return Standard_Complex_Vectors.Vector is
                     43:
                     44:     res : Standard_Complex_Vectors.Vector(e'range);
                     45:     deg : Degrees := new Standard_Natural_Vectors.Vector(e(e'first)'range);
                     46:
                     47:   begin
                     48:     for i in e'range loop
                     49:       for j in e(i)'range loop
                     50:         deg(j) := e(i)(j);
                     51:       end loop;
                     52:       res(i) := Coeff(p,deg);
                     53:     end loop;
                     54:     Clear(deg);
                     55:     return res;
                     56:   end Coeff;
                     57:
                     58:   procedure Test_Eval ( n : in natural;
                     59:                         p : in Standard_Complex_Polynomials.Poly ) is
                     60:
                     61:     ev : constant Standard_Integer_VecVecs.VecVec := Create(p);
                     62:     cf : Standard_Complex_Vectors.Vector(ev'range) := Coeff(p,ev);
                     63:     x : Standard_Complex_Vectors.Vector(1..n);
                     64:     y1,y2 : Complex_Number;
                     65:     ans : character;
                     66:
                     67:   begin
                     68:     put_line("The exponent vectors : "); put(ev);
                     69:     put_line("The coefficients : ");
                     70:     for i in cf'range loop
                     71:       put(cf(i)); new_line;
                     72:     end loop;
                     73:     loop
                     74:       put("Give "); put(n,1); put_line(" complex numbers : ");
                     75:       for i in x'range loop
                     76:         get(x(i));
                     77:       end loop;
                     78:       y1 := Eval(p,x);
                     79:       y2 := Eval(ev,cf,x);
                     80:       put("Eval poly p(x) : "); put(y1); new_line;
                     81:       put("Eval cf*ev^x   : "); put(y2); new_line;
                     82:       put("Do you want more tests ? (y/n) "); get(ans);
                     83:       exit when ans /= 'y';
                     84:     end loop;
                     85:   end Test_Eval;
                     86:
                     87:   procedure Main is
                     88:
                     89:     n : natural;
                     90:     p : Poly;
                     91:
                     92:   begin
                     93:     new_line;
                     94:     put_line("Interactive testing of the operations on exponent vectors.");
                     95:     new_line;
                     96:     Read(n,p);
                     97:     Test_Eval(n,p);
                     98:   end Main;
                     99:
                    100: begin
                    101:   Main;
                    102: end ts_expvec;