Add CX factorization (expert) routines ?GECXX (Reference-LAPACK PR 1161) (#5938)

* Add CX factorization (expert) routines ?GECXX (Reference-LAPACK PR 1161)

Author: martin-frbg

cmake/lapack.cmake

@@ -124,7 +124,7 @@ 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)
+   slatrs3.f strsyl3.f sgelst.f sgedmd.f90 sgedmdq.f90 sgecxx.f)
 
 set(SXLASRC sgesvxx.f sgerfsx.f sla_gerfsx_extended.f sla_geamv.f
    sla_gercond.f sla_gerpvgrw.f ssysvxx.f ssyrfsx.f
@@ -224,7 +224,7 @@ 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)
+   clatrs3.f ctrsyl3.f cgelst.f cgedmd.f90 cgedmdq.f90 cgecxx.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
@@ -317,7 +317,7 @@ 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)
+   dlatrs3.f dtrsyl3.f dgelst.f dgedmd.f90 dgedmdq.f90 dgecxx.f)
 
 set(DXLASRC dgesvxx.f dgerfsx.f dla_gerfsx_extended.f dla_geamv.f
    dla_gercond.f dla_gerpvgrw.f dsysvxx.f dsyrfsx.f
@@ -420,7 +420,7 @@ 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)
+   zlatrs3.f ztrsyl3.f zgelst.f zgedmd.f90 zgedmdq.f90 zgecxx.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
@@ -629,7 +629,7 @@ 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)
+   slatrs3.c strsyl3.c sgelst.c sgedmd.c sgedmdq.c sgecxx.c)
 
 set(SXLASRC sgesvxx.c sgerfsx.c sla_gerfsx_extended.c sla_geamv.c
    sla_gercond.c sla_gerpvgrw.c ssysvxx.c ssyrfsx.c
@@ -728,7 +728,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)
+   clatrs3.c ctrsyl3.c cgelst.c cgedmd.c cgedmdq.c cgecxx.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
@@ -820,7 +820,7 @@ 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)
+   dlatrs3.c dtrsyl3.c dgelst.c dgedmd.c dgedmdq.c dgecxx.c)
 
 set(DXLASRC dgesvxx.c dgerfsx.c dla_gerfsx_extended.c dla_geamv.c
    dla_gercond.c dla_gerpvgrw.c dsysvxx.c dsyrfsx.c
@@ -922,7 +922,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)
+   zgedmd.c zgedmdq.c zgecxx.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

