[Bug fortran/36096] New: F2008 Bessel: Documentation/diagnostic errors

["burnus at gcc dot gnu dot org" <gcc-bugzilla at gcc dot gnu dot org>]

Found at:
http://groups.google.com/group/comp.lang.fortran/browse_thread/thread/e9add97708681397

Please re-check the thread, I have likely missed something.


BESSEL_J0, BESSEL_J1, BESSEL_Y0, and BESSEL_Y1 are elemental, but the manual
says otherwise.http://gcc.gnu.org/onlinedocs/gfortran/Intrinsic-Procedures.html

Note by Dan:
| The j0, j1, y0, y1 functions are elemental.
| The confusion come from the jn and yn functions.
| 
| There are two versions of these.  One is elemental.
| That one is the straight y( :) = jn( n, x( :)) version.
| 
| The transformational version is j/yn( n1, n2, x) version.
| This returns all j/y's between order n1 and order n2 inclusive. 


Furthermore, James Van Buskirk writes:

So gfortran treats the elemental case of BESSEL_J0 elementally
even though its documentation says that X shall be scalar, and
it treats the elemental cases of BESSEL_JN elementally even though
its docs says N and X shall be scalar, and N1723.pdf says the result
of BESSEL_JN(N,X) is scalar.

gfortran doesn't document nor does it implement the transformational
flavor of BESSEL_JN(N1,N2,X), partly because the N1723 documentation
is vague (I think) and partly because transformational intrinsics
are a royal PITB.


The manual says DFLOAT needs an INTEGER actual argument but in fact
it will takes REAL*4, REAL*8, or REAL*10.  LOG_GAMMA was out of
alphabetic order in the manual.  Is it true the gfortran's
implementation of NINT doesn't take a KIND optional argument?
Amusing that there is an INT2 and INT8 but no INT1 nor INT4.  I
suppose you don't want to add extensions like that if you don't
really have to.  Similar with SNGL accepting REAL*8 but not
REAL*10 input.

I didn't see what the default KIND for the result of CEILING was
in the manual.

ERF, ERFC, GAMMA, and LOG_GAMMA returned NaN for REAL*10 inputs in
my test.  EXPONENT returned garbage for REAL*10 and FRACTION
return its input for REAL*10.  SPACING and RRSPACING were hosed in
my test so that I didn't capture any output.

Test is at http://home.comcast.net/~kmbtib/Fortran_stuff/funr1.ZIP
Oops, I forgot to include a README: you can run the tests by
deleting funr1_gfc.out and funr1_gfc.err and then running
funr1_gfc.bat.  If not on Windows, you may have to edit those
*.bat files to make it all go. :)


-- 
           Summary: F2008 Bessel: Documentation/diagnostic errors
           Product: gcc
           Version: 4.4.0
            Status: UNCONFIRMED
          Keywords: documentation, diagnostic
          Severity: normal
          Priority: P3
         Component: fortran
        AssignedTo: unassigned at gcc dot gnu dot org
        ReportedBy: burnus at gcc dot gnu dot org


http://gcc.gnu.org/bugzilla/show_bug.cgi?id=36096