Commit 740d9cb4 authored by agebhard's avatar agebhard

continuing C wrapper ...

parent b8c26736
......@@ -412,9 +412,9 @@ bk.grid <- function(point.obj,
.Package= "baykrig")
ans$bits<-ans$intvec[(16+npr+7*n+2*ntrend+1):(16+npr+7*n+2*ntrend+nz*n+1)]
retval<-list(x=ans$dblvec[1:nx],
y=ans$dblvec[(nx+1):(nx+ny)],
browser()
retval<-list(x=ans$xg,
y=ans$yg,
z=matrix(ans$zg,nx,ny),
var=matrix(ans$varg,nx,ny),
done=matrix(ans$dog, nx, ny),
......
No preview for this file type
leman.bk.grid<-bk.points(leman.grid$x,leman.grid$y,leman.88.pt,"cadpbm",leman.ebp,leman.88.vmsph,maxdist=10,trend=1,method="gqr")
library(baykrig)
leman.bk.grid<-bk.points(leman.grid$x,leman.grid$y,leman.88.pt,"cadpbm",leman.ebp,leman.88.vmsph,maxdist=10,trend=1,method="gqr")
leman.bk.grid$z
bk.grid
leman.bk.grid<-bk.points(leman.grid$x,leman.grid$y,leman.88.pt,"cadpbm",leman.ebp,leman.88.vmsph,maxdist=10,trend=1,method="gqr")
library(baykrig)
leman.bk.grid<-bk.points(leman.grid$x,leman.grid$y,leman.88.pt,"cadpbm",leman.ebp,leman.88.vmsph,maxdist=10,trend=1,method="gqr")
leman.bk.grid$z
>>>>>>> 1.56
ls()
library(baykrig)
library(kdesign)
library(baykrig)
library(baykrig)
library(baykrig)
library(baykrig)
library(baykrig)
library(baykrig)
leman.bk <- bk.grid(point = leman.pt, at = "cadpbm", prior=leman.prior,var.mod.obj = leman.sph, 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")
example(bk.grid)
library(baykrig)
require(rgeostat) if(length(ls(pat="leman.78"))==0)data(leman) if(length(ls(pat="leman.bank"))==0)data(leman.bank) if(length(ls(pat="leman.pt"))==0) leman.pt<-point(leman.88) if(length(ls(pat="leman.pr"))==0) leman.pr<-pair(leman.pt,maxdist=30) if(length(ls(pat="leman.ev"))==0) leman.ev<-est.variogram(leman.pt,leman.pr,"cadpbm") if(length(ls(pat="leman.sph"))==0) leman.sph<-fit.variogram("spherical",leman.ev,0.1,0.2,20,iter=50,plot.it=T) if(length(ls(pat="leman.prior"))==0) example(empirical.prior)
q()
library(baykrig)
library(baykrig)
library(baykrig)
library(baykrig)
leman.bk <- bk.grid(point = leman.pt, at = "cadpbm", prior=leman.prior,var.mod.obj = leman.sph, 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")
library(baykrig)
leman.bk <- bk.grid(point = leman.pt, at = "cadpbm", prior=leman.prior,var.mod.obj = leman.sph, 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 = T,border=leman.bank, duplicate="mean")
ls()
q()
leman.prior.ord
leman.bk.ord <- bk.grid(point = leman.pt, at = "cadpbm", prior=leman.prior.ord,var.mod.obj = leman.sph, 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=0, rsearch = 10, extrap = T,border=leman.bank, duplicate="mean"); save.image()
plot.krige.map(leman.bk.ord)
str(leman.bk.ord$var)
str(leman.bk.ord$z)
leman.bk.ord$z
leman.bk.ord$var
image(leman.bk.ord$z)
contour(leman.bk.ord$z)
leman.bk.ord <- bk.grid(point = leman.pt, at = "cadpbm", prior=leman.prior.ord,var.mod.obj = leman.sph, xsw=min(leman.bank$x),ysw=min(leman.bank$y), xne=max(leman.bank$x), yne=max(leman.bank$y), nx=50, ny=50, trend=0, rsearch = 10, extrap = T,border=leman.bank, duplicate="mean"); save.image()
library(rgeostat);library(baykrig)
leman.bk.ord <- bk.grid(point = leman.pt, at = "cadpbm", prior=leman.prior.ord,var.mod.obj = leman.sph, xsw=min(leman.bank$x),ysw=min(leman.bank$y), xne=max(leman.bank$x), yne=max(leman.bank$y), nx=25, ny=25, trend=0, rsearch = 10, extrap = T,border=leman.bank, duplicate="mean"); save.image()
save.image()
plot.krige.map(leman.bk.ord)
contour(leman.bk.ord$z)
image(leman.bk.ord$z)
image(leman.bk.ord$var)
leman.bk.ord <- bk.grid(point = leman.pt, at = "cadpbm", prior=leman.prior.ord,var.mod.obj = leman.sph, xsw=min(leman.bank$x),ysw=min(leman.bank$y), xne=max(leman.bank$x), yne=max(leman.bank$y), nx=25, ny=25, trend=0, rsearch = 10, extrap = F,border=leman.bank, duplicate="mean"); save.image()
library(rgeostat);library(baykrig);leman.bk.ord <- bk.grid(point = leman.pt, at = "cadpbm", prior=leman.prior.ord,var.mod.obj = leman.sph, xsw=min(leman.bank$x),ysw=min(leman.bank$y), xne=max(leman.bank$x), yne=max(leman.bank$y), nx=25, ny=25, trend=0, rsearch = 10, extrap = F,border=leman.bank, duplicate="mean"); save.image()
library(rgeostat);library(baykrig);leman.bk.ord <- bk.grid(point = leman.pt, at = "cadpbm", prior=leman.prior.ord,var.mod.obj = leman.sph, xsw=min(leman.bank$x),ysw=min(leman.bank$y), xne=max(leman.bank$x), yne=max(leman.bank$y), nx=15, ny=15, trend=0, rsearch = 10, extrap = F,border=leman.bank, duplicate="mean"); save.image()
system("echo $DISPLAY")
x11(display="zeta.home.algepop.net:10.0")
image(leman.bk.ord$z)
lines(leman.bank)
plot.krige.map(leman.bk.ord$z)
plot.krige.map(leman.bk.ord)
str(leman.bk.ord)
library(rgeostat);library(baykrig);leman.bk.ord <- bk.grid(point = leman.pt, at = "cadpbm", prior=leman.prior.ord,var.mod.obj = leman.sph, xsw=min(leman.bank$x),ysw=min(leman.bank$y), xne=max(leman.bank$x), yne=max(leman.bank$y), nx=15, ny=15, trend=0, rsearch = 10, extrap = F,border=leman.bank, duplicate="mean"); save.image()
library(rgeostat);library(baykrig);
leman.bk.ord <- bk.grid(point = leman.pt, at = "cadpbm", prior=leman.prior.ord,var.mod.obj = leman.sph, 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=0, rsearch = 10, extrap = T,border=leman.bank, duplicate="mean");
library(rgeostat);library(baykrig);
leman.bk.ord <- bk.grid(point = leman.pt, at = "cadpbm", prior=leman.prior.ord,var.mod.obj = leman.sph, xsw=min(leman.bank$x),ysw=min(leman.bank$y), xne=max(leman.bank$x), yne=max(leman.bank$y), nx=15, ny=15, trend=0, rsearch = 10, extrap = T,border=leman.bank, duplicate="mean");
This diff is collapsed.
SRCS=\
linpack/dgefa.f \
linpack/dgedi.f \
bkwrp.f \
lsfit.f \
lapack/util/ilaenv.f \
lapack/util/ieeeck.f \
lapack/util/lsame.f \
lapack/double/dpotrf.f \
lapack/double/dtptri.f \
lapack/double/dlasyf.f \
lapack/double/dlacon.f \
lapack/double/dlacpy.f \
lapack/double/dgeqrf.f \
lapack/double/dgerqf.f \
lapack/double/dggqrf.f \
lapack/double/dlaswp.f \
lapack/double/dtrtri.f \
lapack/double/dsycon.f \
lapack/double/dlassq.f \
lapack/double/dormrq.f \
lapack/double/dpotf2.f \
lapack/double/dlatrs.f \
lapack/double/dsytrf.f \
lapack/double/dlarfb.f \
lapack/double/dtrcon.f \
lapack/double/dgeqr2.f \
lapack/double/dgetrs.f \
lapack/double/dpotri.f \
lapack/double/dlauum.f \
lapack/double/dgerq2.f \
lapack/double/drscl.f \
lapack/double/dsysvx.f \
lapack/double/dsytf2.f \
lapack/double/dormr2.f \
lapack/double/dsyrfs.f \
lapack/double/dlauu2.f \
lapack/double/dgesv.f \
lapack/double/dsytrs.f \
lapack/double/dgetf2.f \
lapack/double/dlapy2.f \
lapack/double/dgelqf.f \
lapack/double/dlange.f \
lapack/double/dlascl.f \
lapack/double/dlaset.f \
lapack/double/dormqr.f \
lapack/double/dlantr.f \
lapack/double/dtrti2.f \
lapack/double/dlarft.f \
lapack/double/dlamch.f \
lapack/double/dlarf.f \
lapack/double/dlarfg.f \
lapack/double/dlansy.f \
lapack/double/dlabad.f \
lapack/double/dgelq2.f \
lapack/double/dgetrf.f \
lapack/double/dorml2.f \
lapack/double/dormlq.f \
lapack/double/dggglm.f \
lapack/double/dorm2r.f \
lapack/double/dgels.f \
lapack/double/dlae2.f \
lapack/double/dlaev2.f \
lapack/double/dlanst.f \
lapack/double/dlartg.f \
lapack/double/dlasr.f \
lapack/double/dlasrt.f \
lapack/double/dlatrd.f \
lapack/double/dorg2l.f \
lapack/double/dorg2r.f \
lapack/double/dorgql.f \
lapack/double/dorgqr.f \
lapack/double/dorgtr.f \
lapack/double/dsteqr.f \
lapack/double/dsterf.f \
lapack/double/dsyev.f \
lapack/double/dsytd2.f \
lapack/double/dsytrd.f \
design.f \
bk.f \
bktile.f \
dsysvq.f \
dgelse.f \
slatec/src/j4save.f \
slatec/src/fdump.f \
slatec/src/i1mach.f \
slatec/src/xercnt.f \
slatec/src/xgetua.f \
slatec/src/dpsort.f \
slatec/src/xerhlt.f \
slatec/src/xersve.f \
slatec/src/xerprn.f \
slatec/src/xermsg.f \
bkgrid.f \
preload/xerbla.f \
errmsg.f \
srchnb.f \
dgggle.f \
tools.f \
bkpts.f \
blas/dger.f \
blas/dtrsv.f \
blas/dtpmv.f \
blas/dgemm.f \
blas/dtrsm.f \
blas/dspr.f \
blas/dtrmv.f \
blas/dsymv.f \
blas/dtrmm.f \
blas/dsyrk.f \
blas/dgemv.f \
blas/dsyr.f \
blas/dsyr2.f \
blas/dsyr2k.f \
glsfit.f \
covfn.f \
matpr.c \
bk_grid.c
OBJS=\
linpack/dgefa.o \
linpack/dgedi.o \
bkwrp.o \
lsfit.o \
lapack/util/ilaenv.o \
lapack/util/ieeeck.o \
lapack/util/lsame.o \
lapack/double/dpotrf.o \
lapack/double/dtptri.o \
lapack/double/dlasyf.o \
lapack/double/dlacon.o \
lapack/double/dlacpy.o \
lapack/double/dgeqrf.o \
lapack/double/dgerqf.o \
lapack/double/dggqrf.o \
lapack/double/dlaswp.o \
lapack/double/dtrtri.o \
lapack/double/dsycon.o \
lapack/double/dlassq.o \
lapack/double/dormrq.o \
lapack/double/dpotf2.o \
lapack/double/dlatrs.o \
lapack/double/dsytrf.o \
lapack/double/dlarfb.o \
lapack/double/dtrcon.o \
lapack/double/dgeqr2.o \
lapack/double/dgetrs.o \
lapack/double/dpotri.o \
lapack/double/dlauum.o \
lapack/double/dgerq2.o \
lapack/double/drscl.o \
lapack/double/dsysvx.o \
lapack/double/dsytf2.o \
lapack/double/dormr2.o \
lapack/double/dsyrfs.o \
lapack/double/dlauu2.o \
lapack/double/dgesv.o \
lapack/double/dsytrs.o \
lapack/double/dgetf2.o \
lapack/double/dlapy2.o \
lapack/double/dgelqf.o \
lapack/double/dlange.o \
lapack/double/dlascl.o \
lapack/double/dlaset.o \
lapack/double/dormqr.o \
lapack/double/dlantr.o \
lapack/double/dtrti2.o \
lapack/double/dlarft.o \
lapack/double/dlamch.o \
lapack/double/dlarf.o \
lapack/double/dlarfg.o \
lapack/double/dlansy.o \
lapack/double/dlabad.o \
lapack/double/dgelq2.o \
lapack/double/dgetrf.o \
lapack/double/dorml2.o \
lapack/double/dormlq.o \
lapack/double/dggglm.o \
lapack/double/dorm2r.o \
lapack/double/dgels.o \
lapack/double/dlae2.o \
lapack/double/dlaev2.o \
lapack/double/dlanst.o \
lapack/double/dlartg.o \
lapack/double/dlasr.o \
lapack/double/dlasrt.o \
lapack/double/dlatrd.o \
lapack/double/dorg2l.o \
lapack/double/dorg2r.o \
lapack/double/dorgql.o \
lapack/double/dorgqr.o \
lapack/double/dorgtr.o \
lapack/double/dsteqr.o \
lapack/double/dsterf.o \
lapack/double/dsyev.o \
lapack/double/dsytd2.o \
lapack/double/dsytrd.o \
design.o \
bk.o \
bktile.o \
dsysvq.o \
dgelse.o \
slatec/src/j4save.o \
slatec/src/fdump.o \
slatec/src/i1mach.o \
slatec/src/xercnt.o \
slatec/src/xgetua.o \
slatec/src/dpsort.o \
slatec/src/xerhlt.o \
slatec/src/xersve.o \
slatec/src/xerprn.o \
slatec/src/xermsg.o \
bkgrid.o \
preload/xerbla.o \
errmsg.o \
srchnb.o \
dgggle.o \
tools.o \
bkpts.o \
blas/dger.o \
blas/dtrsv.o \
blas/dtpmv.o \
blas/dgemm.o \
blas/dtrsm.o \
blas/dspr.o \
blas/dtrmv.o \
blas/dsymv.o \
blas/dtrmm.o \
blas/dsyrk.o \
blas/dgemv.o \
blas/dsyr.o \
blas/dsyr2.o \
blas/dsyr2k.o \
glsfit.o \
covfn.o \
matpr.o \
bk_grid.o
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment