| -rw-r--r-- | matrix.lua | 8 |
diff --git a/matrix.lua b/matrix.lua index b28da49e..8f7d788e 100644 --- a/matrix.lua +++ b/matrix.lua @@ -489,17 +489,17 @@ local function vector_op(scalar_op, element_wise, no_inverse) else if ra and rb then local n1, n2 = tonumber(a.size1), tonumber(b.size2) - local c = matrix_new(n1, n2) + local c = matrix_alloc(n1, n2) local NT = gsl.CblasNoTrans - gsl_check(gsl.gsl_blas_dgemm(NT, NT, 1, a, b, 1, c)) + gsl_check(gsl.gsl_blas_dgemm(NT, NT, 1, a, b, 0, c)) return c else if ra then a = mat_complex_of_real(a) end if rb then b = mat_complex_of_real(b) end local n1, n2 = tonumber(a.size1), tonumber(b.size2) - local c = matrix_cnew(n1, n2) + local c = matrix_calloc(n1, n2) local NT = gsl.CblasNoTrans - gsl_check(gsl.gsl_blas_zgemm(NT, NT, 1, a, b, 1, c)) + gsl_check(gsl.gsl_blas_zgemm(NT, NT, 1, a, b, 0, c)) return c end end |