This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
Re: Help: working with large gsl_vector
- From: Linas Vepstas <linas at austin dot ibm dot com>
- To: Toan T Nguyen <ntt at physics dot ucla dot edu>
- Cc: gsl-discuss at sources dot redhat dot com
- Date: Mon, 13 Dec 2004 12:41:00 -0600
- Subject: Re: Help: working with large gsl_vector
- References: <200412102351.20027.ntt@physics.ucla.edu>
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