Add a Bazel-built slicot via f2c
Previously, we were using the host binaries along with slycot, which
didn't support using it at runtime.
We built slicot itself using f2c, and then the normal Bazel C/C++
infrastructure on the resulting C code.
We're also using clapack, which is lapack and BLAS pre-run through f2c
along with various fixups.
The build is horribly noisy and full of warnings and information
messages, but oh well.
Change-Id: If13056ff56d334b7207cd9f0f7c0ce22b999b78e
diff --git a/debian/clapack.BUILD b/debian/clapack.BUILD
new file mode 100644
index 0000000..528b6a7
--- /dev/null
+++ b/debian/clapack.BUILD
@@ -0,0 +1,293 @@
+licenses(['notice'])
+
+load('@//tools/build_rules:fortran.bzl', 'f2c_copts')
+
+genrule(
+ name = 'create_sysdep1',
+ outs = ['extra_includes/sysdep1.h'],
+ srcs = ['F2CLIBS/libf2c/sysdep1.h0'],
+ cmd = 'cp $< $@',
+)
+
+_xerbla_patch = '''
+--- xerbla.c 2018-02-19 19:58:03.685420156 -0500
++++ xerbla.c 2018-02-19 19:59:02.993259128 -0500
+@@ -55,7 +55,7 @@ static integer c__1 = 1;
+ /* .. */
+ /* .. Executable Statements .. */
+
+- printf("** On entry to %6s, parameter number %2i had an illegal value\\n",
++ printf("** On entry to %6s, parameter number %2li had an illegal value\\n",
+\t\tsrname, *info);
+
+
+'''
+genrule(
+ name = 'patch_xerbla',
+ srcs = ['SRC/xerbla.c'],
+ outs = ['patched_xerbla.c'],
+ cmd = '\n'.join([
+ 'cp $< $@',
+ 'patch $@ - <<END',
+ _xerbla_patch,
+ 'END',
+ ]),
+)
+
+_err_patch = '''
+--- err.c 2018-02-19 20:06:40.532033141 -0500
++++ err.c 2018-02-19 20:10:25.907439219 -0500
+@@ -164,8 +164,10 @@ f__fatal(int n, const char *s)
+ \tif (f__curunit) {
+ \t\tfprintf(stderr,"apparent state: unit %d ",
+ \t\t\t(int)(f__curunit-f__units));
+-\t\tfprintf(stderr, f__curunit->ufnm ? "named %s\\n" : "(unnamed)\\n",
+-\t\t\tf__curunit->ufnm);
++\t\tif (f__curunit->ufnm)
++\t\t\tfprintf(stderr, "named %s\\n", f__curunit->ufnm);
++\t\telse
++\t\t\tfprintf(stderr, "(unnamed)\\n");
+ \t\t}
+ \telse
+ \t\tfprintf(stderr,"apparent state: internal I/O\\n");
+'''
+genrule(
+ name = 'patch_err',
+ srcs = ['F2CLIBS/libf2c/err.c'],
+ outs = ['patched_err.c'],
+ cmd = '\n'.join([
+ 'cp $< $@',
+ 'patch $@ - <<END',
+ _err_patch,
+ 'END',
+ ]),
+)
+
+cc_library(
+ name = 'clapack',
+ visibility = ['//visibility:public'],
+ srcs = glob(include = [
+ 'SRC/*.c',
+ 'BLAS/SRC/*.c',
+ ], exclude = [
+ # These are duplicated in SRC (with the ones in SRC a bit more cleaned up).
+ 'BLAS/SRC/xerbla.c',
+ 'BLAS/SRC/xerbla_array.c',
+
+ # We need to use a patched version of this.
+ 'SRC/xerbla.c',
+
+ # Files requiring XBLAS (extended precision), which we don't have.
+ 'SRC/sgesvxx.c',
+ 'SRC/sgerfsx.c',
+ 'SRC/sla_gerfsx_extended.c',
+ 'SRC/sla_geamv.c',
+ 'SRC/sla_gercond.c',
+ 'SRC/sla_gerpvgrw.c',
+ 'SRC/ssysvxx.c',
+ 'SRC/ssyrfsx.c',
+ 'SRC/sla_syrfsx_extended.c',
+ 'SRC/sla_syamv.c',
+ 'SRC/sla_syrcond.c',
+ 'SRC/sla_syrpvgrw.c',
+ 'SRC/sposvxx.c',
+ 'SRC/sporfsx.c',
+ 'SRC/sla_porfsx_extended.c',
+ 'SRC/sla_porcond.c',
+ 'SRC/sla_porpvgrw.c',
+ 'SRC/sgbsvxx.c',
+ 'SRC/sgbrfsx.c',
+ 'SRC/sla_gbrfsx_extended.c',
+ 'SRC/sla_gbamv.c',
+ 'SRC/sla_gbrcond.c',
+ 'SRC/sla_gbrpvgrw.c',
+ 'SRC/sla_lin_berr.c',
+ 'SRC/slarscl2.c',
+ 'SRC/slascl2.c',
+ 'SRC/sla_wwaddw.c',
+ 'SRC/cgesvxx.c',
+ 'SRC/cgerfsx.c',
+ 'SRC/cla_gerfsx_extended.c',
+ 'SRC/cla_geamv.c',
+ 'SRC/cla_gercond_c.c',
+ 'SRC/cla_gercond_x.c',
+ 'SRC/cla_gerpvgrw.c',
+ 'SRC/csysvxx.c',
+ 'SRC/csyrfsx.c',
+ 'SRC/cla_syrfsx_extended.c',
+ 'SRC/cla_syamv.c',
+ 'SRC/cla_syrcond_c.c',
+ 'SRC/cla_syrcond_x.c',
+ 'SRC/cla_syrpvgrw.c',
+ 'SRC/cposvxx.c',
+ 'SRC/cporfsx.c',
+ 'SRC/cla_porfsx_extended.c',
+ 'SRC/cla_porcond_c.c',
+ 'SRC/cla_porcond_x.c',
+ 'SRC/cla_porpvgrw.c',
+ 'SRC/cgbsvxx.c',
+ 'SRC/cgbrfsx.c',
+ 'SRC/cla_gbrfsx_extended.c',
+ 'SRC/cla_gbamv.c',
+ 'SRC/cla_gbrcond_c.c',
+ 'SRC/cla_gbrcond_x.c',
+ 'SRC/cla_gbrpvgrw.c',
+ 'SRC/chesvxx.c',
+ 'SRC/cherfsx.c',
+ 'SRC/cla_herfsx_extended.c',
+ 'SRC/cla_heamv.c',
+ 'SRC/cla_hercond_c.c',
+ 'SRC/cla_hercond_x.c',
+ 'SRC/cla_herpvgrw.c',
+ 'SRC/cla_lin_berr.c',
+ 'SRC/clarscl2.c',
+ 'SRC/clascl2.c',
+ 'SRC/cla_wwaddw.c',
+ 'SRC/dgesvxx.c',
+ 'SRC/dgerfsx.c',
+ 'SRC/dla_gerfsx_extended.c',
+ 'SRC/dla_geamv.c',
+ 'SRC/dla_gercond.c',
+ 'SRC/dla_gerpvgrw.c',
+ 'SRC/dsysvxx.c',
+ 'SRC/dsyrfsx.c',
+ 'SRC/dla_syrfsx_extended.c',
+ 'SRC/dla_syamv.c',
+ 'SRC/dla_syrcond.c',
+ 'SRC/dla_syrpvgrw.c',
+ 'SRC/dposvxx.c',
+ 'SRC/dporfsx.c',
+ 'SRC/dla_porfsx_extended.c',
+ 'SRC/dla_porcond.c',
+ 'SRC/dla_porpvgrw.c',
+ 'SRC/dgbsvxx.c',
+ 'SRC/dgbrfsx.c',
+ 'SRC/dla_gbrfsx_extended.c',
+ 'SRC/dla_gbamv.c',
+ 'SRC/dla_gbrcond.c',
+ 'SRC/dla_gbrpvgrw.c',
+ 'SRC/dla_lin_berr.c',
+ 'SRC/dlarscl2.c',
+ 'SRC/dlascl2.c',
+ 'SRC/dla_wwaddw.c',
+ 'SRC/zgesvxx.c',
+ 'SRC/zgerfsx.c',
+ 'SRC/zla_gerfsx_extended.c',
+ 'SRC/zla_geamv.c',
+ 'SRC/zla_gercond_c.c',
+ 'SRC/zla_gercond_x.c',
+ 'SRC/zla_gerpvgrw.c',
+ 'SRC/zsysvxx.c',
+ 'SRC/zsyrfsx.c',
+ 'SRC/zla_syrfsx_extended.c',
+ 'SRC/zla_syamv.c',
+ 'SRC/zla_syrcond_c.c',
+ 'SRC/zla_syrcond_x.c',
+ 'SRC/zla_syrpvgrw.c',
+ 'SRC/zposvxx.c',
+ 'SRC/zporfsx.c',
+ 'SRC/zla_porfsx_extended.c',
+ 'SRC/zla_porcond_c.c',
+ 'SRC/zla_porcond_x.c',
+ 'SRC/zla_porpvgrw.c',
+ 'SRC/zgbsvxx.c',
+ 'SRC/zgbrfsx.c',
+ 'SRC/zla_gbrfsx_extended.c',
+ 'SRC/zla_gbamv.c',
+ 'SRC/zla_gbrcond_c.c',
+ 'SRC/zla_gbrcond_x.c',
+ 'SRC/zla_gbrpvgrw.c',
+ 'SRC/zhesvxx.c',
+ 'SRC/zherfsx.c',
+ 'SRC/zla_herfsx_extended.c',
+ 'SRC/zla_heamv.c',
+ 'SRC/zla_hercond_c.c',
+ 'SRC/zla_hercond_x.c',
+ 'SRC/zla_herpvgrw.c',
+ 'SRC/zla_lin_berr.c',
+ 'SRC/zlarscl2.c',
+ 'SRC/zlascl2.c',
+ 'SRC/zla_wwaddw.c',
+ ]) + [
+ 'INSTALL/dlamch.c',
+ 'INSTALL/slamch.c',
+
+ 'patched_xerbla.c',
+ 'patched_err.c',
+
+ 'F2CLIBS/libf2c/s_cat.c',
+ 'F2CLIBS/libf2c/d_lg10.c',
+ 'F2CLIBS/libf2c/d_sign.c',
+ 'F2CLIBS/libf2c/i_dnnt.c',
+ 'F2CLIBS/libf2c/pow_di.c',
+ 'F2CLIBS/libf2c/s_copy.c',
+ 'F2CLIBS/libf2c/s_cmp.c',
+ 'F2CLIBS/libf2c/i_nint.c',
+ 'F2CLIBS/libf2c/f77_aloc.c',
+ 'F2CLIBS/libf2c/exit_.c',
+ 'F2CLIBS/libf2c/r_cnjg.c',
+ 'F2CLIBS/libf2c/c_abs.c',
+ 'F2CLIBS/libf2c/r_imag.c',
+ 'F2CLIBS/libf2c/c_div.c',
+ 'F2CLIBS/libf2c/c_exp.c',
+ 'F2CLIBS/libf2c/d_imag.c',
+ 'F2CLIBS/libf2c/r_sign.c',
+ 'F2CLIBS/libf2c/d_cnjg.c',
+ 'F2CLIBS/libf2c/z_abs.c',
+ 'F2CLIBS/libf2c/z_div.c',
+ 'F2CLIBS/libf2c/z_exp.c',
+ 'F2CLIBS/libf2c/z_sqrt.c',
+ 'F2CLIBS/libf2c/pow_dd.c',
+ 'F2CLIBS/libf2c/pow_ri.c',
+ 'F2CLIBS/libf2c/pow_ci.c',
+ 'F2CLIBS/libf2c/pow_ii.c',
+ 'F2CLIBS/libf2c/pow_zi.c',
+ 'F2CLIBS/libf2c/c_sqrt.c',
+ 'F2CLIBS/libf2c/r_lg10.c',
+ 'F2CLIBS/libf2c/i_len.c',
+ 'F2CLIBS/libf2c/cabs.c',
+ 'F2CLIBS/libf2c/sig_die.c',
+ 'F2CLIBS/libf2c/close.c',
+ 'F2CLIBS/libf2c/open.c',
+ 'F2CLIBS/libf2c/endfile.c',
+ 'F2CLIBS/libf2c/util.c',
+ 'F2CLIBS/libf2c/iio.c',
+ 'F2CLIBS/libf2c/fmt.c',
+ 'F2CLIBS/libf2c/rdfmt.c',
+ 'F2CLIBS/libf2c/wrtfmt.c',
+ 'F2CLIBS/libf2c/ctype.c',
+ 'F2CLIBS/libf2c/wref.c',
+ 'F2CLIBS/libf2c/fmtlib.c',
+ 'F2CLIBS/libf2c/lread.c',
+ 'F2CLIBS/libf2c/rsfe.c',
+ 'F2CLIBS/libf2c/sfe.c',
+ 'F2CLIBS/libf2c/dolio.c',
+ 'F2CLIBS/libf2c/wsfe.c',
+
+ 'extra_includes/sysdep1.h',
+ ],
+ hdrs = glob([
+ 'INCLUDE/*.h',
+ 'F2CLIBS/libf2c/*.h',
+ ]),
+ includes = [
+ 'INCLUDE',
+ 'F2CLIBS/libf2c',
+ 'extra_includes',
+ ],
+ copts = f2c_copts + [
+ '-Wno-sign-compare',
+ '-Wno-cast-qual',
+ '-Wno-cast-align',
+ '-Wno-self-assign',
+
+ # Some files don't #include system headers when they should. sysdep1.h
+ # messes with feature test macros, so it always has to come first.
+ '-include', 'sysdep1.h',
+ '-include', 'stdio.h',
+
+ # Don't mangle the names of all the BLAS symbols, because slicot needs to
+ # call them directly.
+ '-DNO_BLAS_WRAP',
+ ],
+)