wolffd@0: function H = structKernelDiag(Psi1, Psi2)
wolffd@0: 
wolffd@0:     H = Psi1' * Psi2;
wolffd@0: 
wolffd@0: end