@@ -30,6 +30,8 @@ set(CSRC
   lapacke_cgebrd_work.c
   lapacke_cgecon.c
   lapacke_cgecon_work.c
+  lapacke_cgecxx.c
+  lapacke_cgecxx_work.c
   lapacke_cgeequ.c
   lapacke_cgeequ_work.c
   lapacke_cgeequb.c
@@ -659,6 +661,8 @@ set(DSRC
   lapacke_dgebrd_work.c
   lapacke_dgecon.c
   lapacke_dgecon_work.c
+  lapacke_dgecxx.c
+  lapacke_dgecxx_work.c
   lapacke_dgeequ.c
   lapacke_dgeequ_work.c
   lapacke_dgeequb.c
@@ -1241,6 +1245,8 @@ set(SSRC
   lapacke_sgebrd_work.c
   lapacke_sgecon.c
   lapacke_sgecon_work.c
+  lapacke_sgecxx.c
+  lapacke_sgecxx_work.c
   lapacke_sgeequ.c
   lapacke_sgeequ_work.c
   lapacke_sgeequb.c
@@ -1817,6 +1823,8 @@ set(ZSRC
   lapacke_zgebrd_work.c
   lapacke_zgecon.c
   lapacke_zgecon_work.c
+  lapacke_zgecxx.c
+  lapacke_zgecxx_work.c
   lapacke_zgeequ.c
   lapacke_zgeequ_work.c
   lapacke_zgeequb.c

exports/gensymbol

@@ -917,6 +917,7 @@ lapackobjs2c="$lapackobjs2c
     clatrs3
     crscl
     ctrsyl3
+    cgecxx
     "
 #    claqz0
 #    claqz1
@@ -932,6 +933,7 @@ lapackobjs2d="$lapackobjs2d
     dlarmm
     dlatrs3
     dtrsyl3
+    dgecxx
     "
 #    dlaqz0
 #    dlaqz1
@@ -947,6 +949,7 @@ lapackobjs2s="$lapackobjs2s
     slarmm
     slatrs3
     strsyl3
+    sgecxx
     "
 
 lapackobjs2z="$lapackobjs2z
@@ -957,6 +960,7 @@ lapackobjs2z="$lapackobjs2z
     zlatrs3
     zrscl
     ztrsyl3
+    zgecxx
     "
 #    zlaqz0
 #    zlaqz1
@@ -1142,6 +1146,8 @@ lapackeobjsc="
     LAPACKE_cgedmd_work
     LAPACKE_cgedmdq
     LAPACKE_cgedmdq_work
+    LAPACKE_cgecxx
+    LAPACKE_cgecxx_work
     LAPACKE_cgeequ
     LAPACKE_cgeequ_work
     LAPACKE_cgeequb
@@ -1813,6 +1819,8 @@ lapackeobjsd="
     LAPACKE_dgedmd_work
     LAPACKE_dgedmdq
     LAPACKE_dgedmdq_work
+    LAPACKE_dgecxx
+    LAPACKE_dgecxx_work
     LAPACKE_dgeequ
     LAPACKE_dgeequ_work
     LAPACKE_dgeequb
@@ -2438,6 +2446,8 @@ lapackeobjss="
     LAPACKE_sgedmd_work
     LAPACKE_sgedmdq
     LAPACKE_sgedmdq_work
+    LAPACKE_sgecxx
+    LAPACKE_sgecxx_work
     LAPACKE_sgeequ
     LAPACKE_sgeequ_work
     LAPACKE_sgeequb
@@ -3059,6 +3069,8 @@ lapackeobjsz="
     LAPACKE_zgedmd_work
     LAPACKE_zgedmdq
     LAPACKE_zgedmdq_work
+    LAPACKE_zgecxx
+    LAPACKE_zgecxx_work
     LAPACKE_zgeequ
     LAPACKE_zgeequ_work
     LAPACKE_zgeequb

exports/gensymbol.pl

@@ -893,7 +893,8 @@
     claqp3rk,
     clatrs3,
     crscl,
-    ctrsyl3
+    ctrsyl3,
+    cgecxx
     );
 #    claqz0
 #    claqz1
@@ -908,7 +909,8 @@
     dlaqp3rk,
     dlarmm,
     dlatrs3,
-    dtrsyl3
+    dtrsyl3,
+    dgecxx
     );
 
 @lapackobjs2s = (@lapackobjs2s,
@@ -918,7 +920,8 @@
     slaqp3rk,
     slarmm,
     slatrs3,
-    strsyl3
+    strsyl3,
+    sgecxx
     );
 
 @lapackobjs2z = (@lapackobjs2z,
@@ -928,7 +931,8 @@
     zlaqp3rk,
     zlatrs3,
     zrscl,
-    ztrsyl3
+    ztrsyl3,
+    zgecxx
     );
 #    zlaqz0
 #    zlaqz1
@@ -1110,6 +1114,8 @@
     LAPACKE_cgedmd_work,
     LAPACKE_cgedmdq,
     LAPACKE_cgedmdq_work,
+    LAPACKE_cgecxx,
+    LAPACKE_cgecxx_work,
     LAPACKE_cgeequ,
     LAPACKE_cgeequ_work,
     LAPACKE_cgeequb,
@@ -1780,6 +1786,8 @@
     LAPACKE_dgedmd_work,
     LAPACKE_dgedmdq,
     LAPACKE_dgedmdq_work,
+    LAPACKE_dgecxx,
+    LAPACKE_dgecxx_work,
     LAPACKE_dgeequ,
     LAPACKE_dgeequ_work,
     LAPACKE_dgeequb,
@@ -2405,6 +2413,8 @@
     LAPACKE_sgedmd_work,
     LAPACKE_sgedmdq,
     LAPACKE_sgedmdq_work,
+    LAPACKE_sgecxx,
+    LAPACKE_sgecxx_work,
     LAPACKE_sgeequ,
     LAPACKE_sgeequ_work,
     LAPACKE_sgeequb,
@@ -3026,6 +3036,8 @@
     LAPACKE_zgedmd_work,
     LAPACKE_zgedmdq,
     LAPACKE_zgedmdq_work,
+    LAPACKE_zgecxx,
+    LAPACKE_zgecxx_work,
     LAPACKE_zgeequ,
     LAPACKE_zgeequ_work,
     LAPACKE_zgeequb,