The survival package uses a generalized cholesky decompostition throughout. If A is a
symmetric matrix A= LDL' where L is lower triangular with 1s on the diagonal, D is
diagonal, and D[i,i] =0 if column i of A is redundant. Being able to read the rank and
dependencies directly off of D is very
Dear Terry,
Why not have the required code in survival and use
importFrom(survival, gchol) in bdsmatrix? Since you maintain both,
that shouldn't be too complicated.
Best regards,
ir. Thierry Onkelinx
Statisticus / Statistician
Vlaamse Overheid / Government of Flanders
INSTITUUT VOOR NATUUR- EN
Hi Martin and Tomas,
Thanks for your reasoned replies. It seems that improving this is going to
take more effort in pinning down exactly what is appropriate than I
anticipated.
Sorry if the intention was to keep the initial discussion of this off the
list, I didn't mean to cause offence by copyi
Please pardon my ignorance, but doesn't OpenBLAS still not always play
nicely with multi-threaded OpenMP? (for example, don't race conditions
sometimes crop up)? If so, it might be nice to have the ability to
temporarily disable multi-threaded OpenMP (effectively:
omp_set_num_threads(1)) for th
Do those issues still arise when OpenBLAS is compiled with USE_OPENMP=1 ?
Keith
> On Jan 9, 2018, at 6:03 PM, Benjamin Tyner wrote:
>
> Please pardon my ignorance, but doesn't OpenBLAS still not always play nicely
> with multi-threaded OpenMP? (for example, don't race conditions sometimes
> c
I didn't do the compile; is there a way to check whether that was used?
If not, I'll inquire with our sysadmin and report back.
In any case, my suggestion was motivated by the fact that some parts of
R use OpenMP while others do not, in the hope that the former could have
their OpenBLAS omelet
Check if libopenblas is linked against libomp or libgomp.
I’d be curious to see any errors that arise when an OpenMP version of OpenBLAS
is linked with R.
Keith
> On Jan 9, 2018, at 11:01 PM, Benjamin Tyner wrote:
>
> I didn't do the compile; is there a way to check whether that was used? I
On Wed, Jan 10, 2018 at 12:04 AM, Keith O'Hara wrote:
>
> Check if libopenblas is linked against libomp or libgomp.
>
> I’d be curious to see any errors that arise when an OpenMP version of
> OpenBLAS is linked with R.
>
> Keith
>
The one time I tried compiling OpenBLAS for Windows 64 with USE O
