def B(value=ones(2, dtype=float), mu=A, tau=tau_B): return mv_normal_like(value, mu, tau)
def B(value = ones(2,dtype=float), mu = A, tau = tau_B): return mv_normal_like(value,mu,tau)
def A(value=ones(2, dtype=float), mu=mu_A, tau=tau_A): return mv_normal_like(value, mu, tau)
def A(value = ones(2,dtype=float), mu=mu_A, tau = tau_A): return mv_normal_like(value,mu,tau)