  We describe an efficient algorithm to compute a large number of baryon-baryon
interactions from $NN$ to $\Xi\Xi$ by means of HAL QCD method, which lays the
groundwork for the nearly physical point lattice QCD calculation with volume
$(96a)^4\approx$($8.2$fm)$^4$. Preliminary results of $\Lambda N$ potential
calculated with quark masses corresponding to
($m_{\pi}$,$m_{K}$)$\approx$(146,525)MeV are presented.
