1function [data,W] = setup_generalwindow(tol)
3% [data,W] = SETUP_GENERALWINDOW(tol)
4% set up fast cheb eval for window
5%
6% See: generalwindow.m
8gam = 0.5; data.gam = 0.5;
9data.tol = tol;
10data.info = {};
12chebApproxInfo.domain = [0,1];
13chebApproxInfo.ninters = max(1,ceil(0.2*log(1/tol))); % *** hack to make sure accurate to tol!;
14chebApproxInfo.nord = 16;
16% table of cheb weights times function values
17data.wei = mktab_wts(@(x) yfun(x,data), chebApproxInfo);
18data.weip = mktab_wts(@(x) ypfun(x,data), chebApproxInfo);
19data.weipp = mktab_wts(@(x) yppfun(x,data), chebApproxInfo);
20data.info = chebApproxInfo;
23% Some more information needed
24theta = log(1/tol); % set up window
25W = ceil(2*theta/(pi*gam));
27end
29function y = yfun(x,data)
30 [y,~,~] = generalwindow(x,data);
31end
32function yp = ypfun(x,data)
33 [~,yp,~] = generalwindow(x,data);
34end
35function ypp = yppfun(x,data)
36 [~,~,ypp] = generalwindow(x,data);
37end