def METHOD_NAME(
    ode_information_operator, shift_vector, transition_matrix, approx_strategy
):
    """Construct measurement models for the IVP with Dirac likelihoods."""
    measmod_initial_condition = randprocs.markov.discrete.LTIGaussian.from_linop(
        transition_matrix=transition_matrix, noise_mean=shift_vector
    )
    if approx_strategy is not None:
        ode_information_operator = approx_strategy(ode_information_operator)
    measmod_ode = ode_information_operator.as_transition()
    return measmod_initial_condition, measmod_ode