with text_io; use text_io;
with Ada.Numerics.Aux;       use Ada.Numerics.Aux;

procedure tlcm is 

   package double_io is new Text_IO.Float_IO (Double);
   use double_io;

  function double_rem (a : Double; b : Double) return Double is
      result : Double;
   begin
      result := a - Double'Floor (a / b) * b;
      return (result);
   end double_rem;

   function double_mod (a : Double; b : Double) return Double is
      result : Double;
   begin
      result := double'floor(a/b);
      result := a - b * result;
      return (result);
   end double_mod;

   function gcd
     (left  :  double;
      right :  double)
      return  double
   is
      min       : double:= double'Min (left, right);
      max       : double:= double'Max (left, right);
      remainder : double;
   begin
      reduce : loop
         if min <= 0.0 then
            return max;
         end if;

         remainder := double_rem(max, min);
         max       := min;
         min       := remainder;
      end loop reduce;
   end gcd;

 function lcm
     (a    : in double;
      b    : in double)
      return Double
   is
      tmp : constant Double :=
        a * b / (gcd (a, b));

   begin
      return Double'Floor (tmp);
   end lcm;

res : double;
begin
	res:=double_mod(5.0, 2.0);
	put("5 mod 2 = "); put(res,3,0); new_line;

	res:=double_mod(11.0, 3.0);
	put("11 mod 3 = "); put(res,3,0); new_line;

	res:=double_rem(5.0, 2.0);
	put("5 rem 2 = "); put(res,3,0); new_line;

	res:=double_rem(6.0, 3.0);
	put("6 rem 3 = "); put(res,3,0); new_line;

	res:=double_rem(11.0, 3.0);
	put("11 rem 3 = "); put(res,3,0); new_line;

	res:=lcm(1.0, 2.0);
	put("lcm(1,2) = "); put(res,3,0); new_line;

	res:=lcm(6.0, 7.0);
	put("lcm(6,7) = "); put(res,3,0); new_line;
end tlcm;