Help: working with large gsl_vector
Linas Vepstas
linas@austin.ibm.com
Mon Dec 13 18:41:00 GMT 2004
On Fri, Dec 10, 2004 at 11:51:19PM -0800, Toan T Nguyen was heard to remark:
>
> Hi,
>
> I'm using the multidimensional minimization procedure in GSL and have problem
> with large vectors. The dimensionality of my vectors is 300000 which means my
> index variable is of long integer type. gsl_vector, however, works only with
> size_t type only. I couldnot even call gsl_vector_alloc(300000) . It says:
double-check. size_t is normally a 32-bit int when compiling for a
32-bit libc, and 64 when compiling for a 64-bit libc. You can do
this with a
printf ("size is=%d\n", sizeof (size_t));
which will print the number of bytes.
> gsl: init_source.c:29: ERROR: vector length n must be positive integer
> Default GSL error handler invoked.
which is not to say there isn't a gsl bug, but you probably should do
more due diligence before offering the diagnosis.
--linas
More information about the Gsl-discuss
mailing list