Revision ade4b694b3f416d6b195ddcd584d6f89ef36ea2f authored by M. Helena Gon\xe7alves on 28 January 2012, 00:00:00 UTC, committed by Gabor Csardi on 28 January 2012, 00:00:00 UTC
1 parent 8cb2c8d
integ1.f
``````      subroutine integ1(logL,bt2,beta2,lpsi,omega,npar,
*x2,y2,theta2,work2,n,li,ls,epsabs,epsrel,key,limit)
implicit double precision (a-h,o-z)
External f1,fpij

DOUBLE PRECISION logL,lpsi,lpsi1,li,ls,sigma,sqr2pi
INTEGER n,npar,y2,key,limit,neval,ier,iord,last,
*m,mpar,y1,k,k1,k2,iaux

DIMENSION x1(5000,10),theta1(5000),
*work1(5000),y1(5000),beta1(10),bt1(10),
*x2(n,npar-1),beta2(npar-1),
*theta2(n),y2(n),work2(n),bt2(npar-1),
*alist(limit),blist(limit),elist(limit),iord(limit),
*rlist(limit)

COMMON/param/x1,theta1,work1,
*y1,beta1,bt1,m,mpar,omega1,lpsi1

DATA sqr2pi/2.506628274631D0/

do 10 k=1,(npar-1)
bt1(k)=bt2(k)
beta1(k)=beta2(k)
10 continue
do 30 k2=1,n
do 40 k=1,(npar-1)
x1(k2,k)=x2(k2,k)
40 continue
y1(k2)=y2(k2)
theta1(k2)=theta2(k2)
work1(k2)=work2(k2)
30 continue
m=n
mpar=npar
omega1=omega
lpsi1=lpsi
sigma=dexp(omega/(2.0d0))
a=li*sigma
b=ls*sigma

CALL dqager(f1,a,b,epsabs,epsrel,key,limit,result,abserr,
*neval,ier,alist,blist,rlist,elist,iord,last,1)
logL=dlog(result/(sqr2pi*sigma))

RETURN
END

``````

Computing file changes ...