Return to ts_mathom.adb CVS log

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

Annotation of OpenXM_contrib/PHC/Ada/Schubert/ts_mathom.adb, Revision 1.1

1.1     ! maekawa     1: with text_io,integer_io;                     use text_io,integer_io;
        !             2: with Standard_Complex_Numbers;               use Standard_Complex_Numbers;
        !             3: with Standard_Complex_Numbers_io;            use Standard_Complex_Numbers_io;
        !             4: with Standard_Complex_Matrices;              use Standard_Complex_Matrices;
        !             5: with Standard_Complex_Matrices_io;           use Standard_Complex_Matrices_io;
        !             6: with Matrix_Homotopies;
        !             7:
        !             8: procedure ts_mathom is
        !             9:
        !            10: -- DESCRIPTION :
        !            11: --   Test of the package Matrix_Homotopies.
        !            12:
        !            13:   procedure Read is
        !            14:
        !            15:     n,m : natural;
        !            16:
        !            17:   begin
        !            18:     new_line;
        !            19:     put_line("Adding a new matrix homotopy");
        !            20:     new_line;
        !            21:     put("Give number of rows : "); get(n);
        !            22:     put("Give number of columns : "); get(m);
        !            23:     declare
        !            24:          start,target : Matrix(1..n,1..m);
        !            25:     begin
        !            26:       put("Give "); put(n,1); put("x"); put(m,1);
        !            27:       put_line(" complex matrix as start : "); get(start);
        !            28:       put("Give "); put(n,1); put("x"); put(m,1);
        !            29:       put_line(" complex matrix as target : "); get(target);
        !            30:       Matrix_Homotopies.Add(start,target);
        !            31:       put_line("Matrix homotopy added...");
        !            32:        end;
        !            33:   end Read;
        !            34:
        !            35:   procedure Modify is
        !            36:
        !            37:     k,n,m : natural;
        !            38:     ans : character;
        !            39:
        !            40:   begin
        !            41:     new_line;
        !            42:     put_line("Modifying or adding a matrix homotopy.");
        !            43:     new_line;
        !            44:     put("Give the number of the map : "); get(k);
        !            45:     put("Give number of rows : "); get(n);
        !            46:     put("Give number of columns : "); get(m);
        !            47:     declare
        !            48:          start,target : Matrix(1..n,1..m);
        !            49:     begin
        !            50:       put("Do you want to modify the start ? (y/n) "); get(ans);
        !            51:       if ans = 'y'
        !            52:        then put("Give "); put(n,1); put("x"); put(m,1);
        !            53:             put_line(" complex matrix as start : "); get(start);
        !            54:             Matrix_Homotopies.Add_Start(k,start);
        !            55:       end if;
        !            56:       put("Do you want to modify the target ? (y/n) "); get(ans);
        !            57:       if ans = 'y'
        !            58:        then put("Give "); put(n,1); put("x"); put(m,1);
        !            59:             put_line(" complex matrix as target : "); get(target);
        !            60:             Matrix_Homotopies.Add_Target(k,target);
        !            61:       end if;
        !            62:       put_line("Matrix homotopy modified...");
        !            63:        end;
        !            64:   end Modify;
        !            65:
        !            66:   procedure Eval is
        !            67:
        !            68:     k : natural;
        !            69:     t : Complex_Number;
        !            70:     ans : character;
        !            71:
        !            72:   begin
        !            73:     new_line;
        !            74:     put_line("Evaluating a matrix homotopy.");
        !            75:     new_line;
        !            76:     loop
        !            77:       put("Give the number of the map : "); get(k);
        !            78:       loop
        !            79:         put("Give the continuation parameter t : "); get(t);
        !            80:         declare
        !            81:           eva : constant Matrix := Matrix_Homotopies.Eval(k,t);
        !            82:         begin
        !            83:           put_line("The evaluated map : "); put(eva);
        !            84:         end;
        !            85:         put("Do you want to evaluate this maps for other values ? (y/n) ");
        !            86:         get(ans);
        !            87:         exit when (ans /= 'y');
        !            88:       end loop;
        !            89:       put("Do you want to evaluate other maps ? (y/n) "); get(ans);
        !            90:       exit when (ans /= 'y');
        !            91:     end loop;
        !            92:   end Eval;
        !            93:
        !            94:   procedure Main is
        !            95:
        !            96:     nb : natural;
        !            97:     ans : character;
        !            98:
        !            99:   begin
        !           100:     new_line;
        !           101:     put_line("Testing the matrix homotopies.");
        !           102:     new_line;
        !           103:     put("Give the number of matrix homotopies : "); get(nb);
        !           104:     Matrix_Homotopies.Init(nb);
        !           105:        loop
        !           106:       new_line;
        !           107:       put_line("Choose one of the following : ");
        !           108:       put_line("  0. Exit this program;");
        !           109:       put_line("  1. Read and Add a new matrix homotopy;");
        !           110:       put_line("  2. Modify or Add a matrix homotopy;");
        !           111:       put_line("  3. Evaluate a matrix homotopy.");
        !           112:       put("Type 0,1,2 or 3 to make your choice : "); get(ans);
        !           113:       exit when (ans = '0');
        !           114:       case ans is
        !           115:         when '1' => Read;
        !           116:         when '2' => Modify;
        !           117:         when '3' => Eval;
        !           118:         when others => null;
        !           119:       end case;
        !           120:     end loop;
        !           121:     Matrix_Homotopies.Clear;
        !           122:   end Main;
        !           123:
        !           124: begin
        !           125:   Main;
        !           126: end ts_mathom;