function D2I2 =smoothie(D2I,H9L,H4,debug=0)
# Usage: D2I2 =smoothie(D2I,H9L,H4)
# Use Laplacians H9L and H4 to improve D2I2 
tmp = (D2I*H9L - H4*D2I);
#tmp(:,sum(H4x)<3)=0;
tmp -= tmp .* (abs(tmp) >1000);
if(debug); print std(nonzeros(tmp)) ; endif
#D2I2 +=tmp/655.36*H4i;
D2I2 = D2I + tmp/655.36;
D2I2 .*= abs(D2I2)>0.000001;
if(debug)
tmp = (H9L * D2I2 - D2I2*H4);
#tmp(:,sum(H4x)<3)=0;
tmp -= tmp .* (abs(tmp) >1000);
std(nonzeros(tmp))
endif
endfunction
