{ kkmap.i kkinit(var v: vector; var u: ltuple); initialize stuff df(var v: vector; var u,fu: ltuple; var dfu: matrix); fu:=F(u) and dfu:=DF(u), for given v findu(var v: vector; var u: ltuple); solve F(u)=0 for given v dunpack1(d1: integer; var u: ltuple; var ve,v1,v2: vector); v2:=(1st part of DU applied to v1) procedure unpack(var v1,v2: vector; var u: ltuple); v2:=U(v1), pass on u from findu(v1,u) uok(var u: ltuple; var v: vector): boolean; check that u is within the allowed range ustar(var v: vector; var u: ltuple); findu and check domains kk(var v1,v2: vector; var u: ltuple); v2:=K(v1), pass on u from findu(v1,u) vvmtprod(var m: matrix; var vv1, vv2: vtuple); vv2:=Sum_j vv1[j]*m[j,.] dkkinit(var u: ltuple; var v,ve: vector; var vv: vtuple); prepare vectors ve and vv[i] for dkk dunpack2(var vv: vtuple; var v1,v2: vector); v2:=(2nd part of DU applied to v1) dkk(d1: integer; var u: ltuple; var vv: vtuple; var ve,vus,hi,ho: vector); ho:=DK(v)hi, given vus Notes: Run kkinit first. The global quantities are procedures from reps.i, scalar.i, vector.i, linear.i const topint= ; strnglen= ; exmax= ; hdim= ; edim= ; vdim= ; lmax=3; up=true; down=false; type longstring=varying[strnglen] of char; hint=array[0..hdim] of integer; hreal=record si: integer; ma: hint; ex: integer; end; scalar=record l: hreal; u: hreal; end; vector=array[0..vdim] of scalar; ltuple=array[0..lmax] of scalar; matrix=array[0..lmax,0..lmax] of scalar; vtuple=array[0..lmax] of vector; var trunc,undeflow,overflow: boolean; e: array[0..edim] of integer; sone,stwo,shalf: scalar; normalize: boolean; beta2,cov,rho,sigma,tau: scalar; ( hiemodel parameters ) vkappa: vector; ( bounds for convolution ) umax,lnormal: ltuple; ------------------------------------------------------------------------------} %include 'kkmap/kkinit'; %include 'kkmap/df'; %include 'kkmap/findu'; %include 'kkmap/dunpack1'; %include 'kkmap/unpack'; %include 'kkmap/uok'; %include 'kkmap/ustar'; %include 'kkmap/kk'; %include 'kkmap/vvmtprod'; %include 'kkmap/dkkinit'; %include 'kkmap/dunpack2'; %include 'kkmap/dkk';