Internal working subroutine for bidiagonal decomposition *************************************************************************/
| 183 | Internal working subroutine for bidiagonal decomposition |
| 184 | *************************************************************************/ |
| 185 | static bool bidiagonalsvddecompositioninternal(ap::real_1d_array& d, |
| 186 | ap::real_1d_array e, |
| 187 | int n, |
| 188 | bool isupper, |
| 189 | bool isfractionalaccuracyrequired, |
| 190 | ap::real_2d_array& u, |
| 191 | int ustart, |
| 192 | int nru, |
| 193 | ap::real_2d_array& c, |
| 194 | int cstart, |
| 195 | int ncc, |
| 196 | ap::real_2d_array& vt, |
| 197 | int vstart, |
| 198 | int ncvt) |
| 199 | { |
| 200 | bool result; |
| 201 | int i; |
| 202 | int idir; |
| 203 | int isub; |
| 204 | int iter; |
| 205 | int j; |
| 206 | int ll; |
| 207 | int lll; |
| 208 | int m; |
| 209 | int maxit; |
| 210 | int oldll; |
| 211 | int oldm; |
| 212 | double abse; |
| 213 | double abss; |
| 214 | double cosl; |
| 215 | double cosr; |
| 216 | double cs; |
| 217 | double eps; |
| 218 | double f; |
| 219 | double g; |
| 220 | double h; |
| 221 | double mu; |
| 222 | double oldcs; |
| 223 | double oldsn; |
| 224 | double r; |
| 225 | double shift; |
| 226 | double sigmn; |
| 227 | double sigmx; |
| 228 | double sinl; |
| 229 | double sinr; |
| 230 | double sll; |
| 231 | double smax; |
| 232 | double smin; |
| 233 | double sminl; |
| 234 | double sminlo; |
| 235 | double sminoa; |
| 236 | double sn; |
| 237 | double thresh; |
| 238 | double tol; |
| 239 | double tolmul; |
| 240 | double unfl; |
| 241 | ap::real_1d_array work0; |
| 242 | ap::real_1d_array work1; |
no test coverage detected