Computing local p-adic height pairings on hyperelliptic curves