Andrew Davison

Thanks, fixed my problems there. On Fri, Jul 26, 2024 at 8:54 AM flexoron ***@***.***> wrote: > Yes, you don't have to install it, just make: > > https://www.gnu.org/software/gsl/ >...

Interesting, if in lu.c / gsl_linalg_LU_decomp() I change: ``` gsl_vector_uint * ipiv = gsl_vector_uint_alloc(minMN); to gsl_vector_uint * ipiv = gsl_vector_uint_calloc(minMN); ``` it works. Seems a bug in GSL 2.8 to...

The C program is simple and is allocating virgin zeroed memory from the OS. In Trealla a lot mallocs & frees have been done. On Fri, Jul 26, 2024 at...

You don't know what else changed that accesses that memory, maybe nothing but ... On Fri, Jul 26, 2024 at 11:37 AM flexoron ***@***.***> wrote: > tpl with GSL v2.7.1...

Added this to that function ``` gsl_vector_uint * ipiv = gsl_vector_uint_alloc(minMN); + for (unsigned i = 0; i < ipiv->size; i++) + printf("*** [%u] %u\n", (unsigned)i, (unsigned)ipiv->data[i]); ``` and it...

LU_decomp_L3() & apply_pivots() can't work if the vector ipiv wasn't zeroed properly, so pivi can be incorrect.

I would suggest submitting a bug report. I did one for GCC about 20 years ago, took about 6 months to get acted on.

The memory can contain anything. If you dump the newly alloc'd vector it contains rubbish, which is useless as it's used without ever being initialized. On Fri, Jul 26, 2024...

I'm talking about this... ``` gsl_vector_uint * ipiv = gsl_vector_uint_alloc(minMN); + for (unsigned i = 0; i < ipiv->size; i++) + printf("*** [%u] %u\n", (unsigned)i, (unsigned)ipiv->data[i]); ``` if the values...

I don't know what you want me to do... gsl_vector_uint_alloc returns unzeroed data (by design) when the following code expected zeroed data.