A density‐matrix divide‐and‐conquer approach for electronic structure calculations of large molecules