def binomial_ll(parameter, sample):
"""
Return the likelihood of the model parameter under sample.
"""
# compile our method body
from qy import (
Value,
type_from_any,
build_engine,
)
q = \
Qy(
return_type = type_from_any(float) ,
argument_types = [type_from_any(float)] * 3,
default_return = Value.from_any(numpy.nan) ,
)
with q.active() as bar:
q.return_(binomial_log_pdf(*q.main_body.argument_values))
engine = build_engine(q.module)
# define the python wrapper
from ctypes import (
CFUNCTYPE,
c_double,
)
prototype = CFUNCTYPE(c_double, c_double, c_double, c_double)
generated = prototype(engine.get_pointer_to_function(q.main))
from qy.support import raise_if_set
def wrapper(parameter, sample):
"""
Return the likelihood of the model parameter under sample.
"""
parameter = numpy.asarray(parameter, Binomial.parameter_dtype)
result = generated(sample, parameter["p"], parameter["n"])
raise_if_set()
return result
# don't let the engine be GCed
wrapper._engine = engine
# replace and compute
global binomial_ll
binomial_ll = wrapper
return binomial_ll(parameter, sample)
def binomial_ll(parameter, sample):
"""
Return the likelihood of the model parameter under sample.
"""
# compile our method body
from qy import (
Value,
type_from_any,
build_engine,
)
q = \
Qy(
return_type = type_from_any(float) ,
argument_types = [type_from_any(float)] * 3,
default_return = Value.from_any(numpy.nan) ,
)
with q.active() as bar:
q.return_(binomial_log_pdf(*q.main_body.argument_values))
engine = build_engine(q.module)
# define the python wrapper
from ctypes import (
CFUNCTYPE,
c_double,
)
prototype = CFUNCTYPE(c_double, c_double, c_double, c_double)
generated = prototype(engine.get_pointer_to_function(q.main))
from qy.support import raise_if_set
def wrapper(parameter, sample):
"""
Return the likelihood of the model parameter under sample.
"""
parameter = numpy.asarray(parameter, Binomial.parameter_dtype)
result = generated(sample, parameter["p"], parameter["n"])
raise_if_set()
return result
# don't let the engine be GCed
wrapper._engine = engine
# replace and compute
global binomial_ll
binomial_ll = wrapper
return binomial_ll(parameter, sample)
def wrapper(*args, **kwargs):
"""
Call a method of the active Qy instance.
"""
return getattr(Qy.get_active(), name)(*args, **kwargs)