Return to ts_jaco.adb CVS log

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

Annotation of OpenXM_contrib/PHC/Ada/Math_Lib/Polynomials/ts_jaco.adb, Revision 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_Complex_Vectors;           use Standard_Complex_Vectors;
        !             6: with Standard_Complex_Vectors_io;        use Standard_Complex_Vectors_io;
        !             7: with Standard_Complex_Polynomials;       use Standard_Complex_Polynomials;
        !             8: with Standard_Complex_Polynomials_io;    use Standard_Complex_Polynomials_io;
        !             9: with Standard_Complex_Poly_Systems;      use Standard_Complex_Poly_Systems;
        !            10: with Standard_Complex_Poly_Systems_io;   use Standard_Complex_Poly_Systems_io;
        !            11: with Standard_Complex_Poly_SysFun;       use Standard_Complex_Poly_SysFun;
        !            12: with Standard_Complex_Jaco_Matrices;     use Standard_Complex_Jaco_Matrices;
        !            13: with Standard_to_Multprec_Convertors;    use Standard_to_Multprec_Convertors;
        !            14: with Multprec_Complex_Polynomials_io;    use Multprec_Complex_Polynomials_io;
        !            15: with Multprec_Complex_Poly_Systems;      use Multprec_Complex_Poly_Systems;
        !            16: with Multprec_Complex_Jaco_Matrices;     use Multprec_Complex_Jaco_Matrices;
        !            17:
        !            18: procedure ts_jaco is
        !            19:
        !            20: -- DESCRIPTION :
        !            21: --   This routine provides basic testing routines for Jacobian matrices.
        !            22:
        !            23:   procedure Test_Standard_Creation is
        !            24:
        !            25:   -- DESCRIPTION :
        !            26:   --   Tests the creation of a standard Jacobian matrix.
        !            27:
        !            28:     lp : Standard_Complex_Poly_Systems.Link_to_Poly_Sys;
        !            29:     n : natural;
        !            30:
        !            31:   begin
        !            32:     new_line;
        !            33:     put_line("Testing the creation of standard Jacobian matrices.");
        !            34:     new_line;
        !            35:     get(lp);
        !            36:     put_line("The system : "); put(lp.all);
        !            37:     n := lp'last;
        !            38:     declare
        !            39:       jm : Standard_Complex_Jaco_Matrices.Jaco_Mat(1..n,1..n) := Create(lp.all);
        !            40:     begin
        !            41:       put_line("The Jacobian matrix : ");
        !            42:       for i in jm'range loop
        !            43:         for j in jm'range loop
        !            44:           put(jm(i,j));
        !            45:         end loop;
        !            46:         new_line;
        !            47:       end loop;
        !            48:     end;
        !            49:   end Test_Standard_Creation;
        !            50:
        !            51:   procedure Test_Multprec_Creation is
        !            52:
        !            53:   -- DESCRIPTION :
        !            54:   --   Tests the creation of a multi-precision Jacobian matrix.
        !            55:
        !            56:     lp : Standard_Complex_Poly_Systems.Link_to_Poly_Sys;
        !            57:     n : natural;
        !            58:
        !            59:   begin
        !            60:     new_line;
        !            61:     put_line("Testing the creation of multi-precision Jacobian matrices.");
        !            62:     new_line;
        !            63:     get(lp);
        !            64:     put_line("The system : "); put(lp.all);
        !            65:     n := lp'last;
        !            66:     declare
        !            67:       mp : Multprec_Complex_Poly_Systems.Poly_Sys(1..n) := Convert(lp.all);
        !            68:       jm : Multprec_Complex_Jaco_Matrices.Jaco_Mat(1..n,1..n) := Create(mp);
        !            69:     begin
        !            70:       put_line("The Jacobian matrix : ");
        !            71:       for i in jm'range loop
        !            72:         for j in jm'range loop
        !            73:           put(jm(i,j));
        !            74:         end loop;
        !            75:         new_line;
        !            76:       end loop;
        !            77:     end;
        !            78:   end Test_Multprec_Creation;
        !            79:
        !            80:   procedure Main is
        !            81:
        !            82:     ans : character;
        !            83:
        !            84:   begin
        !            85:     new_line;
        !            86:     put_line("Interactive testing of the operations on complex polynomials.");
        !            87:     loop
        !            88:       new_line;
        !            89:       put_line("Choose one of the following :                              ");
        !            90:       put_line("  0. Exit this program.                                    ");
        !            91:       put_line("  1. Test creation of standard Jacobian matrices.          ");
        !            92:       put_line("  2. Test creation of multi-precision Jacobian matrices.   ");
        !            93:       put("Type 0,1 or 2 to select : "); get(ans);
        !            94:       exit when (ans = '0');
        !            95:       case ans is
        !            96:         when '1' => Test_Standard_Creation;
        !            97:         when '2' => Test_Multprec_Creation;
        !            98:         when others => null;
        !            99:       end case;
        !           100:     end loop;
        !           101:   end Main;
        !           102:
        !           103: begin
        !           104:   Main;
        !           105: end ts_jaco;