https://github.com/cran/fields
Tip revision: e6e2dec9c9cc857b2226614aaf6c6642000af53c authored by Doug Nychka on 06 February 2009, 00:00:00 UTC
version 5.02
version 5.02
Tip revision: e6e2dec
m2deb.f
subroutine m2deb( nd,x1,n1, x2,n2, par, c,nc,h,work)
implicit double precision (a-h,o-z)
integer n1,n2, ir,nc,nd,ic
real*8 par(nd),x1(n1,nd), x2(n2,nd), c(n2,nc), h(n1,nc)
real*8 work( nc), xtemp1, xtemp2,temp, temp2
c****** work aray must be dimensioned to size n2
c **** loop through columns of output matrix K
c*** outer most loop over columns of x1 and x2 should reduce paging
do 4 ic=1,nc
work(ic)=0
4 continue
do 5 ir= 1, n1
xtemp1= x1(ir,1)
xtemp2= x1(ir,2)
c evaluate row of exp covaraince matrix
do 10 j =1,n2
temp= sqrt((xtemp1- x2(j,1))**2 +(xtemp2- x2(j,2))**2)
temp2= exp( -temp)
do 15 ic= 1,nc
work(ic)= work(ic)+ temp2 * c(j,ic)
15 continue
10 continue
do 25 ic=1,nc
h(ir,ic)= work(ic)
work( ic)=0
25 continue
5 continue
return
end