I have a LTV (linear and time-varying) system. So, $h(\tau, t)$ is the "instantaneous impulse response" at time $t$ such that if the input signal is $x(t) = \delta(t - t_0)$ (an impulse at time $t_0$), the output signal should be $y(t) = h(t - t_0, t_0) = h(\tau, t_0)$ where $\tau = t - t_0$ is the delay/lag from the impulse.
This makes sense to me, intuitively. Obviously, if you input an "impulse at $t_0$", the output signal will be the "impulse response at time $t_0$".
This also has the nice property that if the input is some sort of comb (sum of many impulses with arbitrary complex-valued amplitudes), $$ x(t) = \sum_i A_i \delta(t - t_i) $$ then the output signal is $$ y(t) = \sum_i A_i h(\tau, t_i) $$ which also seems very reasonable to me.
However, this does not agree with the treatment given in most literature that deals with LTV systems where: $$ y(t) = \int_{-\infty}^{+\infty} h(\tau, t) x(t - \tau) d\tau $$
As above, if we have an impulse at $t_0$ such that $x(t) = \delta(t - t_0)$ then $$ \begin{align} y(t) &= \int_{-\infty}^{+\infty} h(\tau, t) \delta(t - \tau - t_0) d\tau \\ &= h(t - t_0, t) \\ &= h(\tau, t) \end{align} $$
This does not match what I expect above!
Am I misunderstanding something? I have a nagging feeling that this is somehow just a matter of definition/convention.
However, I want to make sure I am not accidentally making an implicit assumption of some sort that I don't have to make.
Edit:
Okay, if I define the integral differently, I get what I expected: $$ y(t) = \int_{-\infty}^{+\infty} h(\tau, t - \tau) x(t - \tau) d\tau $$
Is this wrong? How does this differ from the usual treatment of time-varying impulse response?
One immediate problem is that this no longer looks like a "convolution".