I didn't run your code, but is the difference simply a matter of
scaling? There is some discussion of this in the documentation.
On Mon, Jan 11, 2010 at 4:54 AM, Sumit Adhikari
wrote:
> A small correction ... "not returning"
> Sorry for the inconvenience.
>
>
> On Mon, Jan 11, 2010 at 3:22 PM, S
If you want to multiply matrices that the user has typed in, your
program needs to look something like:
1. allocate storage for each of A, B and C matrices (gsl_matrix_alloc
for example)
2. ask the user to type in values, and store them into the matrices A and B
3. call the gsl_cblas_dgemm functio
Are you sure to have thought through whether it is a good idea to port
GSL to the CUDA platform? To me it sounds like a bad idea for two
main reasons:
1) GSL's data types are not opaque, so you can't very well "hide" the
fact that data lives on a GPU. Anyone who tries to access his vector
or mat
I have got the alsa stuff working here. It's first-draft code, and could use
some reorganizing and polishing, but the technique is there and it works.
Using this code implies several new possibilities:
- We can disconnect from alsa without changing anything in csound... we don't
need to kill
On Sat, Nov 11, 2006 at 06:01:25PM -0500, Jordi Gutierrez Hermoso wrote:
> Hello everyone.
>
> I'm writing some simple C++ wrappers for the GSL for personal use, and
> I've already hit a bit of a stump at an early stage. I want to able to
> index my wrapper class with Octave-like syntax. E.g. if A
On Sat, Oct 07, 2006 at 08:30:50PM +0200, Ivan Liu wrote:
> Hi,
>
> I hope to use to declare a gsl_function by pointer, and
> use the integration rountine gsl_integration_qag in the
> following way:
>
>
> gsl_function * Integrand=0;
> Integrand->function = &f;
> Integrand->params = &list;
>
> g
On Fri, Aug 04, 2006 at 02:06:39PM +0100, Wei Cheng wrote:
> Hi,
>
> The cholesky decomposition and LU decomposition in GSL linear algerbra
> section, is only suitable for small matrix, for large matrix LAPACK is
> better as stated on GSL website. I am wondering if there is anyone who has
> experi
On Mon, Aug 07, 2006 at 12:11:30PM +0100, Travers, John C wrote:
> Hi,
>
> I'm about to convert a MATLAB program to C using GSL and FFTW. It is
> essential that I use FFTW as most of the work is done with FFTs. I've
> read through the list archives and browsed around with Google to see if
> anyone
On Fri, Aug 04, 2006 at 02:06:39PM +0100, Wei Cheng wrote:
> Hi,
>
> The cholesky decomposition and LU decomposition in GSL linear algerbra
> section, is only suitable for small matrix, for large matrix LAPACK is
> better as stated on GSL website. I am wondering if there is anyone who has
> experi
On Sat, Jul 15, 2006 at 06:11:54PM +, Mohamed Ali Charmi wrote:
> Is it possible to find all the roots of a function
In general, this is difficult, and the answer is no. Thinking about
your problem in particular might help... here are some quickly-generated
examples of questions that might be
On Sat, Jul 08, 2006 at 08:43:54AM +0200, Petr Ent wrote:
> hi,
> i am pretty new user of gsl, so please sory for this question...
> what is the easiest\fastest way to fit few points with an exponential
> function? i found example of non-linear least squares fitting, but it seems
> bit slow and
On Thu, Jul 06, 2006 at 04:02:06PM +0100, Wei Cheng wrote:
> I know it would be more efficient to do it in BLAS way, but , the reason i
> need the matrix matrix multiplication function to return a matrix result is
> because I have very complicated matrix operations, this matrix product is just
> a
On Thu, Jul 06, 2006 at 03:05:23PM +0100, Wei Cheng wrote:
> I have a matrix-matrix multiplication function based on BLAS matrix-matrix
> mutiplication function,just for the ease of my own use as follow:
>
> //**
>gsl_
On Fri, Jun 30, 2006 at 09:21:58AM +0200, Eric Germaneau wrote:
> I wish to use gsl in oder to perform the derivative of a data set
> contains in an one dimensional array.
> It seems that gsl numerical differentiation estimate the derivative of
> function only.
Indeed, I'm confused... you would l
> On 6/21/06, James Bergstra <[EMAIL PROTECTED]> wrote:
> >
> >On Wed, Jun 21, 2006 at 03:58:50PM -0500, A J wrote:
> >> Could there be an issue using gsl_matrix_fprintf with a gsl_matrix which
> >> doesn't own it's block? I tried doing thi
On Wed, Jun 21, 2006 at 03:58:50PM -0500, A J wrote:
> Could there be an issue using gsl_matrix_fprintf with a gsl_matrix which
> doesn't own it's block? I tried doing this twice:
>
> gsl_matirx_frprintf(f1,X,"%.8g");
>
> In the the first attempt, X owned it's block at the time the function was
>
On Fri, Jun 16, 2006 at 04:23:14PM +0200, Thomas Weber wrote:
> Hi,
>
> Am Donnerstag, den 15.06.2006, 12:10 +0200 schrieb Emiliano Olivieri:
> > double lower,upper;
> > gsl_histogram_get_range (h,gsl_histogram_max_bin (h),&lower,&upper);
> >
> >
> > In this other case, my programs stop with th
deallocates the block itself, but also the data that it
encapsulates. Your matrix buffer was free'd on line 2, so that's why you get a
double free in the else block.
--
James Bergstra
http://www-etud.iro.umontreal.ca/~bergstrj
___
Help-gsl mailing list
Help-gsl@gnu.org
http://lists.gnu.org/mailman/listinfo/help-gsl
hat the todo list has:
- implement revised line search as an alternative implementation of
"directional_minimize.c" ?
- add testing framework for minimization algos
- add rosenbrock to test cases
--
James Bergstra
http://www-etud.iro.umontreal.ca/~bergstrj
_
at facilitate solving A\B. See the
linalg section for docs and examples. There's work going on right now to add QR
factoring to linalg, but I think it's not done yet.
--
James Bergstra
http://www-etud.iro.umontreal.ca/~bergstrj
___
Help-gsl mailing list
Help-gsl@gnu.org
http://lists.gnu.org/mailman/listinfo/help-gsl
a certain physical location in memory,
relative to the first row. This makes many operations fast and easy, but
row-swapping isn't one of them.
--
James Bergstra
http://www-etud.iro.umontreal.ca/~bergstrj
___
Help-gsl mailing list
Help-gsl@gnu.
ip
gsl_matrix_free(B); //B is free'd, but his block isn't
--
James Bergstra
http://www-etud.iro.umontreal.ca/~bergstrj
___
Help-gsl mailing list
Help-gsl@gnu.org
http://lists.gnu.org/mailman/listinfo/help-gsl
f you like, I can set you up with
a project in libmsl (hosted on savannah). It, already has cvs space, autotools
are already configured, etc. Let me know if you like my design idea, and if
you'd like admin status in libmsl. I'm supposed to be doing a masters thesis
too, but it doesn&#x
> i wanna know how to do the sum of vector elements using BLAS.
You actually have to do cblas_ddot with a vector of ones.
This is one of the motivations for supporting the BLAS Technical forum
standard :(
--
James Bergstra
http://www-etud.iro.umontreal.ca/~bergs
A, m)
At the end of this procedure M stores inv(sigma).
--
james bergstra
http://www-etud.iro.umontreal.ca/~bergstrj
___
Help-gsl mailing list
Help-gsl@gnu.org
http://lists.gnu.org/mailman/listinfo/help-gsl
libmsl
package. If you want, check out the implementation of the layer_gaussian
routine:
http://www-etud.iro.umontreal.ca/~bergstrj/msl/html/a4.html#a14
The catch is that layer_gaussian expects that your covariance matrix is already
factored on input, eg. it expects L as i
tion command to read something more
like:
gcc 1.c -o 1.out -L/path/to/libgsl -lgsl -lgslcblas
In passing, these are generic instructions on using gcc, so if you have the same
problem with other libraries, try the same fix.
Good luck
--
james bergstra
ht
ctorA)). How can I do this without
> using loops?
I've written an extension that has exp, and other math.h-ish functions for
vectors and
matrices.
http://www-etud.iro.umontreal.ca/~bergstrj/msl/html/index.html
You can try your luck compiling the libs in cvs, or else you can just steal
ariable values stops changing it enters the my_f
> function succesively more no. of times for each iteration. ... why should it
> be like that ??
Good question. I have almost no idea, but you could look at the references in
the manual, or even just google for conjugate-gradient. Th
posting, I have a more
> complete
> test program which I am trying to get working.
>
> Thanks again.
>
> -dave miller
>
>
> On Tue, 20 Dec 2005 15:59:26 -0600, James Bergstra
> <[EMAIL PROTECTED]> wrote:
>
> >>> >I need to fit a Gaus
> I need to fit a Gaussian function to an array of data points.
What is the nature of your data? Are they just X values or X,Y pairs?
--
james bergstra
http://www-etud.iro.umontreal.ca/~bergstrj
___
Help-gsl mailing list
Help-gsl@gnu.org
h
what you then need to do is
>
> splot "< awk -f addblanks.awk Datafile"
>
> where the matrix data is stored in Datafile as above but without the
> need for blank lines.
>
>
> Hope that helps,
>
> Peter Melchior
>
>
>
x
> gsl_matrices, is it (which would probably solve my above problem)?
>
> Regards
> Thomas
>
>
>
>
> ___
> Help-gsl mailing list
> Help-gsl@gnu.org
> http://lists.gnu.or
aid in my original post, that the function I was hoping to call
was matlab's imagesc() that renders matrices as pixel maps.
Is Gnuplot the right tool for the job?
--
james bergstra
http://www-etud.iro.umontreal.ca/~bergstrj
___
Help-gsl mailing
Dimitris' posting about 3D graphics inspires me to ask a related
question...
What are some good visualization tools for vectors and matrices?
I miss the plotting tools in matlab, what can I do instead?
James
--
james bergstra
http://www-etud.iro.umontreal.ca/~ber
always assumed to be 1, but
other than that, the things are analogous. For some reason, a matrix stride
(from row to row) is in a field called 'tda'.
James Bergstra
On Thu, Dec 08, 2005 at 02:17:26AM +0200, Dimitris Papavasiliou wrote:
> Hi all,
>
> I need to quickly and e
Vancouver, B.C.
> Canada V6T 2A3
> 1-604-222-1047, ext. 6194
> --
>
>
>
> ___
> Help-gsl mailing list
> Help-gsl@gnu.org
> http://lists.gnu.org/mailman/
t; This is somewhat related to Part 1. If I look at the convergence by
> >the value chi squared takes per iteration, a typical [example] picture
> >looks like this:
> >50->30->20.1->20.0->19.99->19.981->19.9809-> >for some 100s of iterations>->19.97->19.5->19.0->18.3->
>
l or privileged. Any unauthorised use or disclosure is
> prohibited. If you have received this e-mail in error, please delete it
> immediately and notify Kane Usher on +61 7 3327 4464.
>
>
> ___
> Help-gsl mailing list
> Help
manage to increase an alredy existing gsl_vector?
> (I mean, without using plain C arrays).
>
> Giulio.
> ___
> Help-gsl mailing list
> Help-gsl@gnu.org
> http://lists.gnu.org/mailman/listinfo/help-gsl
--
james bergstra
http://www-et
* mat = gsl_matrix_alloc(X->size, X->size);
gsl_vector_view diag = gsl_matrix_diagonal(mat);
gsl_matrix_set_all(mat, 0.0); //or whatever number you like
gsl_vector_memcpy(&diag.vector, X);
return mat;
}
--
james bergstra
http://www-etud.iro.umontr
g list
> Help-gsl@gnu.org
> http://lists.gnu.org/mailman/listinfo/help-gsl
--
james bergstra
http://www-etud.iro.umontreal.ca/~bergstrj
___
Help-gsl mailing list
Help-gsl@gnu.org
http://lists.gnu.org/mailman/listinfo/help-gsl
ou Yahoo!?
> >Tired of spam? Yahoo! Mail has the best spam protection around
> >http://mail.yahoo.com
> >
> >
> >___
> >Help-gsl mailing list
> >Help-gsl@gnu.org
> >http://lists.gnu.org/mailman/listinfo/help-gsl
> >
>
> "Home is not where you are born, but where your heart finds peace" -
> Tommy Nordgren, "The dying old crone"
>
>
>
> ___
> Help-gsl mailing list
> Help-gsl@gnu.org
> http://lists.gnu.org/mailman/listinfo/help-gsl
--
james bergstra
http://www-etud.iro.umontreal.ca/~bergstrj
___
Help-gsl mailing list
Help-gsl@gnu.org
http://lists.gnu.org/mailman/listinfo/help-gsl
_
> Help-gsl mailing list
> Help-gsl@gnu.org
> http://lists.gnu.org/mailman/listinfo/help-gsl
--
james bergstra
http://www-etud.iro.umontreal.ca/~bergstrj
___
Help-gsl mailing list
Help-gsl@gnu.org
http://lists.gnu.org/mailman/listinfo/help-gsl
44 matches
Mail list logo
