On Sat, 17 Apr 2021 19:21:12 +0000 Ryan Novosielski <novos...@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 ______________________________________________ R-devel@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-devel
