summaryrefslogtreecommitdiff
path: root/numpy/linalg/umath_linalg.c.src
diff options
context:
space:
mode:
authoralex <argriffi@ncsu.edu>2015-05-21 13:57:46 -0400
committeralex <argriffi@ncsu.edu>2015-05-21 13:57:46 -0400
commitb9f5e85659a79929ccc811498e931608af58c4b7 (patch)
tree70bcc5202a7c41a62405aa4a151dffdba8489d1c /numpy/linalg/umath_linalg.c.src
parentad4aa2555010c73cc0aebb35e09ca7a10e7d956e (diff)
downloadnumpy-b9f5e85659a79929ccc811498e931608af58c4b7.tar.gz
MAINT: more consistent notation in umath_linalg
Diffstat (limited to 'numpy/linalg/umath_linalg.c.src')
-rw-r--r--numpy/linalg/umath_linalg.c.src18
1 files changed, 13 insertions, 5 deletions
diff --git a/numpy/linalg/umath_linalg.c.src b/numpy/linalg/umath_linalg.c.src
index 00f69a03e..60cada325 100644
--- a/numpy/linalg/umath_linalg.c.src
+++ b/numpy/linalg/umath_linalg.c.src
@@ -2464,7 +2464,8 @@ init_@lapack_func@(GESDD_PARAMS_t *params,
fortran_int min_m_n = m<n?m:n;
size_t safe_min_m_n = min_m_n;
size_t s_size = safe_min_m_n * sizeof(@ftyp@);
- fortran_int u_row_count, vt_column_count;
+ fortran_int u_row_count, vt_column_count;
+ size_t safe_u_row_count, safe_vt_column_count;
size_t u_size, vt_size;
fortran_int work_count;
size_t work_size;
@@ -2473,8 +2474,11 @@ init_@lapack_func@(GESDD_PARAMS_t *params,
if (!compute_urows_vtcolumns(jobz, m, n, &u_row_count, &vt_column_count))
goto error;
- u_size = ((size_t)u_row_count) * safe_m * sizeof(@ftyp@);
- vt_size = safe_n * ((size_t)vt_column_count) * sizeof(@ftyp@);
+ safe_u_row_count = u_row_count;
+ safe_vt_column_count = vt_column_count;
+
+ u_size = safe_u_row_count * safe_m * sizeof(@ftyp@);
+ vt_size = safe_n * safe_vt_column_count * sizeof(@ftyp@);
mem_buff = malloc(a_size + s_size + u_size + vt_size + iwork_size);
@@ -2572,6 +2576,7 @@ init_@lapack_func@(GESDD_PARAMS_t *params,
npy_uint8 *mem_buff = NULL, *mem_buff2 = NULL;
npy_uint8 *a,*s, *u, *vt, *work, *rwork, *iwork;
size_t a_size, s_size, u_size, vt_size, work_size, rwork_size, iwork_size;
+ size_t safe_u_row_count, safe_vt_column_count;
fortran_int u_row_count, vt_column_count, work_count;
size_t safe_m = m;
size_t safe_n = n;
@@ -2581,10 +2586,13 @@ init_@lapack_func@(GESDD_PARAMS_t *params,
if (!compute_urows_vtcolumns(jobz, m, n, &u_row_count, &vt_column_count))
goto error;
+ safe_u_row_count = u_row_count;
+ safe_vt_column_count = vt_column_count;
+
a_size = safe_m * safe_n * sizeof(@ftyp@);
s_size = safe_min_m_n * sizeof(@frealtyp@);
- u_size = ((size_t)u_row_count) * safe_m * sizeof(@ftyp@);
- vt_size = safe_n * ((size_t)vt_column_count) * sizeof(@ftyp@);
+ u_size = safe_u_row_count * safe_m * sizeof(@ftyp@);
+ vt_size = safe_n * safe_vt_column_count * sizeof(@ftyp@);
rwork_size = 'N'==jobz?
(7 * safe_min_m_n) :
(5*safe_min_m_n * safe_min_m_n + 5*safe_min_m_n);