A new method for joint dimension reduction of input and output spaces of a function is introduced.Conventional methods focus on reducing either the input or output space, while this coupled approach supports simultaneous reduction of both.The method is suitable for goal-oriented dimension reduction, where input or output quantities of interest are prescribed.Applications include goal-oriented sensor placement and goal-oriented sensitivity analysis, solving combinatorial optimization problems by optimizing gradient-based bounds.