Return to standard_random_matrices.adb CVS log

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

Annotation of OpenXM_contrib/PHC/Ada/Math_Lib/Matrices/standard_random_matrices.adb, Revision 1.1

1.1     ! maekawa     1: with Standard_Random_Numbers;            use Standard_Random_Numbers;
        !             2: with Standard_Complex_Numbers;           use Standard_Complex_Numbers;
        !             3: with Standard_Natural_Vectors;
        !             4: with Standard_Floating_Vectors;
        !             5: with Standard_Complex_Vectors;
        !             6: with Standard_Floating_QR_Decomposition; use Standard_Floating_QR_Decomposition;
        !             7: with Standard_Complex_QR_Decomposition;  use Standard_Complex_QR_Decomposition;
        !             8:
        !             9: package body Standard_Random_Matrices is
        !            10:
        !            11:   function Random_Matrix ( n,m : natural; low,upp : integer )
        !            12:                          return Standard_Integer_Matrices.Matrix is
        !            13:
        !            14:     res : Standard_Integer_Matrices.Matrix(1..n,1..m);
        !            15:
        !            16:   begin
        !            17:     for i in 1..n loop
        !            18:       for j in 1..m loop
        !            19:         res(i,j) := Random(low,upp);
        !            20:       end loop;
        !            21:     end loop;
        !            22:     return res;
        !            23:   end Random_Matrix;
        !            24:
        !            25:   function Random_Matrix ( n,m : natural )
        !            26:                          return Standard_Floating_Matrices.Matrix is
        !            27:
        !            28:     res : Standard_Floating_Matrices.Matrix(1..n,1..m);
        !            29:
        !            30:   begin
        !            31:     for i in 1..n loop
        !            32:       for j in 1..m loop
        !            33:         res(i,j) := Random;
        !            34:       end loop;
        !            35:     end loop;
        !            36:     return res;
        !            37:   end Random_Matrix;
        !            38:
        !            39:   function Orthogonalize ( mat : Standard_Floating_Matrices.Matrix )
        !            40:                          return Standard_Floating_Matrices.Matrix is
        !            41:
        !            42:     n : constant natural := mat'length(1);
        !            43:     m : constant natural := mat'length(2);
        !            44:     res : Standard_Floating_Matrices.Matrix(1..n,1..m);
        !            45:     wrk : Standard_Floating_Matrices.Matrix(1..n,1..m);
        !            46:     bas : Standard_Floating_Matrices.Matrix(1..n,1..n);
        !            47:     qra : Standard_Floating_Vectors.Vector(1..n) := (1..n => 0.0);
        !            48:     pvt : Standard_Natural_Vectors.Vector(1..n) := (1..n => 0);
        !            49:
        !            50:   begin
        !            51:     wrk := mat;
        !            52:     QRD(wrk,qra,pvt,false);
        !            53:     for i in wrk'range(1) loop
        !            54:       for j in wrk'range(2) loop
        !            55:         bas(i,j) := wrk(i,j);
        !            56:       end loop;
        !            57:       for j in m+1..n loop
        !            58:         bas(i,j) := 0.0;
        !            59:          end loop;
        !            60:     end loop;
        !            61:     Basis(bas,mat);
        !            62:     for i in res'range(1) loop
        !            63:       for j in res'range(2) loop
        !            64:         res(i,j) := bas(i,j);
        !            65:       end loop;
        !            66:     end loop;
        !            67:     return res;
        !            68:   end Orthogonalize;
        !            69:
        !            70:   function Random_Orthogonal_Matrix
        !            71:              ( n,m : natural ) return Standard_Floating_Matrices.Matrix is
        !            72:
        !            73:     res : Standard_Floating_Matrices.Matrix(1..n,1..m)
        !            74:         := Orthogonalize(Random_Matrix(n,m));
        !            75:
        !            76:   begin
        !            77:     return res;
        !            78:   end Random_Orthogonal_Matrix;
        !            79:
        !            80:   function Random_Matrix ( n,m : natural )
        !            81:                          return Standard_Complex_Matrices.Matrix is
        !            82:
        !            83:     res : Standard_Complex_Matrices.Matrix(1..n,1..m);
        !            84:
        !            85:   begin
        !            86:     for i in 1..n loop
        !            87:       for j in 1..m loop
        !            88:         res(i,j) := Random;
        !            89:       end loop;
        !            90:     end loop;
        !            91:     return res;
        !            92:   end Random_Matrix;
        !            93:
        !            94:   function Orthogonalize ( mat : Standard_Complex_Matrices.Matrix )
        !            95:                          return Standard_Complex_Matrices.Matrix is
        !            96:
        !            97:     n : constant natural := mat'length(1);
        !            98:     m : constant natural := mat'length(2);
        !            99:     res : Standard_Complex_Matrices.Matrix(1..n,1..m);
        !           100:     wrk : Standard_Complex_Matrices.Matrix(1..n,1..m);
        !           101:     bas : Standard_Complex_Matrices.Matrix(1..n,1..n);
        !           102:     qra : Standard_Complex_Vectors.Vector(1..n) := (1..n => Create(0.0));
        !           103:     pvt : Standard_Natural_Vectors.Vector(1..n) := (1..n => 0);
        !           104:
        !           105:   begin
        !           106:     wrk := mat;
        !           107:     QRD(wrk,qra,pvt,false);
        !           108:     for i in wrk'range(1) loop
        !           109:       for j in wrk'range(2) loop
        !           110:         bas(i,j) := wrk(i,j);
        !           111:       end loop;
        !           112:       for j in m+1..n loop
        !           113:         bas(i,j) := Create(0.0);
        !           114:          end loop;
        !           115:     end loop;
        !           116:     Basis(bas,mat);
        !           117:     for i in res'range(1) loop
        !           118:       for j in res'range(2) loop
        !           119:         res(i,j) := bas(i,j);
        !           120:       end loop;
        !           121:     end loop;
        !           122:     return res;
        !           123:   end Orthogonalize;
        !           124:
        !           125:   function Random_Orthogonal_Matrix
        !           126:                 ( n,m : natural ) return Standard_Complex_Matrices.Matrix is
        !           127:
        !           128:     res : Standard_Complex_Matrices.Matrix(1..n,1..m)
        !           129:         := Orthogonalize(Random_Matrix(n,m));
        !           130:
        !           131:   begin
        !           132:     return res;
        !           133:   end Random_Orthogonal_Matrix;
        !           134:
        !           135: end Standard_Random_Matrices;