Revision 69b0f9dca8eb051f132725ecc679fe1997246e50 authored by Adrian Baddeley on 18 January 2006, 21:47:25 UTC, committed by cran-robot on 18 January 2006, 21:47:25 UTC
1 parent cb2215f
strauss.f
C Output from Public domain Ratfor, version 1.0
subroutine strauss(u,v,ix,x,y,n,par,period,cifval)
implicit double precision(a-h,o-z)
dimension par(3), x(1), y(1), period(2)
logical per
eps = 2.22d-16
zero = 0.d0
per = period(1) .gt. zero
beta = par(1)
gamma = par(2)
r = par(3)
kount = 0
do23000 j = 1,n
if(j .eq. ix)then
continue
else
if(per)then
call dist2(u,v,x(j),y(j),period,d2)
else
d2 = (u-x(j))**2 + (v-y(j))**2
endif
if(d2 .lt. r)then
kount = kount+1
endif
endif
23000 continue
23001 continue
if(gamma .lt. eps )then
if(kount .gt. 0)then
cifval = zero
else
cifval = beta
endif
else
cifval = beta*exp(log(gamma)*dble(kount))
endif
return
end
Computing file changes ...