This is most likely a paradox[1].
We assume that the fundamental reason is that they both treat the boundary conditions BC, the source term S and the thermal conductivity/diffusion k and D as scalar quantities when in fact they are tensors or matrices of second rank having spatial components x, y, z. [2]
for example, S=Sx+Sy+Sz and BC=BCx+BCy+BCz.
1- Spiros Konstantogiannis, Eur J 0f Ph, From the state space of quantum mechanics to position and momentum spaces through a simple relation, September 2020.
2- Abbas, Researchgate, A statistical solution to the Schrödinger equation, IJISRT journal, November 2023.