Return to plane_representations.adb CVS log

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

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

1.1       maekawa     1: with Standard_Floating_Numbers;          use Standard_Floating_Numbers;
                      2: with Standard_Complex_Numbers;           use Standard_Complex_Numbers;
                      3:
                      4: package body Plane_Representations is
                      5:
                      6:   function Localize ( locmap : Standard_Natural_Matrices.Matrix;
                      7:                       plamat : Standard_Complex_Matrices.Matrix )
                      8:                     return Standard_Complex_Matrices.Matrix is
                      9:
                     10:     res : Standard_Complex_Matrices.Matrix(plamat'range(1),plamat'range(2));
                     11:     tol : constant double_float := 10.0**(-10);
                     12:     done_j : boolean;
                     13:
                     14:   begin
                     15:     for j in locmap'range(2) loop
                     16:       done_j := false;
                     17:       for i in locmap'range(1) loop
                     18:         if locmap(i,j) = 1
                     19:          then if AbsVal(plamat(i,j)) > tol
                     20:                then for k in plamat'range(1) loop
                     21:                       if AbsVal(plamat(k,j)) > tol
                     22:                        then res(k,j) := plamat(k,j)/plamat(i,j);
                     23:                        else res(k,j) := plamat(k,j);
                     24:                       end if;
                     25:                     end loop;
                     26:                     res(i,j) := Create(1.0);
                     27:               end if;
                     28:                      done_j := true;
                     29:          else res(i,j) := plamat(i,j);
                     30:         end if;
                     31:         exit when done_j;
                     32:       end loop;
                     33:     end loop;
                     34:     return res;
                     35:   end Localize;
                     36:
                     37:   function Vector_Rep ( plamat : Standard_Complex_Matrices.Matrix )
                     38:                       return Standard_Complex_Vectors.Vector is
                     39:
                     40:     dim : constant natural := plamat'length(1)*plamat'length(2);
                     41:     res : Standard_Complex_Vectors.Vector(1..dim);
                     42:     cnt : natural := 0;
                     43:
                     44:   begin
                     45:     for i in plamat'range(1) loop
                     46:       for j in plamat'range(2) loop
                     47:         cnt := cnt + 1;
                     48:         res(cnt) := plamat(i,j);
                     49:       end loop;
                     50:     end loop;
                     51:     return res;
                     52:   end Vector_Rep;
                     53:
                     54:   function Vector_Rep ( locmap : Standard_Natural_Matrices.Matrix;
                     55:                         plamat : Standard_Complex_Matrices.Matrix )
                     56:                       return Standard_Complex_Vectors.Vector is
                     57:
                     58:     dim : constant natural := plamat'length(1)*plamat'length(2);
                     59:     res : Standard_Complex_Vectors.Vector(1..dim);
                     60:     cnt : natural := 0;
                     61:
                     62:   begin
                     63:     for i in plamat'range(1) loop
                     64:       for j in plamat'range(2) loop
                     65:         if locmap(i,j) = 2
                     66:          then cnt := cnt + 1;
                     67:               res(cnt) := plamat(i,j);
                     68:         end if;
                     69:       end loop;
                     70:     end loop;
                     71:     return res(1..cnt);
                     72:   end Vector_Rep;
                     73:
                     74:   function Matrix_Rep ( locmap : Standard_Natural_Matrices.Matrix;
                     75:                         plavec : Standard_Complex_Vectors.Vector )
                     76:                       return Standard_Complex_Matrices.Matrix is
                     77:
                     78:     res : Standard_Complex_Matrices.Matrix(locmap'range(1),locmap'range(2));
                     79:     cnt : natural := 0;
                     80:
                     81:   begin
                     82:     for i in locmap'range(1) loop
                     83:       for j in locmap'range(2) loop
                     84:         if locmap(i,j) = 0
                     85:          then res(i,j) := Create(0.0);
                     86:          elsif locmap(i,j) = 1
                     87:              then res(i,j) := Create(1.0);
                     88:              else cnt := cnt + 1;
                     89:                   res(i,j) := plavec(cnt);
                     90:         end if;
                     91:       end loop;
                     92:     end loop;
                     93:     return res;
                     94:   end Matrix_Rep;
                     95:
                     96: end Plane_Representations;