[Bug fortran/31930] New: bes[jy]n intrinsics do not follow definition of elemental functions

[dfranke at gcc dot gnu dot org]

F95, 7.1.5, Conformability rules for elemental operations
"An elemental operation is an intrinsic operation or a defined elemental
operation. Two entities are in shape conformance if both are arrays of the same
shape, or one or both are scalars."

Thus, for bes[jy]n the following should be allowed:
 1. n scalar, x scalar
 2. n array,  x scalar
 3. n scalar, x array
 4. n array,  x array  (arrays of same shape)

$> cat besxy.f90
real :: x(2)
integer :: n(2)
x = besjn(n, x)
end

$>gfortran-svn besxy.f90
besxy.f90:3.10:

x = besjn(n, x)
         1
Error: 'ncopies' argument of 'besjn' intrinsic at (1) must be a scalar


In addition, the name of the dummy argument N should not be NCOPIES.


-- 
           Summary: bes[jy]n intrinsics do not follow definition of
                    elemental functions
           Product: gcc
           Version: 4.3.0
            Status: UNCONFIRMED
          Keywords: rejects-valid
          Severity: normal
          Priority: P3
         Component: fortran
        AssignedTo: unassigned at gcc dot gnu dot org
        ReportedBy: dfranke at gcc dot gnu dot org


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