Hi,
Please check the following program (statical linked against libgslcblas.a,
libgsl.a - Dev-C++ 4.9.9.2):
#include <stdio.h>
#include <gsl/gsl_matrix.h>
int
main (void)
{
int i, j;
gsl_matrix * R = gsl_matrix_alloc (3, 3);
gsl_matrix_set_zero(R);
gsl_matrix_set(R, 0, 0, 1);
gsl_matrix_set(R, 0, 1, 1.25);
gsl_matrix_set(R, 1, 1, 1);
gsl_matrix_set(R, 2, 2, 2);
gsl_matrix * QDASH = gsl_matrix_alloc (3, 3);
gsl_matrix_set(QDASH, 0, 0, 1);
gsl_matrix_set(QDASH, 1, 1, 1);
gsl_matrix_set(QDASH, 2, 2, 1);
gsl_matrix_mul_elements(QDASH, R) ;
j=0;
for (i=0; i<3; ++i)
printf("%e %e %e\n", gsl_matrix_get(QDASH, i, j),
gsl_matrix_get(QDASH, i, j+1), gsl_matrix_get(QDASH, i, j+2));
gsl_matrix_free (R);
gsl_matrix_free (QDASH);
// getchar();
return 0;
}
The result of the gsl_matrix_mul_elements should be R as QDASH is matrix
identity - and is not.
Am I doing something wrong or is a bug?
Thanks,
Liviu
Send instant messages to your online friends http://au.messenger.yahoo.com
-------------------------------------------------------------------------
This SF.net email is sponsored by DB2 Express
Download DB2 Express C - the FREE version of DB2 express and take
control of your XML. No limits. Just data. Click to get it now.
http://sourceforge.net/powerbar/db2/
_______________________________________________
GnuWin32-Users mailing list
GnuWin32-Users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/gnuwin32-users
