This is the mail archive of the gsl-discuss@sources.redhat.com mailing list for the GSL project.


Index Nav: [Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav: [Date Prev] [Date Next] [Thread Prev] [Thread Next]
Other format: [Raw text]

Re: Help: working with large gsl_vector


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


Index Nav: [Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav: [Date Prev] [Date Next] [Thread Prev] [Thread Next]