Seol mar théacs é seo: Computing local p-adic height pairings on hyperelliptic curves