From 83519b3773fe59053d3a99a3ae12d186561f18bc Mon Sep 17 00:00:00 2001
From: agebhard <>
Date: Thu, 23 May 2002 08:58:04 +0000
Subject: [PATCH] noch mehr fehler?

---
 R/bk.grid.R           |  4 ++--
 devel/fresh/.Rhistory |  2 +-
 src/glsfit.f          | 15 +++++++++------
 3 files changed, 12 insertions(+), 9 deletions(-)

diff --git a/R/bk.grid.R b/R/bk.grid.R
index c05e94c..12abbdc 100644
--- a/R/bk.grid.R
+++ b/R/bk.grid.R
@@ -210,7 +210,7 @@ bk.grid <- function(point.obj,
     typpr<-prior$info
     typpr[prior$type=="subjective"]<-typpr[prior$type=="subjective"]*(-1)
 
-browser()
+#browser()
     if(FALSE)
     ans<-.C("bk_grid",
 #            xsw     = as.double(xsw),
@@ -338,7 +338,7 @@ if(TRUE)        ans<-.C("bk_grid",
             glsmth  = as.integer(method),
 #           ,.Package= "baykrig"
             )
-browser()
+#browser()
     retval<-list(x=ans$xg,
                  y=ans$yg,
                  z=matrix(ans$zg,nx,ny),
diff --git a/devel/fresh/.Rhistory b/devel/fresh/.Rhistory
index 2b6f7c8..f95f33d 100644
--- a/devel/fresh/.Rhistory
+++ b/devel/fresh/.Rhistory
@@ -29,4 +29,4 @@ library(baykrig)
 library(baykrig)
 leman.bk<- bk.grid(point = leman.88.pt , at = "cadpbm", prior=leman.prior,var.mod.obj = leman.88.vmsph, xsw=min(leman.bank$x),ysw=min(leman.bank$y), xne=max(leman.bank$x), yne=max(leman.bank$y), nx=100, ny=100, trend=1, rsearch = 10, extrap = F,border=leman.bank, duplicate="mean",method="ols")
 library(baykrig)
-leman.bk<- bk.grid(point = leman.88.pt , at = "cadpbm", prior=leman.prior,var.mod.obj = leman.88.vmsph, xsw=min(leman.bank$x),ysw=min(leman.bank$y), xne=max(leman.bank$x), yne=max(leman.bank$y), nx=100, ny=100, trend=1, rsearch = 10, extrap = F,border=leman.bank, duplicate="mean")
+leman.bk<- bk.grid(point = leman.88.pt , at = "cadpbm", prior=leman.prior,var.mod.obj = leman.88.vmsph, xsw=min(leman.bank$x),ysw=min(leman.bank$y), xne=max(leman.bank$x), yne=max(leman.bank$y), nx=5, ny=5, trend=1, rsearch = 10, extrap = F,border=leman.bank, duplicate="mean")
diff --git a/src/glsfit.f b/src/glsfit.f
index 75e4ef5..977d015 100644
--- a/src/glsfit.f
+++ b/src/glsfit.f
@@ -201,11 +201,11 @@ ccc
 c      goto 123
 
 ccccccc FEHLER ccccccccc:
-      CALL DGEMM('N','N',NTREND,NTREND,N,ONE,
-     .            CWORK,LDCWRK,FMAT,LDF,ZERO,COVBTA,LDCVBT)
+c      CALL DGEMM('N','N',NTREND,NTREND,N,ONE,
+c     .            CWORK,LDCWRK,FMAT,LDF,ZERO,COVBTA,LDCVBT)
 
 ccc
-      goto 123
+c      goto 123
       DO 130 I=1,NTREND
          DO 120 J=1,NTREND
             IF( I.EQ.J) THEN
@@ -218,9 +218,11 @@ ccc
 
 
 c     destroys CHLUP with inverse of F'*C**-1*F:
-      CALL DSYSVX( 'N', 'U', NTREND, NTREND, COVBTA, LDCVBT, CWRK2, 
-     $             LDCWK2, IPIV, CWORK, LDCWRK, CHLUP, LDCLUP,
-     $             RCOND, FERR, BERR, WORK, LWORK, IWORK, IERR )
+ccc naechster fehler:
+c      goto 123
+c      CALL DSYSVX( 'N', 'U', NTREND, NTREND, COVBTA, LDCVBT, CWRK2, 
+c     $             LDCWK2, IPIV, CWORK, LDCWRK, CHLUP, LDCLUP,
+c     $             RCOND, FERR, BERR, WORK, LWORK, IWORK, IERR )
       IF ( IERR .NE. 0) THEN         
          IF ( IERR .LT. 0 ) THEN 
             CALL ERRMSG('GLSFIT DSYSVX: argument no IERR wrong',38, 
@@ -242,6 +244,7 @@ c     destroys CHLUP with inverse of F'*C**-1*F:
  140     CONTINUE
  150  CONTINUE
       
+      goto 123
       GO TO 4
 
 c     direct computation of BETA:
-- 
GitLab