We present first results for two-baryon correlation functions, computed using $N_f=2$ flavours of O($a$) improved Wilson quarks, with the aim of explaining potential dibaryon bound states, specifically the H-dibaryon. In particular, we use a GEVP to
isolate the groundstate using two-baryon (hyperon-hyperon) correlation functions $big(langle C_{XY}(t)C_{XY}(0) rangle$, where $XY=LambdaLambda, SigmaSigma, NXi, cdotsbig)$, each of which has an overlap with the H-dibaryon. We employ a `blocking algorithm to handle the large number of contractions, which may easily be extended to N-baryon correlation functions. We also comment on its application to the analysis of single baryon masses ($n$, $Lambda$, $Xi$, $cdots$). This study is performed on an isotropic lattice with $m_pi = 460$ MeV, $m_pi L = 4.7$ and $a = 0.063$ fm.
