function H = structKernelDiag(Psi1, Psi2)

    H = Psi1' * Psi2;

end
