[PATCH] more g77 runtime library functions for gfortran

[Steve Kargl <sgk at troutmask dot apl dot washington dot edu>]

This patch implements the following nonstandard intrinsic
functions included in the g77 runtime library:

  DFLOAT   -- convert an integer to double precision
  [D]BESJ0 -- Bessel function of order 0 and real argument
  [D]BESJ1 -- Bessel function of order 1 and real argument
  [D]BESJN -- Bessel function of order N and real argument
  [D]BESY0 -- Bessel function of order 0 and real argument
  [D]BESY1 -- Bessel function of order 1 and real argument
  [D]BESYN -- Bessel function of order N and real argument

The four files are

  fortran.diff     -- diff for gcc/gcc/fortran
  libgfortran.diff -- diff for gcc/libgfortran
  bessel.c         -- file containing routine;
                      goes in gcc/libgfortran/intrinsics
  bessel.f         -- testsuite program; goes in gfortran.dg

ChangeLog for gcc/gcc/fortran

2004-08-14  Steven G. Kargl  <kargls@comcast.net>
        * check.c (gfc_check_besj0, gfc_check_besjn, gfc_check_dfloat):
          New functions.
        * gfortran.h (GFC_ISYM_BESJ0, GFC_ISYM_BESJ1, GFC_ISYM_BESJN,
          GFC_ISYM_BESY0, GFC_ISYM_BESY1, GFC_ISYM_BESYN,GFC_ISYM_DFLOAT):
          New symbols.
        * intrinsic.c: Add [d]bes{j,y}{0,1,n} to list of intrinsic functions.
        * intrinsic.c: Remove "dfloat" alias, Add "dfloat" to list of
          intrinsic functions.
        * intrinsic.h: Add prototypes.
        * iresolve.c (gfc_resolve_besj0, gfc_resolve_besj1, gfc_resolve_besjn,
          gfc_resolve_besy0, gfc_resolve_besy1, gfc_resolve_besyn): New
          functions.
        * simplify.c (gfc_simplify_dfloat): New function.
        * trans-intrinsics.c: Use new symbols.

ChangeLog for gcc/libgfortran

2004-08-14  Steven G. Kargl  <kargls@comcast.net>
        * intrinsics/bessel.c: New file containing [d]bes{j,y}{0,1,n}.
        * Makefile.am: Add intrinsics/bessel.c.
        * Makefile.in: Regenerate.
        * configure.ac: Add test for {j,y}{0,1,n}[f] in libm.
        * config.h.in: #defines for HAVE_{J,Y}{0,1,n}[F].
        * configure: Regenerate.

-- 
Steve

Attachment: fortran.diff
Description: Text document

Attachment: libgfortran.diff
Description: Text document

Attachment: bessel.c
Description: Text document

Attachment: bessel.f
Description: Text document