The paper describes the slow evolution of two adaptive traits that regulate the interactions between two mutualistic populations (e.g. a flowering plant and its insect pollinator). For frozen values of the traits, the two populations can either coexist or go extinct. The values of the traits for which populations extinction is guaranteed are therefore of no interest from an evolutionary point of view. In other words, the evolutionary dynamics must be studied only in a viable subset of trait space, which is bounded due to the physiological cost of extreme trait values. Thus, evolutionary dynamics experience so-called border collision bifurcations, when a system invariant in trait space hits the border of the viable subset. The unfolding of standard and border collision bifurcations with respect to two parameters of biological interest is presented. The algebraic and boundary-value problems characterizing the border collision bifurcations are described together with some details concerning their computation.