Re: [Help-gsl] Octave-like indexing with the GSL

2006-11-16 Thread Jordi Gutierrez Hermoso
On 16/11/06, Brian Gough <[EMAIL PROTECTED]> wrote: You will have to access the memory directly for that access pattern. The gsl_matrix type is designed so that it can always be passed to a BLAS function, which means that the rows of the matrix must have a "stride" of 1. Well, then I suppose th

Re: [Help-gsl] Gsl-1.8 compiled and running under MSYS-1.10+second order derivative

2006-11-16 Thread Brian Gough
Rene Girard wrote: gcc -O2 -Wall -ansi -I/usr/local/include -c -o t_dn2.o t_dn2.c t_dn2.c: In function ‘g_x’: t_dn2.c:61: warning: dereferencing ‘void *’ pointer t_dn2.c:61: error: void value not ignored as it ought to be > double g_x(double x, void *params) { const double pi = 4.0*atan(1.0

Re: [Help-gsl] Octave-like indexing with the GSL

2006-11-16 Thread Brian Gough
Jordi Gutierrez Hermoso wrote: I was expecting to find just such a method as a matrix view with a clever little hack with the matrix's tda, but I cannot figure out how to do this. Should I just use a number of for-loops with row or column views, or use one of the gsl_matrix view functions for arr

Re: [Help-gsl] Fwd: Bug#398628: GSL has no exception handling

2006-11-16 Thread Brian Gough
Jordi Gutierrez Hermoso wrote: I am forwarding this bug from the Debian Bug Tracking System. Dirk wants to know if you recommend exceptions or not. Sorry, I don't really have any arguments one way or the other. As far as I'm aware there are no problems with -fexceptions apart from the extra s

Re: [Help-gsl] Re: Random number generators and Pthreads

2006-11-16 Thread Joseph Wakeling
Rodney Sparapani wrote: > Here's my 2 cents. Independent parallel streams is a tricky problem. > I'm not an expert in this area so rather than my advice, > take a look at JE Gentle's Random Number Generation and Monte Carlo > Methods, 2nd Ed. He mentions several possible solutions and a free > so

[Help-gsl] Re: Random number generators and Pthreads

2006-11-16 Thread Rodney Sparapani
Martin Jansche wrote: View the problem like this: you need a pool of k random words for your k threads. For debugging purposes, it makes sense if you can initialize that pool to a fixed set of values. For actual simulation runs (don't do this for cryptography), you could initialize the pool in o

Re: [Help-gsl] GSL and MSVC

2006-11-16 Thread Jordi Gutierrez Hermoso
On 16/11/06, Timothy Chan <[EMAIL PROTECTED]> wrote: I am running Microsoft Visual Studio 7.1. I am trying to run the following sample code: [snip] But I get the error: source.obj : error LNK2019: unresolved external symbol _gsl_sf_bessel_J0 referenced in function _main That's a linker error

[Help-gsl] GSL and MSVC

2006-11-16 Thread Timothy Chan
I am running Microsoft Visual Studio 7.1. I am trying to run the following sample code: #include #include int main (void) { double x = 5.0; double y = gsl_sf_bessel_J0 (x); printf ("J0(%g) = %.18e\n", x, y); return 0; } But I get the error: source.obj : error L