Follow up on pr1146

[langou]

Follow up on PR #1146. ( Thanks to @yizeyi18 for the PR.)

While reviewing, @angsch suggested

Shouldn't jobu = 'f' and jobu = 'r' be included here too?

I implemented the suggestion. This ended being easier for me to do a substantial rewrite.

I followed the explanation at

lapack/SRC/cgesvdq.f

Lines 159 to 198 in 034ada1

	*> \param[out] U
	*> \verbatim
	*> U is COMPLEX array, dimension
	*> LDU x M if JOBU = 'A'; see the description of LDU. In this case,
	*> on exit, U contains the M left singular vectors.
	*> LDU x N if JOBU = 'S', 'U', 'R' ; see the description of LDU. In this
	*> case, U contains the leading N or the leading NUMRANK left singular vectors.
	*> LDU x N if JOBU = 'F' ; see the description of LDU. In this case U
	*> contains N x N unitary matrix that can be used to form the left
	*> singular vectors.
	*> If JOBU = 'N', U is not referenced.
	*> \endverbatim
	*>
	*> \param[in] LDU
	*> \verbatim
	*> LDU is INTEGER.
	*> The leading dimension of the array U.
	*> If JOBU = 'A', 'S', 'U', 'R', LDU >= max(1,M).
	*> If JOBU = 'F', LDU >= max(1,N).
	*> Otherwise, LDU >= 1.
	*> \endverbatim
	*>
	*> \param[out] V
	*> \verbatim
	*> V is COMPLEX array, dimension
	*> LDV x N if JOBV = 'A', 'V', 'R' or if JOBA = 'E' .
	*> If JOBV = 'A', or 'V', V contains the N-by-N unitary matrix V**H;
	*> If JOBV = 'R', V contains the first NUMRANK rows of V**H (the right
	*> singular vectors, stored rowwise, of the NUMRANK largest singular values).
	*> If JOBV = 'N' and JOBA = 'E', V is used as a workspace.
	*> If JOBV = 'N', and JOBA.NE.'E', V is not referenced.
	*> \endverbatim
	*>
	*> \param[in] LDV
	*> \verbatim
	*> LDV is INTEGER
	*> The leading dimension of the array V.
	*> If JOBV = 'A', 'V', 'R', or JOBA = 'E', LDV >= max(1,N).
	*> Otherwise, LDV >= 1.
	*> \endverbatim

[mgates3]

Haven't completed review, but this needs fix, I think.

[mgates3]

else if, otherwise for jobu = 'a', the else clause below will set nrows_u = 1, ncols_u = 1.

Same in other precisions.

[langou]

Note: I created a new PR (PR #1221) to address the comment of @mgates3.

[mgates3]

@langou You can just keep pushing changes to the same PR -- even rebase and force push changes.