ADD ?BDSVDMR3 - a MRRR-based bidiagonal SVD with DBDSDC safety fallback (Reference-LAPACK PR 1361)

cmake/lapack.cmake

@@ -134,7 +134,9 @@ set(SLASRC
    ssbev_2stage.f ssbevx_2stage.f ssbevd_2stage.f ssygv_2stage.f
    sgesvdq.f slaorhr_col_getrfnp.f
    slaorhr_col_getrfnp2.f sorgtsqr.f sorgtsqr_row.f sorhr_col.f 
-   slatrs3.f strsyl3.f sgelst.f sgedmd.f90 sgedmdq.f90 sgecxx.f)
+   slatrs3.f strsyl3.f sgelst.f sgedmd.f90 sgedmdq.f90 sgecxx.f
+   sbdsvdmr3.f sbdsvdmr3_work.f sgesvdmr3.f sbdtgk.f slar1v_tgk.f
+   slarrf_tgk.f slarrb_tgk.f slarrv_tgk.f)
 
 set(SXLASRC sgesvxx.f sgerfsx.f sla_gerfsx_extended.f sla_geamv.f
    sla_gercond.f sla_gerpvgrw.f ssysvxx.f ssyrfsx.f
@@ -234,7 +236,8 @@ set(CLASRC
    chbev_2stage.f chbevx_2stage.f chbevd_2stage.f chegv_2stage.f
    cgesvdq.f claunhr_col_getrfnp.f claunhr_col_getrfnp2.f 
    cungtsqr.f cungtsqr_row.f cunhr_col.f 
-   clatrs3.f ctrsyl3.f cgelst.f cgedmd.f90 cgedmdq.f90 cgecxx.f)
+   clatrs3.f ctrsyl3.f cgelst.f cgedmd.f90 cgedmdq.f90 cgecxx.f
+   cgesvdmr3.f)
 
 set(CXLASRC cgesvxx.f cgerfsx.f cla_gerfsx_extended.f cla_geamv.f
    cla_gercond_c.f cla_gercond_x.f cla_gerpvgrw.f
@@ -327,7 +330,9 @@ set(DLASRC
    dsbev_2stage.f dsbevx_2stage.f dsbevd_2stage.f dsygv_2stage.f
    dcombssq.f dgesvdq.f dlaorhr_col_getrfnp.f
    dlaorhr_col_getrfnp2.f dorgtsqr.f dorgtsqr_row.f dorhr_col.f 
-   dlatrs3.f dtrsyl3.f dgelst.f dgedmd.f90 dgedmdq.f90 dgecxx.f)
+   dlatrs3.f dtrsyl3.f dgelst.f dgedmd.f90 dgedmdq.f90 dgecxx.f
+   dbdsvdmr3.f dbdsvdmr3_work.f dgesvdmr3.f dbdtgk.f dlar1v_tgk.f
+   dlarrf_tgk.f dlarrb_tgk.f dlarrv_tgk.f)
 
 set(DXLASRC dgesvxx.f dgerfsx.f dla_gerfsx_extended.f dla_geamv.f
    dla_gercond.f dla_gerpvgrw.f dsysvxx.f dsyrfsx.f
@@ -430,7 +435,8 @@ set(ZLASRC
    zhbev_2stage.f zhbevx_2stage.f zhbevd_2stage.f zhegv_2stage.f
    zgesvdq.f zlaunhr_col_getrfnp.f zlaunhr_col_getrfnp2.f
    zungtsqr.f zungtsqr_row.f zunhr_col.f
-   zlatrs3.f ztrsyl3.f zgelst.f zgedmd.f90 zgedmdq.f90 zgecxx.f)
+   zlatrs3.f ztrsyl3.f zgelst.f zgedmd.f90 zgedmdq.f90 zgecxx.f
+   zgesvdmr3.f)
 
 set(ZXLASRC zgesvxx.f zgerfsx.f zla_gerfsx_extended.f zla_geamv.f
    zla_gercond_c.f zla_gercond_x.f zla_gerpvgrw.f zsysvxx.f zsyrfsx.f
@@ -639,7 +645,9 @@ set(SLASRC
    ssbev_2stage.c ssbevx_2stage.c ssbevd_2stage.c ssygv_2stage.c
    sgesvdq.c slaorhr_col_getrfnp.c
    slaorhr_col_getrfnp2.c sorgtsqr.c sorgtsqr_row.c sorhr_col.c 
-   slatrs3.c strsyl3.c sgelst.c sgedmd.c sgedmdq.c sgecxx.c)
+   slatrs3.c strsyl3.c sgelst.c sgedmd.c sgedmdq.c sgecxx.c
+   sbdsvdmr3.c sbdsvdmr3_work.c sgesvdmr3.c sbdtgk.c slar1v_tgk.c
+   slarrf_tgk.c slarrb_tgk.c slarrv_tgk.c)
 
 set(SXLASRC sgesvxx.c sgerfsx.c sla_gerfsx_extended.c sla_geamv.c
    sla_gercond.c sla_gerpvgrw.c ssysvxx.c ssyrfsx.c
@@ -738,7 +746,7 @@ set(CLASRC
    chbev_2stage.c chbevx_2stage.c chbevd_2stage.c chegv_2stage.c
    cgesvdq.c claunhr_col_getrfnp.c claunhr_col_getrfnp2.c 
    cungtsqr.c cungtsqr_row.c cunhr_col.c 
-   clatrs3.c ctrsyl3.c cgelst.c cgedmd.c cgedmdq.c cgecxx.c)
+   clatrs3.c ctrsyl3.c cgelst.c cgedmd.c cgedmdq.c cgecxx.c cgesvdmr3.c)
 
 set(CXLASRC cgesvxx.c cgerfsx.c cla_gerfsx_extended.c cla_geamv.c
    cla_gercond_c.c cla_gercond_x.c cla_gerpvgrw.c
@@ -830,7 +838,9 @@ set(DLASRC
    dsbev_2stage.c dsbevx_2stage.c dsbevd_2stage.c dsygv_2stage.c
    dcombssq.c dgesvdq.c dlaorhr_col_getrfnp.c
    dlaorhr_col_getrfnp2.c dorgtsqr.c dorgtsqr_row.c dorhr_col.c 
-   dlatrs3.c dtrsyl3.c dgelst.c dgedmd.c dgedmdq.c dgecxx.c)
+   dlatrs3.c dtrsyl3.c dgelst.c dgedmd.c dgedmdq.c dgecxx.c
+   dbdsvdmr3.c dbdsvdmr3_work.c dgesvdmr3.c dbdtgk.c dlar1v_tgk.c
+   dlarrf_tgk.c dlarrb_tgk.c dlarrv_tgk.c)
 
 set(DXLASRC dgesvxx.c dgerfsx.c dla_gerfsx_extended.c dla_geamv.c
    dla_gercond.c dla_gerpvgrw.c dsysvxx.c dsyrfsx.c
@@ -932,7 +942,7 @@ set(ZLASRC
    zhbev_2stage.c zhbevx_2stage.c zhbevd_2stage.c zhegv_2stage.c
    zgesvdq.c zlaunhr_col_getrfnp.c zlaunhr_col_getrfnp2.c
    zungtsqr.c zungtsqr_row.c zunhr_col.c zlatrs3.c ztrsyl3.c zgelst.c
-   zgedmd.c zgedmdq.c zgecxx.c)
+   zgedmd.c zgedmdq.c zgecxx.c zgesvdmr3.c)
 
 set(ZXLASRC zgesvxx.c zgerfsx.c zla_gerfsx_extended.c zla_geamv.c
    zla_gercond_c.c zla_gercond_x.c zla_gerpvgrw.c zsysvxx.c zsyrfsx.c

cmake/lapacke.cmake

@@ -100,6 +100,8 @@ set(CSRC
   lapacke_cgesv_work.c
   lapacke_cgesvd.c
   lapacke_cgesvd_work.c
+  lapacke_cgesvdmr3.c
+  lapacke_cgesvdmr3_work.c
   lapacke_cgesvdq.c
   lapacke_cgesvdq_work.c
   lapacke_cgesvdx.c
@@ -629,6 +631,8 @@ set(DSRC
   lapacke_dbbcsd_work.c
   lapacke_dbdsdc.c
   lapacke_dbdsdc_work.c
+  lapacke_dbdsvdmr3.c
+  lapacke_dbdsvdmr3_work.c
   lapacke_dbdsvdx.c
   lapacke_dbdsvdx_work.c
   lapacke_dbdsqr.c
@@ -731,6 +735,8 @@ set(DSRC
   lapacke_dgesv_work.c
   lapacke_dgesvd.c
   lapacke_dgesvd_work.c
+  lapacke_dgesvdmr3.c
+  lapacke_dgesvdmr3_work.c
   lapacke_dgesvdq.c
   lapacke_dgesvdq_work.c
   lapacke_dgesvdx.c
@@ -1213,6 +1219,8 @@ set(SSRC
   lapacke_sbbcsd_work.c
   lapacke_sbdsdc.c
   lapacke_sbdsdc_work.c
+  lapacke_sbdsvdmr3.c
+  lapacke_sbdsvdmr3_work.c
   lapacke_sbdsvdx.c
   lapacke_sbdsvdx_work.c
   lapacke_sbdsqr.c
@@ -1315,6 +1323,8 @@ set(SSRC
   lapacke_sgesv_work.c
   lapacke_sgesvd.c
   lapacke_sgesvd_work.c
+  lapacke_sgesvdmr3.c
+  lapacke_sgesvdmr3_work.c
   lapacke_sgesvdq.c
   lapacke_sgesvdq_work.c
   lapacke_sgesvdx.c
@@ -1893,6 +1903,8 @@ set(ZSRC
   lapacke_zgesv_work.c
   lapacke_zgesvd.c
   lapacke_zgesvd_work.c
+  lapacke_zgesvdmr3.c
+  lapacke_zgesvdmr3_work.c
   lapacke_zgesvdq.c
   lapacke_zgesvdq_work.c
   lapacke_zgesvdx.c

exports/gensymbol

@@ -924,6 +924,7 @@ lapackobjs2c="$lapackobjs2c
     crscl
     ctrsyl3
     cgecxx
+    cgesvdmr3
     "
 #    claqz0
 #    claqz1
@@ -940,6 +941,14 @@ lapackobjs2d="$lapackobjs2d
     dlatrs3
     dtrsyl3
     dgecxx
+    dbdsvdmr3
+    dbdsvdmr3_work
+    dbdtgk
+    dgesvdmr3
+    dlar1v_tgk
+    dlarrb_tgk
+    dlarrf_tgk
+    dlarrv_tgk
     "
 #    dlaqz0
 #    dlaqz1
@@ -956,6 +965,14 @@ lapackobjs2s="$lapackobjs2s
     slatrs3
     strsyl3
     sgecxx
+    sbdsvdmr3
+    sbdsvdmr3_work
+    sbdtgk
+    sgesvdmr3
+    slar1v_tgk
+    slarrb_tgk
+    slarrf_tgk
+    slarrv_tgk
     "
 
 lapackobjs2z="$lapackobjs2z
@@ -967,6 +984,7 @@ lapackobjs2z="$lapackobjs2z
     zrscl
     ztrsyl3
     zgecxx
+    zgesvdmr3
     "
 #    zlaqz0
 #    zlaqz1
@@ -1750,6 +1768,8 @@ lapackeobjsc="
     LAPACKE_ctz_trans
     LAPACKE_cunhr_col
     LAPACKE_cunhr_col_work
+    LAPACKE_cgesvdmr3
+    LAPACKE_cgesvdmr3_work
 "
 
 lapackeobjsd="
@@ -2377,6 +2397,10 @@ lapackeobjsd="
     LAPACKE_dtrsyl3_work
     LAPACKE_dtz_nancheck
     LAPACKE_dtz_trans
+    LAPACKE_dbdsvdmr3
+    LAPACKE_dbdsvdmr3_work
+    LAPACKE_dgesvdmr3
+    LAPACKE_dgesvdmr3_work
 "
 
 lapackeobjss="
@@ -2996,6 +3020,10 @@ lapackeobjss="
     LAPACKE_strsyl3_work
     LAPACKE_stz_nancheck
     LAPACKE_stz_trans
+    LAPACKE_sbdsvdmr3
+    LAPACKE_sbdsvdmr3_work
+    LAPACKE_sgesvdmr3
+    LAPACKE_sgesvdmr3_work
 "
 
 lapackeobjsz="
@@ -3673,6 +3701,8 @@ lapackeobjsz="
     LAPACKE_ztz_trans
     LAPACKE_zunhr_col
     LAPACKE_zunhr_col_work
+    LAPACKE_zgesvdmr3
+    LAPACKE_zgesvdmr3_work
 "
 ## @(SRCX_OBJ) from `lapack-3.4.1/lapacke/src/Makefile`
 ## Not exported: requires LAPACKE_EXTENDED to be set and depends on the

exports/gensymbol.pl

@@ -897,7 +897,8 @@
     clatrs3,
     crscl,
     ctrsyl3,
-    cgecxx
+    cgecxx,
+    cgesvdmr3
     );
 #    claqz0
 #    claqz1
@@ -913,7 +914,15 @@
     dlarmm,
     dlatrs3,
     dtrsyl3,
-    dgecxx
+    dgecxx,
+    dbdsvdmr3,
+    dbdsvdmr3_work,
+    dbdtgk,
+    dgesvdmr3,
+    dlar1v_tgk,
+    dlarrb_tgk,
+    dlarrf_tgk,
+    dlarrv_tgk
     );
 
 @lapackobjs2s = (@lapackobjs2s,
@@ -924,7 +933,15 @@
     slarmm,
     slatrs3,
     strsyl3,
-    sgecxx
+    sgecxx,
+    sbdsvdmr3,
+    sbdsvdmr3_work,
+    sbdtgk,
+    sgesvdmr3,
+    slar1v_tgk,
+    slarrb_tgk,
+    slarrf_tgk,
+    slarrv_tgk
     );
 
 @lapackobjs2z = (@lapackobjs2z,
@@ -935,7 +952,8 @@
     zlatrs3,
     zrscl,
     ztrsyl3,
-    zgecxx
+    zgecxx,
+    zgesvdmr3
     );
 #    zlaqz0
 #    zlaqz1
@@ -1714,7 +1732,9 @@
     LAPACKE_ctz_nancheck,
     LAPACKE_ctz_trans,
     LAPACKE_cunhr_col,
-    LAPACKE_cunhr_col_work
+    LAPACKE_cunhr_col_work,
+    LAPACKE_cgesvdmr3,
+    LAPACKE_cgesvdmr3_work
 );
 @lapackeobjsd = (
     LAPACKE_dgb_nancheck,
@@ -2341,7 +2361,10 @@
     LAPACKE_dtrsyl3_work,
     LAPACKE_dtz_nancheck,
     LAPACKE_dtz_trans,
-
+    LAPACKE_dbdsvdmr3,
+    LAPACKE_dbdsvdmr3_work,
+    LAPACKE_dgesvdmr3,
+    LAPACKE_dgesvdmr3_work
 );
 @lapackeobjss = (
     LAPACKE_sgb_nancheck,
@@ -2960,7 +2983,10 @@
     LAPACKE_strsyl3_work,
     LAPACKE_stz_nancheck,
     LAPACKE_stz_trans,
-
+    LAPACKE_sbdsvdmr3,
+    LAPACKE_sbdsvdmr3_work,
+    LAPACKE_sgesvdmr3,
+    LAPACKE_sgesvdmr3_work
 );
 @lapackeobjsz = (    
     LAPACKE_zgb_nancheck,
@@ -3517,6 +3543,8 @@
     LAPACKE_ztrsyl3_work,
     LAPACKE_ztz_nancheck,
     LAPACKE_ztz_trans,
+    LAPACKE_zgesvdmr3,
+    LAPACKE_zgesvdmr3_work
 
     ## @(SRCX_OBJ) from `lapack-3.4.1/lapacke/src/Makefile`
     ## Not exported: requires LAPACKE_EXTENDED to be set and depends on the