%
% test of two grid code with full data structure; test3a.m
%
global sigma 
sigma=0;
lmax=10;
itmax=5;
convhist=zeros(itmax,lmax-4);
lv=2^(lmax+1)+lmax-2;
v=zeros(lv,1); f=zeros(lv,1);
vt=v; ft=f;
for l=5:lmax
        ntl=2^(l+1)+l-2; ifu=ntl; ifl=ntl-2^l;
	nl=2^l; nt=nl+1; x=0:nl; x=x'/nl; 
        v(ifl:ifu)=sin(6*pi*x)+sin(17*pi*x);
	f(ifl:ifu)=zeros(nt,1);
	for it = 1:itmax
		[vt(1:ifu),f(1:ifu)]=twog(v(1:ifu),f(1:ifu),2,2);
		convhist(it,l-4)=norm(vt(ifl:ifu),inf)/norm(v(ifl:ifu),inf);
		v(ifl:ifu)=vt(ifl:ifu);
	end
end
convhist