We investigate different types of complex soliton solutions with regard to their stability against linear perturbations. In the Bullough-Dodd scalar field theory we find linearly stable complex PT-symmetric solutions and linearly unstable solutions for which the PT-symmetry is broken. Both types of solutions have real energies. The auxiliary Sturm-Liouville eigenvalue equation in the stability analysis for the PT-symmetric solutions can be solved exactly by supersymmetrically mapping it to an isospectral partner system involving a shifted and scaled inverse cosh-squared potential. We identify exactly one shape mode in form of a bound state solution and scattering states which when used as linear perturbations leave the solutions stable. The auxiliary problem for the solutions with broken PT-symmetry involves a complex shifted and scaled inverse sin-squared potential. The corresponding bound and scattering state solutions have complex eigenvalues, such that when used as linear perturbations for the corresponding soliton solutions lead to their decay or blow up as time evolves.