@@ -222,6 +222,42 @@ extern "C"
222
222
lapack_int* info,
223
223
lapack_int order_len
224
224
);
225
+
226
+
227
+ void sorbdb5_ (
228
+ lapack_int* m1,
229
+ lapack_int* m2,
230
+ lapack_int* n,
231
+ float * x1,
232
+ lapack_int* incx1,
233
+ float * x2,
234
+ lapack_int* incx2,
235
+ float * q1,
236
+ lapack_int* ldq1,
237
+ float * q2,
238
+ lapack_int* ldq2,
239
+ float * work,
240
+ lapack_int* lwork,
241
+ lapack_int* info
242
+ );
243
+
244
+ void dorbdb5_ (
245
+ lapack_int* m1,
246
+ lapack_int* m2,
247
+ lapack_int* n,
248
+ double * x1,
249
+ lapack_int* incx1,
250
+ double * x2,
251
+ lapack_int* incx2,
252
+ double * q1,
253
+ lapack_int* ldq1,
254
+ double * q2,
255
+ lapack_int* ldq2,
256
+ double * work,
257
+ lapack_int* lwork,
258
+ lapack_int* info
259
+ );
260
+
225
261
}
226
262
227
263
@@ -1091,6 +1127,37 @@ inline integer_t xTRCON(
1091
1127
}
1092
1128
1093
1129
1130
+ inline integer_t xUNBDB5 (
1131
+ integer_t m1, integer_t m2, integer_t n,
1132
+ float * x1, integer_t incx1,
1133
+ float * x2, integer_t incx2,
1134
+ float * q1, integer_t ldq1, float * q2, integer_t ldq2,
1135
+ float * work, integer_t lwork)
1136
+ {
1137
+ integer_t info = -1 ;
1138
+ sorbdb5_ (
1139
+ &m1, &m2, &n, x1, &incx1, x2, &incx2, q1, &ldq1, q2, &ldq2,
1140
+ work, &lwork, &info
1141
+ );
1142
+ return info;
1143
+ }
1144
+
1145
+ inline integer_t xUNBDB5 (
1146
+ integer_t m1, integer_t m2, integer_t n,
1147
+ double * x1, integer_t incx1,
1148
+ double * x2, integer_t incx2,
1149
+ double * q1, integer_t ldq1, double * q2, integer_t ldq2,
1150
+ double * work, integer_t lwork)
1151
+ {
1152
+ integer_t info = -1 ;
1153
+ dorbdb5_ (
1154
+ &m1, &m2, &n, x1, &incx1, x2, &incx2, q1, &ldq1, q2, &ldq2,
1155
+ work, &lwork, &info
1156
+ );
1157
+ return info;
1158
+ }
1159
+
1160
+
1094
1161
1095
1162
inline integer_t xUNCSD2BY1 (
1096
1163
char jobu1, char jobu2, char jobv1t,
0 commit comments