import numpy as np

# The logProd function takes a vector of numbers in logspace 
# (i.e., x[i] = log p[i]) and returns the product of those numbers in logspace
# (i.e., logProd(x) = log(product_i p[i]))
from numpy import linalg as LA


c = np.array([[ 1, 2, 3],[-1, 1, 4]])
print c


