Optimize wigner function and marginal function for Bosonic state

[Hugh-888]

Update bosonic wigner function and marginal function
example:

sq = 1
r = torch.tensor(2, dtype=torch.double)
r_prime = (np.cosh(sq)+np.sinh(sq))*r/2
cir = dq.QumodeCircuit(nmode=1, init_state='vac', cutoff=3, backend='bosonic')
cir.cat(wires=0, r=r_prime, theta=0, p=0) # plus cat state
cir.s(0, sq)
cir.homodyne_p(0)
cir.to(torch.double)
re = cir()
s = dq.BosonicState(state=re, nmode=1) # squeezed plus cat state

qrange = 10
prange = 15
marginal_q =s.marginal(wire=0, qrange=qrange)
marginal_p = s.marginal(wire=0, qrange=prange, phi=np.pi/2)
wigner = s.wigner(wire=0, qrange=qrange, prange=prange)

[sansiro77]

GJ