static PyObject *
PyGSL_monte_init(PyGSL_solver *self, PyObject *args)
{
int flag = GSL_EFAILED;
monte_csys * csys;
FUNC_MESS_BEGIN();
assert(PyGSL_solver_check(self));
csys = (monte_csys *)self->c_sys;
assert(csys);
switch(csys->type){
case PyGSL_MONTE_plain:
flag = gsl_monte_plain_init(self->solver);
break;
case PyGSL_MONTE_miser:
flag = gsl_monte_miser_init(self->solver);
break;
case PyGSL_MONTE_vegas:
flag = gsl_monte_vegas_init(self->solver);
break;
default:
DEBUG_MESS(2, "Monte type %d unknown",flag);
PyGSL_ERROR_NULL("Unknown monte type!", GSL_ESANITY);
}
if(PyGSL_ERROR_FLAG(flag) != GSL_SUCCESS){
PyGSL_add_traceback(module, filename, __FUNCTION__, __LINE__);
return NULL;
}
Py_INCREF(Py_None);
FUNC_MESS_END();
return Py_None;
}
gsl_monte_miser_state *
gsl_monte_miser_alloc (size_t dim)
{
gsl_monte_miser_state *s =
(gsl_monte_miser_state *) malloc (sizeof (gsl_monte_miser_state));
if (s == 0)
{
GSL_ERROR_VAL ("failed to allocate space for miser state struct",
GSL_ENOMEM, 0);
}
s->x = (double *) malloc (dim * sizeof (double));
if (s->x == 0)
{
free (s);
GSL_ERROR_VAL ("failed to allocate space for x", GSL_ENOMEM, 0);
}
s->xmid = (double *) malloc (dim * sizeof (double));
if (s->xmid == 0)
{
free (s->x);
free (s);
GSL_ERROR_VAL ("failed to allocate space for xmid", GSL_ENOMEM, 0);
}
s->sigma_l = (double *) malloc (dim * sizeof (double));
if (s->sigma_l == 0)
{
free (s->xmid);
free (s->x);
free (s);
GSL_ERROR_VAL ("failed to allocate space for sigma_l", GSL_ENOMEM, 0);
}
s->sigma_r = (double *) malloc (dim * sizeof (double));
if (s->sigma_r == 0)
{
free (s->sigma_l);
free (s->xmid);
free (s->x);
free (s);
GSL_ERROR_VAL ("failed to allocate space for sigma_r", GSL_ENOMEM, 0);
}
s->fmax_l = (double *) malloc (dim * sizeof (double));
if (s->fmax_l == 0)
{
free (s->sigma_r);
free (s->sigma_l);
free (s->xmid);
free (s->x);
free (s);
GSL_ERROR_VAL ("failed to allocate space for fmax_l", GSL_ENOMEM, 0);
}
s->fmax_r = (double *) malloc (dim * sizeof (double));
if (s->fmax_r == 0)
{
free (s->fmax_l);
free (s->sigma_r);
free (s->sigma_l);
free (s->xmid);
free (s->x);
free (s);
GSL_ERROR_VAL ("failed to allocate space for fmax_r", GSL_ENOMEM, 0);
}
s->fmin_l = (double *) malloc (dim * sizeof (double));
if (s->fmin_l == 0)
{
free (s->fmax_r);
free (s->fmax_l);
free (s->sigma_r);
free (s->sigma_l);
free (s->xmid);
free (s->x);
free (s);
GSL_ERROR_VAL ("failed to allocate space for fmin_l", GSL_ENOMEM, 0);
}
s->fmin_r = (double *) malloc (dim * sizeof (double));
if (s->fmin_r == 0)
{
free (s->fmin_l);
free (s->fmax_r);
free (s->fmax_l);
free (s->sigma_r);
free (s->sigma_l);
free (s->xmid);
free (s->x);
free (s);
GSL_ERROR_VAL ("failed to allocate space for fmin_r", GSL_ENOMEM, 0);
}
s->fsum_l = (double *) malloc (dim * sizeof (double));
if (s->fsum_l == 0)
{
free (s->fmin_r);
free (s->fmin_l);
free (s->fmax_r);
free (s->fmax_l);
free (s->sigma_r);
free (s->sigma_l);
free (s->xmid);
free (s->x);
free (s);
GSL_ERROR_VAL ("failed to allocate space for fsum_l", GSL_ENOMEM, 0);
}
s->fsum_r = (double *) malloc (dim * sizeof (double));
if (s->fsum_r == 0)
{
free (s->fsum_l);
free (s->fmin_r);
free (s->fmin_l);
free (s->fmax_r);
free (s->fmax_l);
free (s->sigma_r);
free (s->sigma_l);
free (s->xmid);
free (s->x);
free (s);
GSL_ERROR_VAL ("failed to allocate space for fsum_r", GSL_ENOMEM, 0);
}
s->fsum2_l = (double *) malloc (dim * sizeof (double));
if (s->fsum2_l == 0)
{
free (s->fsum_r);
free (s->fsum_l);
free (s->fmin_r);
free (s->fmin_l);
free (s->fmax_r);
free (s->fmax_l);
free (s->sigma_r);
free (s->sigma_l);
free (s->xmid);
free (s->x);
free (s);
GSL_ERROR_VAL ("failed to allocate space for fsum2_l", GSL_ENOMEM, 0);
}
s->fsum2_r = (double *) malloc (dim * sizeof (double));
if (s->fsum2_r == 0)
{
free (s->fsum2_l);
free (s->fsum_r);
free (s->fsum_l);
free (s->fmin_r);
free (s->fmin_l);
free (s->fmax_r);
free (s->fmax_l);
free (s->sigma_r);
free (s->sigma_l);
free (s->xmid);
free (s->x);
free (s);
GSL_ERROR_VAL ("failed to allocate space for fsum2_r", GSL_ENOMEM, 0);
}
s->hits_r = (size_t *) malloc (dim * sizeof (size_t));
if (s->hits_r == 0)
{
free (s->fsum2_r);
free (s->fsum2_l);
free (s->fsum_r);
free (s->fsum_l);
free (s->fmin_r);
free (s->fmin_l);
free (s->fmax_r);
free (s->fmax_l);
free (s->sigma_r);
free (s->sigma_l);
free (s->xmid);
free (s->x);
free (s);
GSL_ERROR_VAL ("failed to allocate space for fsum2_r", GSL_ENOMEM, 0);
}
s->hits_l = (size_t *) malloc (dim * sizeof (size_t));
if (s->hits_l == 0)
{
free (s->hits_r);
free (s->fsum2_r);
free (s->fsum2_l);
free (s->fsum_r);
free (s->fsum_l);
free (s->fmin_r);
free (s->fmin_l);
free (s->fmax_r);
free (s->fmax_l);
free (s->sigma_r);
free (s->sigma_l);
free (s->xmid);
free (s->x);
free (s);
GSL_ERROR_VAL ("failed to allocate space for fsum2_r", GSL_ENOMEM, 0);
}
s->dim = dim;
gsl_monte_miser_init (s);
return s;
}
