1 ------------------------------------------------------------------------------
3 -- GNAT COMPILER COMPONENTS --
5 -- S Y S T E M . V A X _ F L O A T _ O P E R A T I O N S --
9 -- Copyright (C) 1997-2000 Free Software Foundation, Inc. --
10 -- (Version for Alpha OpenVMS) --
12 -- GNAT is free software; you can redistribute it and/or modify it under --
13 -- terms of the GNU General Public License as published by the Free Soft- --
14 -- ware Foundation; either version 2, or (at your option) any later ver- --
15 -- sion. GNAT is distributed in the hope that it will be useful, but WITH- --
16 -- OUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY --
17 -- or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License --
18 -- for more details. You should have received a copy of the GNU General --
19 -- Public License distributed with GNAT; see file COPYING. If not, write --
20 -- to the Free Software Foundation, 59 Temple Place - Suite 330, Boston, --
21 -- MA 02111-1307, USA. --
23 -- As a special exception, if other files instantiate generics from this --
24 -- unit, or you link this unit with other files to produce an executable, --
25 -- this unit does not by itself cause the resulting executable to be --
26 -- covered by the GNU General Public License. This exception does not --
27 -- however invalidate any other reasons why the executable file might be --
28 -- covered by the GNU Public License. --
30 -- GNAT was originally developed by the GNAT team at New York University. --
31 -- Extensive contributions were provided by Ada Core Technologies Inc. --
33 ------------------------------------------------------------------------------
35 with System.IO; use System.IO;
36 with System.Machine_Code; use System.Machine_Code;
38 package body System.Vax_Float_Operations is
40 -- Ensure this gets compiled with -O to avoid extra (and possibly
41 -- improper) memory stores.
43 pragma Optimize (Time);
45 -- Declare the functions that do the conversions between floating-point
46 -- formats. Call the operands IEEE float so they get passed in
49 function Cvt_G_T (X : T) return T;
50 function Cvt_T_G (X : T) return T;
51 function Cvt_T_F (X : T) return S;
53 pragma Import (C, Cvt_G_T, "OTS$CVT_FLOAT_G_T");
54 pragma Import (C, Cvt_T_G, "OTS$CVT_FLOAT_T_G");
55 pragma Import (C, Cvt_T_F, "OTS$CVT_FLOAT_T_F");
57 -- In each of the conversion routines that are done with OTS calls,
58 -- we define variables of the corresponding IEEE type so that they are
59 -- passed and kept in the proper register class.
65 function D_To_G (X : D) return G is
70 Asm ("ldg %0,%1", T'Asm_Output ("=f", A), D'Asm_Input ("m", X));
71 Asm ("cvtdg %1,%0", T'Asm_Output ("=f", B), T'Asm_Input ("f", A));
72 Asm ("stg %1,%0", G'Asm_Output ("=m", C), T'Asm_Input ("f", B));
80 function F_To_G (X : F) return G is
85 Asm ("ldf %0,%1", T'Asm_Output ("=f", A), F'Asm_Input ("m", X));
86 Asm ("stg %1,%0", G'Asm_Output ("=m", B), T'Asm_Input ("f", A));
94 function F_To_S (X : F) return S is
99 -- Because converting to a wider FP format is a no-op, we say
100 -- A is 64-bit even though we are loading 32 bits into it.
101 Asm ("ldf %0,%1", T'Asm_Output ("=f", A), F'Asm_Input ("m", X));
103 B := S (Cvt_G_T (A));
111 function G_To_D (X : G) return D is
116 Asm ("ldg %0,%1", T'Asm_Output ("=f", A), G'Asm_Input ("m", X));
117 Asm ("cvtgd %1,%0", T'Asm_Output ("=f", B), T'Asm_Input ("f", A));
118 Asm ("stg %1,%0", D'Asm_Output ("=m", C), T'Asm_Input ("f", B));
126 function G_To_F (X : G) return F is
132 Asm ("ldg %0,%1", T'Asm_Output ("=f", A), G'Asm_Input ("m", X));
133 Asm ("cvtgf %1,%0", S'Asm_Output ("=f", B), T'Asm_Input ("f", A));
134 Asm ("stf %1,%0", F'Asm_Output ("=m", C), S'Asm_Input ("f", B));
142 function G_To_Q (X : G) return Q is
147 Asm ("ldg %0,%1", T'Asm_Output ("=f", A), G'Asm_Input ("m", X));
148 Asm ("cvtgq %1,%0", Q'Asm_Output ("=f", B), T'Asm_Input ("f", A));
156 function G_To_T (X : G) return T is
160 Asm ("ldg %0,%1", T'Asm_Output ("=f", A), G'Asm_Input ("m", X));
169 function F_To_Q (X : F) return Q is
171 return G_To_Q (F_To_G (X));
178 function Q_To_F (X : Q) return F is
183 Asm ("cvtqf %1,%0", S'Asm_Output ("=f", A), Q'Asm_Input ("f", X));
184 Asm ("stf %1,%0", F'Asm_Output ("=m", B), S'Asm_Input ("f", A));
192 function Q_To_G (X : Q) return G is
197 Asm ("cvtqg %1,%0", T'Asm_Output ("=f", A), Q'Asm_Input ("f", X));
198 Asm ("stg %1,%0", G'Asm_Output ("=m", B), T'Asm_Input ("f", A));
206 function S_To_F (X : S) return F is
211 A := Cvt_T_F (T (X));
212 Asm ("stf %1,%0", F'Asm_Output ("=m", B), S'Asm_Input ("f", A));
220 function T_To_D (X : T) return D is
222 return G_To_D (T_To_G (X));
229 function T_To_G (X : T) return G is
235 Asm ("stg %1,%0", G'Asm_Output ("=m", B), T'Asm_Input ("f", A));
243 function Abs_F (X : F) return F is
248 Asm ("ldf %0,%1", S'Asm_Output ("=f", A), F'Asm_Input ("m", X));
249 Asm ("cpys $f31,%1,%0", S'Asm_Output ("=f", B), S'Asm_Input ("f", A));
250 Asm ("stf %1,%0", F'Asm_Output ("=m", C), S'Asm_Input ("f", B));
258 function Abs_G (X : G) return G is
263 Asm ("ldg %0,%1", T'Asm_Output ("=f", A), G'Asm_Input ("m", X));
264 Asm ("cpys $f31,%1,%0", T'Asm_Output ("=f", B), T'Asm_Input ("f", A));
265 Asm ("stg %1,%0", G'Asm_Output ("=m", C), T'Asm_Input ("f", B));
273 function Add_F (X, Y : F) return F is
278 Asm ("ldf %0,%1", S'Asm_Output ("=f", X1), F'Asm_Input ("m", X));
279 Asm ("ldf %0,%1", S'Asm_Output ("=f", Y1), F'Asm_Input ("m", Y));
280 Asm ("addf %1,%2,%0", S'Asm_Output ("=f", R),
281 (S'Asm_Input ("f", X1), S'Asm_Input ("f", Y1)));
282 Asm ("stf %1,%0", F'Asm_Output ("=m", R1), S'Asm_Input ("f", R));
290 function Add_G (X, Y : G) return G is
295 Asm ("ldg %0,%1", T'Asm_Output ("=f", X1), G'Asm_Input ("m", X));
296 Asm ("ldg %0,%1", T'Asm_Output ("=f", Y1), G'Asm_Input ("m", Y));
297 Asm ("addg %1,%2,%0", T'Asm_Output ("=f", R),
298 (T'Asm_Input ("f", X1), T'Asm_Input ("f", Y1)));
299 Asm ("stg %1,%0", G'Asm_Output ("=m", R1), T'Asm_Input ("f", R));
307 procedure Debug_Output_D (Arg : D) is
316 procedure Debug_Output_F (Arg : F) is
325 procedure Debug_Output_G (Arg : G) is
334 Debug_String_Buffer : String (1 .. 32);
335 -- Buffer used by all Debug_String_x routines for returning result
337 function Debug_String_D (Arg : D) return System.Address is
338 Image_String : constant String := D'Image (Arg) & ASCII.NUL;
339 Image_Size : constant Integer := Image_String'Length;
342 Debug_String_Buffer (1 .. Image_Size) := Image_String;
343 return Debug_String_Buffer (1)'Address;
350 function Debug_String_F (Arg : F) return System.Address is
351 Image_String : constant String := F'Image (Arg) & ASCII.NUL;
352 Image_Size : constant Integer := Image_String'Length;
355 Debug_String_Buffer (1 .. Image_Size) := Image_String;
356 return Debug_String_Buffer (1)'Address;
363 function Debug_String_G (Arg : G) return System.Address is
364 Image_String : constant String := G'Image (Arg) & ASCII.NUL;
365 Image_Size : constant Integer := Image_String'Length;
368 Debug_String_Buffer (1 .. Image_Size) := Image_String;
369 return Debug_String_Buffer (1)'Address;
376 function Div_F (X, Y : F) return F is
381 Asm ("ldf %0,%1", S'Asm_Output ("=f", X1), F'Asm_Input ("m", X));
382 Asm ("ldf %0,%1", S'Asm_Output ("=f", Y1), F'Asm_Input ("m", Y));
383 Asm ("divf %1,%2,%0", S'Asm_Output ("=f", R),
384 (S'Asm_Input ("f", X1), S'Asm_Input ("f", Y1)));
385 Asm ("stf %1,%0", F'Asm_Output ("=m", R1), S'Asm_Input ("f", R));
393 function Div_G (X, Y : G) return G is
398 Asm ("ldg %0,%1", T'Asm_Output ("=f", X1), G'Asm_Input ("m", X));
399 Asm ("ldg %0,%1", T'Asm_Output ("=f", Y1), G'Asm_Input ("m", Y));
400 Asm ("divg %1,%2,%0", T'Asm_Output ("=f", R),
401 (T'Asm_Input ("f", X1), T'Asm_Input ("f", Y1)));
402 Asm ("stg %1,%0", G'Asm_Output ("=m", R1), T'Asm_Input ("f", R));
410 function Eq_F (X, Y : F) return Boolean is
414 Asm ("ldf %0,%1", S'Asm_Output ("=f", X1), F'Asm_Input ("m", X));
415 Asm ("ldf %0,%1", S'Asm_Output ("=f", Y1), F'Asm_Input ("m", Y));
416 Asm ("cmpgeq %1,%2,%0", S'Asm_Output ("=f", R),
417 (S'Asm_Input ("f", X1), S'Asm_Input ("f", Y1)));
425 function Eq_G (X, Y : G) return Boolean is
429 Asm ("ldg %0,%1", T'Asm_Output ("=f", X1), G'Asm_Input ("m", X));
430 Asm ("ldg %0,%1", T'Asm_Output ("=f", Y1), G'Asm_Input ("m", Y));
431 Asm ("cmpgeq %1,%2,%0", T'Asm_Output ("=f", R),
432 (T'Asm_Input ("f", X1), T'Asm_Input ("f", Y1)));
440 function Le_F (X, Y : F) return Boolean is
444 Asm ("ldf %0,%1", S'Asm_Output ("=f", X1), F'Asm_Input ("m", X));
445 Asm ("ldf %0,%1", S'Asm_Output ("=f", Y1), F'Asm_Input ("m", Y));
446 Asm ("cmpgle %1,%2,%0", S'Asm_Output ("=f", R),
447 (S'Asm_Input ("f", X1), S'Asm_Input ("f", Y1)));
455 function Le_G (X, Y : G) return Boolean is
459 Asm ("ldg %0,%1", T'Asm_Output ("=f", X1), G'Asm_Input ("m", X));
460 Asm ("ldg %0,%1", T'Asm_Output ("=f", Y1), G'Asm_Input ("m", Y));
461 Asm ("cmpgle %1,%2,%0", T'Asm_Output ("=f", R),
462 (T'Asm_Input ("f", X1), T'Asm_Input ("f", Y1)));
470 function Lt_F (X, Y : F) return Boolean is
474 Asm ("ldf %0,%1", S'Asm_Output ("=f", X1), F'Asm_Input ("m", X));
475 Asm ("ldf %0,%1", S'Asm_Output ("=f", Y1), F'Asm_Input ("m", Y));
476 Asm ("cmpglt %1,%2,%0", S'Asm_Output ("=f", R),
477 (S'Asm_Input ("f", X1), S'Asm_Input ("f", Y1)));
485 function Lt_G (X, Y : G) return Boolean is
489 Asm ("ldg %0,%1", T'Asm_Output ("=f", X1), G'Asm_Input ("m", X));
490 Asm ("ldg %0,%1", T'Asm_Output ("=f", Y1), G'Asm_Input ("m", Y));
491 Asm ("cmpglt %1,%2,%0", T'Asm_Output ("=f", R),
492 (T'Asm_Input ("f", X1), T'Asm_Input ("f", Y1)));
500 function Mul_F (X, Y : F) return F is
505 Asm ("ldf %0,%1", S'Asm_Output ("=f", X1), F'Asm_Input ("m", X));
506 Asm ("ldf %0,%1", S'Asm_Output ("=f", Y1), F'Asm_Input ("m", Y));
507 Asm ("mulf %1,%2,%0", S'Asm_Output ("=f", R),
508 (S'Asm_Input ("f", X1), S'Asm_Input ("f", Y1)));
509 Asm ("stf %1,%0", F'Asm_Output ("=m", R1), S'Asm_Input ("f", R));
517 function Mul_G (X, Y : G) return G is
522 Asm ("ldg %0,%1", T'Asm_Output ("=f", X1), G'Asm_Input ("m", X));
523 Asm ("ldg %0,%1", T'Asm_Output ("=f", Y1), G'Asm_Input ("m", Y));
524 Asm ("mulg %1,%2,%0", T'Asm_Output ("=f", R),
525 (T'Asm_Input ("f", X1), T'Asm_Input ("f", Y1)));
526 Asm ("stg %1,%0", G'Asm_Output ("=m", R1), T'Asm_Input ("f", R));
534 function Neg_F (X : F) return F is
539 Asm ("ldf %0,%1", S'Asm_Output ("=f", A), F'Asm_Input ("m", X));
540 Asm ("cpysn %1,%1,%0", S'Asm_Output ("=f", B), S'Asm_Input ("f", A));
541 Asm ("stf %1,%0", F'Asm_Output ("=m", C), S'Asm_Input ("f", B));
549 function Neg_G (X : G) return G is
554 Asm ("ldg %0,%1", T'Asm_Output ("=f", A), G'Asm_Input ("m", X));
555 Asm ("cpysn %1,%1,%0", T'Asm_Output ("=f", B), T'Asm_Input ("f", A));
556 Asm ("stg %1,%0", G'Asm_Output ("=m", C), T'Asm_Input ("f", B));
564 procedure pd (Arg : D) is
566 Put_Line (D'Image (Arg));
573 procedure pf (Arg : F) is
575 Put_Line (F'Image (Arg));
582 procedure pg (Arg : G) is
584 Put_Line (G'Image (Arg));
591 function Sub_F (X, Y : F) return F is
596 Asm ("ldf %0,%1", S'Asm_Output ("=f", X1), F'Asm_Input ("m", X));
597 Asm ("ldf %0,%1", S'Asm_Output ("=f", Y1), F'Asm_Input ("m", Y));
598 Asm ("subf %1,%2,%0", S'Asm_Output ("=f", R),
599 (S'Asm_Input ("f", X1), S'Asm_Input ("f", Y1)));
600 Asm ("stf %1,%0", F'Asm_Output ("=m", R1), S'Asm_Input ("f", R));
608 function Sub_G (X, Y : G) return G is
613 Asm ("ldg %0,%1", T'Asm_Output ("=f", X1), G'Asm_Input ("m", X));
614 Asm ("ldg %0,%1", T'Asm_Output ("=f", Y1), G'Asm_Input ("m", Y));
615 Asm ("subg %1,%2,%0", T'Asm_Output ("=f", R),
616 (T'Asm_Input ("f", X1), T'Asm_Input ("f", Y1)));
617 Asm ("stg %1,%0", G'Asm_Output ("=m", R1), T'Asm_Input ("f", R));
621 end System.Vax_Float_Operations;