I'm working with a special case of Ichino's triple product formula, which for classical holomorphic newforms $f$, $g$ ,$h$ of weights $k$, $m-k$, $m$ (and central characters satisfying $\chi_f \chi_g = \chi_h$) gives me an identity of the form $$ |\langle fg,h \rangle|^2 = C \cdot L(f\times g\times \bar{h},m-1) $$ (where $m-1$ is the central value for this L-function). This is great in principle, but if you're familiar with Ichino's formula you probably know it requires a bunch of work to actually nail down the constant $C$ correctly (there's lots of normalizations you have to get right, and then you have to deal with some unpleasant local integrals). So I want to be able to computationally check the formula I have in a bunch of cases, to convince myself it's actually right. I'm particularly interested in cases where $g$ and $h$ are CM newforms from Hecke characters of an imaginary quadratic field, and the triple-product $L$-value splits up as a product of two Rankin-Selberg $L$-values (twists of $f$ by Hecke characters).
So, my question is: what's the best way to work with these sort of quantities computationally? I haven't used computer algebra systems for number theory of this sort, and I'm not entirely sure where to start. Googling around I found that Pari and Sage both have algorithms that should (?) be able to compute my $L$-value, and another MO question giving some ways to compute Petersson inner products numerically.
My thought is just to boot up Sage and start messing around with the things suggested in those links, but I figured I'd ask here first to see if those are actually the right tools for the job, or if I'm making things way harder on myself than I need to somehow. (I don't really need a ton of precision - just enough to convince me that my constant isn't off by powers of 2 or by Euler factors or something - but I do want to be able to try it for a bunch of modular forms).