gfortran46: Error: Type of argument 'z' in call to 'dimag' at (1) should be COMPLEX(16), not COMPLEX(8)
Anton Shterenlikht
mexas at bris.ac.uk
Thu Jul 11 08:09:29 UTC 2013
From sgk at troutmask.apl.washington.edu Thu Jul 11 02:17:08 2013
On Thu, Jul 11, 2013 at 12:15:28AM +0100, Anton Shterenlikht wrote:
> These errors are from french/aster on amd64.
> I think I need to rebuld lang/gcc after the
> recent updates to binutils, but anyway,
> which command line options to gfortran46
> would lead to this error?
>
> These routines are GNU extensions specifically for complex(8),
> e.g.: http://gcc.gnu.org/onlinedocs/gcc-4.6.4/gfortran/AIMAG.html#AIMAG
>
> Messages de compilation
> =======================
> /usr/ports/french/aster/work/aster-full-src-10.8.0/instdir/STA10.8/bibfor/algeline/dcargu.f:34.18:
>
> IF (DIMAG(C).GT.0.D0) THEN
> 1
> Error: Type of argument 'z' in call to 'dimag' at (1) should
> be COMPLEX(16), not COMPLEX(8)
Need to see the declaration for C. My guess is that it is declared
as 'complex c' or 'complex*8 c'. In either case, this is the single
precision complex type. The program is using the specific intrinsic
function dimag, which expects a double precision complex argument.
The correct fixes are to properly declare the variables via Fortran's
kind type parameter and to use generic intrinsic functions.
$ cat dcargu.f
FUNCTION DCARGU(C)
IMPLICIT REAL*8 (A-H,O-Z)
C-----------------------------------------------------------------------
C CONFIGURATION MANAGEMENT OF EDF VERSION
C MODIF ALGELINE DATE 14/11/2006 AUTEUR LEBOUVIER F.LEBOUVIER
C ======================================================================
C COPYRIGHT (C) 1991 - 2001 EDF R&D WWW.CODE-ASTER.ORG
C THIS PROGRAM IS FREE SOFTWARE; YOU CAN REDISTRIBUTE IT AND/OR MODIFY
C IT UNDER THE TERMS OF THE GNU GENERAL PUBLIC LICENSE AS PUBLISHED BY
C THE FREE SOFTWARE FOUNDATION; EITHER VERSION 2 OF THE LICENSE, OR
C (AT YOUR OPTION) ANY LATER VERSION.
C
C THIS PROGRAM IS DISTRIBUTED IN THE HOPE THAT IT WILL BE USEFUL, BUT
C WITHOUT ANY WARRANTY; WITHOUT EVEN THE IMPLIED WARRANTY OF
C MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. SEE THE GNU
C GENERAL PUBLIC LICENSE FOR MORE DETAILS.
C
C YOU SHOULD HAVE RECEIVED A COPY OF THE GNU GENERAL PUBLIC LICENSE
C ALONG WITH THIS PROGRAM; IF NOT, WRITE TO EDF R&D CODE_ASTER,
C 1 AVENUE DU GENERAL DE GAULLE, 92141 CLAMART CEDEX, FRANCE.
C ======================================================================
C-----------------------------------------------------------------------
C CALCUL DE L'ARGUMENT D'UN NOMBRE COMPLEXE
C PAR CONVENTION ON AFFECTE 0.D0 SI C = (0.D0,0.D0)
C-----------------------------------------------------------------------
C IN : C : NOMBRE COMPLEXE DONT ON VEUT CALCULER L'ARGUMENT
C-----------------------------------------------------------------------
REAL*8 DCARGU
COMPLEX*16 C
C-----------------------------------------------------------------------
PI = R8PI()
C
IF (DBLE(C).EQ.0.D0) THEN
IF (DIMAG(C).GT.0.D0) THEN
DCARGU = PI/2.D0
ELSE IF (DIMAG(C).LT.0.D0) THEN
DCARGU = -PI/2.D0
ELSE
DCARGU = 0.D0
ENDIF
ELSE IF (DBLE(C).GT.0.D0) THEN
DCARGU = DBLE(ATAN2(DIMAG(C),DBLE(C)))
ELSE IF (DBLE(C).LT.0.D0) THEN
DCARGU = DBLE(ATAN2(DIMAG(C),DBLE(C))) + PI
ENDIF
IF (DCARGU.LT.0.D0) DCARGU = DCARGU + 2.D0*PI
C
END
$
$ gfortran46 -Wall -c dcargu.f
$
So my guess is that in the port build, some other
compiler options are used that cause the failure.
I need to check this.
And I've rebuilt binutils, lang/gcc, math/metis-edf
(with the proposed patch).
Also, french/aster builds fine on ia64 -current.
Anton
More information about the freebsd-fortran
mailing list