gsl_matrix_mul_elements bug?

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
[email protected]
https://lists.sourceforge.net/lists/listinfo/gnuwin32-users