33     double * va = emalloc(m * p * sizeof(double));  in tensor_matmul()
34     double * vb = emalloc(n * p * sizeof(double));  in tensor_matmul()
35     double * vc = emalloc(m * n * sizeof(double));  in tensor_matmul()
120     double * va = emalloc(n * n * sizeof(double));  in tensor_inverse()
121     int * pivots = emalloc(n * sizeof(int));  in tensor_inverse()
183     double * va = emalloc(m * n * sizeof(double));  in tensor_pseudoinverse()
184     double * vu = emalloc(m * m * sizeof(double));  in tensor_pseudoinverse()
185     double * vs = emalloc(k * sizeof(double));  in tensor_pseudoinverse()
186     double * vvt = emalloc(n * n * sizeof(double));  in tensor_pseudoinverse()
187     double * vb = emalloc(n * m * sizeof(double));  in tensor_pseudoinverse()
250     double * va = emalloc(m * n * sizeof(double));  in tensor_ref()
251     int * pivots = emalloc(MIN(m, n) * sizeof(int));  in tensor_ref()
318     double * va = emalloc(n * n * sizeof(double));  in tensor_cholesky()
374     double * va = emalloc(n * n * sizeof(double));  in tensor_lu()
375     int * pivots = emalloc(n * sizeof(int));  in tensor_lu()
472     double * va = emalloc(n * n * sizeof(double));  in tensor_eig()
473     double * wr = emalloc(n * sizeof(double));  in tensor_eig()
474     double * wi = emalloc(n * sizeof(double));  in tensor_eig()
475     double * vr = emalloc(n * n * sizeof(double));  in tensor_eig()
540     double * va = emalloc(n * n * sizeof(double));  in tensor_eig_symmetric()
541     double * wr = emalloc(n * sizeof(double));  in tensor_eig_symmetric()
606     double * va = emalloc(m * n * sizeof(double));  in tensor_svd()
607     double * vu = emalloc(m * m * sizeof(double));  in tensor_svd()
608     double * vs = emalloc(k * sizeof(double));  in tensor_svd()
609     double * vvt = emalloc(n * n * sizeof(double));  in tensor_svd()