function B = ihaar_col(A);

% This function inverts the sums and differences column-wise for the Haar
% Transform. The result is then placed back into the input matrix.

n = length(A);
for j = 1:n
num = 1;
for i = 1:(n/2)
B(num,j) = (A(i,j) + A(i+(n/2),j))/sqrt(2);
B(num+1,j) = (A(i,j) - A(i+(n/2),j))/sqrt(2);
num = num + 2;
end
end