import continuous_component_model_test as ccmt

hyper_map = dict()
hyper_map["mu"] = 1.0
hyper_map["s"] = 1.0
hyper_map["nu"] = 1.0
hyper_map["r"] = 1.0

ccmt.set_string_double_map(hyper_map)
component_model = ccmt.p_ContinuousComponentModel(hyper_map)
print component_model.calc_marginal_logp()
component_model.insert_element(1.0)
print component_model.calc_marginal_logp()
component_model.remove_element(1.0)
print component_model.calc_marginal_logp()

import continuous_component_model_test as ccmt

hyper_map = dict()
hyper_map["mu"] = 1.0
hyper_map["s"] = 1.0
hyper_map["nu"] = 1.0
hyper_map["r"] = 1.0

ccmt.set_string_double_map(hyper_map)
component_model = ccmt.p_ContinuousComponentModel(hyper_map)
print component_model.calc_marginal_logp()
component_model.insert_element(1.0)
print component_model.calc_marginal_logp()
component_model.remove_element(1.0)
print component_model.calc_marginal_logp()

p_matrix = ccmt.p_matrix(5,5)
print p_matrix

test boost matrix class