Matrizenmultiplikation in ADA
2008-11-28 13:52:29

   function "*"   (m1 , m2 : Matrix) 
      return Matrix 
   is
      res : Matrix (m1'Range (1), m2'Range (2));
      Hilfsvar : Float := 0.0;
      Spalten : Natural := 1;
   begin 
      if 
        m1'First (2) = m2'First (1)
         and m1'Last (2) = m2'Last (2)
      then
         for K in m1'Range (2) loop
            for I in m1'Range (1) loop
               Hilfsvar := 0.0;
               for J in m2'Range (2) loop
                  Hilfsvar = Hilfsvar + m1(I,J) * m2(J,I);
               end loop;
            res (K, I) := Hilfsvar;
            end loop;
         end loop;
      else
         raise Dimension_dont_match;
      end if;
   end "*";