blob: 7fe39189442c52e3dbbbfaa2bff6ef7641867938 [file] [log] [blame]
(* Copyright (C) 2009 Free Software Foundation, Inc. *)
(* This file is part of GNU Modula-2.
GNU Modula-2 is free software; you can redistribute it and/or modify it under
the terms of the GNU General Public License as published by the Free
Software Foundation; either version 2, or (at your option) any later
version.
GNU Modula-2 is distributed in the hope that it will be useful, but WITHOUT ANY
WARRANTY; without even the implied warranty of MERCHANTABILITY or
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
for more details.
You should have received a copy of the GNU General Public License along
with gm2; see the file COPYING. If not, write to the Free Software
Foundation, 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA. *)
MODULE arith5 ;
FROM libc IMPORT printf ;
FROM M2RTS IMPORT Halt ;
CONST
i = CMPLX (0.0, 1.0) ;
one = CMPLX (1.0, 0.0) ;
two = CMPLX (2.0, 0.0) ;
zero = CMPLX (0.0, 0.0) ;
three= CMPLX (3.0, 0.0) ;
foo = CMPLX (4.0, 1.0) ;
PROCEDURE Assert (b: BOOLEAN; l: CARDINAL) ;
BEGIN
IF NOT b
THEN
Halt ('failure', __FILE__, __FUNCTION__, l)
END
END Assert ;
PROCEDURE getOne () : LONGCOMPLEX ;
BEGIN
RETURN one
END getOne ;
VAR
a, b, c, d: LONGCOMPLEX ;
BEGIN
a := i ;
b := a*getOne() ;
c := i*i ;
d := -one ;
Assert(c=-one, __LINE__) ;
Assert(c=d, __LINE__) ;
printf("RE(foo) = %Lg\n", RE(foo)) ;
printf("IM(foo) = %Lg\n", IM(foo)) ;
printf("RE(i*i) = %Lg\n", RE(i*i)) ;
printf("IM(i*i) = %Lg\n", IM(i*i)) ;
printf("RE(-one) = %Lg\n", RE(-one)) ;
printf("IM(-one) = %Lg\n", IM(-one)) ;
Assert(i*i=-one, __LINE__)
END arith5.