function f=D4etoile(x);

%                     f=D4etoile(x);
%
% Action de l'opérateur adjoint de D4 sur un vecteur ligne x. 

N=length(x);
G=zeros(1,8);
G(1)= -.010597;
G(2)= -.032883; 
G(3)= .030841; 
G(4)= .187035; 
G(5)= -.027984; 
G(6)= -.630881; 
G(7)= .714846; 
G(8)= -.230378; 
ff=zeros(1,2*N+2);
for k=9:2*N,
	if (mod(k,2) == 1)
		ff(k)=G(1)*x((k-1)/2)+G(3)*x((k-3)/2)+G(5)*x((k-5)/2)+G(7)*x((k-7)/2);
	else
		ff(k)=G(2)*x((k-2)/2)+G(4)*x((k-4)/2)+G(6)*x((k-6)/2)+G(8)*x((k-8)/2);
	end;
end;
tt=1:2*N-7;
f(tt)=ff(tt+9);