Return to ts_brackpols.adb CVS log

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

Annotation of OpenXM_contrib/PHC/Ada/Schubert/ts_brackpols.adb, Revision 1.1.1.1

1.1       maekawa     1: with text_io,integer_io;               use text_io,integer_io;
                      2: with Complex_Numbers;                  use Complex_Numbers;
                      3: with Brackets,Brackets_io;             use Brackets,Brackets_io;
                      4: with Bracket_Monomials;                use Bracket_Monomials;
                      5: with Bracket_Monomials_io;             use Bracket_Monomials_io;
                      6: with Bracket_Polynomials;              use Bracket_Polynomials;
                      7: with Bracket_Polynomials_io;           use Bracket_Polynomials_io;
                      8:
                      9: procedure ts_brackpols is
                     10:
                     11:   procedure Compare ( p : in Bracket_Polynomial; t : in Bracket_Term ) is
                     12:
                     13:   -- DESCRIPTION :
                     14:   --   Compares every monomial in p with the term t.
                     15:
                     16:     procedure Compare ( tt : in Bracket_Term; continue : out boolean ) is
                     17:     begin
                     18:       if tt < t
                     19:        then put(tt); put(" < "); put(t);
                     20:        elsif tt > t
                     21:            then put(tt); put(" > "); put(t);
                     22:            else put(tt); put(" = "); put(t);
                     23:       end if;
                     24:       new_line;
                     25:       if tt > t
                     26:        then put(tt); put(" > "); put(t);
                     27:        elsif tt < t
                     28:            then put(tt); put(" < "); put(t);
                     29:            else put(tt); put(" = "); put(t);
                     30:       end if;
                     31:       new_line;
                     32:       continue := true;
                     33:     end Compare;
                     34:     procedure Compare_Terms is new Enumerate_Terms(Compare);
                     35:
                     36:   begin
                     37:     Compare_Terms(p);
                     38:   end Compare;
                     39:
                     40:   procedure Interactive_Read ( d : in natural; t : out Bracket_Term ) is
                     41:
                     42:     m : natural;
                     43:     b : Bracket(1..d);
                     44:     s,c : integer;
                     45:     bm : Bracket_Monomial;
                     46:
                     47:   begin
                     48:     c := 1;
                     49:     put("Give the number of brackets in the monomial : "); get(m);
                     50:     for i in 1..m loop
                     51:       put("Give "); put(d,1); put(" numbers for the ");
                     52:       put(i,1); put("th bracket : "); get(b,s);
                     53:       c := s*c;
                     54:       Multiply(bm,b);
                     55:     end loop;
                     56:     put_line("The bracket monomial : "); put(bm); new_line;
                     57:     if c > 0
                     58:      then t.coeff := CMPLX(1.0);
                     59:      else t.coeff := -CMPLX(1.0);
                     60:     end if;
                     61:     t.monom := bm;
                     62:   end Interactive_Read;
                     63:
                     64:   procedure Interactive_Read ( d : in natural; p : out Bracket_Polynomial ) is
                     65:
                     66:     m : natural;
                     67:     bp : Bracket_Polynomial;
                     68:
                     69:   begin
                     70:     put_line("Reading a bracket polynomial.");
                     71:     put("Give the number of monomials : "); get(m);
                     72:     for i in 1..m loop
                     73:       declare
                     74:         bt : Bracket_Term;
                     75:       begin
                     76:         put("Reading bracket monomial "); put(i,1); new_line;
                     77:         Interactive_Read(d,bt);
                     78:         put_line("The bracket term : "); put(bt);
                     79:         Compare(bp,bt);
                     80:         Add(bp,bt);
                     81:         put_line("The bracket polynomial : "); put(bp);
                     82:       end;
                     83:     end loop;
                     84:     p := bp;
                     85:   end Interactive_Read;
                     86:
                     87:   procedure Main is
                     88:
                     89:     d : natural;
                     90:     bp : Bracket_Polynomial;
                     91:
                     92:   begin
                     93:     put("Give the number of entries in the brackets : "); get(d);
                     94:     Interactive_Read(d,bp);
                     95:     put_line("The bracket polynomial : "); put(bp);
                     96:   end Main;
                     97:
                     98: begin
                     99:   Main;
                    100: end ts_brackpols;