We have derived an expression for the global stress in inhomogeneous complex copolymer liquids. We apply the principle of virtual work to the free energy as defined in the dynamic mean-field density functional method. This method automatically provides the full stress tensor (deviatoric and isotropic parts) and hence an equation of state for inhomogeneous compressible copolymer melts. The excluded volume interactions and cohesive interactions between chains have been explicitly taken into account. Therefore the expressions for the stress and thermodynamic pressure have a wide range of validity. The connectivity of the chains is automatically accounted for and the fret. energy adapts very well to changes in the molecule properties. In the limiting case of homogeneous systems it simplifies to known results. In order to study rheological properties of copolymer melts and npT-ensemble simulations, the pressure and stress components have to be calculated at any given moment in time. We show how the pressure and stress can be numerically evaluated during simulations using a Green propagator algorithm, instead of having to calculate the time dependent configuration distribution function explicitly from a Smoluchowski equation. We provide illustrative numerical results that indicate how the pressure changes during microphase separation.