| author | Francesco Abbate <francesco.abbate@gmail.com> | 2013年08月17日 19:31:45 +0200 |
|---|---|---|
| committer | Francesco Abbate <francesco.abbate@gmail.com> | 2013年08月17日 19:44:08 +0200 |
| commit | 903757db1d6224de0e37ad4c6002e1897f7b30f7 (patch) | |
| tree | c934c4d6c7e7bd50eb3293ef415fb7f3fbcf6a86 /matrix.lua | |
| parent | f746cadb80a64ae34e2a00063eaf963e80c768a5 (diff) | |
| download | gsl-shell-903757db1d6224de0e37ad4c6002e1897f7b30f7.tar.gz | |
| -rw-r--r-- | matrix.lua | 7 |
diff --git a/matrix.lua b/matrix.lua index 1e7a86fc..42519243 100644 --- a/matrix.lua +++ b/matrix.lua @@ -93,7 +93,7 @@ local function matrix_new(n1, n2, f) end end else - ffi.fill(m.data, n1 * n2 * ffi.sizeof('double'), 0) + gsl.gsl_matrix_set_zero(m) end return m end @@ -110,7 +110,7 @@ local function matrix_cnew(n1, n2, f) end end else - ffi.fill(m.data, 2 * n1 * n2 * ffi.sizeof('double'), 0) + gsl.gsl_matrix_complex_set_zero(m) end return m end @@ -616,8 +616,7 @@ ffi.metatype(gsl_complex, complex_mt) local function matrix_new_unit(n) local m = matrix_alloc(n, n) - for k = 0, n*n - 1 do m.data[k] = 0 end - for k = 0, n-1 do m.data[k*(n+1)] = 1 end + gsl.gsl_matrix_set_identity(m) return m end |