I'm reading about boost::math::interpolators::cubic_hermite and the constructor takes three arguments: x values, y values and slopes (dy/dx). Can this library be used when the slope at a point is vertical? If so, how do you specify that? I tried std::numeric_limits<double>::infinity() but that doesn't work. Cory