Return to ts_cheby.adb CVS log

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

Annotation of OpenXM_contrib/PHC/Ada/Schubert/ts_cheby.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_Floating_Vectors;          use Standard_Floating_Vectors;
                      5: with Standard_Floating_Vectors_io;       use Standard_Floating_Vectors_io;
                      6: with Chebychev_Polynomials;              use Chebychev_Polynomials;
                      7:
                      8: procedure ts_cheby is
                      9:
                     10:   procedure Test_Create is
                     11:
                     12:     k : natural;
                     13:
                     14:   begin
                     15:     new_line;
                     16:     put_line("Testing the creation of Chebychev polynomials.");
                     17:     new_line;
                     18:     put("Give the degree k : "); get(k);
                     19:     for i in 0..k loop
                     20:       declare
                     21:         pi : constant Vector := Create(i);
                     22:       begin
                     23:         put("The coefficients of the "); put(i,1);
                     24:         put_line("-degree Chebychev polynomial : ");
                     25:         put_line(pi);
                     26:       end;
                     27:     end loop;
                     28:   end Test_Create;
                     29:
                     30:   procedure Test_Diff is
                     31:
                     32:     k : natural;
                     33:     ans : character;
                     34:
                     35:   begin
                     36:     new_line;
                     37:     put_line("Testing the differentiation of Chebychev polynomials.");
                     38:     loop
                     39:       new_line;
                     40:       put("Give the degree k : "); get(k);
                     41:       declare
                     42:         pk : constant Vector := Create(k);
                     43:         dp : constant Vector := Diff(pk);
                     44:       begin
                     45:         put("The coefficients of the "); put(k,1);
                     46:         put_line("-degree Chebychev polynomial : ");
                     47:         put_line(pk);
                     48:         put_line("The coefficients of its derivative : ");
                     49:         put_line(dp);
                     50:         for i in 1..k loop
                     51:           declare
                     52:             dpi : constant Vector := Diff(pk,i);
                     53:           begin
                     54:             put("The coefficients after deriving "); put(i,1);
                     55:             put_line(" times : "); put_line(dpi);
                     56:           end;
                     57:         end loop;
                     58:       end;
                     59:       put("Do you want more tests (y/n) "); get(ans);
                     60:       exit when (ans /= 'y');
                     61:     end loop;
                     62:   end Test_Diff;
                     63:
                     64:   procedure Test_Int is
                     65:
                     66:     k : natural;
                     67:     ans : character;
                     68:
                     69:   begin
                     70:     new_line;
                     71:     put_line("Testing the Antidifferentiation of Chebychev polynomials.");
                     72:     loop
                     73:       new_line;
                     74:       put("Give the degree k : "); get(k);
                     75:       declare
                     76:         pk : constant Vector := Create(k);
                     77:         dp : constant Vector := Int(pk);
                     78:       begin
                     79:         put("The coefficients of the "); put(k,1);
                     80:         put_line("-degree Chebychev polynomial : ");
                     81:         put_line(pk);
                     82:         put_line("The coefficients of its antiderivative : ");
                     83:         put_line(dp);
                     84:         for i in 1..k loop
                     85:           declare
                     86:             dpi : constant Vector := Int(pk,i);
                     87:           begin
                     88:             put("The coefficients after antideriving "); put(i,1);
                     89:             put_line(" times : "); put_line(dpi);
                     90:           end;
                     91:         end loop;
                     92:       end;
                     93:       put("Do you want more tests (y/n) "); get(ans);
                     94:       exit when (ans /= 'y');
                     95:     end loop;
                     96:   end Test_Int;
                     97:
                     98:   procedure Test_Eval is
                     99:
                    100:     k : natural;
                    101:     ans : character;
                    102:
                    103:   begin
                    104:     new_line;
                    105:     put_line("Testing the evaluation of Chebychev polynomials.");
                    106:     loop
                    107:       new_line;
                    108:       put("Give the degree k : "); get(k);
                    109:       declare
                    110:         p : constant Vector := Create(k);
                    111:         x : double_float;
                    112:       begin
                    113:         loop
                    114:           put("Give x : "); get(x);
                    115:           put("p(x) : "); put(Eval(p,x)); new_line;
                    116:           put("COS-ARCCOS-eval : "); put(Eval(k,x)); new_line;
                    117:           put("Do you want more evaluations ? (y/n) "); get(ans);
                    118:           exit when (ans /= 'y');
                    119:         end loop;
                    120:       end;
                    121:       put("Do you want other polynomials to evaluate ? (y/n) "); get(ans);
                    122:       exit when (ans /= 'y');
                    123:     end loop;
                    124:   end Test_Eval;
                    125:
                    126:   procedure Main is
                    127:
                    128:     ans : character;
                    129:
                    130:   begin
                    131:     new_line;
                    132:     put_line("Testing the manipulation of Chebychev polynomials.");
                    133:     loop
                    134:       new_line;
                    135:       put_line("Choose one of the following :");
                    136:       put_line("  0. Exit this program.");
                    137:       put_line("  1. Creation of Chebychev polynomials.");
                    138:       put_line("  2. Differentiation of Chebychev polynomials.");
                    139:       put_line("  3. Antidifferentiation of Chebychev polynomials.");
                    140:       put_line("  4. Evaluation of Chebychev polynomials.");
                    141:       put("Type 0,1,2,3 or 4 to make your choice : "); get(ans);
                    142:       exit when (ans = '0');
                    143:       case ans is
                    144:         when '1' => Test_Create;
                    145:         when '2' => Test_Diff;
                    146:         when '3' => Test_Int;
                    147:         when '4' => Test_Eval;
                    148:         when others => null;
                    149:       end case;
                    150:     end loop;
                    151:   end Main;
                    152:
                    153: begin
                    154:   Main;
                    155: end ts_cheby;