% Returns the skew of A
function S = skewPart(A, tol)
    S = 1/2*(A - A');
end
