[R] solve() function freezes CLI in GNU R 3.6.3

[Rasmus Liland]

On 2020-05-09 11:40 -0400, J C Nash wrote:
> 
> > solve(D)
>      [,1] [,2]
> [1,] -2.0  1.0
> [2,]  1.5 -0.5
> > D %*% solve(D)
>      [,1]         [,2]
> [1,]    1 1.110223e-16
> [2,]    0 1.000000e+00
> >

Dear list,

I get another solution on my Linux i7-7500U 
laptop, but the same solution on my FreeBSD 
E3-1240Lv5 machine with a really old R 
version (without BLAS) ...

How is this possible?

> D %*% solve(D)
             [,1] [,2]
[1,] 1.000000e+00    0
[2,] 8.881784e-16    1
> sessionInfo()
R version 3.6.3 (2020-02-29)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Arch Linux

Matrix products: default
BLAS:   /usr/lib/libopenblasp-r0.3.9.so
LAPACK: /usr/lib/liblapack.so.3.9.0

locale:
 [1] LC_CTYPE=en_GB.UTF-8       LC_NUMERIC=C              
 [3] LC_TIME=en_DK.UTF-8        LC_COLLATE=en_GB.UTF-8    
 [5] LC_MONETARY=nb_NO.UTF-8    LC_MESSAGES=en_GB.UTF-8   
 [7] LC_PAPER=nb_NO.UTF-8       LC_NAME=C                 
 [9] LC_ADDRESS=C               LC_TELEPHONE=C            
[11] LC_MEASUREMENT=nb_NO.UTF-8 LC_IDENTIFICATION=C       

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

loaded via a namespace (and not attached):
[1] compiler_3.6.3
> 

>From the machine running FreeBSD:

> D %*% solve(D)
     [,1]         [,2]
[1,]    1 1.110223e-16
[2,]    0 1.000000e+00
> sessionInfo()
R version 3.5.2 (2018-12-20)
Platform: amd64-portbld-freebsd12.0 (64-bit)
Running under: FreeBSD hmm 12.0-RELEASE FreeBSD 12.0-RELEASE r341666 GENERIC  amd64

Matrix products: default
LAPACK: /usr/local/lib/R/lib/libRlapack.so

locale:
[1] C/en_GB.UTF-8/en_GB.UTF-8/C/en_GB.UTF-8/en_GB.UTF-8

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base

loaded via a namespace (and not attached):
[1] compiler_3.5.2
>

Best,
Rasmus