|
From: | Patrick Alken |
Subject: | Re: [Bug-gsl] gsl_linalg_LU_invert |
Date: | Fri, 20 Apr 2012 13:46:47 -0600 |
User-agent: | Mozilla/5.0 (X11; Linux i686 on x86_64; rv:11.0) Gecko/20120327 Thunderbird/11.0.1 |
On 04/20/2012 12:03 PM, Guardincerri, Elena wrote:
my m matrix is defined as: double ma[]= {1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 1.0}; gsl_matrix_complex_view m = gsl_matrix_complex_view_array(ma, 3, 3);
Looks like a buffer overflow problem. Complex matrices have the form: [ r, i, r, i, r, i ...] so your matrix array is a factor 2 too small
[Prev in Thread] | Current Thread | [Next in Thread] |