Skip to content
GitLab
Menu
Projects
Groups
Snippets
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in
Toggle navigation
Menu
Open sidebar
Gebhardt, Albrecht
akima
Commits
e0adeb2f
Commit
e0adeb2f
authored
Jan 01, 2002
by
agebhard
Browse files
T/F changed to TRUE/FALSE
parent
c4fdbfe2
Changes
5
Hide whitespace changes
Inline
Side-by-side
DESCRIPTION
View file @
e0adeb2f
Package: akima
Version: 0.3-
2
Version: 0.3-
3
Title: Interpolation of irregularly spaced data
Author: Fortran code by H. Akima
R port by Albrecht Gebhardt <albrecht.gebhardt@uni-klu.ac.at>
...
...
R/interp.new.R
View file @
e0adeb2f
"interp.new"
<-
function
(
x
,
y
,
z
,
xo
=
seq
(
min
(
x
),
max
(
x
),
length
=
40
),
yo
=
seq
(
min
(
y
),
max
(
y
),
length
=
40
),
linear
=
F
,
yo
=
seq
(
min
(
y
),
max
(
y
),
length
=
40
),
linear
=
F
ALSE
,
ncp
=
NULL
,
extrap
=
FALSE
,
duplicate
=
"error"
,
dupfun
=
NULL
)
{
...
...
@@ -50,7 +50,7 @@
n
<-
length
(
x
)
}
else
{
ord
<-
(
hist
(
i
,
plot
=
F
,
freq
=
T
,
breaks
=
seq
(
0.5
,
max
(
i
)
+0.5
,
1
))
$
counts
==
1
)
ord
<-
(
hist
(
i
,
plot
=
F
ALSE
,
freq
=
T
RUE
,
breaks
=
seq
(
0.5
,
max
(
i
)
+0.5
,
1
))
$
counts
==
1
)
x
<-
x
[
ord
]
y
<-
y
[
ord
]
z
<-
z
[
ord
]
...
...
@@ -63,7 +63,7 @@
zo
<-
matrix
(
0
,
nx
,
ny
)
storage.mode
(
zo
)
<-
"double"
miss
<-
!
extrap
#if not extrapolating use missing values
extrap
<-
matrix
(
T
,
nx
,
ny
)
extrap
<-
matrix
(
T
RUE
,
nx
,
ny
)
if
(
!
is.null
(
ncp
)){
if
(
extrap
&
ncp
==
0
)
warning
(
"Cannot extrapolate with linear option"
)
...
...
R/interp.old.R
View file @
e0adeb2f
...
...
@@ -39,7 +39,7 @@
n
<-
length
(
x
)
}
else
{
ord
<-
(
hist
(
i
,
plot
=
F
,
freq
=
T
,
breaks
=
seq
(
0.5
,
max
(
i
)
+0.5
,
1
))
$
counts
==
1
)
ord
<-
(
hist
(
i
,
plot
=
F
ALSE
,
freq
=
T
RUE
,
breaks
=
seq
(
0.5
,
max
(
i
)
+0.5
,
1
))
$
counts
==
1
)
x
<-
x
[
ord
]
y
<-
y
[
ord
]
z
<-
z
[
ord
]
...
...
R/interpp.new.R
View file @
e0adeb2f
...
...
@@ -44,7 +44,7 @@
n
<-
length
(
x
)
}
else
{
ord
<-
(
hist
(
i
,
plot
=
F
,
freq
=
T
,
breaks
=
seq
(
0.5
,
max
(
i
)
+0.5
,
1
))
$
counts
==
1
)
ord
<-
(
hist
(
i
,
plot
=
F
ALSE
,
freq
=
T
RUE
,
breaks
=
seq
(
0.5
,
max
(
i
)
+0.5
,
1
))
$
counts
==
1
)
x
<-
x
[
ord
]
y
<-
y
[
ord
]
z
<-
z
[
ord
]
...
...
@@ -57,7 +57,7 @@
zo
<-
rep
(
0
,
np
)
storage.mode
(
zo
)
<-
"double"
miss
<-
!
extrap
#if not extrapolating use missing values
extrap
<-
seq
(
T
,
np
)
extrap
<-
seq
(
T
RUE
,
np
)
if
(
extrap
&
ncp
==
0
)
warning
(
"Cannot extrapolate with linear option"
)
ans
<-
.Fortran
(
"sdbi3p"
,
...
...
R/interpp.old.R
View file @
e0adeb2f
...
...
@@ -44,7 +44,7 @@
n
<-
length
(
x
)
}
else
{
ord
<-
(
hist
(
i
,
plot
=
F
,
freq
=
T
,
breaks
=
seq
(
0.5
,
max
(
i
)
+0.5
,
1
))
$
counts
==
1
)
ord
<-
(
hist
(
i
,
plot
=
F
ALSE
,
freq
=
T
RUE
,
breaks
=
seq
(
0.5
,
max
(
i
)
+0.5
,
1
))
$
counts
==
1
)
x
<-
x
[
ord
]
y
<-
y
[
ord
]
z
<-
z
[
ord
]
...
...
Write
Preview
Markdown
is supported
0%
Try again
or
attach a new file
.
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment