The implementation of a local correlation (LC) treatment of multireference (MR) configuration interaction approaches within the COLUMBUS program system is reported. The LC treatment is based on the weak pairs approximation of Sæbø and Pulay (Ann. Rev. Phys. Chem. 1993, 44, 213) and a geometrical analysis of Walter et al. (Chem. Phys. Lett. 2001, 346, 177). The removal of simultaneous single excitations out of the weak pairs is based on the reference doubly occupied space only, leading to a straightforward program implementation and a conceptual simplicity in terms of well-defined localized orbitals. Reductions of up to an order of magnitude in the configuration space expansion and in computer time for the Davidson diagonalization step are found. The selection of the active and the virtual orbital spaces is not affected by this procedure. This treatment is successfully applied to the singlet biradical heptazethrene and its different acceptor-donor substituents: 4,12-dicyanoheptazethrene, 4,12-diaminoheptazethrene, and 4-Amino-12-cyanoheptazethrene. Simultaneous insertion of pairs of donor and acceptor groups increases the biradical character; for push-pull substitution, this effect is significantly smaller. In addition, results obtained from spin-corrected unrestricted density functional theory calculations are supported by our MR calculations.