2004-07-15 Toon Moene <toon@moene.indiv.nluug.nl>

	* libf2c: Removed.
	* gcc/gccbug.in: Updated because of libf2c removal.
	* maintainer-scripts/gcc_release: Ditto.


git-svn-id: svn+ssh://gcc.gnu.org/svn/gcc/trunk@84759 138bc75d-0d04-0410-961f-82ee72b054a4

132 files changed, 0 insertions, 3076 deletions

diff --git a/libf2c/libF77/F77_aloc.c b/libf2c/libF77/F77_aloc.c
deleted file mode 100644
index b2861029198..00000000000
--- a/libf2c/libF77/F77_aloc.c
+++ /dev/null
@@ -1,24 +0,0 @@
-#include "f2c.h"
-#undef abs
-#undef min
-#undef max
-#include <stdio.h>
-
-static integer memfailure = 3;
-
-#include <stdlib.h>
-extern void G77_exit_0 (integer *);
-
-char *
-F77_aloc (integer Len, char *whence)
-{
-  char *rv;
-  unsigned int uLen = (unsigned int) Len;	/* for K&R C */
-
-  if (!(rv = (char *) malloc (uLen)))
-    {
-      fprintf (stderr, "malloc(%u) failure in %s\n", uLen, whence);
-      G77_exit_0 (&memfailure);
-    }
-  return rv;
-}
diff --git a/libf2c/libF77/Makefile.in b/libf2c/libF77/Makefile.in
deleted file mode 100644
index 07af6eb2773..00000000000
--- a/libf2c/libF77/Makefile.in
+++ /dev/null
@@ -1,255 +0,0 @@
-# Makefile for GNU F77 compiler runtime.
-# Copyright 1990 - 1994 by AT&T Bell Laboratories and Bellcore (see the
-# file `Notice').
-# Portions of this file Copyright (C) 1995-1998, 2001 Free Software Foundation, Inc.
-#   Contributed by Dave Love (d.love@dl.ac.uk).
-#
-#This file is part of GNU Fortran.
-#
-#GNU Fortran is free software; you can redistribute it and/or modify
-#it under the terms of the GNU General Public License as published by
-#the Free Software Foundation; either version 2, or (at your option)
-#any later version.
-#
-#GNU Fortran is distributed in the hope that it will be useful,
-#but WITHOUT ANY WARRANTY; without even the implied warranty of
-#MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
-#GNU General Public License for more details.
-#
-#You should have received a copy of the GNU General Public License
-#along with GNU Fortran; see the file COPYING.  If not, write to
-#the Free Software Foundation, 59 Temple Place - Suite 330, Boston, MA
-#02111-1307, USA.
-
-#### Start of system configuration section. ####
-
-# $(srcdir) must be set to the g77 runtime libF77 source directory.
-srcdir = @srcdir@
-VPATH = @srcdir@
-
-# configure sets this to all the -D options appropriate for the
-# configuration.
-DEFS = @DEFS@
-
-F2C_H_DIR = @srcdir@/..
-G2C_H_DIR = ..
-CC = @CC@
-CFLAGS = @CFLAGS@
-CPPFLAGS = @CPPFLAGS@
-AR = @AR@
-ARFLAGS = rc
-RANLIB = @RANLIB@
-@SET_MAKE@
-
-SHELL = @SHELL@
-
-#### End of system configuration section. ####
-
-ALL_CFLAGS = -I. -I$(srcdir) -I$(G2C_H_DIR) -I$(F2C_H_DIR) $(CPPFLAGS) $(DEFS) $(WARN_CFLAGS) $(CFLAGS)
-
-.SUFFIXES:
-.SUFFIXES: .c .lo
-
-.c.lo:
-	@LIBTOOL@ --mode=compile $(CC) -c -DSkip_f2c_Undefs $(ALL_CFLAGS) $<
-
-MISC =	F77_aloc.lo VersionF.lo s_rnge.lo abort_.lo getarg_.lo iargc_.lo\
-	getenv_.lo signal_.lo s_stop.lo s_paus.lo system_.lo cabs.lo\
-	derf_.lo derfc_.lo erf_.lo erfc_.lo sig_die.lo exit_.lo setarg.lo setsig.lo
-POW =	pow_ci.lo pow_dd.lo pow_di.lo pow_hh.lo pow_ii.lo  pow_ri.lo pow_zi.lo pow_zz.lo \
-	pow_qq.lo
-CX =	c_abs.lo c_cos.lo c_div.lo c_exp.lo c_log.lo c_sin.lo c_sqrt.lo
-DCX =	z_abs.lo z_cos.lo z_div.lo z_exp.lo z_log.lo z_sin.lo z_sqrt.lo
-REAL =	r_abs.lo r_acos.lo r_asin.lo r_atan.lo r_atn2.lo r_cnjg.lo r_cos.lo\
-	r_cosh.lo r_dim.lo r_exp.lo r_imag.lo r_int.lo\
-	r_lg10.lo r_log.lo r_mod.lo r_nint.lo r_sign.lo\
-	r_sin.lo r_sinh.lo r_sqrt.lo r_tan.lo r_tanh.lo
-DBL =	d_abs.lo d_acos.lo d_asin.lo d_atan.lo d_atn2.lo\
-	d_cnjg.lo d_cos.lo d_cosh.lo d_dim.lo d_exp.lo\
-	d_imag.lo d_int.lo d_lg10.lo d_log.lo d_mod.lo\
-	d_nint.lo d_prod.lo d_sign.lo d_sin.lo d_sinh.lo\
-	d_sqrt.lo d_tan.lo d_tanh.lo
-INT =	i_abs.lo i_dim.lo i_dnnt.lo i_indx.lo i_len.lo i_mod.lo i_nint.lo i_sign.lo
-HALF =	h_abs.lo h_dim.lo h_dnnt.lo h_indx.lo h_len.lo h_mod.lo  h_nint.lo h_sign.lo
-CMP =	l_ge.lo l_gt.lo l_le.lo l_lt.lo hl_ge.lo hl_gt.lo hl_le.lo hl_lt.lo
-EFL =	ef1asc_.lo ef1cmc_.lo
-CHAR =	s_cat.lo s_cmp.lo s_copy.lo
-F90BIT = lbitbits.lo lbitshft.lo qbitbits.lo qbitshft.lo
-
-OBJS = $(MISC) $(POW) $(CX) $(DCX) $(REAL) $(DBL) $(INT) \
-	 $(HALF) $(CMP) $(EFL) $(CHAR) $(F90BIT)
-
-all: ../s-libf77 ../libfrtbegin.a
-
-../s-libf77: $(OBJS)
-	-rm -f $@.T $@
-	objs='$(OBJS)'; for name in $$objs; do \
-	  echo libF77/$${name} >> $@.T; done
-	mv $@.T $@
-
-Makefile: Makefile.in config.status
-	$(SHELL) config.status
-
-config.status: configure
-	rm -f config.cache
-	CONFIG_SITE=no-such-file CC='$(CC)' AR='$(AR)' CFLAGS='$(CFLAGS)' \
-	  CPPFLAGS='$(CPPFLAGS)' $(SHELL) config.status --recheck
-
-${srcdir}/configure: configure.in
-	rm -f config.cache
-	cd ${srcdir} && autoconf
-
-VersionF.lo: Version.c
-	@LIBTOOL@ --mode=compile $(CC) -c $(ALL_CFLAGS) $(srcdir)/Version.c -o $@
-
-frtbegin.o : main.c
-	$(CC) -c $(ALL_CFLAGS) $(srcdir)/main.c -o $@
-
-../libfrtbegin.a: frtbegin.o
-	-rm -f $@
-	$(AR) $(ARFLAGS) $@ frtbegin.o
-	$(RANLIB) $@
-
-F77_aloc.lo: F77_aloc.c
-s_rnge.lo: s_rnge.c
-abort_.lo: abort_.c
-getarg_.lo: getarg_.c
-iargc_.lo: iargc_.c
-getenv_.lo: getenv_.c
-signal_.lo: signal_.c
-s_stop.lo: s_stop.c
-s_paus.lo: s_paus.c
-system_.lo: system_.c
-cabs.lo: cabs.c
-derf_.lo: derf_.c
-derfc_.lo: derfc_.c
-erf_.lo: erf_.c
-erfc_.lo: erfc_.c
-sig_die.lo: sig_die.c
-exit_.lo: exit_.c
-setarg.lo: setarg.c
-setsig.lo: setsig.c
-pow_ci.lo: pow_ci.c
-pow_dd.lo: pow_dd.c
-pow_di.lo: pow_di.c
-pow_hh.lo: pow_hh.c
-pow_ii.lo: pow_ii.c
-pow_ri.lo:  pow_ri.c
-pow_zi.lo: pow_zi.c
-pow_zz.lo: pow_zz.c
-pow_qq.lo:  pow_qq.c
-c_abs.lo: c_abs.c
-c_cos.lo: c_cos.c
-c_div.lo: c_div.c
-c_exp.lo: c_exp.c
-c_log.lo: c_log.c
-c_sin.lo: c_sin.c
-c_sqrt.lo: c_sqrt.c
-z_abs.lo: z_abs.c
-z_cos.lo: z_cos.c
-z_div.lo: z_div.c
-z_exp.lo: z_exp.c
-z_log.lo: z_log.c
-	@LIBTOOL@ --mode=compile $(CC) -c -DSkip_f2c_Undefs $(ALL_CFLAGS) -ffloat-store $(srcdir)/z_log.c
-z_sin.lo: z_sin.c
-z_sqrt.lo: z_sqrt.c
-r_abs.lo: r_abs.c
-r_acos.lo: r_acos.c
-r_asin.lo: r_asin.c
-r_atan.lo: r_atan.c
-r_atn2.lo: r_atn2.c
-r_cnjg.lo: r_cnjg.c
-r_cos.lo: r_cos.c
-r_cosh.lo: r_cosh.c
-r_dim.lo: r_dim.c
-r_exp.lo: r_exp.c
-r_imag.lo: r_imag.c
-r_int.lo: r_int.c
-r_lg10.lo: r_lg10.c
-r_log.lo: r_log.c
-r_mod.lo: r_mod.c
-r_nint.lo: r_nint.c
-r_sign.lo: r_sign.c
-r_sin.lo: r_sin.c
-r_sinh.lo: r_sinh.c
-r_sqrt.lo: r_sqrt.c
-r_tan.lo: r_tan.c
-r_tanh.lo: r_tanh.c
-d_abs.lo: d_abs.c
-d_acos.lo: d_acos.c
-d_asin.lo: d_asin.c
-d_atan.lo: d_atan.c
-d_atn2.lo: d_atn2.c
-d_cnjg.lo: d_cnjg.c
-d_cos.lo: d_cos.c
-d_cosh.lo: d_cosh.c
-d_dim.lo: d_dim.c
-d_exp.lo: d_exp.c
-d_imag.lo: d_imag.c
-d_int.lo: d_int.c
-d_lg10.lo: d_lg10.c
-d_log.lo: d_log.c
-d_mod.lo: d_mod.c
-d_nint.lo: d_nint.c
-d_prod.lo: d_prod.c
-d_sign.lo: d_sign.c
-d_sin.lo: d_sin.c
-d_sinh.lo: d_sinh.c
-d_sqrt.lo: d_sqrt.c
-d_tan.lo: d_tan.c
-d_tanh.lo: d_tanh.c
-i_abs.lo: i_abs.c
-i_dim.lo: i_dim.c
-i_dnnt.lo: i_dnnt.c
-i_indx.lo: i_indx.c
-i_len.lo: i_len.c
-i_mod.lo: i_mod.c
-i_nint.lo: i_nint.c
-i_sign.lo: i_sign.c
-h_abs.lo:  h_abs.c
-h_dim.lo: h_dim.c
-h_dnnt.lo: h_dnnt.c
-h_indx.lo: h_indx.c
-h_len.lo: h_len.c
-h_mod.lo: h_mod.c
-h_nint.lo:  h_nint.c
-h_sign.lo: h_sign.c
-l_ge.lo: l_ge.c
-l_gt.lo: l_gt.c
-l_le.lo: l_le.c
-l_lt.lo: l_lt.c
-hl_ge.lo: hl_ge.c
-hl_gt.lo: hl_gt.c
-hl_le.lo: hl_le.c
-hl_lt.lo: hl_lt.c
-ef1asc_.lo: ef1asc_.c
-ef1cmc_.lo: ef1cmc_.c
-s_cat.lo: s_cat.c
-s_cmp.lo: s_cmp.c
-s_copy.lo: s_copy.c
-lbitbits.lo: lbitbits.c
-lbitshft.lo: lbitshft.c
-qbitbits.lo: qbitbits.c
-qbitshft.lo: qbitshft.c
-
-# Not quite all these actually do depend on f2c.h...
-$(MISC) $(POW) $(CX) $(DCX) $(REAL) $(DBL) $(INT) $(HALF) $(CMP) $(EFL) \
-  $(CHAR) $(F90BIT): $(F2C_H_DIR)/f2c.h $(G2C_H_DIR)/g2c.h
-
-check install uninstall install-strip dist installcheck installdirs:
-
-mostlyclean:
-	rm -f *.o *.lo ../libfrtbegin.a
-	rm -rf .libs
-
-clean: mostlyclean
-	rm -f config.log
-	rm -f ../s-libf77
-
-distclean: clean
-	rm -f config.cache config.status Makefile ../s-libf77 configure
-
-maintainer-clean:
-
-.PHONY: mostlyclean clean distclean maintainer-clean all check uninstall \
-        install-strip dist installcheck installdirs archive
diff --git a/libf2c/libF77/Notice b/libf2c/libF77/Notice
deleted file mode 100644
index 261b719bc57..00000000000
--- a/libf2c/libF77/Notice
+++ /dev/null
@@ -1,23 +0,0 @@
-/****************************************************************
-Copyright 1990 - 1997 by AT&T, Lucent Technologies and Bellcore.
-
-Permission to use, copy, modify, and distribute this software
-and its documentation for any purpose and without fee is hereby
-granted, provided that the above copyright notice appear in all
-copies and that both that the copyright notice and this
-permission notice and warranty disclaimer appear in supporting
-documentation, and that the names of AT&T, Bell Laboratories,
-Lucent or Bellcore or any of their entities not be used in
-advertising or publicity pertaining to distribution of the
-software without specific, written prior permission.
-
-AT&T, Lucent and Bellcore disclaim all warranties with regard to
-this software, including all implied warranties of
-merchantability and fitness.  In no event shall AT&T, Lucent or
-Bellcore be liable for any special, indirect or consequential
-damages or any damages whatsoever resulting from loss of use,
-data or profits, whether in an action of contract, negligence or
-other tortious action, arising out of or in connection with the
-use or performance of this software.
-****************************************************************/
-
diff --git a/libf2c/libF77/README.netlib b/libf2c/libF77/README.netlib
deleted file mode 100644
index 5e532ee0990..00000000000
--- a/libf2c/libF77/README.netlib
+++ /dev/null
@@ -1,112 +0,0 @@
-If your compiler does not recognize ANSI C headers,
-compile with KR_headers defined:  either add -DKR_headers
-to the definition of CFLAGS in the makefile, or insert
-
-#define KR_headers
-
-at the top of f2c.h , cabs.c , main.c , and sig_die.c .
-
-Under MS-DOS, compile s_paus.c with -DMSDOS.
-
-If you have a really ancient K&R C compiler that does not understand
-void, add -Dvoid=int to the definition of CFLAGS in the makefile.
-
-If you use a C++ compiler, first create a local f2c.h by appending
-f2ch.add to the usual f2c.h, e.g., by issuing the command
-	make f2c.h
-which assumes f2c.h is installed in /usr/include .
-
-If your system lacks onexit() and you are not using an ANSI C
-compiler, then you should compile main.c, s_paus.c, s_stop.c, and
-sig_die.c with NO_ONEXIT defined.  See the comments about onexit in
-the makefile.
-
-If your system has a double drem() function such that drem(a,b)
-is the IEEE remainder function (with double a, b), then you may
-wish to compile r_mod.c and d_mod.c with IEEE_drem defined.
-On some systems, you may also need to compile with -Ddrem=remainder .
-
-To check for transmission errors, issue the command
-	make check
-This assumes you have the xsum program whose source, xsum.c,
-is distributed as part of "all from f2c/src".  If you do not
-have xsum, you can obtain xsum.c by sending the following E-mail
-message to netlib@netlib.bell-labs.com
-	send xsum.c from f2c/src
-
-The makefile assumes you have installed f2c.h in a standard
-place (and does not cause recompilation when f2c.h is changed);
-f2c.h comes with "all from f2c" (the source for f2c) and is
-available separately ("f2c.h from f2c").
-
-Most of the routines in libF77 are support routines for Fortran
-intrinsic functions or for operations that f2c chooses not
-to do "in line".  There are a few exceptions, summarized below --
-functions and subroutines that appear to your program as ordinary
-external Fortran routines.
-
-1.	CALL ABORT prints a message and causes a core dump.
-
-2.	ERF(r) and DERF(d) and the REAL and DOUBLE PRECISION
-	error functions (with x REAL and d DOUBLE PRECISION);
-	DERF must be declared DOUBLE PRECISION in your program.
-	Both ERF and DERF assume your C library provides the
-	underlying erf() function (which not all systems do).
-
-3.	ERFC(r) and DERFC(d) are the complementary error functions:
-	ERFC(r) = 1 - ERF(r) and DERFC(d) = 1.d0 - DERFC(d)
-	(except that their results may be more accurate than
-	explicitly evaluating the above formulae would give).
-	Again, ERFC and r are REAL, and DERFC and d are DOUBLE
-	PRECISION (and must be declared as such in your program),
-	and ERFC and DERFC rely on your system's erfc().
-
-4.	CALL GETARG(n,s), where n is an INTEGER and s is a CHARACTER
-	variable, sets s to the n-th command-line argument (or to
-	all blanks if there are fewer than n command-line arguments);
-	CALL GETARG(0,s) sets s to the name of the program (on systems
-	that support this feature).  See IARGC below.
-
-5.	CALL GETENV(name, value), where name and value are of type
-	CHARACTER, sets value to the environment value, $name, of
-	name (or to blanks if $name has not been set).
-
-6.	NARGS = IARGC() sets NARGS to the number of command-line
-	arguments (an INTEGER value).
-
-7.	CALL SIGNAL(n,func), where n is an INTEGER and func is an
-	EXTERNAL procedure, arranges for func to be invoked when
-	signal n occurs (on systems where this makes sense).
-
-8.	CALL SYSTEM(cmd), where cmd is of type CHARACTER, passes
-	cmd to the system's command processor (on systems where
-	this can be done).
-
-The makefile does not attempt to compile pow_qq.c, qbitbits.c,
-and qbitshft.c, which are meant for use with INTEGER*8.  To use
-INTEGER*8, you must modify f2c.h to declare longint and ulongint
-appropriately; then add pow_qq.o to the POW = line in the makefile,
-and add " qbitbits.o qbitshft.o" to the makefile's F90BIT = line.
-
-Following Fortran 90, s_cat.c and s_copy.c allow the target of a
-(character string) assignment to be appear on its right-hand, at
-the cost of some extra overhead for all run-time concatenations.
-If you prefer the  extra efficiency that comes with the Fortran 77
-requirement that the left-hand side of a character assignment not
-be involved in the right-hand side, compile s_cat.c and s_copy.c
-with -DNO_OVERWRITE .
-
-If your system lacks a ranlib command, you don't need it.
-Either comment out the makefile's ranlib invocation, or install
-a harmless "ranlib" command somewhere in your PATH, such as the
-one-line shell script
-
-	exit 0
-
-or (on some systems)
-
-	exec /usr/bin/ar lts $1 >/dev/null
-
-If your compiler complains about the signal calls in main.c, s_paus.c,
-and signal_.c, you may need to adjust signal1.h suitably.  See the
-comments in signal1.h.
diff --git a/libf2c/libF77/Version.c b/libf2c/libF77/Version.c
deleted file mode 100644
index a27b0df7696..00000000000
--- a/libf2c/libF77/Version.c
+++ /dev/null
@@ -1,94 +0,0 @@
-const char __LIBF77_VERSION__[] = "@(#) LIBF77 VERSION 20000929\n";
-extern const char __LIBI77_VERSION__[];
-extern const char __LIBU77_VERSION__[];
-
-/*
-2.00	11 June 1980.  File version.c added to library.
-2.01	31 May 1988.  s_paus() flushes stderr; names of hl_* fixed
-	[ d]erf[c ] added
-	 8 Aug. 1989: #ifdefs for f2c -i2 added to s_cat.c
-	29 Nov. 1989: s_cmp returns long (for f2c)
-	30 Nov. 1989: arg types from f2c.h
-	12 Dec. 1989: s_rnge allows long names
-	19 Dec. 1989: getenv_ allows unsorted environment
-	28 Mar. 1990: add exit(0) to end of main()
-	 2 Oct. 1990: test signal(...) == SIG_IGN rather than & 01 in main
-	17 Oct. 1990: abort() calls changed to sig_die(...,1)
-	22 Oct. 1990: separate sig_die from main
-	25 Apr. 1991: minor, theoretically invisible tweaks to s_cat, sig_die
-	31 May  1991: make system_ return status
-	18 Dec. 1991: change long to ftnlen (for -i2) many places
-	28 Feb. 1992: repair z_sqrt.c (scribbled on input, gave wrong answer)
-	18 July 1992: for n < 0, repair handling of 0**n in pow_[dr]i.c
-			and m**n in pow_hh.c and pow_ii.c;
-			catch SIGTRAP in main() for error msg before abort
-	23 July 1992: switch to ANSI prototypes unless KR_headers is #defined
-	23 Oct. 1992: fix botch in signal_.c (erroneous deref of 2nd arg);
-			change Cabs to f__cabs.
-	12 March 1993: various tweaks for C++
-	 2 June 1994: adjust so abnormal terminations invoke f_exit just once
-	16 Sept. 1994: s_cmp: treat characters as unsigned in comparisons.
-	19 Sept. 1994: s_paus: flush after end of PAUSE; add -DMSDOS
-	12 Jan. 1995:	pow_[dhiqrz][hiq]: adjust x**i to work on machines
-			that sign-extend right shifts when i is the most
-			negative integer.
-	26 Jan. 1995: adjust s_cat.c, s_copy.c to permit the left-hand side
-			of character assignments to appear on the right-hand
-			side (unless compiled with -DNO_OVERWRITE).
-	27 Jan. 1995: minor tweak to s_copy.c: copy forward whenever
-			possible (for better cache behavior).
-	30 May 1995:  added subroutine exit(rc) integer rc. Version not changed.
-	29 Aug. 1995: add F77_aloc.c; use it in s_cat.c and system_.c.
-	6 Sept. 1995: fix return type of system_ under -DKR_headers.
-	19 Dec. 1995: s_cat.c: fix bug when 2nd or later arg overlaps lhs.
-	19 Mar. 1996: s_cat.c: supply missing break after overlap detection.
-	13 May 1996:  add [lq]bitbits.c and [lq]bitshft.c (f90 bit intrinsics).
-	19 June 1996: add casts to unsigned in [lq]bitshft.c.
-	26 Feb. 1997: adjust functions with a complex output argument
-			to permit aliasing it with input arguments.
-			(For now, at least, this is just for possible
-			benefit of g77.)
-	4 April 1997: [cz]_div.c: tweaks invisible on most systems (that may
-			affect systems using gratuitous extra precision).
-	19 Sept. 1997: [de]time_.c (Unix systems only): change return
-			type to double.
-	2 May 1999:	getenv_.c: omit environ in favor of getenv().
-			c_cos.c, c_exp.c, c_sin.c, d_cnjg.c, r_cnjg.c,
-			z_cos.c, z_exp.c, z_log.c, z_sin.c: cope fully with
-			overlapping arguments caused by equivalence.
-	3 May 1999:	"invisible" tweaks to omit compiler warnings in
-			abort_.c, ef1asc_.c, s_rnge.c, s_stop.c.
-	7 Sept. 1999: [cz]_div.c: arrange for compilation under
-			-DIEEE_COMPLEX_DIVIDE to make these routines
-			avoid calling sig_die when the denominator
-			vanishes; instead, they return pairs of NaNs
-			or Infinities, depending whether the numerator
-			also vanishes or not.  VERSION not changed.
-	15 Nov. 1999: s_rnge.c: add casts for the case of
-			sizeof(ftnint) == sizeof(int) < sizeof(long).
-	10 March 2000: z_log.c: improve accuracy of Real(log(z)) for, e.g.,
-			z near (+-1,eps) with |eps| small.  For the old
-			evaluation, compile with -DPre20000310 .
-	20 April 2000: s_cat.c: tweak argument types to accord with
-			calls by f2c when ftnint and ftnlen are of
-			different sizes (different numbers of bits).
-	4 July 2000: adjustments to permit compilation by C++ compilers;
-			VERSION string remains unchanged. NOT APPLIED FOR G77.
-	29 Sept. 2000: dtime_.c, etime_.c: use floating-point divide.
-			dtime_.d, erf_.c, erfc_.c, etime.c: for use with
-			"f2c -R", compile with -DREAL=float.
-*/
-
-#include <stdio.h>
-
-void
-g77__fvers__ ()
-{
-  fputs ("GNU Fortran library.\n", stderr);
-#if defined __GNUC__ && defined __VERSION__
-  fprintf (stderr, "Compiled by GCC %s\n", __VERSION__);
-#endif
-  fputs (__LIBF77_VERSION__, stderr);
-  fputs (__LIBI77_VERSION__, stderr);
-  fputs (__LIBU77_VERSION__, stderr);
-}
diff --git a/libf2c/libF77/abort_.c b/libf2c/libF77/abort_.c
deleted file mode 100644
index 761bc3bd6e9..00000000000
--- a/libf2c/libF77/abort_.c
+++ /dev/null
@@ -1,11 +0,0 @@
-#include <stdio.h>
-#include "f2c.h"
-
-extern void sig_die (char *, int);
-
-int
-G77_abort_0 (void)
-{
-  sig_die ("Fortran abort routine called", 1);
-  return 0;			/* not reached */
-}
diff --git a/libf2c/libF77/c_abs.c b/libf2c/libF77/c_abs.c
deleted file mode 100644
index 3fc4d7c0a6a..00000000000
--- a/libf2c/libF77/c_abs.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-extern double f__cabs (double, double);
-
-double
-c_abs (complex * z)
-{
-  return (f__cabs (z->r, z->i));
-}
diff --git a/libf2c/libF77/c_cos.c b/libf2c/libF77/c_cos.c
deleted file mode 100644
index 59659754978..00000000000
--- a/libf2c/libF77/c_cos.c
+++ /dev/null
@@ -1,12 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include "math.h"
-
-void
-c_cos (complex * r, complex * z)
-{
-  double zi = z->i, zr = z->r;
-  r->r = cos (zr) * cosh (zi);
-  r->i = -sin (zr) * sinh (zi);
-}
diff --git a/libf2c/libF77/c_div.c b/libf2c/libF77/c_div.c
deleted file mode 100644
index b5ede0eb801..00000000000
--- a/libf2c/libF77/c_div.c
+++ /dev/null
@@ -1,43 +0,0 @@
-#include "f2c.h"
-
-extern void sig_die (char *, int);
-void
-c_div (complex * c, complex * a, complex * b)
-{
-  double ratio, den;
-  double abr, abi, cr;
-
-  if ((abr = b->r) < 0.)
-    abr = -abr;
-  if ((abi = b->i) < 0.)
-    abi = -abi;
-  if (abr <= abi)
-    {
-      if (abi == 0)
-	{
-#ifdef IEEE_COMPLEX_DIVIDE
-	  float af, bf;
-	  af = bf = abr;
-	  if (a->i != 0 || a->r != 0)
-	    af = 1.;
-	  c->i = c->r = af / bf;
-	  return;
-#else
-	  sig_die ("complex division by zero", 1);
-#endif
-	}
-      ratio = (double) b->r / b->i;
-      den = b->i * (1 + ratio * ratio);
-      cr = (a->r * ratio + a->i) / den;
-      c->i = (a->i * ratio - a->r) / den;
-    }
-
-  else
-    {
-      ratio = (double) b->i / b->r;
-      den = b->r * (1 + ratio * ratio);
-      cr = (a->r + a->i * ratio) / den;
-      c->i = (a->i - a->r * ratio) / den;
-    }
-  c->r = cr;
-}
diff --git a/libf2c/libF77/c_exp.c b/libf2c/libF77/c_exp.c
deleted file mode 100644
index 56a8695420c..00000000000
--- a/libf2c/libF77/c_exp.c
+++ /dev/null
@@ -1,14 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include "math.h"
-
-void
-c_exp (complex * r, complex * z)
-{
-  double expx, zi = z->i;
-
-  expx = exp (z->r);
-  r->r = expx * cos (zi);
-  r->i = expx * sin (zi);
-}
diff --git a/libf2c/libF77/c_log.c b/libf2c/libF77/c_log.c
deleted file mode 100644
index 7d5b9513167..00000000000
--- a/libf2c/libF77/c_log.c
+++ /dev/null
@@ -1,13 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include "math.h"
-extern double f__cabs (double, double);
-
-void
-c_log (complex * r, complex * z)
-{
-  double zi, zr;
-  r->i = atan2 (zi = z->i, zr = z->r);
-  r->r = log (f__cabs (zr, zi));
-}
diff --git a/libf2c/libF77/c_sin.c b/libf2c/libF77/c_sin.c
deleted file mode 100644
index 44bce9f7739..00000000000
--- a/libf2c/libF77/c_sin.c
+++ /dev/null
@@ -1,12 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include "math.h"
-
-void
-c_sin (complex * r, complex * z)
-{
-  double zi = z->i, zr = z->r;
-  r->r = sin (zr) * cosh (zi);
-  r->i = cos (zr) * sinh (zi);
-}
diff --git a/libf2c/libF77/c_sqrt.c b/libf2c/libF77/c_sqrt.c
deleted file mode 100644
index 81b72fc372e..00000000000
--- a/libf2c/libF77/c_sqrt.c
+++ /dev/null
@@ -1,30 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include "math.h"
-extern double f__cabs (double, double);
-
-void
-c_sqrt (complex * r, complex * z)
-{
-  double mag, t;
-  double zi = z->i, zr = z->r;
-
-  if ((mag = f__cabs (zr, zi)) == 0.)
-    r->r = r->i = 0.;
-  else if (zr > 0)
-    {
-      r->r = t = sqrt (0.5 * (mag + zr));
-      t = zi / t;
-      r->i = 0.5 * t;
-    }
-  else
-    {
-      t = sqrt (0.5 * (mag - zr));
-      if (zi < 0)
-	t = -t;
-      r->i = t;
-      t = zi / t;
-      r->r = 0.5 * t;
-    }
-}
diff --git a/libf2c/libF77/cabs.c b/libf2c/libF77/cabs.c
deleted file mode 100644
index 5d2142e1825..00000000000
--- a/libf2c/libF77/cabs.c
+++ /dev/null
@@ -1,24 +0,0 @@
-#undef abs
-#include <math.h>
-double
-f__cabs (double real, double imag)
-{
-  double temp;
-
-  if (real < 0)
-    real = -real;
-  if (imag < 0)
-    imag = -imag;
-  if (imag > real)
-    {
-      temp = real;
-      real = imag;
-      imag = temp;
-    }
-  if ((real + imag) == real)
-    return (real);
-
-  temp = imag / real;
-  temp = real * sqrt (1.0 + temp * temp);	/*overflow!! */
-  return (temp);
-}
diff --git a/libf2c/libF77/configure.in b/libf2c/libF77/configure.in
deleted file mode 100644
index b9f0f222fd7..00000000000
--- a/libf2c/libF77/configure.in
+++ /dev/null
@@ -1,111 +0,0 @@
-# Process this file with autoconf to produce a configure script.
-#   Copyright (C) 1995, 1997, 1998, 2001 Free Software Foundation, Inc.
-#   Contributed by Dave Love (d.love@dl.ac.uk).
-#
-#This file is part of GNU Fortran.
-#
-#GNU Fortran is free software; you can redistribute it and/or modify
-#it under the terms of the GNU General Public License as published by
-#the Free Software Foundation; either version 2, or (at your option)
-#any later version.
-#
-#GNU Fortran is distributed in the hope that it will be useful,
-#but WITHOUT ANY WARRANTY; without even the implied warranty of
-#MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
-#GNU General Public License for more details.
-#
-#You should have received a copy of the GNU General Public License
-#along with GNU Fortran; see the file COPYING.  If not, write to
-#the Free Software Foundation, 59 Temple Place - Suite 330, Boston, MA
-#02111-1307, USA.
-
-AC_PREREQ(2.13)
-AC_INIT(getarg_.c)
-
-dnl Checks for programs.
-
-dnl FIXME AC_PROG_CC wants CC to be able to link things, but it may
-dnl not be able to.
-define([AC_PROG_CC_WORKS],[])
-
-# For g77 we'll set CC to point at the built gcc, but this will get it into
-# the makefiles
-AC_PROG_CC
-
-LIBTOOL='$(SHELL) ../libtool'
-AC_SUBST(LIBTOOL)
-
-test "$AR" || AR=ar
-AC_SUBST(AR)
-if test "$RANLIB"; then :
-  AC_SUBST(RANLIB)
-else
-  AC_PROG_RANLIB
-fi
-AC_PROG_MAKE_SET
-
-dnl Checks for libraries.
-
-dnl Checks for header files.
-# Sanity check for the cross-compilation case:
-AC_CHECK_HEADER(stdio.h,:,
-  [AC_MSG_ERROR([Can't find stdio.h.
-You must have a usable C system for the target already installed, at least
-including headers and, preferably, the library, before you can configure
-the G77 runtime system.  If necessary, install gcc now with \`LANGUAGES=c',
-then the target library, then build with \`LANGUAGES=f77'.])])
-
-AC_HEADER_STDC
-AC_MSG_CHECKING(for posix)
-AC_CACHE_VAL(g77_cv_header_posix,
-  AC_EGREP_CPP(yes,
-  [#include <sys/types.h>
-#include <unistd.h>
-#ifdef _POSIX_VERSION
-  yes
-#endif
-], 
-  g77_cv_header_posix=yes,
-  g77_cv_header_posix=no))
-AC_MSG_RESULT($g77_cv_header_posix)
-
-# We can rely on the GNU library being posix-ish.  I guess checking the 
-# header isn't actually like checking the functions, though...
-AC_MSG_CHECKING(for GNU library)
-AC_CACHE_VAL(g77_cv_lib_gnu,
-  AC_EGREP_CPP(yes,
-    [#include <stdio.h>
-#ifdef __GNU_LIBRARY__
-  yes
-#endif
-], 
-  g77_cv_lib_gnu=yes, g77_cv_lib_gnu=no))
-AC_MSG_RESULT($g77_cv_lib_gnu)
-
-dnl Checks for library functions.
-AC_TYPE_SIGNAL
-# we'll get atexit by default
-if test $ac_cv_header_stdc != yes; then
-AC_CHECK_FUNC(atexit,
-  AC_DEFINE(onexit,atexit),dnl just in case
-  [AC_DEFINE(NO_ONEXIT)
-  AC_CHECK_FUNC(onexit,,
-    [AC_CHECK_FUNC(on_exit,
-     AC_DEFINE(onexit,on_exit),)])])
-else true
-fi
-
-dnl perhaps should check also for remainder
-dnl Unfortunately, the message implies we're just checking for -lm...
-AC_CHECK_LIB(m,drem,AC_DEFINE(IEEE_drem))
-
-AC_DEFINE(Skip_f2c_Undefs)
-AC_DEFINE(IEEE_COMPLEX_DIVIDE)
-
-AC_OUTPUT(Makefile)
-
-dnl Local Variables:
-dnl comment-start: "dnl "
-dnl comment-end: ""
-dnl comment-start-skip: "\\bdnl\\b\\s *"
-dnl End:
diff --git a/libf2c/libF77/d_abs.c b/libf2c/libF77/d_abs.c
deleted file mode 100644
index a43a5c7f0a8..00000000000
--- a/libf2c/libF77/d_abs.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-double
-d_abs (doublereal * x)
-{
-  if (*x >= 0)
-    return (*x);
-  return (-*x);
-}
diff --git a/libf2c/libF77/d_acos.c b/libf2c/libF77/d_acos.c
deleted file mode 100644
index 41c4f172c65..00000000000
--- a/libf2c/libF77/d_acos.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include <math.h>
-double
-d_acos (doublereal * x)
-{
-  return (acos (*x));
-}
diff --git a/libf2c/libF77/d_asin.c b/libf2c/libF77/d_asin.c
deleted file mode 100644
index 65603891501..00000000000
--- a/libf2c/libF77/d_asin.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include <math.h>
-double
-d_asin (doublereal * x)
-{
-  return (asin (*x));
-}
diff --git a/libf2c/libF77/d_atan.c b/libf2c/libF77/d_atan.c
deleted file mode 100644
index e25fa2e4177..00000000000
--- a/libf2c/libF77/d_atan.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include <math.h>
-double
-d_atan (doublereal * x)
-{
-  return (atan (*x));
-}
diff --git a/libf2c/libF77/d_atn2.c b/libf2c/libF77/d_atn2.c
deleted file mode 100644
index e0b217863a5..00000000000
--- a/libf2c/libF77/d_atn2.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include <math.h>
-double
-d_atn2 (doublereal * x, doublereal * y)
-{
-  return (atan2 (*x, *y));
-}
diff --git a/libf2c/libF77/d_cnjg.c b/libf2c/libF77/d_cnjg.c
deleted file mode 100644
index bc46ae69598..00000000000
--- a/libf2c/libF77/d_cnjg.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-void
-d_cnjg (doublecomplex * r, doublecomplex * z)
-{
-  doublereal zi = z->i;
-  r->r = z->r;
-  r->i = -zi;
-}
diff --git a/libf2c/libF77/d_cos.c b/libf2c/libF77/d_cos.c
deleted file mode 100644
index 010db6b5421..00000000000
--- a/libf2c/libF77/d_cos.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include <math.h>
-double
-d_cos (doublereal * x)
-{
-  return (cos (*x));
-}
diff --git a/libf2c/libF77/d_cosh.c b/libf2c/libF77/d_cosh.c
deleted file mode 100644
index 00938bda13e..00000000000
--- a/libf2c/libF77/d_cosh.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include <math.h>
-double
-d_cosh (doublereal * x)
-{
-  return (cosh (*x));
-}
diff --git a/libf2c/libF77/d_dim.c b/libf2c/libF77/d_dim.c
deleted file mode 100644
index a4274ca6859..00000000000
--- a/libf2c/libF77/d_dim.c
+++ /dev/null
@@ -1,7 +0,0 @@
-#include "f2c.h"
-
-double
-d_dim (doublereal * a, doublereal * b)
-{
-  return (*a > *b ? *a - *b : 0);
-}
diff --git a/libf2c/libF77/d_exp.c b/libf2c/libF77/d_exp.c
deleted file mode 100644
index 7b4f3e529c4..00000000000
--- a/libf2c/libF77/d_exp.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include <math.h>
-double
-d_exp (doublereal * x)
-{
-  return (exp (*x));
-}
diff --git a/libf2c/libF77/d_imag.c b/libf2c/libF77/d_imag.c
deleted file mode 100644
index cc937644b75..00000000000
--- a/libf2c/libF77/d_imag.c
+++ /dev/null
@@ -1,7 +0,0 @@
-#include "f2c.h"
-
-double
-d_imag (doublecomplex * z)
-{
-  return (z->i);
-}
diff --git a/libf2c/libF77/d_int.c b/libf2c/libF77/d_int.c
deleted file mode 100644
index f7ab8b0cb6a..00000000000
--- a/libf2c/libF77/d_int.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include <math.h>
-double
-d_int (doublereal * x)
-{
-  return ((*x > 0) ? floor (*x) : -floor (-*x));
-}
diff --git a/libf2c/libF77/d_lg10.c b/libf2c/libF77/d_lg10.c
deleted file mode 100644
index d0f86e554d0..00000000000
--- a/libf2c/libF77/d_lg10.c
+++ /dev/null
@@ -1,11 +0,0 @@
-#include "f2c.h"
-
-#define log10e 0.43429448190325182765
-
-#undef abs
-#include <math.h>
-double
-d_lg10 (doublereal * x)
-{
-  return (log10e * log (*x));
-}
diff --git a/libf2c/libF77/d_log.c b/libf2c/libF77/d_log.c
deleted file mode 100644
index 95dc767fb2a..00000000000
--- a/libf2c/libF77/d_log.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include <math.h>
-double
-d_log (doublereal * x)
-{
-  return (log (*x));
-}
diff --git a/libf2c/libF77/d_mod.c b/libf2c/libF77/d_mod.c
deleted file mode 100644
index 15cedefadd1..00000000000
--- a/libf2c/libF77/d_mod.c
+++ /dev/null
@@ -1,33 +0,0 @@
-#include "f2c.h"
-
-#ifdef IEEE_drem
-double drem (double, double);
-#else
-#undef abs
-#include <math.h>
-#endif
-double
-d_mod (doublereal * x, doublereal * y)
-{
-#ifdef IEEE_drem
-  double xa, ya, z;
-  if ((ya = *y) < 0.)
-    ya = -ya;
-  z = drem (xa = *x, ya);
-  if (xa > 0)
-    {
-      if (z < 0)
-	z += ya;
-    }
-  else if (z > 0)
-    z -= ya;
-  return z;
-#else
-  double quotient;
-  if ((quotient = *x / *y) >= 0)
-    quotient = floor (quotient);
-  else
-    quotient = -floor (-quotient);
-  return (*x - (*y) * quotient);
-#endif
-}
diff --git a/libf2c/libF77/d_nint.c b/libf2c/libF77/d_nint.c
deleted file mode 100644
index 8be42758644..00000000000
--- a/libf2c/libF77/d_nint.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include <math.h>
-double
-d_nint (doublereal * x)
-{
-  return ((*x) >= 0 ? floor (*x + .5) : -floor (.5 - *x));
-}
diff --git a/libf2c/libF77/d_prod.c b/libf2c/libF77/d_prod.c
deleted file mode 100644
index 11fe2c7f68e..00000000000
--- a/libf2c/libF77/d_prod.c
+++ /dev/null
@@ -1,7 +0,0 @@
-#include "f2c.h"
-
-double
-d_prod (real * x, real * y)
-{
-  return ((*x) * (*y));
-}
diff --git a/libf2c/libF77/d_sign.c b/libf2c/libF77/d_sign.c
deleted file mode 100644
index da8d24ba7a1..00000000000
--- a/libf2c/libF77/d_sign.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-double
-d_sign (doublereal * a, doublereal * b)
-{
-  double x;
-  x = (*a >= 0 ? *a : -*a);
-  return (*b >= 0 ? x : -x);
-}
diff --git a/libf2c/libF77/d_sin.c b/libf2c/libF77/d_sin.c
deleted file mode 100644
index 24b37a439e5..00000000000
--- a/libf2c/libF77/d_sin.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include <math.h>
-double
-d_sin (doublereal * x)
-{
-  return (sin (*x));
-}
diff --git a/libf2c/libF77/d_sinh.c b/libf2c/libF77/d_sinh.c
deleted file mode 100644
index dc9dc439073..00000000000
--- a/libf2c/libF77/d_sinh.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include <math.h>
-double
-d_sinh (doublereal * x)
-{
-  return (sinh (*x));
-}
diff --git a/libf2c/libF77/d_sqrt.c b/libf2c/libF77/d_sqrt.c
deleted file mode 100644
index 0a458823318..00000000000
--- a/libf2c/libF77/d_sqrt.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include <math.h>
-double
-d_sqrt (doublereal * x)
-{
-  return (sqrt (*x));
-}
diff --git a/libf2c/libF77/d_tan.c b/libf2c/libF77/d_tan.c
deleted file mode 100644
index 370c1b85fef..00000000000
--- a/libf2c/libF77/d_tan.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include <math.h>
-double
-d_tan (doublereal * x)
-{
-  return (tan (*x));
-}
diff --git a/libf2c/libF77/d_tanh.c b/libf2c/libF77/d_tanh.c
deleted file mode 100644
index df81ea0df22..00000000000
--- a/libf2c/libF77/d_tanh.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include <math.h>
-double
-d_tanh (doublereal * x)
-{
-  return (tanh (*x));
-}
diff --git a/libf2c/libF77/derf_.c b/libf2c/libF77/derf_.c
deleted file mode 100644
index b78fde0e9b4..00000000000
--- a/libf2c/libF77/derf_.c
+++ /dev/null
@@ -1,8 +0,0 @@
-#include "f2c.h"
-
-extern double erf (double);
-double
-G77_derf_0 (doublereal * x)
-{
-  return (erf (*x));
-}
diff --git a/libf2c/libF77/derfc_.c b/libf2c/libF77/derfc_.c
deleted file mode 100644
index 78e8e889c38..00000000000
--- a/libf2c/libF77/derfc_.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-extern double erfc (double);
-
-double
-G77_derfc_0 (doublereal * x)
-{
-  return (erfc (*x));
-}
diff --git a/libf2c/libF77/dtime_.c b/libf2c/libF77/dtime_.c
deleted file mode 100644
index e2ea1c6b5d6..00000000000
--- a/libf2c/libF77/dtime_.c
+++ /dev/null
@@ -1,49 +0,0 @@
-#include "time.h"
-
-#ifdef MSDOS
-#undef USE_CLOCK
-#define USE_CLOCK
-#endif
-
-#ifndef USE_CLOCK
-#define _INCLUDE_POSIX_SOURCE	/* for HP-UX */
-#define _INCLUDE_XOPEN_SOURCE	/* for HP-UX */
-#include "sys/types.h"
-#include "sys/times.h"
-#endif
-
-#undef Hz
-#ifdef CLK_TCK
-#define Hz CLK_TCK
-#else
-#ifdef HZ
-#define Hz HZ
-#else
-#define Hz 60
-#endif
-#endif
-
-double
-dtime_ (float *tarray)
-{
-#ifdef USE_CLOCK
-#ifndef CLOCKS_PER_SECOND
-#define CLOCKS_PER_SECOND Hz
-#endif
-  static double t0;
-  double t = clock ();
-  tarray[1] = 0;
-  tarray[0] = (t - t0) / CLOCKS_PER_SECOND;
-  t0 = t;
-  return tarray[0];
-#else
-  struct tms t;
-  static struct tms t0;
-
-  times (&t);
-  tarray[0] = (double) (t.tms_utime - t0.tms_utime) / Hz;
-  tarray[1] = (double) (t.tms_stime - t0.tms_stime) / Hz;
-  t0 = t;
-  return tarray[0] + tarray[1];
-#endif
-}
diff --git a/libf2c/libF77/ef1asc_.c b/libf2c/libF77/ef1asc_.c
deleted file mode 100644
index d9bea344e4f..00000000000
--- a/libf2c/libF77/ef1asc_.c
+++ /dev/null
@@ -1,15 +0,0 @@
-/* EFL support routine to copy string b to string a */
-
-#include "f2c.h"
-
-
-#define M	( (long) (sizeof(long) - 1) )
-#define EVEN(x)	( ( (x)+ M) & (~M) )
-
-extern void s_copy (char *, char *, ftnlen, ftnlen);
-int
-G77_ef1asc_0 (ftnint * a, ftnlen * la, ftnint * b, ftnlen * lb)
-{
-  s_copy ((char *) a, (char *) b, EVEN (*la), *lb);
-  return 0;			/* ignored return value */
-}
diff --git a/libf2c/libF77/ef1cmc_.c b/libf2c/libF77/ef1cmc_.c
deleted file mode 100644
index 2e102fb6b5d..00000000000
--- a/libf2c/libF77/ef1cmc_.c
+++ /dev/null
@@ -1,10 +0,0 @@
-/* EFL support routine to compare two character strings */
-
-#include "f2c.h"
-
-extern integer s_cmp (char *, char *, ftnlen, ftnlen);
-integer
-G77_ef1cmc_0 (ftnint * a, ftnlen * la, ftnint * b, ftnlen * lb)
-{
-  return (s_cmp ((char *) a, (char *) b, *la, *lb));
-}
diff --git a/libf2c/libF77/erf_.c b/libf2c/libF77/erf_.c
deleted file mode 100644
index fadbfaf038a..00000000000
--- a/libf2c/libF77/erf_.c
+++ /dev/null
@@ -1,8 +0,0 @@
-#include "f2c.h"
-
-extern double erf (double);
-double
-G77_erf_0 (real * x)
-{
-  return (erf (*x));
-}
diff --git a/libf2c/libF77/erfc_.c b/libf2c/libF77/erfc_.c
deleted file mode 100644
index 7f3ff8a5dbd..00000000000
--- a/libf2c/libF77/erfc_.c
+++ /dev/null
@@ -1,8 +0,0 @@
-#include "f2c.h"
-
-extern double erfc (double);
-double
-G77_erfc_0 (real * x)
-{
-  return (erfc (*x));
-}
diff --git a/libf2c/libF77/etime_.c b/libf2c/libF77/etime_.c
deleted file mode 100644
index cc64612c8d9..00000000000
--- a/libf2c/libF77/etime_.c
+++ /dev/null
@@ -1,43 +0,0 @@
-#include "time.h"
-
-#ifdef MSDOS
-#undef USE_CLOCK
-#define USE_CLOCK
-#endif
-
-#ifndef USE_CLOCK
-#define _INCLUDE_POSIX_SOURCE	/* for HP-UX */
-#define _INCLUDE_XOPEN_SOURCE	/* for HP-UX */
-#include "sys/types.h"
-#include "sys/times.h"
-#endif
-
-#undef Hz
-#ifdef CLK_TCK
-#define Hz CLK_TCK
-#else
-#ifdef HZ
-#define Hz HZ
-#else
-#define Hz 60
-#endif
-#endif
-
-double
-etime_ (float *tarray)
-{
-#ifdef USE_CLOCK
-#ifndef CLOCKS_PER_SECOND
-#define CLOCKS_PER_SECOND Hz
-#endif
-  double t = clock ();
-  tarray[1] = 0;
-  return tarray[0] = t / CLOCKS_PER_SECOND;
-#else
-  struct tms t;
-
-  times (&t);
-  return (tarray[0] = (double) t.tms_utime / Hz)
-    + (tarray[1] = (double) t.tms_stime / Hz);
-#endif
-}
diff --git a/libf2c/libF77/exit_.c b/libf2c/libF77/exit_.c
deleted file mode 100644
index adf3d859274..00000000000
--- a/libf2c/libF77/exit_.c
+++ /dev/null
@@ -1,25 +0,0 @@
-/* This gives the effect of
-
-	subroutine exit(rc)
-	integer*4 rc
-	stop
-	end
-
- * with the added side effect of supplying rc as the program's exit code.
- */
-
-#include "f2c.h"
-#undef abs
-#undef min
-#undef max
-#include <stdlib.h>
-extern void f_exit (void);
-
-void
-G77_exit_0 (integer * rc)
-{
-#ifdef NO_ONEXIT
-  f_exit ();
-#endif
-  exit (*rc);
-}
diff --git a/libf2c/libF77/f2ch.add b/libf2c/libF77/f2ch.add
deleted file mode 100644
index 04b13e8de5d..00000000000
--- a/libf2c/libF77/f2ch.add
+++ /dev/null
@@ -1,163 +0,0 @@
-/* If you are using a C++ compiler, append the following to f2c.h
-   for compiling libF77 and libI77. */
-
-#ifdef __cplusplus
-extern "C"
-{
-  extern int abort_ (void);
-  extern double c_abs (complex *);
-  extern void c_cos (complex *, complex *);
-  extern void c_div (complex *, complex *, complex *);
-  extern void c_exp (complex *, complex *);
-  extern void c_log (complex *, complex *);
-  extern void c_sin (complex *, complex *);
-  extern void c_sqrt (complex *, complex *);
-  extern double d_abs (double *);
-  extern double d_acos (double *);
-  extern double d_asin (double *);
-  extern double d_atan (double *);
-  extern double d_atn2 (double *, double *);
-  extern void d_cnjg (doublecomplex *, doublecomplex *);
-  extern double d_cos (double *);
-  extern double d_cosh (double *);
-  extern double d_dim (double *, double *);
-  extern double d_exp (double *);
-  extern double d_imag (doublecomplex *);
-  extern double d_int (double *);
-  extern double d_lg10 (double *);
-  extern double d_log (double *);
-  extern double d_mod (double *, double *);
-  extern double d_nint (double *);
-  extern double d_prod (float *, float *);
-  extern double d_sign (double *, double *);
-  extern double d_sin (double *);
-  extern double d_sinh (double *);
-  extern double d_sqrt (double *);
-  extern double d_tan (double *);
-  extern double d_tanh (double *);
-  extern double derf_ (double *);
-  extern double derfc_ (double *);
-  extern integer do_fio (ftnint *, char *, ftnlen);
-  extern integer do_lio (ftnint *, ftnint *, char *, ftnlen);
-  extern integer do_uio (ftnint *, char *, ftnlen);
-  extern integer e_rdfe (void);
-  extern integer e_rdue (void);
-  extern integer e_rsfe (void);
-  extern integer e_rsfi (void);
-  extern integer e_rsle (void);
-  extern integer e_rsli (void);
-  extern integer e_rsue (void);
-  extern integer e_wdfe (void);
-  extern integer e_wdue (void);
-  extern integer e_wsfe (void);
-  extern integer e_wsfi (void);
-  extern integer e_wsle (void);
-  extern integer e_wsli (void);
-  extern integer e_wsue (void);
-  extern int ef1asc_ (ftnint *, ftnlen *, ftnint *, ftnlen *);
-  extern integer ef1cmc_ (ftnint *, ftnlen *, ftnint *, ftnlen *);
-  extern double erf (double);
-  extern double erf_ (float *);
-  extern double erfc (double);
-  extern double erfc_ (float *);
-  extern integer f_back (alist *);
-  extern integer f_clos (cllist *);
-  extern integer f_end (alist *);
-  extern void f_exit (void);
-  extern integer f_inqu (inlist *);
-  extern integer f_open (olist *);
-  extern integer f_rew (alist *);
-  extern int flush_ (void);
-  extern void getarg_ (integer *, char *, ftnlen);
-  extern void getenv_ (char *, char *, ftnlen, ftnlen);
-  extern short h_abs (short *);
-  extern short h_dim (short *, short *);
-  extern short h_dnnt (double *);
-  extern short h_indx (char *, char *, ftnlen, ftnlen);
-  extern short h_len (char *, ftnlen);
-  extern short h_mod (short *, short *);
-  extern short h_nint (float *);
-  extern short h_sign (short *, short *);
-  extern short hl_ge (char *, char *, ftnlen, ftnlen);
-  extern short hl_gt (char *, char *, ftnlen, ftnlen);
-  extern short hl_le (char *, char *, ftnlen, ftnlen);
-  extern short hl_lt (char *, char *, ftnlen, ftnlen);
-  extern integer i_abs (integer *);
-  extern integer i_dim (integer *, integer *);
-  extern integer i_dnnt (double *);
-  extern integer i_indx (char *, char *, ftnlen, ftnlen);
-  extern integer i_len (char *, ftnlen);
-  extern integer i_mod (integer *, integer *);
-  extern integer i_nint (float *);
-  extern integer i_sign (integer *, integer *);
-  extern integer iargc_ (void);
-  extern ftnlen l_ge (char *, char *, ftnlen, ftnlen);
-  extern ftnlen l_gt (char *, char *, ftnlen, ftnlen);
-  extern ftnlen l_le (char *, char *, ftnlen, ftnlen);
-  extern ftnlen l_lt (char *, char *, ftnlen, ftnlen);
-  extern void pow_ci (complex *, complex *, integer *);
-  extern double pow_dd (double *, double *);
-  extern double pow_di (double *, integer *);
-  extern short pow_hh (short *, shortint *);
-  extern integer pow_ii (integer *, integer *);
-  extern double pow_ri (float *, integer *);
-  extern void pow_zi (doublecomplex *, doublecomplex *, integer *);
-  extern void pow_zz (doublecomplex *, doublecomplex *, doublecomplex *);
-  extern double r_abs (float *);
-  extern double r_acos (float *);
-  extern double r_asin (float *);
-  extern double r_atan (float *);
-  extern double r_atn2 (float *, float *);
-  extern void r_cnjg (complex *, complex *);
-  extern double r_cos (float *);
-  extern double r_cosh (float *);
-  extern double r_dim (float *, float *);
-  extern double r_exp (float *);
-  extern double r_imag (complex *);
-  extern double r_int (float *);
-  extern double r_lg10 (float *);
-  extern double r_log (float *);
-  extern double r_mod (float *, float *);
-  extern double r_nint (float *);
-  extern double r_sign (float *, float *);
-  extern double r_sin (float *);
-  extern double r_sinh (float *);
-  extern double r_sqrt (float *);
-  extern double r_tan (float *);
-  extern double r_tanh (float *);
-  extern void s_cat (char *, char **, integer *, integer *, ftnlen);
-  extern integer s_cmp (char *, char *, ftnlen, ftnlen);
-  extern void s_copy (char *, char *, ftnlen, ftnlen);
-  extern int s_paus (char *, ftnlen);
-  extern integer s_rdfe (cilist *);
-  extern integer s_rdue (cilist *);
-  extern integer s_rnge (char *, integer, char *, integer);
-  extern integer s_rsfe (cilist *);
-  extern integer s_rsfi (icilist *);
-  extern integer s_rsle (cilist *);
-  extern integer s_rsli (icilist *);
-  extern integer s_rsne (cilist *);
-  extern integer s_rsni (icilist *);
-  extern integer s_rsue (cilist *);
-  extern int s_stop (char *, ftnlen);
-  extern integer s_wdfe (cilist *);
-  extern integer s_wdue (cilist *);
-  extern integer s_wsfe (cilist *);
-  extern integer s_wsfi (icilist *);
-  extern integer s_wsle (cilist *);
-  extern integer s_wsli (icilist *);
-  extern integer s_wsne (cilist *);
-  extern integer s_wsni (icilist *);
-  extern integer s_wsue (cilist *);
-  extern void sig_die (char *, int);
-  extern integer signal_ (integer *, void (*)(int));
-  extern integer system_ (char *, ftnlen);
-  extern double z_abs (doublecomplex *);
-  extern void z_cos (doublecomplex *, doublecomplex *);
-  extern void z_div (doublecomplex *, doublecomplex *, doublecomplex *);
-  extern void z_exp (doublecomplex *, doublecomplex *);
-  extern void z_log (doublecomplex *, doublecomplex *);
-  extern void z_sin (doublecomplex *, doublecomplex *);
-  extern void z_sqrt (doublecomplex *, doublecomplex *);
-}
-#endif
diff --git a/libf2c/libF77/getarg_.c b/libf2c/libF77/getarg_.c
deleted file mode 100644
index b35043b4a48..00000000000
--- a/libf2c/libF77/getarg_.c
+++ /dev/null
@@ -1,25 +0,0 @@
-#include "f2c.h"
-
-/*
- * subroutine getarg(k, c)
- * returns the kth unix command argument in fortran character
- * variable argument c
-*/
-
-void
-G77_getarg_0 (ftnint * n, register char *s, ftnlen ls)
-{
-  extern int f__xargc;
-  extern char **f__xargv;
-  register char *t;
-  register int i;
-
-  if (*n >= 0 && *n < f__xargc)
-    t = f__xargv[*n];
-  else
-    t = "";
-  for (i = 0; i < ls && *t != '\0'; ++i)
-    *s++ = *t++;
-  for (; i < ls; ++i)
-    *s++ = ' ';
-}
diff --git a/libf2c/libF77/getenv_.c b/libf2c/libF77/getenv_.c
deleted file mode 100644
index 4067b8ce214..00000000000
--- a/libf2c/libF77/getenv_.c
+++ /dev/null
@@ -1,49 +0,0 @@
-#include "f2c.h"
-#undef abs
-#include <stdlib.h>
-#include <string.h>
-extern char *F77_aloc (ftnlen, char *);
-
-/*
- * getenv - f77 subroutine to return environment variables
- *
- * called by:
- *	call getenv (ENV_NAME, char_var)
- * where:
- *	ENV_NAME is the name of an environment variable
- *	char_var is a character variable which will receive
- *		the current value of ENV_NAME, or all blanks
- *		if ENV_NAME is not defined
- */
-
-void
-G77_getenv_0 (char *fname, char *value, ftnlen flen, ftnlen vlen)
-{
-  char buf[256], *ep, *fp;
-  integer i;
-
-  if (flen <= 0)
-    goto add_blanks;
-  for (i = 0; i < (integer) sizeof (buf); i++)
-    {
-      if (i == flen || (buf[i] = fname[i]) == ' ')
-	{
-	  buf[i] = 0;
-	  ep = getenv (buf);
-	  goto have_ep;
-	}
-    }
-  while (i < flen && fname[i] != ' ')
-    i++;
-  strncpy (fp = F77_aloc (i + 1, "getenv_"), fname, (int) i);
-  fp[i] = 0;
-  ep = getenv (fp);
-  free (fp);
-have_ep:
-  if (ep)
-    while (*ep && vlen-- > 0)
-      *value++ = *ep++;
-add_blanks:
-  while (vlen-- > 0)
-    *value++ = ' ';
-}
diff --git a/libf2c/libF77/h_abs.c b/libf2c/libF77/h_abs.c
deleted file mode 100644
index 9db19ca34e9..00000000000
--- a/libf2c/libF77/h_abs.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-shortint
-h_abs (shortint * x)
-{
-  if (*x >= 0)
-    return (*x);
-  return (-*x);
-}
diff --git a/libf2c/libF77/h_dim.c b/libf2c/libF77/h_dim.c
deleted file mode 100644
index 1519478747d..00000000000
--- a/libf2c/libF77/h_dim.c
+++ /dev/null
@@ -1,7 +0,0 @@
-#include "f2c.h"
-
-shortint
-h_dim (shortint * a, shortint * b)
-{
-  return (*a > *b ? *a - *b : 0);
-}
diff --git a/libf2c/libF77/h_dnnt.c b/libf2c/libF77/h_dnnt.c
deleted file mode 100644
index 46c83bbd28e..00000000000
--- a/libf2c/libF77/h_dnnt.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include <math.h>
-shortint
-h_dnnt (doublereal * x)
-{
-  return (shortint) (*x >= 0. ? floor (*x + .5) : -floor (.5 - *x));
-}
diff --git a/libf2c/libF77/h_indx.c b/libf2c/libF77/h_indx.c
deleted file mode 100644
index 2353b2b143a..00000000000
--- a/libf2c/libF77/h_indx.c
+++ /dev/null
@@ -1,23 +0,0 @@
-#include "f2c.h"
-
-shortint
-h_indx (char *a, char *b, ftnlen la, ftnlen lb)
-{
-  ftnlen i, n;
-  char *s, *t, *bend;
-
-  n = la - lb + 1;
-  bend = b + lb;
-
-  for (i = 0; i < n; ++i)
-    {
-      s = a + i;
-      t = b;
-      while (t < bend)
-	if (*s++ != *t++)
-	  goto no;
-      return ((shortint) i + 1);
-    no:;
-    }
-  return (0);
-}
diff --git a/libf2c/libF77/h_len.c b/libf2c/libF77/h_len.c
deleted file mode 100644
index 506016e67d9..00000000000
--- a/libf2c/libF77/h_len.c
+++ /dev/null
@@ -1,7 +0,0 @@
-#include "f2c.h"
-
-shortint
-h_len (char *s __attribute__ ((__unused__)), ftnlen n)
-{
-  return (n);
-}
diff --git a/libf2c/libF77/h_mod.c b/libf2c/libF77/h_mod.c
deleted file mode 100644
index c04e0df823b..00000000000
--- a/libf2c/libF77/h_mod.c
+++ /dev/null
@@ -1,7 +0,0 @@
-#include "f2c.h"
-
-shortint
-h_mod (short *a, short *b)
-{
-  return (*a % *b);
-}
diff --git a/libf2c/libF77/h_nint.c b/libf2c/libF77/h_nint.c
deleted file mode 100644
index a8c366a41be..00000000000
--- a/libf2c/libF77/h_nint.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include <math.h>
-shortint
-h_nint (real * x)
-{
-  return (shortint) (*x >= 0 ? floor (*x + .5) : -floor (.5 - *x));
-}
diff --git a/libf2c/libF77/h_sign.c b/libf2c/libF77/h_sign.c
deleted file mode 100644
index 70402325be2..00000000000
--- a/libf2c/libF77/h_sign.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-shortint
-h_sign (shortint * a, shortint * b)
-{
-  shortint x;
-  x = (*a >= 0 ? *a : -*a);
-  return (*b >= 0 ? x : -x);
-}
diff --git a/libf2c/libF77/hl_ge.c b/libf2c/libF77/hl_ge.c
deleted file mode 100644
index 988686d8d1d..00000000000
--- a/libf2c/libF77/hl_ge.c
+++ /dev/null
@@ -1,8 +0,0 @@
-#include "f2c.h"
-
-extern integer s_cmp (char *, char *, ftnlen, ftnlen);
-shortlogical
-hl_ge (char *a, char *b, ftnlen la, ftnlen lb)
-{
-  return (s_cmp (a, b, la, lb) >= 0);
-}
diff --git a/libf2c/libF77/hl_gt.c b/libf2c/libF77/hl_gt.c
deleted file mode 100644
index 0024ca7a0d2..00000000000
--- a/libf2c/libF77/hl_gt.c
+++ /dev/null
@@ -1,8 +0,0 @@
-#include "f2c.h"
-
-extern integer s_cmp (char *, char *, ftnlen, ftnlen);
-shortlogical
-hl_gt (char *a, char *b, ftnlen la, ftnlen lb)
-{
-  return (s_cmp (a, b, la, lb) > 0);
-}
diff --git a/libf2c/libF77/hl_le.c b/libf2c/libF77/hl_le.c
deleted file mode 100644
index 76aa3e12da7..00000000000
--- a/libf2c/libF77/hl_le.c
+++ /dev/null
@@ -1,8 +0,0 @@
-#include "f2c.h"
-
-extern integer s_cmp (char *, char *, ftnlen, ftnlen);
-shortlogical
-hl_le (char *a, char *b, ftnlen la, ftnlen lb)
-{
-  return (s_cmp (a, b, la, lb) <= 0);
-}
diff --git a/libf2c/libF77/hl_lt.c b/libf2c/libF77/hl_lt.c
deleted file mode 100644
index 68a47fa98eb..00000000000
--- a/libf2c/libF77/hl_lt.c
+++ /dev/null
@@ -1,8 +0,0 @@
-#include "f2c.h"
-
-extern integer s_cmp (char *, char *, ftnlen, ftnlen);
-shortlogical
-hl_lt (char *a, char *b, ftnlen la, ftnlen lb)
-{
-  return (s_cmp (a, b, la, lb) < 0);
-}
diff --git a/libf2c/libF77/i_abs.c b/libf2c/libF77/i_abs.c
deleted file mode 100644
index 2ed183a6225..00000000000
--- a/libf2c/libF77/i_abs.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-integer
-i_abs (integer * x)
-{
-  if (*x >= 0)
-    return (*x);
-  return (-*x);
-}
diff --git a/libf2c/libF77/i_dim.c b/libf2c/libF77/i_dim.c
deleted file mode 100644
index 66ef7c99fb4..00000000000
--- a/libf2c/libF77/i_dim.c
+++ /dev/null
@@ -1,7 +0,0 @@
-#include "f2c.h"
-
-integer
-i_dim (integer * a, integer * b)
-{
-  return (*a > *b ? *a - *b : 0);
-}
diff --git a/libf2c/libF77/i_dnnt.c b/libf2c/libF77/i_dnnt.c
deleted file mode 100644
index 7a3783d29e2..00000000000
--- a/libf2c/libF77/i_dnnt.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include <math.h>
-integer
-i_dnnt (doublereal * x)
-{
-  return (integer) (*x >= 0. ? floor (*x + .5) : -floor (.5 - *x));
-}
diff --git a/libf2c/libF77/i_indx.c b/libf2c/libF77/i_indx.c
deleted file mode 100644
index 5b8e13693ba..00000000000
--- a/libf2c/libF77/i_indx.c
+++ /dev/null
@@ -1,23 +0,0 @@
-#include "f2c.h"
-
-integer
-i_indx (char *a, char *b, ftnlen la, ftnlen lb)
-{
-  ftnlen i, n;
-  char *s, *t, *bend;
-
-  n = la - lb + 1;
-  bend = b + lb;
-
-  for (i = 0; i < n; ++i)
-    {
-      s = a + i;
-      t = b;
-      while (t < bend)
-	if (*s++ != *t++)
-	  goto no;
-      return (i + 1);
-    no:;
-    }
-  return (0);
-}
diff --git a/libf2c/libF77/i_len.c b/libf2c/libF77/i_len.c
deleted file mode 100644
index 2d5a3a4381a..00000000000
--- a/libf2c/libF77/i_len.c
+++ /dev/null
@@ -1,7 +0,0 @@
-#include "f2c.h"
-
-integer
-i_len (char *s __attribute__ ((__unused__)), ftnlen n)
-{
-  return (n);
-}
diff --git a/libf2c/libF77/i_mod.c b/libf2c/libF77/i_mod.c
deleted file mode 100644
index 7ed7b391c19..00000000000
--- a/libf2c/libF77/i_mod.c
+++ /dev/null
@@ -1,7 +0,0 @@
-#include "f2c.h"
-
-integer
-i_mod (integer * a, integer * b)
-{
-  return (*a % *b);
-}
diff --git a/libf2c/libF77/i_nint.c b/libf2c/libF77/i_nint.c
deleted file mode 100644
index c4eaff48d33..00000000000
--- a/libf2c/libF77/i_nint.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include <math.h>
-integer
-i_nint (real * x)
-{
-  return (integer) (*x >= 0 ? floor (*x + .5) : -floor (.5 - *x));
-}
diff --git a/libf2c/libF77/i_sign.c b/libf2c/libF77/i_sign.c
deleted file mode 100644
index cf090086d00..00000000000
--- a/libf2c/libF77/i_sign.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-integer
-i_sign (integer * a, integer * b)
-{
-  integer x;
-  x = (*a >= 0 ? *a : -*a);
-  return (*b >= 0 ? x : -x);
-}
diff --git a/libf2c/libF77/iargc_.c b/libf2c/libF77/iargc_.c
deleted file mode 100644
index c3165709a53..00000000000
--- a/libf2c/libF77/iargc_.c
+++ /dev/null
@@ -1,8 +0,0 @@
-#include "f2c.h"
-
-ftnint
-G77_iargc_0 (void)
-{
-  extern int f__xargc;
-  return (f__xargc - 1);
-}
diff --git a/libf2c/libF77/l_ge.c b/libf2c/libF77/l_ge.c
deleted file mode 100644
index 78af8d04f5e..00000000000
--- a/libf2c/libF77/l_ge.c
+++ /dev/null
@@ -1,8 +0,0 @@
-#include "f2c.h"
-
-extern integer s_cmp (char *, char *, ftnlen, ftnlen);
-logical
-l_ge (char *a, char *b, ftnlen la, ftnlen lb)
-{
-  return (s_cmp (a, b, la, lb) >= 0);
-}
diff --git a/libf2c/libF77/l_gt.c b/libf2c/libF77/l_gt.c
deleted file mode 100644
index be7e4894719..00000000000
--- a/libf2c/libF77/l_gt.c
+++ /dev/null
@@ -1,8 +0,0 @@
-#include "f2c.h"
-
-extern integer s_cmp (char *, char *, ftnlen, ftnlen);
-logical
-l_gt (char *a, char *b, ftnlen la, ftnlen lb)
-{
-  return (s_cmp (a, b, la, lb) > 0);
-}
diff --git a/libf2c/libF77/l_le.c b/libf2c/libF77/l_le.c
deleted file mode 100644
index d2886fb7d48..00000000000
--- a/libf2c/libF77/l_le.c
+++ /dev/null
@@ -1,8 +0,0 @@
-#include "f2c.h"
-
-extern integer s_cmp (char *, char *, ftnlen, ftnlen);
-logical
-l_le (char *a, char *b, ftnlen la, ftnlen lb)
-{
-  return (s_cmp (a, b, la, lb) <= 0);
-}
diff --git a/libf2c/libF77/l_lt.c b/libf2c/libF77/l_lt.c
deleted file mode 100644
index ff151f890a2..00000000000
--- a/libf2c/libF77/l_lt.c
+++ /dev/null
@@ -1,8 +0,0 @@
-#include "f2c.h"
-
-extern integer s_cmp (char *, char *, ftnlen, ftnlen);
-logical
-l_lt (char *a, char *b, ftnlen la, ftnlen lb)
-{
-  return (s_cmp (a, b, la, lb) < 0);
-}
diff --git a/libf2c/libF77/lbitbits.c b/libf2c/libF77/lbitbits.c
deleted file mode 100644
index 3b28ae97aba..00000000000
--- a/libf2c/libF77/lbitbits.c
+++ /dev/null
@@ -1,58 +0,0 @@
-#include "f2c.h"
-
-#ifndef LONGBITS
-#define LONGBITS 32
-#endif
-
-integer
-lbit_bits (integer a, integer b, integer len)
-{
-  /* Assume 2's complement arithmetic */
-
-  unsigned long x, y;
-
-  x = (unsigned long) a;
-  y = (unsigned long) -1L;
-  x >>= b;
-  y <<= len;
-  return (integer) (x & ~y);
-}
-
-integer
-lbit_cshift (integer a, integer b, integer len)
-{
-  unsigned long x, y, z;
-
-  x = (unsigned long) a;
-  if (len <= 0)
-    {
-      if (len == 0)
-	return 0;
-      goto full_len;
-    }
-  if (len >= LONGBITS)
-    {
-    full_len:
-      if (b >= 0)
-	{
-	  b %= LONGBITS;
-	  return (integer) (x << b | x >> (LONGBITS - b));
-	}
-      b = -b;
-      b %= LONGBITS;
-      return (integer) (x << (LONGBITS - b) | x >> b);
-    }
-  y = z = (unsigned long) -1;
-  y <<= len;
-  z &= ~y;
-  y &= x;
-  x &= z;
-  if (b >= 0)
-    {
-      b %= len;
-      return (integer) (y | (z & (x << b | x >> (len - b))));
-    }
-  b = -b;
-  b %= len;
-  return (integer) (y | (z & (x >> b | x << (len - b))));
-}
diff --git a/libf2c/libF77/lbitshft.c b/libf2c/libF77/lbitshft.c
deleted file mode 100644
index bfbb7c01833..00000000000
--- a/libf2c/libF77/lbitshft.c
+++ /dev/null
@@ -1,7 +0,0 @@
-#include "f2c.h"
-
-integer
-lbit_shift (integer a, integer b)
-{
-  return b >= 0 ? a << b : (integer) ((uinteger) a >> -b);
-}
diff --git a/libf2c/libF77/main.c b/libf2c/libF77/main.c
deleted file mode 100644
index a3955cbc1a7..00000000000
--- a/libf2c/libF77/main.c
+++ /dev/null
@@ -1,35 +0,0 @@
-/* STARTUP PROCEDURE FOR UNIX FORTRAN PROGRAMS */
-
-#include <stdio.h>
-#include "signal1.h"
-
-#include <stdlib.h>
-
-extern void f_exit (void);
-#ifndef NO_ONEXIT
-#define ONEXIT atexit
-extern int atexit (void (*)(void));
-#endif
-
-extern void f_init (void);
-extern int MAIN__ (void);
-extern void f_setarg (int, char **);
-extern void f_setsig (void);
-
-int
-main (int argc, char **argv)
-{
-  f_setarg (argc, argv);
-  f_setsig ();
-  f_init ();
-#ifndef NO_ONEXIT
-  ONEXIT (f_exit);
-#endif
-  MAIN__ ();
-#ifdef NO_ONEXIT
-  f_exit ();
-#endif
-  exit (0);			/* exit(0) rather than return(0) to bypass Cray bug */
-  return 0;			/* For compilers that complain of missing return values; */
-  /* others will complain that this is unreachable code. */
-}
diff --git a/libf2c/libF77/makefile.netlib b/libf2c/libF77/makefile.netlib
deleted file mode 100644
index 230ca7e9f93..00000000000
--- a/libf2c/libF77/makefile.netlib
+++ /dev/null
@@ -1,103 +0,0 @@
-.SUFFIXES: .c .o
-CC = cc
-SHELL = /bin/sh
-CFLAGS = -O
-
-# If your system lacks onexit() and you are not using an
-# ANSI C compiler, then you should add -DNO_ONEXIT to CFLAGS,
-# e.g., by changing the above "CFLAGS =" line to
-# CFLAGS = -O -DNO_ONEXIT
-
-# On at least some Sun systems, it is more appropriate to change the
-# "CFLAGS =" line to
-# CFLAGS = -O -Donexit=on_exit
-
-# compile, then strip unnecessary symbols
-.c.o:
-	$(CC) -c -DSkip_f2c_Undefs $(CFLAGS) $*.c
-	ld -r -x -o $*.xxx $*.o
-	mv $*.xxx $*.o
-## Under Solaris (and other systems that do not understand ld -x),
-## omit -x in the ld line above.
-## If your system does not have the ld command, comment out
-## or remove both the ld and mv lines above.
-
-MISC =	F77_aloc.o Version.o main.o s_rnge.o abort_.o getarg_.o iargc_.o \
-	getenv_.o signal_.o s_stop.o s_paus.o system_.o cabs.o\
-	derf_.o derfc_.o erf_.o erfc_.o sig_die.o exit_.o
-POW =	pow_ci.o pow_dd.o pow_di.o pow_hh.o pow_ii.o  pow_ri.o pow_zi.o pow_zz.o
-CX =	c_abs.o c_cos.o c_div.o c_exp.o c_log.o c_sin.o c_sqrt.o
-DCX =	z_abs.o z_cos.o z_div.o z_exp.o z_log.o z_sin.o z_sqrt.o
-REAL =	r_abs.o r_acos.o r_asin.o r_atan.o r_atn2.o r_cnjg.o r_cos.o\
-	r_cosh.o r_dim.o r_exp.o r_imag.o r_int.o\
-	r_lg10.o r_log.o r_mod.o r_nint.o r_sign.o\
-	r_sin.o r_sinh.o r_sqrt.o r_tan.o r_tanh.o
-DBL =	d_abs.o d_acos.o d_asin.o d_atan.o d_atn2.o\
-	d_cnjg.o d_cos.o d_cosh.o d_dim.o d_exp.o\
-	d_imag.o d_int.o d_lg10.o d_log.o d_mod.o\
-	d_nint.o d_prod.o d_sign.o d_sin.o d_sinh.o\
-	d_sqrt.o d_tan.o d_tanh.o
-INT =	i_abs.o i_dim.o i_dnnt.o i_indx.o i_len.o i_mod.o i_nint.o i_sign.o
-HALF =	h_abs.o h_dim.o h_dnnt.o h_indx.o h_len.o h_mod.o  h_nint.o h_sign.o
-CMP =	l_ge.o l_gt.o l_le.o l_lt.o hl_ge.o hl_gt.o hl_le.o hl_lt.o
-EFL =	ef1asc_.o ef1cmc_.o
-CHAR =	F77_aloc.o s_cat.o s_cmp.o s_copy.o
-F90BIT = lbitbits.o lbitshft.o
-QINT =	pow_qq.o qbitbits.o qbitshft.o
-TIME =	dtime_.o etime_.o
-
-all: signal1.h libF77.a
-
-# You may need to adjust signal1.h suitably for your system...
-signal1.h: signal1.h0
-	cp signal1.h0 signal1.h
-
-# If you get an error compiling dtime_.c or etime_.c, try adding
-# -DUSE_CLOCK to the CFLAGS assignment above; if that does not work,
-# omit $(TIME) from the dependency list for libF77.a below.
-
-# For INTEGER*8 support (which requires system-dependent adjustments to
-# f2c.h), add $(QINT) to the libf2c.a dependency list below...
-
-libF77.a : $(MISC) $(POW) $(CX) $(DCX) $(REAL) $(DBL) $(INT) \
-	$(HALF) $(CMP) $(EFL) $(CHAR) $(F90BIT) $(TIME)
-	ar r libF77.a $?
-	-ranlib libF77.a
-
-### If your system lacks ranlib, you don't need it; see README.
-
-Version.o: Version.c
-	$(CC) -c Version.c
-
-# To compile with C++, first "make f2c.h"
-f2c.h: f2ch.add
-	cat /usr/include/f2c.h f2ch.add >f2c.h
-
-install:	libF77.a
-	mv libF77.a /usr/lib
-	ranlib /usr/lib/libF77.a
-
-clean:
-	rm -f libF77.a *.o
-
-check:
-	xsum F77_aloc.c Notice README Version.c abort_.c c_abs.c c_cos.c \
-	c_div.c c_exp.c c_log.c c_sin.c c_sqrt.c cabs.c d_abs.c d_acos.c \
-	d_asin.c d_atan.c d_atn2.c d_cnjg.c d_cos.c d_cosh.c d_dim.c \
-	d_exp.c d_imag.c d_int.c d_lg10.c d_log.c d_mod.c d_nint.c \
-	d_prod.c d_sign.c d_sin.c d_sinh.c d_sqrt.c d_tan.c d_tanh.c \
-	derf_.c derfc_.c dtime_.c \
-	ef1asc_.c ef1cmc_.c erf_.c erfc_.c etime_.c exit_.c f2ch.add \
-	getarg_.c getenv_.c h_abs.c h_dim.c h_dnnt.c h_indx.c h_len.c \
-	h_mod.c h_nint.c h_sign.c hl_ge.c hl_gt.c hl_le.c hl_lt.c \
-	i_abs.c i_dim.c i_dnnt.c i_indx.c i_len.c i_mod.c i_nint.c \
-	i_sign.c iargc_.c l_ge.c l_gt.c l_le.c l_lt.c lbitbits.c lbitshft.c \
-	main.c makefile pow_ci.c pow_dd.c pow_di.c pow_hh.c pow_ii.c \
-	pow_qq.c pow_ri.c pow_zi.c pow_zz.c qbitbits.c qbitshft.c \
-	r_abs.c r_acos.c r_asin.c r_atan.c r_atn2.c \
-	r_cnjg.c r_cos.c r_cosh.c r_dim.c r_exp.c r_imag.c r_int.c r_lg10.c \
-	r_log.c r_mod.c r_nint.c r_sign.c r_sin.c r_sinh.c r_sqrt.c \
-	r_tan.c r_tanh.c s_cat.c s_cmp.c s_copy.c \
-	s_paus.c s_rnge.c s_stop.c sig_die.c signal1.h0 signal_.c system_.c \
-	z_abs.c z_cos.c z_div.c z_exp.c z_log.c z_sin.c z_sqrt.c >zap
-	cmp zap libF77.xsum && rm zap || diff libF77.xsum zap
diff --git a/libf2c/libF77/pow_ci.c b/libf2c/libF77/pow_ci.c
deleted file mode 100644
index 1df3eb34b00..00000000000
--- a/libf2c/libF77/pow_ci.c
+++ /dev/null
@@ -1,16 +0,0 @@
-#include "f2c.h"
-
-extern void pow_zi (doublecomplex *, doublecomplex *, integer *);
-void
-pow_ci (complex * p, complex * a, integer * b)	/* p = a**b  */
-{
-  doublecomplex p1, a1;
-
-  a1.r = a->r;
-  a1.i = a->i;
-
-  pow_zi (&p1, &a1, b);
-
-  p->r = p1.r;
-  p->i = p1.i;
-}
diff --git a/libf2c/libF77/pow_dd.c b/libf2c/libF77/pow_dd.c
deleted file mode 100644
index 0ab208ebf62..00000000000
--- a/libf2c/libF77/pow_dd.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include <math.h>
-double
-pow_dd (doublereal * ap, doublereal * bp)
-{
-  return (pow (*ap, *bp));
-}
diff --git a/libf2c/libF77/pow_di.c b/libf2c/libF77/pow_di.c
deleted file mode 100644
index d2298a0a1f0..00000000000
--- a/libf2c/libF77/pow_di.c
+++ /dev/null
@@ -1,32 +0,0 @@
-#include "f2c.h"
-
-double
-pow_di (doublereal * ap, integer * bp)
-{
-  double pow, x;
-  integer n;
-  unsigned long u;
-
-  pow = 1;
-  x = *ap;
-  n = *bp;
-
-  if (n != 0)
-    {
-      if (n < 0)
-	{
-	  n = -n;
-	  x = 1 / x;
-	}
-      for (u = n;;)
-	{
-	  if (u & 01)
-	    pow *= x;
-	  if (u >>= 1)
-	    x *= x;
-	  else
-	    break;
-	}
-    }
-  return (pow);
-}
diff --git a/libf2c/libF77/pow_hh.c b/libf2c/libF77/pow_hh.c
deleted file mode 100644
index 3379d8a7f7c..00000000000
--- a/libf2c/libF77/pow_hh.c
+++ /dev/null
@@ -1,31 +0,0 @@
-#include "f2c.h"
-
-shortint
-pow_hh (shortint * ap, shortint * bp)
-{
-  shortint pow, x, n;
-  unsigned u;
-
-  x = *ap;
-  n = *bp;
-
-  if (n <= 0)
-    {
-      if (n == 0 || x == 1)
-	return 1;
-      if (x != -1)
-	return x == 0 ? 1 / x : 0;
-      n = -n;
-    }
-  u = n;
-  for (pow = 1;;)
-    {
-      if (u & 01)
-	pow *= x;
-      if (u >>= 1)
-	x *= x;
-      else
-	break;
-    }
-  return (pow);
-}
diff --git a/libf2c/libF77/pow_ii.c b/libf2c/libF77/pow_ii.c
deleted file mode 100644
index 92347453fb0..00000000000
--- a/libf2c/libF77/pow_ii.c
+++ /dev/null
@@ -1,31 +0,0 @@
-#include "f2c.h"
-
-integer
-pow_ii (integer * ap, integer * bp)
-{
-  integer pow, x, n;
-  unsigned long u;
-
-  x = *ap;
-  n = *bp;
-
-  if (n <= 0)
-    {
-      if (n == 0 || x == 1)
-	return 1;
-      if (x != -1)
-	return x == 0 ? 1 / x : 0;
-      n = -n;
-    }
-  u = n;
-  for (pow = 1;;)
-    {
-      if (u & 01)
-	pow *= x;
-      if (u >>= 1)
-	x *= x;
-      else
-	break;
-    }
-  return (pow);
-}
diff --git a/libf2c/libF77/pow_qq.c b/libf2c/libF77/pow_qq.c
deleted file mode 100644
index 0cec5ca7d8c..00000000000
--- a/libf2c/libF77/pow_qq.c
+++ /dev/null
@@ -1,31 +0,0 @@
-#include "f2c.h"
-
-longint
-pow_qq (longint * ap, longint * bp)
-{
-  longint pow, x, n;
-  unsigned long long u;		/* system-dependent */
-
-  x = *ap;
-  n = *bp;
-
-  if (n <= 0)
-    {
-      if (n == 0 || x == 1)
-	return 1;
-      if (x != -1)
-	return x == 0 ? 1 / x : 0;
-      n = -n;
-    }
-  u = n;
-  for (pow = 1;;)
-    {
-      if (u & 01)
-	pow *= x;
-      if (u >>= 1)
-	x *= x;
-      else
-	break;
-    }
-  return (pow);
-}
diff --git a/libf2c/libF77/pow_ri.c b/libf2c/libF77/pow_ri.c
deleted file mode 100644
index 792db0c2018..00000000000
--- a/libf2c/libF77/pow_ri.c
+++ /dev/null
@@ -1,32 +0,0 @@
-#include "f2c.h"
-
-double
-pow_ri (real * ap, integer * bp)
-{
-  double pow, x;
-  integer n;
-  unsigned long u;
-
-  pow = 1;
-  x = *ap;
-  n = *bp;
-
-  if (n != 0)
-    {
-      if (n < 0)
-	{
-	  n = -n;
-	  x = 1 / x;
-	}
-      for (u = n;;)
-	{
-	  if (u & 01)
-	    pow *= x;
-	  if (u >>= 1)
-	    x *= x;
-	  else
-	    break;
-	}
-    }
-  return (pow);
-}
diff --git a/libf2c/libF77/pow_zi.c b/libf2c/libF77/pow_zi.c
deleted file mode 100644
index 214db3d7a0f..00000000000
--- a/libf2c/libF77/pow_zi.c
+++ /dev/null
@@ -1,50 +0,0 @@
-#include "f2c.h"
-
-extern void z_div (doublecomplex *, doublecomplex *, doublecomplex *);
-void
-pow_zi (doublecomplex * p, doublecomplex * a, integer * b)	/* p = a**b  */
-{
-  integer n;
-  unsigned long u;
-  double t;
-  doublecomplex q, x;
-  static doublecomplex one = { 1.0, 0.0 };
-
-  n = *b;
-  q.r = 1;
-  q.i = 0;
-
-  if (n == 0)
-    goto done;
-  if (n < 0)
-    {
-      n = -n;
-      z_div (&x, &one, a);
-    }
-  else
-    {
-      x.r = a->r;
-      x.i = a->i;
-    }
-
-  for (u = n;;)
-    {
-      if (u & 01)
-	{
-	  t = q.r * x.r - q.i * x.i;
-	  q.i = q.r * x.i + q.i * x.r;
-	  q.r = t;
-	}
-      if (u >>= 1)
-	{
-	  t = x.r * x.r - x.i * x.i;
-	  x.i = 2 * x.r * x.i;
-	  x.r = t;
-	}
-      else
-	break;
-    }
-done:
-  p->i = q.i;
-  p->r = q.r;
-}
diff --git a/libf2c/libF77/pow_zz.c b/libf2c/libF77/pow_zz.c
deleted file mode 100644
index d5cfbf3d931..00000000000
--- a/libf2c/libF77/pow_zz.c
+++ /dev/null
@@ -1,25 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include <math.h>
-extern double f__cabs (double, double);
-void
-pow_zz (doublecomplex * r, doublecomplex * a, doublecomplex * b)
-{
-  double logr, logi, x, y;
-
-  if (a->r == 0.0 && a->i == 0.0)
-    {
-	/* Algorithm below doesn't cope.  */
-        r->r = r->i = 0.0;
-        return;
-    }
-  logr = log (f__cabs (a->r, a->i));
-  logi = atan2 (a->i, a->r);
-
-  x = exp (logr * b->r - logi * b->i);
-  y = logr * b->i + logi * b->r;
-
-  r->r = x * cos (y);
-  r->i = x * sin (y);
-}
diff --git a/libf2c/libF77/qbitbits.c b/libf2c/libF77/qbitbits.c
deleted file mode 100644
index f72858e7a33..00000000000
--- a/libf2c/libF77/qbitbits.c
+++ /dev/null
@@ -1,62 +0,0 @@
-#include "f2c.h"
-
-#ifndef LONGBITS
-#define LONGBITS 32
-#endif
-
-#ifndef LONG8BITS
-#define LONG8BITS (2*LONGBITS)
-#endif
-
-integer
-qbit_bits (longint a, integer b, integer len)
-{
-  /* Assume 2's complement arithmetic */
-
-  ulongint x, y;
-
-  x = (ulongint) a;
-  y = (ulongint) - 1L;
-  x >>= b;
-  y <<= len;
-  return (longint) (x & y);
-}
-
-longint
-qbit_cshift (longint a, integer b, integer len)
-{
-  ulongint x, y, z;
-
-  x = (ulongint) a;
-  if (len <= 0)
-    {
-      if (len == 0)
-	return 0;
-      goto full_len;
-    }
-  if (len >= LONG8BITS)
-    {
-    full_len:
-      if (b >= 0)
-	{
-	  b %= LONG8BITS;
-	  return (longint) (x << b | x >> (LONG8BITS - b));
-	}
-      b = -b;
-      b %= LONG8BITS;
-      return (longint) (x << (LONG8BITS - b) | x >> b);
-    }
-  y = z = (unsigned long) -1;
-  y <<= len;
-  z &= ~y;
-  y &= x;
-  x &= z;
-  if (b >= 0)
-    {
-      b %= len;
-      return (longint) (y | (z & (x << b | x >> (len - b))));
-    }
-  b = -b;
-  b %= len;
-  return (longint) (y | (z & (x >> b | x << (len - b))));
-}
diff --git a/libf2c/libF77/qbitshft.c b/libf2c/libF77/qbitshft.c
deleted file mode 100644
index ce740edde19..00000000000
--- a/libf2c/libF77/qbitshft.c
+++ /dev/null
@@ -1,7 +0,0 @@
-#include "f2c.h"
-
-longint
-qbit_shift (longint a, integer b)
-{
-  return b >= 0 ? a << b : (longint) ((ulongint) a >> -b);
-}
diff --git a/libf2c/libF77/r_abs.c b/libf2c/libF77/r_abs.c
deleted file mode 100644
index 6f62724ef17..00000000000
--- a/libf2c/libF77/r_abs.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-double
-r_abs (real * x)
-{
-  if (*x >= 0)
-    return (*x);
-  return (-*x);
-}
diff --git a/libf2c/libF77/r_acos.c b/libf2c/libF77/r_acos.c
deleted file mode 100644
index d761cfdc5c8..00000000000
--- a/libf2c/libF77/r_acos.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include <math.h>
-double
-r_acos (real * x)
-{
-  return (acos (*x));
-}
diff --git a/libf2c/libF77/r_asin.c b/libf2c/libF77/r_asin.c
deleted file mode 100644
index b8c73c786af..00000000000
--- a/libf2c/libF77/r_asin.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include <math.h>
-double
-r_asin (real * x)
-{
-  return (asin (*x));
-}
diff --git a/libf2c/libF77/r_atan.c b/libf2c/libF77/r_atan.c
deleted file mode 100644
index 33a6589fda9..00000000000
--- a/libf2c/libF77/r_atan.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include <math.h>
-double
-r_atan (real * x)
-{
-  return (atan (*x));
-}
diff --git a/libf2c/libF77/r_atn2.c b/libf2c/libF77/r_atn2.c
deleted file mode 100644
index 076d8743b1e..00000000000
--- a/libf2c/libF77/r_atn2.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include <math.h>
-double
-r_atn2 (real * x, real * y)
-{
-  return (atan2 (*x, *y));
-}
diff --git a/libf2c/libF77/r_cnjg.c b/libf2c/libF77/r_cnjg.c
deleted file mode 100644
index 5f849291b20..00000000000
--- a/libf2c/libF77/r_cnjg.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-void
-r_cnjg (complex * r, complex * z)
-{
-  real zi = z->i;
-  r->r = z->r;
-  r->i = -zi;
-}
diff --git a/libf2c/libF77/r_cos.c b/libf2c/libF77/r_cos.c
deleted file mode 100644
index ed556e81e87..00000000000
--- a/libf2c/libF77/r_cos.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include <math.h>
-double
-r_cos (real * x)
-{
-  return (cos (*x));
-}
diff --git a/libf2c/libF77/r_cosh.c b/libf2c/libF77/r_cosh.c
deleted file mode 100644
index b22e0cf0a3d..00000000000
--- a/libf2c/libF77/r_cosh.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include <math.h>
-double
-r_cosh (real * x)
-{
-  return (cosh (*x));
-}
diff --git a/libf2c/libF77/r_dim.c b/libf2c/libF77/r_dim.c
deleted file mode 100644
index 48d2fc7f767..00000000000
--- a/libf2c/libF77/r_dim.c
+++ /dev/null
@@ -1,7 +0,0 @@
-#include "f2c.h"
-
-double
-r_dim (real * a, real * b)
-{
-  return (*a > *b ? *a - *b : 0);
-}
diff --git a/libf2c/libF77/r_exp.c b/libf2c/libF77/r_exp.c
deleted file mode 100644
index 7c1ceea5895..00000000000
--- a/libf2c/libF77/r_exp.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include <math.h>
-double
-r_exp (real * x)
-{
-  return (exp (*x));
-}
diff --git a/libf2c/libF77/r_imag.c b/libf2c/libF77/r_imag.c
deleted file mode 100644
index 784abc8434b..00000000000
--- a/libf2c/libF77/r_imag.c
+++ /dev/null
@@ -1,7 +0,0 @@
-#include "f2c.h"
-
-double
-r_imag (complex * z)
-{
-  return (z->i);
-}
diff --git a/libf2c/libF77/r_int.c b/libf2c/libF77/r_int.c
deleted file mode 100644
index 3c1b28ea4c7..00000000000
--- a/libf2c/libF77/r_int.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include <math.h>
-double
-r_int (real * x)
-{
-  return ((*x > 0) ? floor (*x) : -floor (-*x));
-}
diff --git a/libf2c/libF77/r_lg10.c b/libf2c/libF77/r_lg10.c
deleted file mode 100644
index 563e73c4d8c..00000000000
--- a/libf2c/libF77/r_lg10.c
+++ /dev/null
@@ -1,11 +0,0 @@
-#include "f2c.h"
-
-#define log10e 0.43429448190325182765
-
-#undef abs
-#include <math.h>
-double
-r_lg10 (real * x)
-{
-  return (log10e * log (*x));
-}
diff --git a/libf2c/libF77/r_log.c b/libf2c/libF77/r_log.c
deleted file mode 100644
index eaaecc836f2..00000000000
--- a/libf2c/libF77/r_log.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include <math.h>
-double
-r_log (real * x)
-{
-  return (log (*x));
-}
diff --git a/libf2c/libF77/r_mod.c b/libf2c/libF77/r_mod.c
deleted file mode 100644
index 9518d66ce99..00000000000
--- a/libf2c/libF77/r_mod.c
+++ /dev/null
@@ -1,33 +0,0 @@
-#include "f2c.h"
-
-#ifdef IEEE_drem
-double drem (double, double);
-#else
-#undef abs
-#include <math.h>
-#endif
-double
-r_mod (real * x, real * y)
-{
-#ifdef IEEE_drem
-  double xa, ya, z;
-  if ((ya = *y) < 0.)
-    ya = -ya;
-  z = drem (xa = *x, ya);
-  if (xa > 0)
-    {
-      if (z < 0)
-	z += ya;
-    }
-  else if (z > 0)
-    z -= ya;
-  return z;
-#else
-  double quotient;
-  if ((quotient = (double) *x / *y) >= 0)
-    quotient = floor (quotient);
-  else
-    quotient = -floor (-quotient);
-  return (*x - (*y) * quotient);
-#endif
-}
diff --git a/libf2c/libF77/r_nint.c b/libf2c/libF77/r_nint.c
deleted file mode 100644
index f2713d588b1..00000000000
--- a/libf2c/libF77/r_nint.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include <math.h>
-double
-r_nint (real * x)
-{
-  return ((*x) >= 0 ? floor (*x + .5) : -floor (.5 - *x));
-}
diff --git a/libf2c/libF77/r_sign.c b/libf2c/libF77/r_sign.c
deleted file mode 100644
index f53c6bf9220..00000000000
--- a/libf2c/libF77/r_sign.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-double
-r_sign (real * a, real * b)
-{
-  double x;
-  x = (*a >= 0 ? *a : -*a);
-  return (*b >= 0 ? x : -x);
-}
diff --git a/libf2c/libF77/r_sin.c b/libf2c/libF77/r_sin.c
deleted file mode 100644
index 5a5ef136a8c..00000000000
--- a/libf2c/libF77/r_sin.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include <math.h>
-double
-r_sin (real * x)
-{
-  return (sin (*x));
-}
diff --git a/libf2c/libF77/r_sinh.c b/libf2c/libF77/r_sinh.c
deleted file mode 100644
index 723c7ab29c3..00000000000
--- a/libf2c/libF77/r_sinh.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include <math.h>
-double
-r_sinh (real * x)
-{
-  return (sinh (*x));
-}
diff --git a/libf2c/libF77/r_sqrt.c b/libf2c/libF77/r_sqrt.c
deleted file mode 100644
index ed832ba70e1..00000000000
--- a/libf2c/libF77/r_sqrt.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include <math.h>
-double
-r_sqrt (real * x)
-{
-  return (sqrt (*x));
-}
diff --git a/libf2c/libF77/r_tan.c b/libf2c/libF77/r_tan.c
deleted file mode 100644
index 4ef913e6522..00000000000
--- a/libf2c/libF77/r_tan.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include <math.h>
-double
-r_tan (real * x)
-{
-  return (tan (*x));
-}
diff --git a/libf2c/libF77/r_tanh.c b/libf2c/libF77/r_tanh.c
deleted file mode 100644
index 6f2552a9f97..00000000000
--- a/libf2c/libF77/r_tanh.c
+++ /dev/null
@@ -1,9 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include <math.h>
-double
-r_tanh (real * x)
-{
-  return (tanh (*x));
-}
diff --git a/libf2c/libF77/s_cat.c b/libf2c/libF77/s_cat.c
deleted file mode 100644
index 4e8da1b3fdb..00000000000
--- a/libf2c/libF77/s_cat.c
+++ /dev/null
@@ -1,70 +0,0 @@
-/* Unless compiled with -DNO_OVERWRITE, this variant of s_cat allows the
- * target of a concatenation to appear on its right-hand side (contrary
- * to the Fortran 77 Standard, but in accordance with Fortran 90).
- */
-
-#include "f2c.h"
-#ifndef NO_OVERWRITE
-#include <stdio.h>
-#undef abs
-#undef min
-#undef max
-#include <stdlib.h>
-extern char *F77_aloc (ftnlen, char *);
-#include <string.h>
-#endif /* NO_OVERWRITE */
-
-void
-s_cat (char *lp, char *rpp[], ftnint rnp[], ftnint * np, ftnlen ll)
-{
-  ftnlen i, nc;
-  char *rp;
-  ftnlen n = *np;
-#ifndef NO_OVERWRITE
-  ftnlen L, m;
-  char *lp0, *lp1;
-
-  lp0 = 0;
-  lp1 = lp;
-  L = ll;
-  i = 0;
-  while (i < n)
-    {
-      rp = rpp[i];
-      m = rnp[i++];
-      if (rp >= lp1 || rp + m <= lp)
-	{
-	  if ((L -= m) <= 0)
-	    {
-	      n = i;
-	      break;
-	    }
-	  lp1 += m;
-	  continue;
-	}
-      lp0 = lp;
-      lp = lp1 = F77_aloc (L = ll, "s_cat");
-      break;
-    }
-  lp1 = lp;
-#endif /* NO_OVERWRITE */
-  for (i = 0; i < n; ++i)
-    {
-      nc = ll;
-      if (rnp[i] < nc)
-	nc = rnp[i];
-      ll -= nc;
-      rp = rpp[i];
-      while (--nc >= 0)
-	*lp++ = *rp++;
-    }
-  while (--ll >= 0)
-    *lp++ = ' ';
-#ifndef NO_OVERWRITE
-  if (lp0)
-    {
-      memcpy (lp0, lp1, L);
-      free (lp1);
-    }
-#endif
-}
diff --git a/libf2c/libF77/s_cmp.c b/libf2c/libF77/s_cmp.c
deleted file mode 100644
index 5b43c9edb9c..00000000000
--- a/libf2c/libF77/s_cmp.c
+++ /dev/null
@@ -1,49 +0,0 @@
-#include "f2c.h"
-
-/* compare two strings */
-
-integer
-s_cmp (char *a0, char *b0, ftnlen la, ftnlen lb)
-{
-  register unsigned char *a, *aend, *b, *bend;
-  a = (unsigned char *) a0;
-  b = (unsigned char *) b0;
-  aend = a + la;
-  bend = b + lb;
-
-  if (la <= lb)
-    {
-      while (a < aend)
-	if (*a != *b)
-	  return (*a - *b);
-	else
-	  {
-	    ++a;
-	    ++b;
-	  }
-
-      while (b < bend)
-	if (*b != ' ')
-	  return (' ' - *b);
-	else
-	  ++b;
-    }
-
-  else
-    {
-      while (b < bend)
-	if (*a == *b)
-	  {
-	    ++a;
-	    ++b;
-	  }
-	else
-	  return (*a - *b);
-      while (a < aend)
-	if (*a != ' ')
-	  return (*a - ' ');
-	else
-	  ++a;
-    }
-  return (0);
-}
diff --git a/libf2c/libF77/s_copy.c b/libf2c/libF77/s_copy.c
deleted file mode 100644
index a91071eccab..00000000000
--- a/libf2c/libF77/s_copy.c
+++ /dev/null
@@ -1,50 +0,0 @@
-/* Unless compiled with -DNO_OVERWRITE, this variant of s_copy allows the
- * target of an assignment to appear on its right-hand side (contrary
- * to the Fortran 77 Standard, but in accordance with Fortran 90),
- * as in  a(2:5) = a(4:7) .
- */
-
-#include "f2c.h"
-
-/* assign strings:  a = b */
-
-void
-s_copy (register char *a, register char *b, ftnlen la, ftnlen lb)
-{
-  register char *aend, *bend;
-
-  aend = a + la;
-
-  if (la <= lb)
-#ifndef NO_OVERWRITE
-    if (a <= b || a >= b + la)
-#endif
-      while (a < aend)
-	*a++ = *b++;
-#ifndef NO_OVERWRITE
-    else
-      for (b += la; a < aend;)
-	*--aend = *--b;
-#endif
-
-  else
-    {
-      bend = b + lb;
-#ifndef NO_OVERWRITE
-      if (a <= b || a >= bend)
-#endif
-	while (b < bend)
-	  *a++ = *b++;
-#ifndef NO_OVERWRITE
-      else
-	{
-	  a += lb;
-	  while (b < bend)
-	    *--a = *--bend;
-	  a += lb;
-	}
-#endif
-      while (a < aend)
-	*a++ = ' ';
-    }
-}
diff --git a/libf2c/libF77/s_paus.c b/libf2c/libF77/s_paus.c
deleted file mode 100644
index 79691366e6f..00000000000
--- a/libf2c/libF77/s_paus.c
+++ /dev/null
@@ -1,71 +0,0 @@
-#include <stdio.h>
-#include "f2c.h"
-#define PAUSESIG 15
-
-#include "signal1.h"
-#undef abs
-#undef min
-#undef max
-#include <stdlib.h>
-extern int getpid (void), isatty (int), pause (void);
-
-extern void f_exit (void);
-
-static void
-waitpause (Sigarg)
-{
-  Use_Sigarg;
-  return;
-}
-
-static void
-s_1paus (FILE * fin)
-{
-  fprintf (stderr,
-	   "To resume execution, type go.  Other input will terminate the job.\n");
-  fflush (stderr);
-  if (getc (fin) != 'g' || getc (fin) != 'o' || getc (fin) != '\n')
-    {
-      fprintf (stderr, "STOP\n");
-#ifdef NO_ONEXIT
-      f_exit ();
-#endif
-      exit (0);
-    }
-}
-
-int
-s_paus (char *s, ftnlen n)
-{
-  fprintf (stderr, "PAUSE ");
-  if (n > 0)
-    fprintf (stderr, " %.*s", (int) n, s);
-  fprintf (stderr, " statement executed\n");
-  if (isatty (fileno (stdin)))
-    s_1paus (stdin);
-  else
-    {
-#if (defined (MSDOS) && !defined (GO32)) || defined (_WIN32)
-      FILE *fin;
-      fin = fopen ("con", "r");
-      if (!fin)
-	{
-	  fprintf (stderr, "s_paus: can't open con!\n");
-	  fflush (stderr);
-	  exit (1);
-	}
-      s_1paus (fin);
-      fclose (fin);
-#else
-      fprintf (stderr,
-	       "To resume execution, execute a   kill -%d %d   command\n",
-	       PAUSESIG, getpid ());
-      signal1 (PAUSESIG, waitpause);
-      fflush (stderr);
-      pause ();
-#endif
-    }
-  fprintf (stderr, "Execution resumes after PAUSE.\n");
-  fflush (stderr);
-  return 0;			/* NOT REACHED */
-}
diff --git a/libf2c/libF77/s_rnge.c b/libf2c/libF77/s_rnge.c
deleted file mode 100644
index 6c054c707f5..00000000000
--- a/libf2c/libF77/s_rnge.c
+++ /dev/null
@@ -1,22 +0,0 @@
-#include <stdio.h>
-#include "f2c.h"
-
-/* called when a subscript is out of range */
-
-extern void sig_die (char *, int);
-integer
-s_rnge (char *varn, ftnint offset, char *procn, ftnint line)
-{
-  register int i;
-
-  fprintf (stderr, "Subscript out of range on file line %ld, procedure ",
-	   (long) line);
-  while ((i = *procn) && i != '_' && i != ' ')
-    putc (*procn++, stderr);
-  fprintf (stderr, ".\nAttempt to access the %ld-th element of variable ",
-	   (long) offset + 1);
-  while ((i = *varn) && i != ' ')
-    putc (*varn++, stderr);
-  sig_die (".", 1);
-  return 0;			/* not reached */
-}
diff --git a/libf2c/libF77/s_stop.c b/libf2c/libF77/s_stop.c
deleted file mode 100644
index 571416be7bb..00000000000
--- a/libf2c/libF77/s_stop.c
+++ /dev/null
@@ -1,32 +0,0 @@
-#include <stdio.h>
-#include "f2c.h"
-
-#undef abs
-#undef min
-#undef max
-#include <stdlib.h>
-void f_exit (void);
-
-int
-s_stop (char *s, ftnlen n)
-{
-  int i;
-
-  if (n > 0)
-    {
-      fprintf (stderr, "STOP ");
-      for (i = 0; i < n; ++i)
-	putc (*s++, stderr);
-      fprintf (stderr, " statement executed\n");
-    }
-#ifdef NO_ONEXIT
-  f_exit ();
-#endif
-  exit (0);
-
-/* We cannot avoid (useless) compiler diagnostics here:		*/
-/* some compilers complain if there is no return statement,	*/
-/* and others complain that this one cannot be reached.		*/
-
-  return 0;			/* NOT REACHED */
-}
diff --git a/libf2c/libF77/setarg.c b/libf2c/libF77/setarg.c
deleted file mode 100644
index 49515746f48..00000000000
--- a/libf2c/libF77/setarg.c
+++ /dev/null
@@ -1,14 +0,0 @@
-/* Set up the global argc/argv info for use by getarg_, iargc_, and
-   g77's inlined intrinsic equivalents.  */
-
-#include <stdlib.h>
-
-int f__xargc;
-char **f__xargv;
-
-void
-f_setarg (int argc, char **argv)
-{
-  f__xargc = argc;
-  f__xargv = argv;
-}
diff --git a/libf2c/libF77/setsig.c b/libf2c/libF77/setsig.c
deleted file mode 100644
index 96826be94fd..00000000000
--- a/libf2c/libF77/setsig.c
+++ /dev/null
@@ -1,86 +0,0 @@
-/* Set up the signal behavior.  */
-
-#include <stdio.h>
-#include "signal1.h"
-
-#ifndef SIGIOT
-#ifdef SIGABRT
-#define SIGIOT SIGABRT
-#endif
-#endif
-
-#include <stdlib.h>
-
-extern void sig_die (char *, int);
-
-static void
-sigfdie (Sigarg)
-{
-  Use_Sigarg;
-  sig_die ("Floating Exception", 1);
-}
-
-
-static void
-sigidie (Sigarg)
-{
-  Use_Sigarg;
-  sig_die ("IOT Trap", 1);
-}
-
-#ifdef SIGQUIT
-static void
-sigqdie (Sigarg)
-{
-  Use_Sigarg;
-  sig_die ("Quit signal", 1);
-}
-#endif
-
-
-static void
-sigindie (Sigarg)
-{
-  Use_Sigarg;
-  sig_die ("Interrupt", 0);
-}
-
-static void
-sigtdie (Sigarg)
-{
-  Use_Sigarg;
-  sig_die ("Killed", 0);
-}
-
-#ifdef SIGTRAP
-static void
-sigtrdie (Sigarg)
-{
-  Use_Sigarg;
-  sig_die ("Trace trap", 1);
-}
-#endif
-
-
-void
-f_setsig ()
-{
-  signal1 (SIGFPE, sigfdie);	/* ignore underflow, enable overflow */
-#ifdef SIGIOT
-  signal1 (SIGIOT, sigidie);
-#endif
-#ifdef SIGTRAP
-  signal1 (SIGTRAP, sigtrdie);
-#endif
-#ifdef SIGQUIT
-  if (signal1 (SIGQUIT, sigqdie) == SIG_IGN)
-    signal1 (SIGQUIT, SIG_IGN);
-#endif
-  if (signal1 (SIGINT, sigindie) == SIG_IGN)
-    signal1 (SIGINT, SIG_IGN);
-  signal1 (SIGTERM, sigtdie);
-
-#ifdef pdp11
-  ldfps (01200);		/* detect overflow as an exception */
-#endif
-}
diff --git a/libf2c/libF77/sig_die.c b/libf2c/libF77/sig_die.c
deleted file mode 100644
index 622462e248b..00000000000
--- a/libf2c/libF77/sig_die.c
+++ /dev/null
@@ -1,37 +0,0 @@
-#include <stdio.h>
-#include <signal.h>
-
-#ifndef SIGIOT
-#ifdef SIGABRT
-#define SIGIOT SIGABRT
-#endif
-#endif
-
-#include <stdlib.h>
-extern void f_exit (void);
-
-void
-sig_die (register char *s, int kill)
-{
-  /* print error message, then clear buffers */
-  fprintf (stderr, "%s\n", s);
-
-  if (kill)
-    {
-      fflush (stderr);
-      f_exit ();
-      fflush (stderr);
-      /* now get a core */
-#ifdef SIGIOT
-      signal (SIGIOT, SIG_DFL);
-#endif
-      abort ();
-    }
-  else
-    {
-#ifdef NO_ONEXIT
-      f_exit ();
-#endif
-      exit (1);
-    }
-}
diff --git a/libf2c/libF77/signal1.h b/libf2c/libF77/signal1.h
deleted file mode 100644
index b559211e8e4..00000000000
--- a/libf2c/libF77/signal1.h
+++ /dev/null
@@ -1,5 +0,0 @@
-/* The g77 implementation of libf2c directly includes signal1.h0,
-   instead of copying it to signal1.h, since that seems easier to
-   cope with at this point.  */
-
-#include "signal1.h0"
diff --git a/libf2c/libF77/signal1.h0 b/libf2c/libF77/signal1.h0
deleted file mode 100644
index 0e2fcf2c1e6..00000000000
--- a/libf2c/libF77/signal1.h0
+++ /dev/null
@@ -1,26 +0,0 @@
-/* You may need to adjust the definition of signal1 to supply a */
-/* cast to the correct argument type.  This detail is system- and */
-/* compiler-dependent.   The #define below assumes signal.h declares */
-/* type SIG_PF for the signal function's second argument. */
-
-/* For some C++ compilers, "#define Sigarg_t ..." may be appropriate. */
-
-#include <signal.h>
-
-#ifndef Sigret_t
-#define Sigret_t void
-#endif
-#ifndef Sigarg_t
-#define Sigarg_t int
-#endif /*Sigarg_t */
-
-#ifdef USE_SIG_PF		/* compile with -DUSE_SIG_PF under IRIX */
-#define sig_pf SIG_PF
-#else
-typedef Sigret_t (*sig_pf) (Sigarg_t);
-#endif
-
-#define signal1(a,b) signal(a,(sig_pf)b)
-
-#define Sigarg int n
-#define Use_Sigarg n = n	/* shut up compiler warning */
diff --git a/libf2c/libF77/signal_.c b/libf2c/libF77/signal_.c
deleted file mode 100644
index f67831c3c43..00000000000
--- a/libf2c/libF77/signal_.c
+++ /dev/null
@@ -1,11 +0,0 @@
-#include "f2c.h"
-#include "signal1.h"
-
-void *
-G77_signal_0 (integer * sigp, sig_pf proc)
-{
-  int sig;
-  sig = (int) *sigp;
-
-  return (void *) signal (sig, proc);
-}
diff --git a/libf2c/libF77/system_.c b/libf2c/libF77/system_.c
deleted file mode 100644
index dbbd0bea9ec..00000000000
--- a/libf2c/libF77/system_.c
+++ /dev/null
@@ -1,28 +0,0 @@
-/* f77 interface to system routine */
-
-#include "f2c.h"
-
-#undef abs
-#undef min
-#undef max
-#include <stdlib.h>
-extern char *F77_aloc (ftnlen, char *);
-
-integer
-G77_system_0 (register char *s, ftnlen n)
-{
-  char buff0[256], *buff;
-  register char *bp, *blast;
-  integer rv;
-
-  buff = bp = n < (ftnlen) sizeof (buff0) ? buff0 : F77_aloc (n + 1, "system_");
-  blast = bp + n;
-
-  while (bp < blast && *s)
-    *bp++ = *s++;
-  *bp = 0;
-  rv = system (buff);
-  if (buff != buff0)
-    free (buff);
-  return rv;
-}
diff --git a/libf2c/libF77/z_abs.c b/libf2c/libF77/z_abs.c
deleted file mode 100644
index 2419c0eac49..00000000000
--- a/libf2c/libF77/z_abs.c
+++ /dev/null
@@ -1,8 +0,0 @@
-#include "f2c.h"
-
-double f__cabs (double, double);
-double
-z_abs (doublecomplex * z)
-{
-  return (f__cabs (z->r, z->i));
-}
diff --git a/libf2c/libF77/z_cos.c b/libf2c/libF77/z_cos.c
deleted file mode 100644
index 0f4cd71bd5a..00000000000
--- a/libf2c/libF77/z_cos.c
+++ /dev/null
@@ -1,11 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include "math.h"
-void
-z_cos (doublecomplex * r, doublecomplex * z)
-{
-  double zi = z->i, zr = z->r;
-  r->r = cos (zr) * cosh (zi);
-  r->i = -sin (zr) * sinh (zi);
-}
diff --git a/libf2c/libF77/z_div.c b/libf2c/libF77/z_div.c
deleted file mode 100644
index a5fc527bdc2..00000000000
--- a/libf2c/libF77/z_div.c
+++ /dev/null
@@ -1,41 +0,0 @@
-#include "f2c.h"
-
-extern void sig_die (char *, int);
-void
-z_div (doublecomplex * c, doublecomplex * a, doublecomplex * b)
-{
-  double ratio, den;
-  double abr, abi, cr;
-
-  if ((abr = b->r) < 0.)
-    abr = -abr;
-  if ((abi = b->i) < 0.)
-    abi = -abi;
-  if (abr <= abi)
-    {
-      if (abi == 0)
-	{
-#ifdef IEEE_COMPLEX_DIVIDE
-	  if (a->i != 0 || a->r != 0)
-	    abi = 1.;
-	  c->i = c->r = abi / abr;
-	  return;
-#else
-	  sig_die ("complex division by zero", 1);
-#endif
-	}
-      ratio = b->r / b->i;
-      den = b->i * (1 + ratio * ratio);
-      cr = (a->r * ratio + a->i) / den;
-      c->i = (a->i * ratio - a->r) / den;
-    }
-
-  else
-    {
-      ratio = b->i / b->r;
-      den = b->r * (1 + ratio * ratio);
-      cr = (a->r + a->i * ratio) / den;
-      c->i = (a->i - a->r * ratio) / den;
-    }
-  c->r = cr;
-}
diff --git a/libf2c/libF77/z_exp.c b/libf2c/libF77/z_exp.c
deleted file mode 100644
index 16f51e74a45..00000000000
--- a/libf2c/libF77/z_exp.c
+++ /dev/null
@@ -1,13 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include "math.h"
-void
-z_exp (doublecomplex * r, doublecomplex * z)
-{
-  double expx, zi = z->i;
-
-  expx = exp (z->r);
-  r->r = expx * cos (zi);
-  r->i = expx * sin (zi);
-}
diff --git a/libf2c/libF77/z_log.c b/libf2c/libF77/z_log.c
deleted file mode 100644
index f56b12ed7e3..00000000000
--- a/libf2c/libF77/z_log.c
+++ /dev/null
@@ -1,59 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include "math.h"
-extern double f__cabs (double, double);
-void
-z_log (doublecomplex * r, doublecomplex * z)
-{
-  double s, s0, t, t2, u, v;
-  double zi = z->i, zr = z->r;
-
-  r->i = atan2 (zi, zr);
-#ifdef Pre20000310
-  r->r = log (f__cabs (zr, zi));
-#else
-  if (zi < 0)
-    zi = -zi;
-  if (zr < 0)
-    zr = -zr;
-  if (zr < zi)
-    {
-      t = zi;
-      zi = zr;
-      zr = t;
-    }
-  t = zi / zr;
-  s = zr * sqrt (1 + t * t);
-  /* now s = f__cabs(zi,zr), and zr = |zr| >= |zi| = zi */
-  if ((t = s - 1) < 0)
-    t = -t;
-  if (t > .01)
-    r->r = log (s);
-  else
-    {
-
-#ifdef Comment
-
-      log (1 + x) = x - x ^ 2 / 2 + x ^ 3 / 3 - x ^ 4 / 4 + -...
-	= x (1 - x / 2 + x ^ 2 / 3 - +...)
-	[sqrt (y ^ 2 + z ^ 2) - 1] *[sqrt (y ^ 2 + z ^ 2) + 1] =
-	y ^ 2 + z ^ 2 - 1, so sqrt (y ^ 2 + z ^ 2) - 1 =
-	(y ^ 2 + z ^ 2 - 1) /[sqrt (y ^ 2 + z ^ 2) + 1]
-#endif /*Comment */
-	t = ((zr * zr - 1.) + zi * zi) / (s + 1);
-      t2 = t * t;
-      s = 1. - 0.5 * t;
-      u = v = 1;
-      do
-	{
-	  s0 = s;
-	  u *= t2;
-	  v += 2;
-	  s += u / v - t * u / (v + 1);
-	}
-      while (s > s0);
-      r->r = s * t;
-    }
-#endif
-}
diff --git a/libf2c/libF77/z_sin.c b/libf2c/libF77/z_sin.c
deleted file mode 100644
index 8cb44cf53ca..00000000000
--- a/libf2c/libF77/z_sin.c
+++ /dev/null
@@ -1,11 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include "math.h"
-void
-z_sin (doublecomplex * r, doublecomplex * z)
-{
-  double zi = z->i, zr = z->r;
-  r->r = sin (zr) * cosh (zi);
-  r->i = cos (zr) * sinh (zi);
-}
diff --git a/libf2c/libF77/z_sqrt.c b/libf2c/libF77/z_sqrt.c
deleted file mode 100644
index 954c2fa1889..00000000000
--- a/libf2c/libF77/z_sqrt.c
+++ /dev/null
@@ -1,25 +0,0 @@
-#include "f2c.h"
-
-#undef abs
-#include "math.h"
-extern double f__cabs (double, double);
-void
-z_sqrt (doublecomplex * r, doublecomplex * z)
-{
-  double mag, zi = z->i, zr = z->r;
-
-  if ((mag = f__cabs (zr, zi)) == 0.)
-    r->r = r->i = 0.;
-  else if (zr > 0)
-    {
-      r->r = sqrt (0.5 * (mag + zr));
-      r->i = zi / r->r / 2;
-    }
-  else
-    {
-      r->i = sqrt (0.5 * (mag - zr));
-      if (zi < 0)
-	r->i = -r->i;
-      r->r = zi / r->i / 2;
-    }
-}