Evaluation of solvent effects in isotropic and anisotropic dielectrics, and in ionic solutions, with a unified integral equation method: theoretical bases, implementation and numerical applications