diff --git a/academic/arpack/README b/academic/arpack/README
new file mode 100644
index 0000000000..300d65039c
--- /dev/null
+++ b/academic/arpack/README
@@ -0,0 +1,8 @@
+arpack (Fortran library for solving large scale eigenvalue problems)
+
+ARPACK is a collection of Fortran77 subroutines designed to solve
+large scale eigenvalue problems.
+
+This requires blas and lapack libraries.
+
+arpack is an optional dependency for the octave package.
diff --git a/academic/arpack/arpack-2.1-redhat.patch b/academic/arpack/arpack-2.1-redhat.patch
new file mode 100644
index 0000000000..945d0157f5
--- /dev/null
+++ b/academic/arpack/arpack-2.1-redhat.patch
@@ -0,0 +1,45 @@
+--- ARPACK/ARmake.inc~	1996-04-26 21:12:57.000000000 +0200
++++ ARPACK/ARmake.inc	2003-07-19 11:57:10.000000000 +0200
+@@ -25,7 +25,7 @@
+ # | other than your home directory.      |
+ # %--------------------------------------%
+ #
+-home = $(HOME)/ARPACK
++home = $(PWD)/.
+ #
+ #  %--------------------------------------%
+ #  | The platform identifier to suffix to |
+@@ -43,7 +43,7 @@
+ UTILdir      = $(home)/UTIL
+ SRCdir       = $(home)/SRC
+ #
+-DIRS        = $(BLASdir) $(LAPACKdir) $(UTILdir) $(SRCdir)
++#DIRS        = $(BLASdir) $(LAPACKdir) $(UTILdir) $(SRCdir)
+ #
+ # %-------------------------------------------------------------------%
+ # | Comment out the previous line and uncomment the following         |
+@@ -51,13 +51,13 @@
+ # | NOTE: ARPACK assumes the use of LAPACK version 2 codes.           |
+ # %-------------------------------------------------------------------%
+ #
+-#DIRS         = $(UTILdir) $(SRCdir)
++DIRS         = $(UTILdir) $(SRCdir)
+ #
+ # %---------------------------------------------------%
+ # | The name of the libraries to be created/linked to |
+ # %---------------------------------------------------%
+ #
+-ARPACKLIB  = $(home)/libarpack_$(PLAT).a
++ARPACKLIB  = $(home)/libarpack.a
+ LAPACKLIB = 
+ BLASLIB = 
+ #
+@@ -112,7 +112,7 @@
+ LN      = ln
+ LNFLAGS = -s
+ 
+-MAKE    = /bin/make
++MAKE    = make
+ 
+ RM      = rm
+ RMFLAGS = -f
diff --git a/academic/arpack/arpack-etime.patch b/academic/arpack/arpack-etime.patch
new file mode 100644
index 0000000000..5a6e9862ec
--- /dev/null
+++ b/academic/arpack/arpack-etime.patch
@@ -0,0 +1,11 @@
+diff -up ARPACK/UTIL/second.f.etime ARPACK/UTIL/second.f
+--- ARPACK/UTIL/second.f.etime	2008-09-24 22:40:53.000000000 +0200
++++ ARPACK/UTIL/second.f	2008-09-24 22:42:03.000000000 +0200
+@@ -21,7 +21,6 @@
+ *     ..
+ *     .. External Functions ..
+       REAL               ETIME
+-      EXTERNAL           ETIME
+ *     ..
+ *     .. Executable Statements ..
+ *
diff --git a/academic/arpack/arpack-second-bug.patch b/academic/arpack/arpack-second-bug.patch
new file mode 100644
index 0000000000..99998b26d2
--- /dev/null
+++ b/academic/arpack/arpack-second-bug.patch
@@ -0,0 +1,1829 @@
+diff -up ARPACK/UTIL/second.f.sb ARPACK/UTIL/second.f
+--- ARPACK/UTIL/second.f.sb	1995-10-02 19:56:44.000000000 +0100
++++ ARPACK/UTIL/second.f	2007-10-24 02:50:57.000000000 +0200
+@@ -1,4 +1,4 @@
+-      SUBROUTINE SECOND( T )
++      SUBROUTINE secnd2( T )
+ *
+       REAL       T
+ *
+diff -up ARPACK/SRC/znaupd.f.sb ARPACK/SRC/znaupd.f
+--- ARPACK/SRC/znaupd.f.sb	2002-08-15 07:50:57.000000000 +0200
++++ ARPACK/SRC/znaupd.f	2007-10-24 02:50:57.000000000 +0200
+@@ -451,7 +451,7 @@ c        | & message level for debugging
+ c        %-------------------------------%
+ c
+          call zstatn
+-         call second (t0)
++         call secnd2 (t0)
+          msglvl = mcaupd
+ c
+ c        %----------------%
+@@ -611,7 +611,7 @@ c
+      &               '_naupd: Associated Ritz estimates')
+       end if
+ c
+-      call second (t1)
++      call secnd2 (t1)
+       tcaupd = t1 - t0
+ c
+       if (msglvl .gt. 0) then
+diff -up ARPACK/SRC/dgetv0.f.sb ARPACK/SRC/dgetv0.f
+--- ARPACK/SRC/dgetv0.f.sb	1999-04-08 00:42:47.000000000 +0200
++++ ARPACK/SRC/dgetv0.f	2007-10-24 02:50:57.000000000 +0200
+@@ -214,7 +214,7 @@ c        | Initialize timing statistics 
+ c        | & message level for debugging |
+ c        %-------------------------------%
+ c
+-         call second (t0)
++         call secnd2 (t0)
+          msglvl = mgetv0
+ c 
+          ierr   = 0
+@@ -241,7 +241,7 @@ c        | Force the starting vector int
+ c        | the generalized problem when B is possibly (singular).   |
+ c        %----------------------------------------------------------%
+ c
+-         call second (t2)
++         call secnd2 (t2)
+          if (bmat .eq. 'G') then
+             nopx = nopx + 1
+             ipntr(1) = 1
+@@ -265,7 +265,7 @@ c
+       if (orth)  go to 40
+ c 
+       if (bmat .eq. 'G') then
+-         call second (t3)
++         call secnd2 (t3)
+          tmvopx = tmvopx + (t3 - t2)
+       end if
+ c 
+@@ -274,7 +274,7 @@ c     | Starting vector is now in the ra
+ c     | Compute B-norm of starting vector.                   |
+ c     %------------------------------------------------------%
+ c
+-      call second (t2)
++      call secnd2 (t2)
+       first = .TRUE.
+       if (bmat .eq. 'G') then
+          nbx = nbx + 1
+@@ -290,7 +290,7 @@ c 
+    20 continue
+ c
+       if (bmat .eq. 'G') then
+-         call second (t3)
++         call secnd2 (t3)
+          tmvbx = tmvbx + (t3 - t2)
+       end if
+ c 
+@@ -333,7 +333,7 @@ c     %---------------------------------
+ c     | Compute the B-norm of the orthogonalized starting vector |
+ c     %----------------------------------------------------------%
+ c
+-      call second (t2)
++      call secnd2 (t2)
+       if (bmat .eq. 'G') then
+          nbx = nbx + 1
+          call dcopy (n, resid, 1, workd(n+1), 1)
+@@ -348,7 +348,7 @@ c 
+    40 continue
+ c
+       if (bmat .eq. 'G') then
+-         call second (t3)
++         call secnd2 (t3)
+          tmvbx = tmvbx + (t3 - t2)
+       end if
+ c 
+@@ -406,7 +406,7 @@ c
+       end if
+       ido = 99
+ c 
+-      call second (t1)
++      call secnd2 (t1)
+       tgetv0 = tgetv0 + (t1 - t0)
+ c 
+  9000 continue
+diff -up ARPACK/SRC/ssgets.f.sb ARPACK/SRC/ssgets.f
+--- ARPACK/SRC/ssgets.f.sb	1996-04-20 18:27:30.000000000 +0200
++++ ARPACK/SRC/ssgets.f	2007-10-24 02:50:57.000000000 +0200
+@@ -148,7 +148,7 @@ c     | Initialize timing statistics  |
+ c     | & message level for debugging |
+ c     %-------------------------------%
+ c
+-      call second (t0)
++      call secnd2 (t0)
+       msglvl = msgets
+ c 
+       if (which .eq. 'BE') then
+@@ -198,7 +198,7 @@ c     
+          call scopy (np, ritz, 1, shifts, 1)
+       end if
+ c 
+-      call second (t1)
++      call secnd2 (t1)
+       tsgets = tsgets + (t1 - t0)
+ c
+       if (msglvl .gt. 0) then
+diff -up ARPACK/SRC/snaup2.f.sb ARPACK/SRC/snaup2.f
+--- ARPACK/SRC/snaup2.f.sb	2002-08-23 07:31:49.000000000 +0200
++++ ARPACK/SRC/snaup2.f	2007-10-24 02:50:57.000000000 +0200
+@@ -258,7 +258,7 @@ c     %-----------------------%
+ c
+       if (ido .eq. 0) then
+ c 
+-         call second (t0)
++         call secnd2 (t0)
+ c 
+          msglvl = mnaup2
+ c 
+@@ -758,7 +758,7 @@ c        | the first step of the next ca
+ c        %---------------------------------------------%
+ c
+          cnorm = .true.
+-         call second (t2)
++         call secnd2 (t2)
+          if (bmat .eq. 'G') then
+             nbx = nbx + 1
+             call scopy (n, resid, 1, workd(n+1), 1)
+@@ -783,7 +783,7 @@ c        | WORKD(1:N) := B*RESID        
+ c        %----------------------------------%
+ c
+          if (bmat .eq. 'G') then
+-            call second (t3)
++            call secnd2 (t3)
+             tmvbx = tmvbx + (t3 - t2)
+          end if
+ c 
+@@ -822,7 +822,7 @@ c     %------------%
+ c     | Error Exit |
+ c     %------------%
+ c
+-      call second (t1)
++      call secnd2 (t1)
+       tnaup2 = t1 - t0
+ c     
+  9000 continue
+diff -up ARPACK/SRC/dsgets.f.sb ARPACK/SRC/dsgets.f
+--- ARPACK/SRC/dsgets.f.sb	1996-04-20 18:27:14.000000000 +0200
++++ ARPACK/SRC/dsgets.f	2007-10-24 02:50:57.000000000 +0200
+@@ -148,7 +148,7 @@ c     | Initialize timing statistics  |
+ c     | & message level for debugging |
+ c     %-------------------------------%
+ c
+-      call second (t0)
++      call secnd2 (t0)
+       msglvl = msgets
+ c 
+       if (which .eq. 'BE') then
+@@ -198,7 +198,7 @@ c     
+          call dcopy (np, ritz, 1, shifts, 1)
+       end if
+ c 
+-      call second (t1)
++      call secnd2 (t1)
+       tsgets = tsgets + (t1 - t0)
+ c
+       if (msglvl .gt. 0) then
+diff -up ARPACK/SRC/dsconv.f.sb ARPACK/SRC/dsconv.f
+--- ARPACK/SRC/dsconv.f.sb	1996-04-20 18:27:11.000000000 +0200
++++ ARPACK/SRC/dsconv.f	2007-10-24 02:50:57.000000000 +0200
+@@ -106,7 +106,7 @@ c     %-----------------------%
+ c     | Executable Statements |
+ c     %-----------------------%
+ c
+-      call second (t0)
++      call secnd2 (t0)
+ c
+       eps23 = dlamch('Epsilon-Machine') 
+       eps23 = eps23**(2.0D+0 / 3.0D+0)
+@@ -126,7 +126,7 @@ c
+ c
+    10 continue
+ c 
+-      call second (t1)
++      call secnd2 (t1)
+       tsconv = tsconv + (t1 - t0)
+ c 
+       return
+diff -up ARPACK/SRC/dneigh.f.sb ARPACK/SRC/dneigh.f
+--- ARPACK/SRC/dneigh.f.sb	1996-04-20 18:27:06.000000000 +0200
++++ ARPACK/SRC/dneigh.f	2007-10-24 02:50:57.000000000 +0200
+@@ -170,7 +170,7 @@ c     | Initialize timing statistics  |
+ c     | & message level for debugging |
+ c     %-------------------------------%
+ c
+-      call second (t0)
++      call secnd2 (t0)
+       msglvl = mneigh
+ c 
+       if (msglvl .gt. 2) then
+@@ -301,7 +301,7 @@ c
+      &              '_neigh: Ritz estimates for the eigenvalues of H')
+       end if
+ c
+-      call second (t1)
++      call secnd2 (t1)
+       tneigh = tneigh + (t1 - t0)
+ c
+  9000 continue
+diff -up ARPACK/SRC/dsaitr.f.sb ARPACK/SRC/dsaitr.f
+--- ARPACK/SRC/dsaitr.f.sb	1996-08-28 16:21:43.000000000 +0200
++++ ARPACK/SRC/dsaitr.f	2007-10-24 02:50:57.000000000 +0200
+@@ -300,7 +300,7 @@ c        | Initialize timing statistics 
+ c        | & message level for debugging |
+ c        %-------------------------------%
+ c
+-         call second (t0)
++         call secnd2 (t0)
+          msglvl = msaitr
+ c 
+ c        %------------------------------%
+@@ -420,7 +420,7 @@ c              | which spans OP and exit
+ c              %------------------------------------------------%
+ c
+                info = j - 1
+-               call second (t1)
++               call secnd2 (t1)
+                tsaitr = tsaitr + (t1 - t0)
+                ido = 99
+                go to 9000
+@@ -460,7 +460,7 @@ c        %------------------------------
+ c
+          step3 = .true.
+          nopx  = nopx + 1
+-         call second (t2)
++         call secnd2 (t2)
+          call dcopy (n, v(1,j), 1, workd(ivj), 1)
+          ipntr(1) = ivj
+          ipntr(2) = irj
+@@ -479,7 +479,7 @@ c        | Back from reverse communicati
+ c        | WORKD(IRJ:IRJ+N-1) := OP*v_{j}.   |
+ c        %-----------------------------------%
+ c
+-         call second (t3)
++         call secnd2 (t3)
+          tmvopx = tmvopx + (t3 - t2)
+ c 
+          step3 = .false.
+@@ -500,7 +500,7 @@ c        | assumed to have A*v_{j}.     
+ c        %-------------------------------------------%
+ c
+          if (mode .eq. 2) go to 65
+-         call second (t2)
++         call secnd2 (t2)
+          if (bmat .eq. 'G') then
+             nbx = nbx + 1
+             step4 = .true.
+@@ -524,7 +524,7 @@ c        | WORKD(IPJ:IPJ+N-1) := B*OP*v_
+ c        %-----------------------------------%
+ c
+          if (bmat .eq. 'G') then
+-            call second (t3)
++            call secnd2 (t3)
+             tmvbx = tmvbx + (t3 - t2)
+          end if 
+ c
+@@ -592,12 +592,12 @@ c
+          else
+             h(j,1) = rnorm
+          end if
+-         call second (t4)
++         call secnd2 (t4)
+ c 
+          orth1 = .true.
+          iter  = 0
+ c 
+-         call second (t2)
++         call secnd2 (t2)
+          if (bmat .eq. 'G') then
+             nbx = nbx + 1
+             call dcopy (n, resid, 1, workd(irj), 1)
+@@ -621,7 +621,7 @@ c        | WORKD(IPJ:IPJ+N-1) := B*r_{j}
+ c        %---------------------------------------------------%
+ c
+          if (bmat .eq. 'G') then
+-            call second (t3)
++            call secnd2 (t3)
+             tmvbx = tmvbx + (t3 - t2)
+          end if
+ c 
+@@ -695,7 +695,7 @@ c
+          h(j,2) = h(j,2) + workd(irj + j - 1)
+ c 
+          orth2 = .true.
+-         call second (t2)
++         call secnd2 (t2)
+          if (bmat .eq. 'G') then
+             nbx = nbx + 1
+             call dcopy (n, resid, 1, workd(irj), 1)
+@@ -719,7 +719,7 @@ c        | Back from reverse communicati
+ c        %---------------------------------------------------%
+ c
+          if (bmat .eq. 'G') then
+-            call second (t3)
++            call secnd2 (t3)
+             tmvbx = tmvbx + (t3 - t2)
+          end if
+ c
+@@ -791,7 +791,7 @@ c 
+          rstart = .false.
+          orth2  = .false.
+ c 
+-         call second (t5)
++         call secnd2 (t5)
+          titref = titref + (t5 - t4)
+ c 
+ c        %----------------------------------------------------------%
+@@ -815,7 +815,7 @@ c        %------------------------------
+ c
+          j = j + 1
+          if (j .gt. k+np) then
+-            call second (t1)
++            call secnd2 (t1)
+             tsaitr = tsaitr + (t1 - t0)
+             ido = 99
+ c
+diff -up ARPACK/SRC/dnaup2.f.sb ARPACK/SRC/dnaup2.f
+--- ARPACK/SRC/dnaup2.f.sb	2002-08-23 07:31:50.000000000 +0200
++++ ARPACK/SRC/dnaup2.f	2007-10-24 02:50:57.000000000 +0200
+@@ -258,7 +258,7 @@ c     %-----------------------%
+ c
+       if (ido .eq. 0) then
+ c 
+-         call second (t0)
++         call secnd2 (t0)
+ c 
+          msglvl = mnaup2
+ c 
+@@ -758,7 +758,7 @@ c        | the first step of the next ca
+ c        %---------------------------------------------%
+ c
+          cnorm = .true.
+-         call second (t2)
++         call secnd2 (t2)
+          if (bmat .eq. 'G') then
+             nbx = nbx + 1
+             call dcopy (n, resid, 1, workd(n+1), 1)
+@@ -783,7 +783,7 @@ c        | WORKD(1:N) := B*RESID        
+ c        %----------------------------------%
+ c
+          if (bmat .eq. 'G') then
+-            call second (t3)
++            call secnd2 (t3)
+             tmvbx = tmvbx + (t3 - t2)
+          end if
+ c 
+@@ -822,7 +822,7 @@ c     %------------%
+ c     | Error Exit |
+ c     %------------%
+ c
+-      call second (t1)
++      call secnd2 (t1)
+       tnaup2 = t1 - t0
+ c     
+  9000 continue
+diff -up ARPACK/SRC/cneigh.f.sb ARPACK/SRC/cneigh.f
+--- ARPACK/SRC/cneigh.f.sb	1996-04-20 19:15:53.000000000 +0200
++++ ARPACK/SRC/cneigh.f	2007-10-24 02:50:57.000000000 +0200
+@@ -171,7 +171,7 @@ c     | Initialize timing statistics  |
+ c     | & message level for debugging |
+ c     %-------------------------------%
+ c
+-      call second (t0)
++      call secnd2 (t0)
+       msglvl = mceigh
+ c 
+       if (msglvl .gt. 2) then
+@@ -244,7 +244,7 @@ c
+      &              '_neigh: Ritz estimates for the eigenvalues of H')
+       end if
+ c
+-      call second(t1)
++      call secnd2(t1)
+       tceigh = tceigh + (t1 - t0)
+ c
+  9000 continue
+diff -up ARPACK/SRC/zngets.f.sb ARPACK/SRC/zngets.f
+--- ARPACK/SRC/zngets.f.sb	1996-04-20 19:16:02.000000000 +0200
++++ ARPACK/SRC/zngets.f	2007-10-24 02:50:57.000000000 +0200
+@@ -137,7 +137,7 @@ c     | Initialize timing statistics  |
+ c     | & message level for debugging |
+ c     %-------------------------------%
+ c 
+-      call second (t0)
++      call secnd2 (t0)
+       msglvl = mcgets
+ c 
+       call zsortc (which, .true., kev+np, ritz, bounds)
+@@ -157,7 +157,7 @@ c     
+ c
+       end if
+ c     
+-      call second (t1)
++      call secnd2 (t1)
+       tcgets = tcgets + (t1 - t0)
+ c
+       if (msglvl .gt. 0) then
+diff -up ARPACK/SRC/sseigt.f.sb ARPACK/SRC/sseigt.f
+--- ARPACK/SRC/sseigt.f.sb	1996-08-27 07:29:04.000000000 +0200
++++ ARPACK/SRC/sseigt.f	2007-10-24 02:50:57.000000000 +0200
+@@ -138,7 +138,7 @@ c     | Initialize timing statistics  |
+ c     | & message level for debugging |
+ c     %-------------------------------% 
+ c
+-      call second (t0)
++      call secnd2 (t0)
+       msglvl = mseigt
+ c
+       if (msglvl .gt. 0) then
+@@ -168,7 +168,7 @@ c
+          bounds(k) = rnorm*abs(bounds(k))
+    30 continue
+ c 
+-      call second (t1)
++      call secnd2 (t1)
+       tseigt = tseigt + (t1 - t0)
+ c
+  9000 continue
+diff -up ARPACK/SRC/dnapps.f.sb ARPACK/SRC/dnapps.f
+--- ARPACK/SRC/dnapps.f.sb	1998-05-20 16:58:56.000000000 +0200
++++ ARPACK/SRC/dnapps.f	2007-10-24 02:50:57.000000000 +0200
+@@ -237,7 +237,7 @@ c     | Initialize timing statistics  |
+ c     | & message level for debugging |
+ c     %-------------------------------%
+ c
+-      call second (t0)
++      call secnd2 (t0)
+       msglvl = mnapps
+       kplusp = kev + np 
+ c 
+@@ -635,7 +635,7 @@ c
+       end if
+ c 
+  9000 continue
+-      call second (t1)
++      call secnd2 (t1)
+       tnapps = tnapps + (t1 - t0)
+ c 
+       return
+diff -up ARPACK/SRC/dsapps.f.sb ARPACK/SRC/dsapps.f
+--- ARPACK/SRC/dsapps.f.sb	1998-05-20 16:58:59.000000000 +0200
++++ ARPACK/SRC/dsapps.f	2007-10-24 02:50:57.000000000 +0200
+@@ -213,7 +213,7 @@ c     | Initialize timing statistics  |
+ c     | & message level for debugging |
+ c     %-------------------------------%
+ c
+-      call second (t0)
++      call secnd2 (t0)
+       msglvl = msapps
+ c 
+       kplusp = kev + np 
+@@ -503,7 +503,7 @@ c
+          end if
+       end if
+ c
+-      call second (t1)
++      call secnd2 (t1)
+       tsapps = tsapps + (t1 - t0)
+ c 
+  9000 continue 
+diff -up ARPACK/SRC/ssaup2.f.sb ARPACK/SRC/ssaup2.f
+--- ARPACK/SRC/ssaup2.f.sb	1998-05-20 16:58:59.000000000 +0200
++++ ARPACK/SRC/ssaup2.f	2007-10-24 02:50:57.000000000 +0200
+@@ -262,7 +262,7 @@ c        | Initialize timing statistics 
+ c        | & message level for debugging |
+ c        %-------------------------------%
+ c
+-         call second (t0)
++         call secnd2 (t0)
+          msglvl = msaup2
+ c
+ c        %---------------------------------%
+@@ -770,7 +770,7 @@ c        | the first step of the next ca
+ c        %---------------------------------------------%
+ c
+          cnorm = .true.
+-         call second (t2)
++         call secnd2 (t2)
+          if (bmat .eq. 'G') then
+             nbx = nbx + 1
+             call scopy (n, resid, 1, workd(n+1), 1)
+@@ -795,7 +795,7 @@ c        | WORKD(1:N) := B*RESID        
+ c        %----------------------------------%
+ c
+          if (bmat .eq. 'G') then
+-            call second (t3)
++            call secnd2 (t3)
+             tmvbx = tmvbx + (t3 - t2)
+          end if
+ c 
+@@ -837,7 +837,7 @@ c     %------------%
+ c     | Error exit |
+ c     %------------%
+ c
+-      call second (t1)
++      call secnd2 (t1)
+       tsaup2 = t1 - t0
+ c 
+  9000 continue
+diff -up ARPACK/SRC/ssaitr.f.sb ARPACK/SRC/ssaitr.f
+--- ARPACK/SRC/ssaitr.f.sb	1996-08-28 16:21:43.000000000 +0200
++++ ARPACK/SRC/ssaitr.f	2007-10-24 02:50:57.000000000 +0200
+@@ -300,7 +300,7 @@ c        | Initialize timing statistics 
+ c        | & message level for debugging |
+ c        %-------------------------------%
+ c
+-         call second (t0)
++         call secnd2 (t0)
+          msglvl = msaitr
+ c 
+ c        %------------------------------%
+@@ -420,7 +420,7 @@ c              | which spans OP and exit
+ c              %------------------------------------------------%
+ c
+                info = j - 1
+-               call second (t1)
++               call secnd2 (t1)
+                tsaitr = tsaitr + (t1 - t0)
+                ido = 99
+                go to 9000
+@@ -460,7 +460,7 @@ c        %------------------------------
+ c
+          step3 = .true.
+          nopx  = nopx + 1
+-         call second (t2)
++         call secnd2 (t2)
+          call scopy (n, v(1,j), 1, workd(ivj), 1)
+          ipntr(1) = ivj
+          ipntr(2) = irj
+@@ -479,7 +479,7 @@ c        | Back from reverse communicati
+ c        | WORKD(IRJ:IRJ+N-1) := OP*v_{j}.   |
+ c        %-----------------------------------%
+ c
+-         call second (t3)
++         call secnd2 (t3)
+          tmvopx = tmvopx + (t3 - t2)
+ c 
+          step3 = .false.
+@@ -500,7 +500,7 @@ c        | assumed to have A*v_{j}.     
+ c        %-------------------------------------------%
+ c
+          if (mode .eq. 2) go to 65
+-         call second (t2)
++         call secnd2 (t2)
+          if (bmat .eq. 'G') then
+             nbx = nbx + 1
+             step4 = .true.
+@@ -524,7 +524,7 @@ c        | WORKD(IPJ:IPJ+N-1) := B*OP*v_
+ c        %-----------------------------------%
+ c
+          if (bmat .eq. 'G') then
+-            call second (t3)
++            call secnd2 (t3)
+             tmvbx = tmvbx + (t3 - t2)
+          end if 
+ c
+@@ -592,12 +592,12 @@ c
+          else
+             h(j,1) = rnorm
+          end if
+-         call second (t4)
++         call secnd2 (t4)
+ c 
+          orth1 = .true.
+          iter  = 0
+ c 
+-         call second (t2)
++         call secnd2 (t2)
+          if (bmat .eq. 'G') then
+             nbx = nbx + 1
+             call scopy (n, resid, 1, workd(irj), 1)
+@@ -621,7 +621,7 @@ c        | WORKD(IPJ:IPJ+N-1) := B*r_{j}
+ c        %---------------------------------------------------%
+ c
+          if (bmat .eq. 'G') then
+-            call second (t3)
++            call secnd2 (t3)
+             tmvbx = tmvbx + (t3 - t2)
+          end if
+ c 
+@@ -695,7 +695,7 @@ c
+          h(j,2) = h(j,2) + workd(irj + j - 1)
+ c 
+          orth2 = .true.
+-         call second (t2)
++         call secnd2 (t2)
+          if (bmat .eq. 'G') then
+             nbx = nbx + 1
+             call scopy (n, resid, 1, workd(irj), 1)
+@@ -719,7 +719,7 @@ c        | Back from reverse communicati
+ c        %---------------------------------------------------%
+ c
+          if (bmat .eq. 'G') then
+-            call second (t3)
++            call secnd2 (t3)
+             tmvbx = tmvbx + (t3 - t2)
+          end if
+ c
+@@ -791,7 +791,7 @@ c 
+          rstart = .false.
+          orth2  = .false.
+ c 
+-         call second (t5)
++         call secnd2 (t5)
+          titref = titref + (t5 - t4)
+ c 
+ c        %----------------------------------------------------------%
+@@ -815,7 +815,7 @@ c        %------------------------------
+ c
+          j = j + 1
+          if (j .gt. k+np) then
+-            call second (t1)
++            call secnd2 (t1)
+             tsaitr = tsaitr + (t1 - t0)
+             ido = 99
+ c
+diff -up ARPACK/SRC/snapps.f.sb ARPACK/SRC/snapps.f
+--- ARPACK/SRC/snapps.f.sb	1998-05-20 16:58:55.000000000 +0200
++++ ARPACK/SRC/snapps.f	2007-10-24 02:50:57.000000000 +0200
+@@ -237,7 +237,7 @@ c     | Initialize timing statistics  |
+ c     | & message level for debugging |
+ c     %-------------------------------%
+ c
+-      call second (t0)
++      call secnd2 (t0)
+       msglvl = mnapps
+       kplusp = kev + np 
+ c 
+@@ -635,7 +635,7 @@ c
+       end if
+ c 
+  9000 continue
+-      call second (t1)
++      call secnd2 (t1)
+       tnapps = tnapps + (t1 - t0)
+ c 
+       return
+diff -up ARPACK/SRC/dnaitr.f.sb ARPACK/SRC/dnaitr.f
+--- ARPACK/SRC/dnaitr.f.sb	1996-08-27 18:09:00.000000000 +0200
++++ ARPACK/SRC/dnaitr.f	2007-10-24 02:50:57.000000000 +0200
+@@ -319,7 +319,7 @@ c        | Initialize timing statistics 
+ c        | & message level for debugging |
+ c        %-------------------------------%
+ c
+-         call second (t0)
++         call secnd2 (t0)
+          msglvl = mnaitr
+ c 
+ c        %------------------------------%
+@@ -430,7 +430,7 @@ c              | which spans OP and exit
+ c              %------------------------------------------------%
+ c
+                info = j - 1
+-               call second (t1)
++               call secnd2 (t1)
+                tnaitr = tnaitr + (t1 - t0)
+                ido = 99
+                go to 9000
+@@ -470,7 +470,7 @@ c        %------------------------------
+ c
+          step3 = .true.
+          nopx  = nopx + 1
+-         call second (t2)
++         call secnd2 (t2)
+          call dcopy (n, v(1,j), 1, workd(ivj), 1)
+          ipntr(1) = ivj
+          ipntr(2) = irj
+@@ -490,7 +490,7 @@ c        | WORKD(IRJ:IRJ+N-1) := OP*v_{j
+ c        | if step3 = .true.                |
+ c        %----------------------------------%
+ c
+-         call second (t3)
++         call secnd2 (t3)
+          tmvopx = tmvopx + (t3 - t2)
+  
+          step3 = .false.
+@@ -506,7 +506,7 @@ c        | STEP 4:  Finish extending the
+ c        |          factorization to length j.   |
+ c        %---------------------------------------%
+ c
+-         call second (t2)
++         call secnd2 (t2)
+          if (bmat .eq. 'G') then
+             nbx = nbx + 1
+             step4 = .true.
+@@ -531,7 +531,7 @@ c        | if step4 = .true.            
+ c        %----------------------------------%
+ c
+          if (bmat .eq. 'G') then
+-            call second (t3)
++            call secnd2 (t3)
+             tmvbx = tmvbx + (t3 - t2)
+          end if
+ c 
+@@ -576,11 +576,11 @@ c
+ c
+          if (j .gt. 1) h(j,j-1) = betaj
+ c
+-         call second (t4)
++         call secnd2 (t4)
+ c 
+          orth1 = .true.
+ c
+-         call second (t2)
++         call secnd2 (t2)
+          if (bmat .eq. 'G') then
+             nbx = nbx + 1
+             call dcopy (n, resid, 1, workd(irj), 1)
+@@ -604,7 +604,7 @@ c        | WORKD(IPJ:IPJ+N-1) := B*r_{j}
+ c        %---------------------------------------------------%
+ c
+          if (bmat .eq. 'G') then
+-            call second (t3)
++            call secnd2 (t3)
+             tmvbx = tmvbx + (t3 - t2)
+          end if
+ c 
+@@ -681,7 +681,7 @@ c
+          call daxpy (j, one, workd(irj), 1, h(1,j), 1)
+ c 
+          orth2 = .true.
+-         call second (t2)
++         call secnd2 (t2)
+          if (bmat .eq. 'G') then
+             nbx = nbx + 1
+             call dcopy (n, resid, 1, workd(irj), 1)
+@@ -705,7 +705,7 @@ c        | Back from reverse communicati
+ c        %---------------------------------------------------%
+ c
+          if (bmat .eq. 'G') then
+-            call second (t3)
++            call secnd2 (t3)
+             tmvbx = tmvbx + (t3 - t2)
+          end if
+ c
+@@ -783,7 +783,7 @@ c 
+          rstart = .false.
+          orth2  = .false.
+ c 
+-         call second (t5)
++         call secnd2 (t5)
+          titref = titref + (t5 - t4)
+ c 
+ c        %------------------------------------%
+@@ -792,7 +792,7 @@ c        %------------------------------
+ c
+          j = j + 1
+          if (j .gt. k+np) then
+-            call second (t1)
++            call secnd2 (t1)
+             tnaitr = tnaitr + (t1 - t0)
+             ido = 99
+             do 110 i = max(1,k), k+np-1
+diff -up ARPACK/SRC/sneigh.f.sb ARPACK/SRC/sneigh.f
+--- ARPACK/SRC/sneigh.f.sb	1996-04-20 18:27:22.000000000 +0200
++++ ARPACK/SRC/sneigh.f	2007-10-24 02:50:57.000000000 +0200
+@@ -170,7 +170,7 @@ c     | Initialize timing statistics  |
+ c     | & message level for debugging |
+ c     %-------------------------------%
+ c
+-      call second (t0)
++      call secnd2 (t0)
+       msglvl = mneigh
+ c 
+       if (msglvl .gt. 2) then
+@@ -301,7 +301,7 @@ c
+      &              '_neigh: Ritz estimates for the eigenvalues of H')
+       end if
+ c
+-      call second (t1)
++      call secnd2 (t1)
+       tneigh = tneigh + (t1 - t0)
+ c
+  9000 continue
+diff -up ARPACK/SRC/dsaup2.f.sb ARPACK/SRC/dsaup2.f
+--- ARPACK/SRC/dsaup2.f.sb	1998-05-20 16:59:00.000000000 +0200
++++ ARPACK/SRC/dsaup2.f	2007-10-24 02:50:57.000000000 +0200
+@@ -262,7 +262,7 @@ c        | Initialize timing statistics 
+ c        | & message level for debugging |
+ c        %-------------------------------%
+ c
+-         call second (t0)
++         call secnd2 (t0)
+          msglvl = msaup2
+ c
+ c        %---------------------------------%
+@@ -770,7 +770,7 @@ c        | the first step of the next ca
+ c        %---------------------------------------------%
+ c
+          cnorm = .true.
+-         call second (t2)
++         call secnd2 (t2)
+          if (bmat .eq. 'G') then
+             nbx = nbx + 1
+             call dcopy (n, resid, 1, workd(n+1), 1)
+@@ -795,7 +795,7 @@ c        | WORKD(1:N) := B*RESID        
+ c        %----------------------------------%
+ c
+          if (bmat .eq. 'G') then
+-            call second (t3)
++            call secnd2 (t3)
+             tmvbx = tmvbx + (t3 - t2)
+          end if
+ c 
+@@ -837,7 +837,7 @@ c     %------------%
+ c     | Error exit |
+ c     %------------%
+ c
+-      call second (t1)
++      call secnd2 (t1)
+       tsaup2 = t1 - t0
+ c 
+  9000 continue
+diff -up ARPACK/SRC/snaupd.f.sb ARPACK/SRC/snaupd.f
+--- ARPACK/SRC/snaupd.f.sb	2002-08-23 07:29:34.000000000 +0200
++++ ARPACK/SRC/snaupd.f	2007-10-24 02:50:57.000000000 +0200
+@@ -476,7 +476,7 @@ c        | & message level for debugging
+ c        %-------------------------------%
+ c
+          call sstatn
+-         call second (t0)
++         call secnd2 (t0)
+          msglvl = mnaupd
+ c
+ c        %----------------%
+@@ -640,7 +640,7 @@ c
+      &               '_naupd: Associated Ritz estimates')
+       end if
+ c
+-      call second (t1)
++      call secnd2 (t1)
+       tnaupd = t1 - t0
+ c
+       if (msglvl .gt. 0) then
+diff -up ARPACK/SRC/znapps.f.sb ARPACK/SRC/znapps.f
+--- ARPACK/SRC/znapps.f.sb	1998-05-20 16:45:03.000000000 +0200
++++ ARPACK/SRC/znapps.f	2007-10-24 02:50:57.000000000 +0200
+@@ -240,7 +240,7 @@ c     | Initialize timing statistics  |
+ c     | & message level for debugging |
+ c     %-------------------------------%
+ c
+-      call second (t0)
++      call secnd2 (t0)
+       msglvl = mcapps
+ c 
+       kplusp = kev + np 
+@@ -495,7 +495,7 @@ c
+       end if
+ c
+  9000 continue
+-      call second (t1)
++      call secnd2 (t1)
+       tcapps = tcapps + (t1 - t0)
+ c 
+       return
+diff -up ARPACK/SRC/znaitr.f.sb ARPACK/SRC/znaitr.f
+--- ARPACK/SRC/znaitr.f.sb	1996-08-27 18:12:57.000000000 +0200
++++ ARPACK/SRC/znaitr.f	2007-10-24 02:50:57.000000000 +0200
+@@ -326,7 +326,7 @@ c        | Initialize timing statistics 
+ c        | & message level for debugging |
+ c        %-------------------------------%
+ c
+-         call second (t0)
++         call secnd2 (t0)
+          msglvl = mcaitr
+ c 
+ c        %------------------------------%
+@@ -437,7 +437,7 @@ c              | which spans OP and exit
+ c              %------------------------------------------------%
+ c
+                info = j - 1
+-               call second (t1)
++               call secnd2 (t1)
+                tcaitr = tcaitr + (t1 - t0)
+                ido = 99
+                go to 9000
+@@ -477,7 +477,7 @@ c        %------------------------------
+ c
+          step3 = .true.
+          nopx  = nopx + 1
+-         call second (t2)
++         call secnd2 (t2)
+          call zcopy (n, v(1,j), 1, workd(ivj), 1)
+          ipntr(1) = ivj
+          ipntr(2) = irj
+@@ -497,7 +497,7 @@ c        | WORKD(IRJ:IRJ+N-1) := OP*v_{j
+ c        | if step3 = .true.                |
+ c        %----------------------------------%
+ c
+-         call second (t3)
++         call secnd2 (t3)
+          tmvopx = tmvopx + (t3 - t2)
+  
+          step3 = .false.
+@@ -513,7 +513,7 @@ c        | STEP 4:  Finish extending the
+ c        |          factorization to length j.   |
+ c        %---------------------------------------%
+ c
+-         call second (t2)
++         call secnd2 (t2)
+          if (bmat .eq. 'G') then
+             nbx = nbx + 1
+             step4 = .true.
+@@ -538,7 +538,7 @@ c        | if step4 = .true.            
+ c        %----------------------------------%
+ c
+          if (bmat .eq. 'G') then
+-            call second (t3)
++            call secnd2 (t3)
+             tmvbx = tmvbx + (t3 - t2)
+          end if
+ c 
+@@ -583,11 +583,11 @@ c
+ c
+          if (j .gt. 1) h(j,j-1) = dcmplx(betaj, rzero)
+ c
+-         call second (t4)
++         call secnd2 (t4)
+ c 
+          orth1 = .true.
+ c 
+-         call second (t2)
++         call secnd2 (t2)
+          if (bmat .eq. 'G') then
+             nbx = nbx + 1
+             call zcopy (n, resid, 1, workd(irj), 1)
+@@ -611,7 +611,7 @@ c        | WORKD(IPJ:IPJ+N-1) := B*r_{j}
+ c        %---------------------------------------------------%
+ c
+          if (bmat .eq. 'G') then
+-            call second (t3)
++            call secnd2 (t3)
+             tmvbx = tmvbx + (t3 - t2)
+          end if
+ c 
+@@ -689,7 +689,7 @@ c
+          call zaxpy (j, one, workd(irj), 1, h(1,j), 1)
+ c 
+          orth2 = .true.
+-         call second (t2)
++         call secnd2 (t2)
+          if (bmat .eq. 'G') then
+             nbx = nbx + 1
+             call zcopy (n, resid, 1, workd(irj), 1)
+@@ -713,7 +713,7 @@ c        | Back from reverse communicati
+ c        %---------------------------------------------------%
+ c
+          if (bmat .eq. 'G') then
+-            call second (t3)
++            call secnd2 (t3)
+             tmvbx = tmvbx + (t3 - t2)
+          end if 
+ c
+@@ -791,7 +791,7 @@ c 
+          rstart = .false.
+          orth2  = .false.
+ c 
+-         call second (t5)
++         call secnd2 (t5)
+          titref = titref + (t5 - t4)
+ c 
+ c        %------------------------------------%
+@@ -800,7 +800,7 @@ c        %------------------------------
+ c
+          j = j + 1
+          if (j .gt. k+np) then
+-            call second (t1)
++            call secnd2 (t1)
+             tcaitr = tcaitr + (t1 - t0)
+             ido = 99
+             do 110 i = max(1,k), k+np-1
+diff -up ARPACK/SRC/zgetv0.f.sb ARPACK/SRC/zgetv0.f
+--- ARPACK/SRC/zgetv0.f.sb	1999-04-08 00:46:45.000000000 +0200
++++ ARPACK/SRC/zgetv0.f	2007-10-24 02:50:57.000000000 +0200
+@@ -211,7 +211,7 @@ c        | Initialize timing statistics 
+ c        | & message level for debugging |
+ c        %-------------------------------%
+ c
+-         call second (t0)
++         call secnd2 (t0)
+          msglvl = mgetv0
+ c 
+          ierr   = 0
+@@ -238,7 +238,7 @@ c        | Force the starting vector int
+ c        | the generalized problem when B is possibly (singular).   |
+ c        %----------------------------------------------------------%
+ c
+-         call second (t2)
++         call secnd2 (t2)
+          if (bmat .eq. 'G') then
+             nopx = nopx + 1
+             ipntr(1) = 1
+@@ -261,7 +261,7 @@ c     %---------------------------------
+ c
+       if (orth)  go to 40
+ c 
+-      call second (t3)
++      call secnd2 (t3)
+       tmvopx = tmvopx + (t3 - t2)
+ c 
+ c     %------------------------------------------------------%
+@@ -269,7 +269,7 @@ c     | Starting vector is now in the ra
+ c     | Compute B-norm of starting vector.                   |
+ c     %------------------------------------------------------%
+ c
+-      call second (t2)
++      call secnd2 (t2)
+       first = .TRUE.
+       if (bmat .eq. 'G') then
+          nbx = nbx + 1
+@@ -285,7 +285,7 @@ c 
+    20 continue
+ c
+       if (bmat .eq. 'G') then
+-         call second (t3)
++         call secnd2 (t3)
+          tmvbx = tmvbx + (t3 - t2)
+       end if
+ c 
+@@ -328,7 +328,7 @@ c     %---------------------------------
+ c     | Compute the B-norm of the orthogonalized starting vector |
+ c     %----------------------------------------------------------%
+ c
+-      call second (t2)
++      call secnd2 (t2)
+       if (bmat .eq. 'G') then
+          nbx = nbx + 1
+          call zcopy (n, resid, 1, workd(n+1), 1)
+@@ -343,7 +343,7 @@ c 
+    40 continue
+ c
+       if (bmat .eq. 'G') then
+-         call second (t3)
++         call secnd2 (t3)
+          tmvbx = tmvbx + (t3 - t2)
+       end if
+ c 
+@@ -401,7 +401,7 @@ c
+       end if
+       ido = 99
+ c 
+-      call second (t1)
++      call secnd2 (t1)
+       tgetv0 = tgetv0 + (t1 - t0)
+ c 
+  9000 continue
+diff -up ARPACK/SRC/zneigh.f.sb ARPACK/SRC/zneigh.f
+--- ARPACK/SRC/zneigh.f.sb	1996-04-20 19:16:01.000000000 +0200
++++ ARPACK/SRC/zneigh.f	2007-10-24 02:50:57.000000000 +0200
+@@ -171,7 +171,7 @@ c     | Initialize timing statistics  |
+ c     | & message level for debugging |
+ c     %-------------------------------%
+ c
+-      call second (t0)
++      call secnd2 (t0)
+       msglvl = mceigh
+ c 
+       if (msglvl .gt. 2) then
+@@ -244,7 +244,7 @@ c
+      &              '_neigh: Ritz estimates for the eigenvalues of H')
+       end if
+ c
+-      call second(t1)
++      call secnd2(t1)
+       tceigh = tceigh + (t1 - t0)
+ c
+  9000 continue
+diff -up ARPACK/SRC/cnaupd.f.sb ARPACK/SRC/cnaupd.f
+--- ARPACK/SRC/cnaupd.f.sb	2002-08-15 07:50:57.000000000 +0200
++++ ARPACK/SRC/cnaupd.f	2007-10-24 02:50:57.000000000 +0200
+@@ -451,7 +451,7 @@ c        | & message level for debugging
+ c        %-------------------------------%
+ c
+          call cstatn
+-         call second (t0)
++         call secnd2 (t0)
+          msglvl = mcaupd
+ c
+ c        %----------------%
+@@ -611,7 +611,7 @@ c
+      &               '_naupd: Associated Ritz estimates')
+       end if
+ c
+-      call second (t1)
++      call secnd2 (t1)
+       tcaupd = t1 - t0
+ c
+       if (msglvl .gt. 0) then
+diff -up ARPACK/SRC/cgetv0.f.sb ARPACK/SRC/cgetv0.f
+--- ARPACK/SRC/cgetv0.f.sb	1999-04-08 00:46:44.000000000 +0200
++++ ARPACK/SRC/cgetv0.f	2007-10-24 02:50:57.000000000 +0200
+@@ -211,7 +211,7 @@ c        | Initialize timing statistics 
+ c        | & message level for debugging |
+ c        %-------------------------------%
+ c
+-         call second (t0)
++         call secnd2 (t0)
+          msglvl = mgetv0
+ c 
+          ierr   = 0
+@@ -238,7 +238,7 @@ c        | Force the starting vector int
+ c        | the generalized problem when B is possibly (singular).   |
+ c        %----------------------------------------------------------%
+ c
+-         call second (t2)
++         call secnd2 (t2)
+          if (bmat .eq. 'G') then
+             nopx = nopx + 1
+             ipntr(1) = 1
+@@ -261,7 +261,7 @@ c     %---------------------------------
+ c
+       if (orth)  go to 40
+ c 
+-      call second (t3)
++      call secnd2 (t3)
+       tmvopx = tmvopx + (t3 - t2)
+ c 
+ c     %------------------------------------------------------%
+@@ -269,7 +269,7 @@ c     | Starting vector is now in the ra
+ c     | Compute B-norm of starting vector.                   |
+ c     %------------------------------------------------------%
+ c
+-      call second (t2)
++      call secnd2 (t2)
+       first = .TRUE.
+       if (bmat .eq. 'G') then
+          nbx = nbx + 1
+@@ -285,7 +285,7 @@ c 
+    20 continue
+ c
+       if (bmat .eq. 'G') then
+-         call second (t3)
++         call secnd2 (t3)
+          tmvbx = tmvbx + (t3 - t2)
+       end if
+ c 
+@@ -328,7 +328,7 @@ c     %---------------------------------
+ c     | Compute the B-norm of the orthogonalized starting vector |
+ c     %----------------------------------------------------------%
+ c
+-      call second (t2)
++      call secnd2 (t2)
+       if (bmat .eq. 'G') then
+          nbx = nbx + 1
+          call ccopy (n, resid, 1, workd(n+1), 1)
+@@ -343,7 +343,7 @@ c 
+    40 continue
+ c
+       if (bmat .eq. 'G') then
+-         call second (t3)
++         call secnd2 (t3)
+          tmvbx = tmvbx + (t3 - t2)
+       end if
+ c 
+@@ -401,7 +401,7 @@ c
+       end if
+       ido = 99
+ c 
+-      call second (t1)
++      call secnd2 (t1)
+       tgetv0 = tgetv0 + (t1 - t0)
+ c 
+  9000 continue
+diff -up ARPACK/SRC/sgetv0.f.sb ARPACK/SRC/sgetv0.f
+--- ARPACK/SRC/sgetv0.f.sb	1999-04-08 00:42:46.000000000 +0200
++++ ARPACK/SRC/sgetv0.f	2007-10-24 02:50:57.000000000 +0200
+@@ -214,7 +214,7 @@ c        | Initialize timing statistics 
+ c        | & message level for debugging |
+ c        %-------------------------------%
+ c
+-         call second (t0)
++         call secnd2 (t0)
+          msglvl = mgetv0
+ c 
+          ierr   = 0
+@@ -241,7 +241,7 @@ c        | Force the starting vector int
+ c        | the generalized problem when B is possibly (singular).   |
+ c        %----------------------------------------------------------%
+ c
+-         call second (t2)
++         call secnd2 (t2)
+          if (bmat .eq. 'G') then
+             nopx = nopx + 1
+             ipntr(1) = 1
+@@ -265,7 +265,7 @@ c
+       if (orth)  go to 40
+ c 
+       if (bmat .eq. 'G') then
+-         call second (t3)
++         call secnd2 (t3)
+          tmvopx = tmvopx + (t3 - t2)
+       end if
+ c 
+@@ -274,7 +274,7 @@ c     | Starting vector is now in the ra
+ c     | Compute B-norm of starting vector.                   |
+ c     %------------------------------------------------------%
+ c
+-      call second (t2)
++      call secnd2 (t2)
+       first = .TRUE.
+       if (bmat .eq. 'G') then
+          nbx = nbx + 1
+@@ -290,7 +290,7 @@ c 
+    20 continue
+ c
+       if (bmat .eq. 'G') then
+-         call second (t3)
++         call secnd2 (t3)
+          tmvbx = tmvbx + (t3 - t2)
+       end if
+ c 
+@@ -333,7 +333,7 @@ c     %---------------------------------
+ c     | Compute the B-norm of the orthogonalized starting vector |
+ c     %----------------------------------------------------------%
+ c
+-      call second (t2)
++      call secnd2 (t2)
+       if (bmat .eq. 'G') then
+          nbx = nbx + 1
+          call scopy (n, resid, 1, workd(n+1), 1)
+@@ -348,7 +348,7 @@ c 
+    40 continue
+ c
+       if (bmat .eq. 'G') then
+-         call second (t3)
++         call secnd2 (t3)
+          tmvbx = tmvbx + (t3 - t2)
+       end if
+ c 
+@@ -406,7 +406,7 @@ c
+       end if
+       ido = 99
+ c 
+-      call second (t1)
++      call secnd2 (t1)
+       tgetv0 = tgetv0 + (t1 - t0)
+ c 
+  9000 continue
+diff -up ARPACK/SRC/dnaupd.f.sb ARPACK/SRC/dnaupd.f
+--- ARPACK/SRC/dnaupd.f.sb	2002-08-23 07:29:34.000000000 +0200
++++ ARPACK/SRC/dnaupd.f	2007-10-24 02:50:57.000000000 +0200
+@@ -476,7 +476,7 @@ c        | & message level for debugging
+ c        %-------------------------------%
+ c
+          call dstatn
+-         call second (t0)
++         call secnd2 (t0)
+          msglvl = mnaupd
+ c
+ c        %----------------%
+@@ -640,7 +640,7 @@ c
+      &               '_naupd: Associated Ritz estimates')
+       end if
+ c
+-      call second (t1)
++      call secnd2 (t1)
+       tnaupd = t1 - t0
+ c
+       if (msglvl .gt. 0) then
+diff -up ARPACK/SRC/cnaitr.f.sb ARPACK/SRC/cnaitr.f
+--- ARPACK/SRC/cnaitr.f.sb	1996-08-27 18:12:56.000000000 +0200
++++ ARPACK/SRC/cnaitr.f	2007-10-24 02:50:57.000000000 +0200
+@@ -326,7 +326,7 @@ c        | Initialize timing statistics 
+ c        | & message level for debugging |
+ c        %-------------------------------%
+ c
+-         call second (t0)
++         call secnd2 (t0)
+          msglvl = mcaitr
+ c 
+ c        %------------------------------%
+@@ -437,7 +437,7 @@ c              | which spans OP and exit
+ c              %------------------------------------------------%
+ c
+                info = j - 1
+-               call second (t1)
++               call secnd2 (t1)
+                tcaitr = tcaitr + (t1 - t0)
+                ido = 99
+                go to 9000
+@@ -477,7 +477,7 @@ c        %------------------------------
+ c
+          step3 = .true.
+          nopx  = nopx + 1
+-         call second (t2)
++         call secnd2 (t2)
+          call ccopy (n, v(1,j), 1, workd(ivj), 1)
+          ipntr(1) = ivj
+          ipntr(2) = irj
+@@ -497,7 +497,7 @@ c        | WORKD(IRJ:IRJ+N-1) := OP*v_{j
+ c        | if step3 = .true.                |
+ c        %----------------------------------%
+ c
+-         call second (t3)
++         call secnd2 (t3)
+          tmvopx = tmvopx + (t3 - t2)
+  
+          step3 = .false.
+@@ -513,7 +513,7 @@ c        | STEP 4:  Finish extending the
+ c        |          factorization to length j.   |
+ c        %---------------------------------------%
+ c
+-         call second (t2)
++         call secnd2 (t2)
+          if (bmat .eq. 'G') then
+             nbx = nbx + 1
+             step4 = .true.
+@@ -538,7 +538,7 @@ c        | if step4 = .true.            
+ c        %----------------------------------%
+ c
+          if (bmat .eq. 'G') then
+-            call second (t3)
++            call secnd2 (t3)
+             tmvbx = tmvbx + (t3 - t2)
+          end if
+ c 
+@@ -583,11 +583,11 @@ c
+ c
+          if (j .gt. 1) h(j,j-1) = cmplx(betaj, rzero)
+ c
+-         call second (t4)
++         call secnd2 (t4)
+ c 
+          orth1 = .true.
+ c 
+-         call second (t2)
++         call secnd2 (t2)
+          if (bmat .eq. 'G') then
+             nbx = nbx + 1
+             call ccopy (n, resid, 1, workd(irj), 1)
+@@ -611,7 +611,7 @@ c        | WORKD(IPJ:IPJ+N-1) := B*r_{j}
+ c        %---------------------------------------------------%
+ c
+          if (bmat .eq. 'G') then
+-            call second (t3)
++            call secnd2 (t3)
+             tmvbx = tmvbx + (t3 - t2)
+          end if
+ c 
+@@ -689,7 +689,7 @@ c
+          call caxpy (j, one, workd(irj), 1, h(1,j), 1)
+ c 
+          orth2 = .true.
+-         call second (t2)
++         call secnd2 (t2)
+          if (bmat .eq. 'G') then
+             nbx = nbx + 1
+             call ccopy (n, resid, 1, workd(irj), 1)
+@@ -713,7 +713,7 @@ c        | Back from reverse communicati
+ c        %---------------------------------------------------%
+ c
+          if (bmat .eq. 'G') then
+-            call second (t3)
++            call secnd2 (t3)
+             tmvbx = tmvbx + (t3 - t2)
+          end if 
+ c
+@@ -791,7 +791,7 @@ c 
+          rstart = .false.
+          orth2  = .false.
+ c 
+-         call second (t5)
++         call secnd2 (t5)
+          titref = titref + (t5 - t4)
+ c 
+ c        %------------------------------------%
+@@ -800,7 +800,7 @@ c        %------------------------------
+ c
+          j = j + 1
+          if (j .gt. k+np) then
+-            call second (t1)
++            call secnd2 (t1)
+             tcaitr = tcaitr + (t1 - t0)
+             ido = 99
+             do 110 i = max(1,k), k+np-1
+diff -up ARPACK/SRC/dsaupd.f.sb ARPACK/SRC/dsaupd.f
+--- ARPACK/SRC/dsaupd.f.sb	2001-04-11 00:52:40.000000000 +0200
++++ ARPACK/SRC/dsaupd.f	2007-10-24 02:50:57.000000000 +0200
+@@ -478,7 +478,7 @@ c        | & message level for debugging
+ c        %-------------------------------%
+ c
+          call dstats 
+-         call second (t0)
++         call secnd2 (t0)
+          msglvl = msaupd
+ c
+          ierr   = 0
+@@ -638,7 +638,7 @@ c
+      &               '_saupd: corresponding error bounds')
+       end if 
+ c
+-      call second (t1)
++      call secnd2 (t1)
+       tsaupd = t1 - t0
+ c 
+       if (msglvl .gt. 0) then
+diff -up ARPACK/SRC/znaup2.f.sb ARPACK/SRC/znaup2.f
+--- ARPACK/SRC/znaup2.f.sb	2000-06-01 22:28:53.000000000 +0200
++++ ARPACK/SRC/znaup2.f	2007-10-24 02:50:57.000000000 +0200
+@@ -264,7 +264,7 @@ c     %-----------------------%
+ c
+       if (ido .eq. 0) then
+ c 
+-         call second (t0)
++         call secnd2 (t0)
+ c 
+          msglvl = mcaup2
+ c 
+@@ -724,7 +724,7 @@ c        | the first step of the next ca
+ c        %---------------------------------------------%
+ c
+          cnorm = .true.
+-         call second (t2)
++         call secnd2 (t2)
+          if (bmat .eq. 'G') then
+             nbx = nbx + 1
+             call zcopy  (n, resid, 1, workd(n+1), 1)
+@@ -749,7 +749,7 @@ c        | WORKD(1:N) := B*RESID        
+ c        %----------------------------------%
+ c
+          if (bmat .eq. 'G') then
+-            call second (t3)
++            call secnd2 (t3)
+             tmvbx = tmvbx + (t3 - t2)
+          end if
+ c 
+@@ -788,7 +788,7 @@ c     %------------%
+ c     | Error Exit |
+ c     %------------%
+ c
+-      call second (t1)
++      call secnd2 (t1)
+       tcaup2 = t1 - t0
+ c     
+  9000 continue
+diff -up ARPACK/SRC/snaitr.f.sb ARPACK/SRC/snaitr.f
+--- ARPACK/SRC/snaitr.f.sb	1996-08-27 18:09:01.000000000 +0200
++++ ARPACK/SRC/snaitr.f	2007-10-24 02:50:57.000000000 +0200
+@@ -319,7 +319,7 @@ c        | Initialize timing statistics 
+ c        | & message level for debugging |
+ c        %-------------------------------%
+ c
+-         call second (t0)
++         call secnd2 (t0)
+          msglvl = mnaitr
+ c 
+ c        %------------------------------%
+@@ -430,7 +430,7 @@ c              | which spans OP and exit
+ c              %------------------------------------------------%
+ c
+                info = j - 1
+-               call second (t1)
++               call secnd2 (t1)
+                tnaitr = tnaitr + (t1 - t0)
+                ido = 99
+                go to 9000
+@@ -470,7 +470,7 @@ c        %------------------------------
+ c
+          step3 = .true.
+          nopx  = nopx + 1
+-         call second (t2)
++         call secnd2 (t2)
+          call scopy (n, v(1,j), 1, workd(ivj), 1)
+          ipntr(1) = ivj
+          ipntr(2) = irj
+@@ -490,7 +490,7 @@ c        | WORKD(IRJ:IRJ+N-1) := OP*v_{j
+ c        | if step3 = .true.                |
+ c        %----------------------------------%
+ c
+-         call second (t3)
++         call secnd2 (t3)
+          tmvopx = tmvopx + (t3 - t2)
+  
+          step3 = .false.
+@@ -506,7 +506,7 @@ c        | STEP 4:  Finish extending the
+ c        |          factorization to length j.   |
+ c        %---------------------------------------%
+ c
+-         call second (t2)
++         call secnd2 (t2)
+          if (bmat .eq. 'G') then
+             nbx = nbx + 1
+             step4 = .true.
+@@ -531,7 +531,7 @@ c        | if step4 = .true.            
+ c        %----------------------------------%
+ c
+          if (bmat .eq. 'G') then
+-            call second (t3)
++            call secnd2 (t3)
+             tmvbx = tmvbx + (t3 - t2)
+          end if
+ c 
+@@ -576,11 +576,11 @@ c
+ c
+          if (j .gt. 1) h(j,j-1) = betaj
+ c
+-         call second (t4)
++         call secnd2 (t4)
+ c 
+          orth1 = .true.
+ c
+-         call second (t2)
++         call secnd2 (t2)
+          if (bmat .eq. 'G') then
+             nbx = nbx + 1
+             call scopy (n, resid, 1, workd(irj), 1)
+@@ -604,7 +604,7 @@ c        | WORKD(IPJ:IPJ+N-1) := B*r_{j}
+ c        %---------------------------------------------------%
+ c
+          if (bmat .eq. 'G') then
+-            call second (t3)
++            call secnd2 (t3)
+             tmvbx = tmvbx + (t3 - t2)
+          end if
+ c 
+@@ -681,7 +681,7 @@ c
+          call saxpy (j, one, workd(irj), 1, h(1,j), 1)
+ c 
+          orth2 = .true.
+-         call second (t2)
++         call secnd2 (t2)
+          if (bmat .eq. 'G') then
+             nbx = nbx + 1
+             call scopy (n, resid, 1, workd(irj), 1)
+@@ -705,7 +705,7 @@ c        | Back from reverse communicati
+ c        %---------------------------------------------------%
+ c
+          if (bmat .eq. 'G') then
+-            call second (t3)
++            call secnd2 (t3)
+             tmvbx = tmvbx + (t3 - t2)
+          end if
+ c
+@@ -783,7 +783,7 @@ c 
+          rstart = .false.
+          orth2  = .false.
+ c 
+-         call second (t5)
++         call secnd2 (t5)
+          titref = titref + (t5 - t4)
+ c 
+ c        %------------------------------------%
+@@ -792,7 +792,7 @@ c        %------------------------------
+ c
+          j = j + 1
+          if (j .gt. k+np) then
+-            call second (t1)
++            call secnd2 (t1)
+             tnaitr = tnaitr + (t1 - t0)
+             ido = 99
+             do 110 i = max(1,k), k+np-1
+diff -up ARPACK/SRC/cnapps.f.sb ARPACK/SRC/cnapps.f
+--- ARPACK/SRC/cnapps.f.sb	1998-05-20 16:45:03.000000000 +0200
++++ ARPACK/SRC/cnapps.f	2007-10-24 02:50:57.000000000 +0200
+@@ -240,7 +240,7 @@ c     | Initialize timing statistics  |
+ c     | & message level for debugging |
+ c     %-------------------------------%
+ c
+-      call second (t0)
++      call secnd2 (t0)
+       msglvl = mcapps
+ c 
+       kplusp = kev + np 
+@@ -495,7 +495,7 @@ c
+       end if
+ c
+  9000 continue
+-      call second (t1)
++      call secnd2 (t1)
+       tcapps = tcapps + (t1 - t0)
+ c 
+       return
+diff -up ARPACK/SRC/dnconv.f.sb ARPACK/SRC/dnconv.f
+--- ARPACK/SRC/dnconv.f.sb	1996-04-20 18:27:05.000000000 +0200
++++ ARPACK/SRC/dnconv.f	2007-10-24 02:50:57.000000000 +0200
+@@ -119,7 +119,7 @@ c     |                                 
+ c     | for some appropriate choice of norm.                        |
+ c     %-------------------------------------------------------------%
+ c
+-      call second (t0)
++      call secnd2 (t0)
+ c
+ c     %---------------------------------%
+ c     | Get machine dependent constant. |
+@@ -134,7 +134,7 @@ c
+          if (bounds(i) .le. tol*temp)   nconv = nconv + 1
+    20 continue
+ c 
+-      call second (t1)
++      call secnd2 (t1)
+       tnconv = tnconv + (t1 - t0)
+ c 
+       return
+diff -up ARPACK/SRC/dngets.f.sb ARPACK/SRC/dngets.f
+--- ARPACK/SRC/dngets.f.sb	1996-04-20 18:27:07.000000000 +0200
++++ ARPACK/SRC/dngets.f	2007-10-24 02:50:57.000000000 +0200
+@@ -152,7 +152,7 @@ c     | Initialize timing statistics  |
+ c     | & message level for debugging |
+ c     %-------------------------------%
+ c 
+-      call second (t0)
++      call secnd2 (t0)
+       msglvl = mngets
+ c 
+ c     %----------------------------------------------------%
+@@ -208,7 +208,7 @@ c     
+          call dsortc ( 'SR', .true., np, bounds, ritzr, ritzi )
+       end if
+ c     
+-      call second (t1)
++      call secnd2 (t1)
+       tngets = tngets + (t1 - t0)
+ c
+       if (msglvl .gt. 0) then
+diff -up ARPACK/SRC/ssconv.f.sb ARPACK/SRC/ssconv.f
+--- ARPACK/SRC/ssconv.f.sb	1996-04-20 18:27:27.000000000 +0200
++++ ARPACK/SRC/ssconv.f	2007-10-24 02:50:57.000000000 +0200
+@@ -106,7 +106,7 @@ c     %-----------------------%
+ c     | Executable Statements |
+ c     %-----------------------%
+ c
+-      call second (t0)
++      call secnd2 (t0)
+ c
+       eps23 = slamch('Epsilon-Machine') 
+       eps23 = eps23**(2.0E+0 / 3.0E+0)
+@@ -126,7 +126,7 @@ c
+ c
+    10 continue
+ c 
+-      call second (t1)
++      call secnd2 (t1)
+       tsconv = tsconv + (t1 - t0)
+ c 
+       return
+diff -up ARPACK/SRC/ssapps.f.sb ARPACK/SRC/ssapps.f
+--- ARPACK/SRC/ssapps.f.sb	1998-05-20 16:58:58.000000000 +0200
++++ ARPACK/SRC/ssapps.f	2007-10-24 02:50:57.000000000 +0200
+@@ -213,7 +213,7 @@ c     | Initialize timing statistics  |
+ c     | & message level for debugging |
+ c     %-------------------------------%
+ c
+-      call second (t0)
++      call secnd2 (t0)
+       msglvl = msapps
+ c 
+       kplusp = kev + np 
+@@ -503,7 +503,7 @@ c
+          end if
+       end if
+ c
+-      call second (t1)
++      call secnd2 (t1)
+       tsapps = tsapps + (t1 - t0)
+ c 
+  9000 continue 
+diff -up ARPACK/SRC/ssaupd.f.sb ARPACK/SRC/ssaupd.f
+--- ARPACK/SRC/ssaupd.f.sb	2001-04-11 00:52:40.000000000 +0200
++++ ARPACK/SRC/ssaupd.f	2007-10-24 02:50:57.000000000 +0200
+@@ -478,7 +478,7 @@ c        | & message level for debugging
+ c        %-------------------------------%
+ c
+          call sstats
+-         call second (t0)
++         call secnd2 (t0)
+          msglvl = msaupd
+ c
+          ierr   = 0
+@@ -638,7 +638,7 @@ c
+      &               '_saupd: corresponding error bounds')
+       end if 
+ c
+-      call second (t1)
++      call secnd2 (t1)
+       tsaupd = t1 - t0
+ c 
+       if (msglvl .gt. 0) then
+diff -up ARPACK/SRC/snconv.f.sb ARPACK/SRC/snconv.f
+--- ARPACK/SRC/snconv.f.sb	1996-04-20 18:27:21.000000000 +0200
++++ ARPACK/SRC/snconv.f	2007-10-24 02:50:57.000000000 +0200
+@@ -119,7 +119,7 @@ c     |                                 
+ c     | for some appropriate choice of norm.                        |
+ c     %-------------------------------------------------------------%
+ c
+-      call second (t0)
++      call secnd2 (t0)
+ c
+ c     %---------------------------------%
+ c     | Get machine dependent constant. |
+@@ -134,7 +134,7 @@ c
+          if (bounds(i) .le. tol*temp)   nconv = nconv + 1
+    20 continue
+ c 
+-      call second (t1)
++      call secnd2 (t1)
+       tnconv = tnconv + (t1 - t0)
+ c 
+       return
+diff -up ARPACK/SRC/cnaup2.f.sb ARPACK/SRC/cnaup2.f
+--- ARPACK/SRC/cnaup2.f.sb	2000-06-01 22:28:53.000000000 +0200
++++ ARPACK/SRC/cnaup2.f	2007-10-24 02:50:57.000000000 +0200
+@@ -264,7 +264,7 @@ c     %-----------------------%
+ c
+       if (ido .eq. 0) then
+ c 
+-         call second (t0)
++         call secnd2 (t0)
+ c 
+          msglvl = mcaup2
+ c 
+@@ -724,7 +724,7 @@ c        | the first step of the next ca
+ c        %---------------------------------------------%
+ c
+          cnorm = .true.
+-         call second (t2)
++         call secnd2 (t2)
+          if (bmat .eq. 'G') then
+             nbx = nbx + 1
+             call ccopy (n, resid, 1, workd(n+1), 1)
+@@ -749,7 +749,7 @@ c        | WORKD(1:N) := B*RESID        
+ c        %----------------------------------%
+ c
+          if (bmat .eq. 'G') then
+-            call second (t3)
++            call secnd2 (t3)
+             tmvbx = tmvbx + (t3 - t2)
+          end if
+ c 
+@@ -788,7 +788,7 @@ c     %------------%
+ c     | Error Exit |
+ c     %------------%
+ c
+-      call second (t1)
++      call secnd2 (t1)
+       tcaup2 = t1 - t0
+ c     
+  9000 continue
+diff -up ARPACK/SRC/dseigt.f.sb ARPACK/SRC/dseigt.f
+--- ARPACK/SRC/dseigt.f.sb	1996-08-27 07:29:04.000000000 +0200
++++ ARPACK/SRC/dseigt.f	2007-10-24 02:50:57.000000000 +0200
+@@ -138,7 +138,7 @@ c     | Initialize timing statistics  |
+ c     | & message level for debugging |
+ c     %-------------------------------% 
+ c
+-      call second (t0)
++      call secnd2 (t0)
+       msglvl = mseigt
+ c
+       if (msglvl .gt. 0) then
+@@ -168,7 +168,7 @@ c
+          bounds(k) = rnorm*abs(bounds(k))
+    30 continue
+ c 
+-      call second (t1)
++      call secnd2 (t1)
+       tseigt = tseigt + (t1 - t0)
+ c
+  9000 continue
+diff -up ARPACK/SRC/sngets.f.sb ARPACK/SRC/sngets.f
+--- ARPACK/SRC/sngets.f.sb	1996-04-20 18:27:24.000000000 +0200
++++ ARPACK/SRC/sngets.f	2007-10-24 02:50:57.000000000 +0200
+@@ -152,7 +152,7 @@ c     | Initialize timing statistics  |
+ c     | & message level for debugging |
+ c     %-------------------------------%
+ c 
+-      call second (t0)
++      call secnd2 (t0)
+       msglvl = mngets
+ c 
+ c     %----------------------------------------------------%
+@@ -208,7 +208,7 @@ c     
+          call ssortc ( 'SR', .true., np, bounds, ritzr, ritzi )
+       end if
+ c     
+-      call second (t1)
++      call secnd2 (t1)
+       tngets = tngets + (t1 - t0)
+ c
+       if (msglvl .gt. 0) then
+diff -up ARPACK/SRC/cngets.f.sb ARPACK/SRC/cngets.f
+--- ARPACK/SRC/cngets.f.sb	1996-04-20 19:15:56.000000000 +0200
++++ ARPACK/SRC/cngets.f	2007-10-24 02:50:57.000000000 +0200
+@@ -137,7 +137,7 @@ c     | Initialize timing statistics  |
+ c     | & message level for debugging |
+ c     %-------------------------------%
+ c 
+-      call second (t0)
++      call secnd2 (t0)
+       msglvl = mcgets
+ c 
+       call csortc (which, .true., kev+np, ritz, bounds)
+@@ -157,7 +157,7 @@ c     
+ c
+       end if
+ c     
+-      call second (t1)
++      call secnd2 (t1)
+       tcgets = tcgets + (t1 - t0)
+ c
+       if (msglvl .gt. 0) then
diff --git a/academic/arpack/arpack.SlackBuild b/academic/arpack/arpack.SlackBuild
new file mode 100644
index 0000000000..c976fdc587
--- /dev/null
+++ b/academic/arpack/arpack.SlackBuild
@@ -0,0 +1,101 @@
+#!/bin/sh
+
+# Slackware build script for arpack
+
+# Written by Pablo Santamaria (pablosantamaria@gmail.com)
+
+PRGNAM=arpack
+VERSION=${VERSION:-2.1}
+BUILD=${BUILD:-1}
+TAG=${TAG:-_SBo}
+
+# Automatically determine the architecture we're building on:
+if [ -z "$ARCH" ]; then
+  case "$( uname -m )" in
+    i?86) ARCH=i486 ;;
+    arm*) ARCH=arm ;;
+    # Unless $ARCH is already set, use uname -m for all other archs:
+       *) ARCH=$( uname -m ) ;;
+  esac
+fi
+
+CWD=$(pwd)
+TMP=${TMP:-/tmp/SBo}
+PKG=$TMP/package-$PRGNAM
+OUTPUT=${OUTPUT:-/tmp}
+
+if [ "$ARCH" = "i486" ]; then
+  SLKCFLAGS="-O2 -march=i486 -mtune=i686"
+  LIBDIRSUFFIX=""
+elif [ "$ARCH" = "i686" ]; then
+  SLKCFLAGS="-O2 -march=i686 -mtune=i686"
+  LIBDIRSUFFIX=""
+elif [ "$ARCH" = "x86_64" ]; then
+  SLKCFLAGS="-O2 -fPIC"
+  LIBDIRSUFFIX="64"
+else
+  SLKCFLAGS="-O2"
+  LIBDIRSUFFIX=""
+fi
+
+MAJOR=$(echo $VERSION | cut -d. -f1)
+
+set -e
+
+rm -rf $PKG
+mkdir -p $TMP $PKG $OUTPUT
+cd $TMP
+rm -rf ARPACK
+tar xvf $CWD/${PRGNAM}96.tar.gz
+tar xvf $CWD/patch.tar.gz
+
+patch -p0 < $CWD/arpack-2.1-redhat.patch
+patch -p0 < $CWD/arpack-etime.patch
+patch -p0 < $CWD/arpack-second-bug.patch
+
+cd ARPACK
+chown -R root:root .
+find . \
+ \( -perm 777 -o -perm 775 -o -perm 711 -o -perm 555 -o -perm 511 \) \
+ -exec chmod 755 {} \; -o \
+ \( -perm 666 -o -perm 664 -o -perm 600 -o -perm 444 -o -perm 440 -o -perm 400 \) \
+ -exec chmod 644 {} \;
+
+mkdir static shared
+cd shared
+for dir in ../SRC ../UTIL; do
+  make -f $dir/Makefile VPATH=$dir srcdir=$dir \
+       FC=gfortran FFLAGS="-std=legacy $SLKCFLAGS -fPIC" \
+       single double complex complex16
+done
+gcc -shared -lblas -llapack -Wl,-soname,libarpack.so.$MAJOR -o libarpack.so.$VERSION *.o
+cd ..
+cd static
+for dir in ../SRC ../UTIL; do
+  make -f $dir/Makefile VPATH=$dir srcdir=$dir \
+       FC=gfortran FFLAGS="-std=legacy $SLKCFLAGS" \
+       single double complex complex16
+done
+ar rv libarpack.a *.o
+ranlib libarpack.a
+cd ..
+
+mkdir -p $PKG/usr/lib${LIBDIRSUFFIX}
+install -p -m644 static/libarpack.a $PKG/usr/lib${LIBDIRSUFFIX}
+install -p -m755 shared/libarpack.so.$VERSION $PKG/usr/lib${LIBDIRSUFFIX}
+( cd $PKG/usr/lib${LIBDIRSUFFIX}
+  ln -s libarpack.so.$VERSION libarpack.so.$MAJOR
+  ln -s libarpack.so.$MAJOR libarpack.so )
+
+find $PKG | xargs file | grep -e "executable" -e "shared object" | grep ELF \
+  | cut -f 1 -d : | xargs strip --strip-unneeded 2> /dev/null || true
+
+mkdir -p $PKG/usr/doc/$PRGNAM-$VERSION
+cp -a DOCUMENTS/* $PKG/usr/doc/$PRGNAM-$VERSION
+cat $CWD/$PRGNAM.SlackBuild > $PKG/usr/doc/$PRGNAM-$VERSION/$PRGNAM.SlackBuild
+
+mkdir -p $PKG/install
+cat $CWD/slack-desc > $PKG/install/slack-desc
+
+cd $PKG
+/sbin/makepkg -l y -c n $OUTPUT/$PRGNAM-$VERSION-$ARCH-$BUILD$TAG.${PKGTYPE:-tgz}
diff --git a/academic/arpack/arpack.info b/academic/arpack/arpack.info
new file mode 100644
index 0000000000..0eb33f272f
--- /dev/null
+++ b/academic/arpack/arpack.info
@@ -0,0 +1,10 @@
+PRGNAM="arpack"
+VERSION="2.1"
+HOMEPAGE="http://www.caam.rice.edu/software/ARPACK/"
+DOWNLOAD="http://www.caam.rice.edu/software/ARPACK/SRC/arpack96.tar.gz http://www.caam.rice.edu/software/ARPACK/SRC/patch.tar.gz"
+MD5SUM="fffaa970198b285676f4156cebc8626e 14830d758f195f272b8594a493501fa2"
+DOWNLOAD_x86_64=""
+MD5SUM_x86_64=""
+MAINTAINER="Pablo Santamaria"
+EMAIL="pablosantamaria@gmail.com"
+APPROVED="dsomero"
diff --git a/academic/arpack/slack-desc b/academic/arpack/slack-desc
new file mode 100644
index 0000000000..4d381e6550
--- /dev/null
+++ b/academic/arpack/slack-desc
@@ -0,0 +1,19 @@
+# HOW TO EDIT THIS FILE:
+# The "handy ruler" below makes it easier to edit a package description.  Line
+# up the first '|' above the ':' following the base package name, and the '|'
+# on the right side marks the last column you can put a character in.  You must
+# make exactly 11 lines for the formatting to be correct.  It's also
+# customary to leave one space after the ':' except on otherwise blank lines.
+
+      |-----handy-ruler------------------------------------------------------|
+arpack: arpack (Fortran library for solving large scale eigenvalue problems)
+arpack:
+arpack: ARPACK is a collection of Fortran77 subroutines designed to solve
+arpack: large scale eigenvalue problems.
+arpack:
+arpack:
+arpack:
+arpack:
+arpack:
+arpack:
+arpack: