% Returns the symmetric of A
function S = symmetricPart(A, tol)
    S = 1/2*(A + A');
end
