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
