Commit 760972cc authored by agebhard's avatar agebhard
Browse files

Fortran interface to C entry points added

parent 8fb00b82
...@@ -79,7 +79,7 @@ plot(leman.bk0) ...@@ -79,7 +79,7 @@ plot(leman.bk0)
# compare with # compare with
require(rgeostat) require(rgeostat)
leman.krg<-krige.grid.c(min(leman.bank$x),min(leman.bank$y), leman.krg<-krige.grid(min(leman.bank$x),min(leman.bank$y),
max(leman.bank$x), max(leman.bank$y), nx=100, ny=100, point = max(leman.bank$x), max(leman.bank$y), nx=100, ny=100, point =
leman.88.pt, at = "cadpbm", var.mod.obj = leman.88.vmsph, leman.88.pt, at = "cadpbm", var.mod.obj = leman.88.vmsph,
rsearch = 7.5,extrap = F,border=leman.bank) rsearch = 7.5,extrap = F,border=leman.bank)
......
...@@ -419,7 +419,7 @@ c LAMBDA = (C0VEC+FWORK*PHIPR*F0WORK)*(COVMAT+FWORK*PHIPR*FWORK')^-1 ...@@ -419,7 +419,7 @@ c LAMBDA = (C0VEC+FWORK*PHIPR*F0WORK)*(COVMAT+FWORK*PHIPR*FWORK')^-1
c name="kwork\0" c name="kwork\0"
c call matpr(name,kwork,ns,ns,LDkwrk,dbglvl) c call matpr(name,kwork,ns,ns,LDkwrk,dbglvl)
c name="rhswork\0" c name="rhswork\0"
call matpr(name,rhswork,ns,NDO,LDkWrk,dbglvl) c call matpr(name,rhswork,ns,NDO,LDkWrk,dbglvl)
C BETTER USE DSYSVX ? C BETTER USE DSYSVX ?
CALL DGESV(NS,NDO,KWORK,LDKWRK,IPIV,RHSWORK,LDKWRK,INFO) CALL DGESV(NS,NDO,KWORK,LDKWRK,IPIV,RHSWORK,LDKWRK,INFO)
......
#include "bk_grid.h" #include "bk_grid.h"
void F77_NAME(bk_grid)(double *xsw,
double *ysw,
double *xne,
double *yne,
double *angle,
int *nx,
int *ny,
double *dx,
double *dy,
double *xg,
double *yg,
double *zg,
double *varg,
int *dog,
double *lon,
double *lat,
double *z,
int *extrap,
int *n,
int *covtype,
double *covpar,
int *trend,
int *ntrend,
double *mupr,
int *ldmpr,
double *phipr,
int *ldphpr,
double *lonpr,
double *latpr,
int *npr,
int *typpr,
double *rsearch,
int *nsearch,
int *nsmin,
int *nsmax,
int *lwork,
int *mode,
double *lambda,
double *lambd0,
int *searchnb,
int *ierr,
int *glsmth){
/* simple Fortran wrapper */
bk_grid(xsw,
ysw,
xne,
yne,
angle,
nx,
ny,
dx,
dy,
xg,
yg,
zg,
varg,
dog,
lon,
lat,
z,
extrap,
n,
covtype,
covpar,
trend,
ntrend,
mupr,
ldmpr,
phipr,
ldphpr,
lonpr,
latpr,
npr,
typpr,
rsearch,
nsearch,
nsmin,
nsmax,
lwork,
mode,
lambda,
lambd0,
searchnb,
ierr,
glsmth);
}
void bk_grid(double *xsw, void bk_grid(double *xsw,
double *ysw, double *ysw,
double *xne, double *xne,
......
...@@ -45,6 +45,49 @@ void bk_grid(double *xsw, ...@@ -45,6 +45,49 @@ void bk_grid(double *xsw,
int *searchnb, int *searchnb,
int *ierr, int *ierr,
int *glsmth); int *glsmth);
void F77_NAME(bk_grid)(double *xsw,
double *ysw,
double *xne,
double *yne,
double *angle,
int *nx,
int *ny,
double *dx,
double *dy,
double *xg,
double *yg,
double *zg,
double *varg,
int *dog,
double *lon,
double *lat,
double *z,
int *extrap,
int *n,
int *covtype,
double *covpar,
int *trend,
int *ntrend,
double *mupr,
int *ldmpr,
double *phipr,
int *ldphpr,
double *lonpr,
double *latpr,
int *npr,
int *typpr,
double *rsearch,
int *nsearch,
int *nsmin,
int *nsmax,
int *lwork,
int *mode,
double *lambda,
double *lambd0,
int *searchnb,
int *ierr,
int *glsmth);
void F77_NAME(bkgrid)(double *xsw, void F77_NAME(bkgrid)(double *xsw,
double *ysw, double *ysw,
......
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