hi: below almost gets you there except for the endpoints. see filter for
more information.
d=0.5
L=20
x=seq(20, by=1, length.out=20)
temp <- as.numeric(filter(x,filter=c(d,(1-2*d),d),sides=2))
print(temp)
On Sun, Nov 9, 2008 at 10:57 PM, stephen sefick wrote:
#Is there a way to vectorize the for loop
#maybe a fancy indexing trick?
#thanks
d=0.5
L=20
x=seq(20, by=1, length.out=20)
reflecting <- function(pre, d, L){
r=L-1
x=rep(0, L)
for(j in 2:r){
x[j]=((1-(2*d))*pre[j])+(d*pre[(j+1)])+(d*pre[(j-1)])
}
x[1]=((1-d)*pre[1])+(d*pre[2])
x[L]=((1-d)*pre[L])+(d*pre[(L-1)])
y=x
}
f = reflecting(x, d, L)
--
Stephen Sefick
Research Scientist
Southeastern Natural Sciences Academy
Let's not spend our time and resources thinking about things that are
so little or so large that all they really do for us is puff us up and
make us feel like gods. We are mammals, and have not exhausted the
annoying little problems of being mammals.
-K. Mullis
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
