Mercurial > forge

changeset 4815:9f7194a5867e octave-forge

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
remove unused variabled & externals
author highegg
date Tue, 29 Apr 2008 05:56:25 +0000
parents cdc85bc99433
children 97cf5c1a758a
files main/octgpr/src/dsdacc.f main/octgpr/src/infgpr.f main/octgpr/src/nldgpr.f main/octgpr/src/nllgpr.f main/octgpr/src/optdrv.f main/octgpr/src/stheta.f
diffstat 6 files changed, 7 insertions(+), 8 deletions(-) [+]
line wrap: on
line diff
--- a/main/octgpr/src/dsdacc.f	Mon Apr 28 17:01:39 2008 +0000
+++ b/main/octgpr/src/dsdacc.f	Tue Apr 29 05:56:25 2008 +0000
@@ -30,7 +30,6 @@
       integer n
       double precision x(*),y(*),x2,xy
       integer i
-      external xerbla
       xy = 0
       x2 = 0
       do i = 1,n
--- a/main/octgpr/src/infgpr.f	Mon Apr 28 17:01:39 2008 +0000
+++ b/main/octgpr/src/infgpr.f	Tue Apr 29 05:56:25 2008 +0000
@@ -67,7 +67,7 @@
       external corr
       external dwdis2,dsdacc,dtpsv,dcopy,xerbla
       real*8 tmp,dwdis2
-      integer i,j,k,info
+      integer i,k,info
 c argument checks
       info = 0
       if (ndim < 0) then
--- a/main/octgpr/src/nldgpr.f	Mon Apr 28 17:01:39 2008 +0000
+++ b/main/octgpr/src/nldgpr.f	Tue Apr 29 05:56:25 2008 +0000
@@ -39,10 +39,10 @@
 c               info = 0 no problem
 c               info = 1 invalid correlation matrix from nllgpr (dtrtri
 c                        failure)
-      integer ndim,nx,nlin,info
+      integer ndim,nx,info
       real*8 X(ndim,nx),theta(ndim),nu,var
       real*8 R(nx,0:nx),dtheta(ndim),dnu(2)
-      external dtrsv,dtrtri,dsyr,xerbla
+      external dtrsv,dpotri,dsyr,xerbla
       integer i,j,k
       real*8 tmp,tmp2,m2
 
@@ -61,7 +61,7 @@
 c form L' \ m = R \ ones(nx,1)
       call dtrsv('L','T','N',nx,R(1,1),nx,R(1,0),1)
 
-c invert the triangular factor
+c turn the factorization into inversion
       call dpotri('L',nx,R(1,1),nx,info)
       if (info /= 0) goto 501
 
--- a/main/octgpr/src/nllgpr.f	Mon Apr 28 17:01:39 2008 +0000
+++ b/main/octgpr/src/nllgpr.f	Tue Apr 29 05:56:25 2008 +0000
@@ -65,7 +65,7 @@
       external corr
       external dwdis2,dcopy,daxpy,dscal,dtrsv,dtrsm,xerbla,
      +dpotrf,dtrtrs,dgeqr2,dorm2r
-      integer i,j,k,nl1
+      integer i,j,nl1
       real*8 sums,sum,dwdis2
 
 c argument checks
--- a/main/octgpr/src/optdrv.f	Mon Apr 28 17:01:39 2008 +0000
+++ b/main/octgpr/src/optdrv.f	Tue Apr 29 05:56:25 2008 +0000
@@ -72,7 +72,7 @@
       logical l2nu
       double precision VM(*),CP(13)
       double precision stp(ndim+1),grd(ndim+1),Zg(ndim),Zba(ndim),
-     +snlo,snup,eps,ss,relr
+     +snlo,snup,eps,relr
       double precision dlamch,dnrm2,ddot
       external dlamch,dspmid,trstp,dnrm2,ddot
       integer iB,iba,iW,iZ,i
--- a/main/octgpr/src/stheta.f	Mon Apr 28 17:01:39 2008 +0000
+++ b/main/octgpr/src/stheta.f	Tue Apr 29 05:56:25 2008 +0000
@@ -35,7 +35,7 @@
       integer ndim,nx
       real*8 X(ndim,nx),theta(ndim)
       real*8 mX(ndim)
-      integer i,j,k
+      integer i,k
       do k = 1,ndim
         mX(k) = 0
       end do