Return to ts_sagbi.adb CVS log

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

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

1.1       maekawa     1: with text_io,integer_io;                 use text_io,integer_io;
                      2: with Communications_with_User;           use Communications_with_User;
                      3: with Symbol_Table;
                      4: with Standard_Complex_Numbers_io;        use Standard_Complex_Numbers_io;
                      5: with Standard_Natural_Vectors_io;        use Standard_Natural_Vectors_io;
                      6: with Standard_Complex_Vectors;           use Standard_Complex_Vectors;
                      7: with Standard_Complex_Vectors_io;        use Standard_Complex_Vectors_io;
                      8: with Standard_Floating_Matrices;         use Standard_Floating_Matrices;
                      9: with Standard_Floating_Matrices_io;      use Standard_Floating_Matrices_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 Matrix_Indeterminates;              use Matrix_Indeterminates;
                     14: with Bracket_Expansions;                 use Bracket_Expansions;
                     15: with SAGBI_Homotopies;                   use SAGBI_Homotopies;
                     16: with Driver_for_SAGBI_Homotopies;
                     17:
                     18: procedure ts_sagbi is
                     19:
                     20:   procedure Write ( p : in Poly ) is
                     21:
                     22:        procedure Write_Term ( t : in Term; cont : out boolean ) is
                     23:
                     24:     begin
                     25:       put("t.cf : "); put(t.cf); new_line;
                     26:          put("t.dg : ");
                     27:       if t.dg /= null
                     28:           then put(t.dg.all); new_line;
                     29:       end if;
                     30:       cont := true;
                     31:     end Write_Term;
                     32:     procedure Write_Terms is new Visiting_Iterator(Write_Term);
                     33:
                     34:   begin
                     35:        Write_Terms(p);
                     36:   end Write;
                     37:
                     38:   procedure Write_Coeff ( p : in Poly ) is
                     39:
                     40:     cff : constant Vector := Coeff(p);
                     41:
                     42:   begin
                     43:     put_line("The coefficient vector : ");
                     44:     put_line(cff);
                     45:   end Write_Coeff;
                     46:
                     47:   procedure Test_SAGBI_Homotopy ( file : in file_type; n,d : in natural ) is
                     48:
                     49:     p,l,p1 : Poly;
                     50:     mat : Matrix(1..n,1..d);
                     51:     ans : character;
                     52:
                     53:   begin
                     54:     loop
                     55:       Matrix_Indeterminates.Initialize_Symbols(n,d);
                     56:       p := Lifted_Localized_Laplace_Expansion(n,d);
                     57:       Symbol_Table.Replace((n-d)*d+1,"t    ");
                     58:       new_line;
                     59:       put_line("The generic polynomial : "); put(p); new_line;
                     60:       Write_Coeff(p);
                     61:       p1 := Lifted_Localized_Laplace_Expansion(n,d);
                     62:       put_line("The new generic polynomial : "); put(p1); new_line;
                     63:       Write(p1);
                     64:       put("Give a floating-point "); put(n,1); put("x"); put(d,1);
                     65:       put_line("-matrix : "); get(mat);
                     66:       l := Intersection_Condition(mat,p);
                     67:       put_line("The specific polynomial : "); put(l); new_line;
                     68:       Write_Coeff(l);
                     69:       declare
                     70:         cff : constant Vector := Intersection_Coefficients(mat,Coeff(p));
                     71:       begin
                     72:         put_line("The computed coefficient vector : ");
                     73:         put_line(cff);
                     74:       end;
                     75:       Matrix_Indeterminates.Clear_Symbols;
                     76:       put("Do you want more tests ? (y/n) "); get(ans);
                     77:       exit when ans /= 'y';
                     78:     end loop;
                     79:   end Test_SAGBI_Homotopy;
                     80:
                     81:   procedure Main is
                     82:
                     83:     n,d : natural;
                     84:     file : file_type;
                     85:
                     86:   begin
                     87:     new_line;
                     88:     put_line("SAGBI Homotopies to intersect planes in projective space.");
                     89:     new_line;
                     90:     put_line("Reading the name of the output file.");
                     91:     Read_Name_and_Create_File(file);
                     92:     new_line;
                     93:     put("Give number of elements to choose from : "); get(n);
                     94:     put("Give the number of entries in bracket : "); get(d);
                     95:     put(file,"SAGBI Homotopies for n = "); put(file,n,1);
                     96:     put(file," and d = "); put(file,d,1); new_line(file);
                     97:     Test_SAGBI_Homotopy(file,n,d);
                     98:    -- Driver_for_SAGBI_Homotopies(file,n,d);
                     99:   end Main;
                    100:
                    101: begin
                    102:   Main;
                    103: end ts_sagbi;