[Rd] R 4.0.1-4.0.5 built with Intel Composer 19.0-19.1.1, errors in "make check" on CentOS 7.7-7.9

[Ivan Krylov]

On Sat, 17 Apr 2021 19:21:12 +0000
Ryan Novosielski <novosirj using rutgers.edu> wrote:

> I tried actual Pastebin and it told me that the
> reg-tests-1d.Rout.fail was offensive. https://paste.debian.net says
> it’s too large. Let’s give that another whack; both are here:
> http://www.rnovosielski.ftml.net/r-project/

This link works fine for me, thanks!

For some reason, this build of R fails a test for LAPACK accidentally
ignoring NAs when computing Frobenius norms of matrices:

https://github.com/wch/r-source/commit/db10ee5237b1f9db83a693903c4293650a43244a
https://github.com/wch/r-source/commit/2f546cf778ae3bae8ef2e82c613658c72098a528

Does the following program print NaN on the last line of its output on
your machine?

program testdlange
use ieee_arithmetic, only: ieee_value, ieee_quiet_nan

intrinsic transpose

interface
    double precision function dlange (NORM, M, N, A, LDA, WORK)
        character :: NORM
        integer :: M
        integer :: N
        double precision, dimension( lda, * ) :: A
        integer :: LDA
        double precision, dimension( * ) :: WORK
    end function
end interface

double precision, dimension(2,3) :: A
double precision :: norm

    A = reshape([0, 0, 0, 0, 0, -1], shape(A))
    print '(3f5.1)', transpose(A)

    norm = dlange('F', size(A, 1), size(A, 2), A, size(A, 1), [0D0])
    print *, norm

    A(1,2) = ieee_value(A(1,2), ieee_quiet_nan)
    print '(3f5.1)', transpose(A)

    norm = dlange('F', size(A, 1), size(A, 2), A, size(A, 1), [0D0])
    print *, norm

end program

-- 
Best regards,
Ivan