diff --git a/mex/sources/bytecode/SparseMatrix.cc b/mex/sources/bytecode/SparseMatrix.cc index 9fc01320b07ae9042a497b10aa01af39bb7892cc..1c8da9fab40bd0a0eb8ab2132b6ba44e18dc0cff 100644 --- a/mex/sources/bytecode/SparseMatrix.cc +++ b/mex/sources/bytecode/SparseMatrix.cc @@ -1796,13 +1796,12 @@ mxArray * dynSparseMatrix::mult_SAT_B(const mxArray *A_m, const mxArray *B_m) { size_t n_A = mxGetN(A_m); - size_t m_A = mxGetM(A_m); mwIndex *A_i = mxGetIr(A_m); mwIndex *A_j = mxGetJc(A_m); double *A_d = mxGetPr(A_m); size_t n_B = mxGetN(B_m); double *B_d = mxGetPr(B_m); - mxArray *C_m = mxCreateDoubleMatrix(m_A, n_B, mxREAL); + mxArray *C_m = mxCreateDoubleMatrix(n_A, n_B, mxREAL); double *C_d = mxGetPr(C_m); for (int j = 0; j < static_cast<int>(n_B); j++) for (unsigned int i = 0; i < n_A; i++) @@ -1823,14 +1822,13 @@ mxArray * dynSparseMatrix::Sparse_mult_SAT_B(const mxArray *A_m, const mxArray *B_m) { size_t n_A = mxGetN(A_m); - size_t m_A = mxGetM(A_m); mwIndex *A_i = mxGetIr(A_m); mwIndex *A_j = mxGetJc(A_m); double *A_d = mxGetPr(A_m); size_t n_B = mxGetN(B_m); size_t m_B = mxGetM(B_m); double *B_d = mxGetPr(B_m); - mxArray *C_m = mxCreateSparse(m_A, n_B, m_A*n_B, mxREAL); + mxArray *C_m = mxCreateSparse(n_A, n_B, n_A*n_B, mxREAL); mwIndex *C_i = mxGetIr(C_m); mwIndex *C_j = mxGetJc(C_m); double *C_d = mxGetPr(C_m); @@ -1868,7 +1866,6 @@ mxArray * dynSparseMatrix::Sparse_mult_SAT_SB(const mxArray *A_m, const mxArray *B_m) { size_t n_A = mxGetN(A_m); - size_t m_A = mxGetM(A_m); mwIndex *A_i = mxGetIr(A_m); mwIndex *A_j = mxGetJc(A_m); double *A_d = mxGetPr(A_m); @@ -1876,7 +1873,7 @@ dynSparseMatrix::Sparse_mult_SAT_SB(const mxArray *A_m, const mxArray *B_m) mwIndex *B_i = mxGetIr(B_m); mwIndex *B_j = mxGetJc(B_m); double *B_d = mxGetPr(B_m); - mxArray *C_m = mxCreateSparse(m_A, n_B, m_A*n_B, mxREAL); + mxArray *C_m = mxCreateSparse(n_A, n_B, n_A*n_B, mxREAL); mwIndex *C_i = mxGetIr(C_m); mwIndex *C_j = mxGetJc(C_m); double *C_d = mxGetPr(C_m);