void
scew_element_list_del(scew_element* element, XML_Char const* name)
{
unsigned int i = 0;
unsigned int count = 0;
scew_element** list = NULL;
if ((element == NULL) || (name == NULL))
{
return;
}
list = scew_element_list(element, name, &count);
if (list == NULL)
{
return;
}
for (i = 0; i < count; i++)
{
scew_element_free(list[i]);
}
scew_element_list_free(list);
}
local_data->nherds_detected_by_prodtype =
RPT_new_reporting ("num-units-detected-by-production-type", NULL, RPT_group, RPT_never, FALSE);
local_data->cumul_nherds_detected =
RPT_new_reporting ("cumulative-num-units-detected", NULL, RPT_integer, RPT_never, TRUE);
local_data->cumul_nherds_detected_by_prodtype =
RPT_new_reporting ("cumulative-num-units-detected-by-production-type", NULL, RPT_group,
RPT_never, TRUE);
g_ptr_array_add (m->outputs, local_data->detections);
g_ptr_array_add (m->outputs, local_data->day_1st_detection);
g_ptr_array_add (m->outputs, local_data->nherds_detected);
g_ptr_array_add (m->outputs, local_data->nherds_detected_by_prodtype);
g_ptr_array_add (m->outputs, local_data->cumul_nherds_detected);
g_ptr_array_add (m->outputs, local_data->cumul_nherds_detected_by_prodtype);
/* Set the reporting frequency for the output variables. */
ee = scew_element_list (params, "output", &noutputs);
#if INFO
g_log (G_LOG_DOMAIN, G_LOG_LEVEL_INFO, "%i output variables", noutputs);
#endif
for (i = 0; i < noutputs; i++)
{
e = ee[i];
variable_name = scew_element_contents (scew_element_by_name (e, "variable-name"));
/* Do the outputs include a variable with this name? */
for (j = 0; j < m->outputs->len; j++)
{
output = (RPT_reporting_t *) g_ptr_array_index (m->outputs, j);
if (strcmp (output->name, variable_name) == 0)
break;
}
if (j == m->outputs->len)
/**
* Adds a set of parameters to a contact spread model.
*/
void
set_params (struct naadsm_model_t_ *self, PAR_parameter_t * params)
{
local_data_t *local_data = (local_data_t *) (self->model_data);
unsigned int nzones = ZON_zone_list_length (local_data->zones);
unsigned int nprod_types = local_data->production_types->len;
scew_element *e, **ee;
gboolean success;
gboolean has_destruction_cost_params = FALSE;
gboolean has_vaccination_cost_params = FALSE;
gboolean has_surveillance_cost_param = FALSE;
destruction_cost_data_t destruction_cost_params = {};
vaccination_cost_data_t vaccination_cost_params = {};
double surveillance_cost_param = 0;
gboolean *production_type = NULL;
gboolean *zone = NULL;
unsigned int i, j;
unsigned int noutputs;
RPT_reporting_t *output;
const XML_Char *variable_name;
#if DEBUG
g_log (G_LOG_DOMAIN, G_LOG_LEVEL_DEBUG, "----- ENTER set_params (%s)", MODEL_NAME);
#endif
/* Make sure the right XML subtree was sent. */
g_assert (strcmp (scew_element_name (params), MODEL_NAME) == 0);
#if DEBUG
g_log (G_LOG_DOMAIN, G_LOG_LEVEL_DEBUG, "setting production types");
#endif
/* Destruction Cost Parameters */
e = scew_element_by_name (params, "appraisal");
if (e != NULL)
{
destruction_cost_params.appraisal = PAR_get_money (e, &success);
if (success)
{
has_destruction_cost_params = TRUE;
}
else
{
g_warning ("%s: setting per-unit appraisal cost to 0", MODEL_NAME);
destruction_cost_params.appraisal = 0;
}
}
else
{
g_warning ("%s: per-unit appraisal cost missing, setting to 0", MODEL_NAME);
destruction_cost_params.appraisal = 0;
}
e = scew_element_by_name (params, "euthanasia");
if (e != NULL)
{
destruction_cost_params.euthanasia = PAR_get_money (e, &success);
if (success)
{
has_destruction_cost_params = TRUE;
}
else
{
g_warning ("%s: setting per-animal euthanasia cost to 0", MODEL_NAME);
destruction_cost_params.euthanasia = 0;
}
}
else
{
g_warning ("%s: per-animal euthanasia cost missing, setting to 0", MODEL_NAME);
destruction_cost_params.euthanasia = 0;
}
e = scew_element_by_name (params, "indemnification");
if (e != NULL)
{
destruction_cost_params.indemnification = PAR_get_money (e, &success);
if (success)
{
has_destruction_cost_params = TRUE;
}
else
{
g_warning ("%s: setting per-animal indemnification cost to 0", MODEL_NAME);
destruction_cost_params.indemnification = 0;
}
}
else
{
g_warning ("%s: per-animal indemnification cost missing, setting to 0", MODEL_NAME);
destruction_cost_params.indemnification = 0;
}
e = scew_element_by_name (params, "carcass-disposal");
if (e != NULL)
{
destruction_cost_params.carcass_disposal = PAR_get_money (e, &success);
if (success)
{
has_destruction_cost_params = TRUE;
}
else
{
g_warning ("%s: setting per-animal carcass disposal cost to 0", MODEL_NAME);
destruction_cost_params.carcass_disposal = 0;
}
}
else
{
g_warning ("%s: per-animal carcass disposal cost missing, setting to 0", MODEL_NAME);
destruction_cost_params.carcass_disposal = 0;
}
e = scew_element_by_name (params, "cleaning-disinfecting");
if (e != NULL)
{
destruction_cost_params.cleaning_disinfecting = PAR_get_money (e, &success);
if (success)
{
has_destruction_cost_params = TRUE;
}
else
{
g_warning ("%s: setting per-unit cleaning and disinfecting cost to 0", MODEL_NAME);
destruction_cost_params.cleaning_disinfecting = 0;
}
}
else
{
g_warning ("%s: per-unit cleaning and disinfecting cost missing, setting to 0", MODEL_NAME);
destruction_cost_params.cleaning_disinfecting = 0;
}
/* Vaccination Cost Parameters */
e = scew_element_by_name (params, "vaccination-fixed");
if (e != NULL)
{
vaccination_cost_params.vaccination_fixed = PAR_get_money (e, &success);
if (success)
{
has_vaccination_cost_params = TRUE;
}
else
{
g_warning ("%s: setting per-unit vaccination cost to 0", MODEL_NAME);
vaccination_cost_params.vaccination_fixed = 0;
}
}
else
{
g_warning ("%s: per-unit vaccination cost missing, setting to 0", MODEL_NAME);
vaccination_cost_params.vaccination_fixed = 0;
}
e = scew_element_by_name (params, "vaccination");
if (e != NULL)
{
vaccination_cost_params.vaccination = PAR_get_money (e, &success);
if (success)
{
has_vaccination_cost_params = TRUE;
}
else
{
g_warning ("%s: setting per-animal vaccination cost to 0", MODEL_NAME);
vaccination_cost_params.vaccination = 0;
}
}
else
{
g_warning ("%s: per-animal vaccination cost missing, setting to 0", MODEL_NAME);
vaccination_cost_params.vaccination = 0;
}
e = scew_element_by_name (params, "baseline-vaccination-capacity");
if (e != NULL)
{
vaccination_cost_params.baseline_capacity = (unsigned int) PAR_get_unitless (e, &success);
if (success)
{
has_vaccination_cost_params = TRUE;
}
else
{
g_warning ("%s: setting baseline vaccination capacity to 1,000,000", MODEL_NAME);
vaccination_cost_params.baseline_capacity = 1000000;
}
}
else
{
g_warning ("%s: baseline vaccination capacity missing, setting to 1,000,000", MODEL_NAME);
vaccination_cost_params.baseline_capacity = 1000000;
}
e = scew_element_by_name (params, "additional-vaccination");
if (e != NULL)
{
vaccination_cost_params.extra_vaccination = PAR_get_money (e, &success);
if (success)
{
has_vaccination_cost_params = TRUE;
}
else
{
g_warning ("%s: setting additional per-animal vaccination cost to 0", MODEL_NAME);
vaccination_cost_params.extra_vaccination = 0;
}
}
else
{
g_warning ("%s: additional per-animal vaccination cost missing, setting to 0", MODEL_NAME);
vaccination_cost_params.extra_vaccination = 0;
}
/* No vaccinations have been performed yet. */
vaccination_cost_params.capacity_used = 0;
/* Surveillance Cost Parameters */
e = scew_element_by_name (params, "surveillance");
if (e != NULL)
{
surveillance_cost_param = PAR_get_money (e, &success);
if (success)
{
has_surveillance_cost_param = TRUE;
}
else
{
g_warning ("%s: setting per-animal zone surveillance cost to 0", MODEL_NAME);
surveillance_cost_param = 0;
}
}
else
{
g_warning ("%s: per-animal zone surveillance cost missing, setting to 0", MODEL_NAME);
surveillance_cost_param = 0;
}
/* Set the reporting frequency for the output variables. */
ee = scew_element_list (params, "output", &noutputs);
#if DEBUG
g_debug ("%i output variables", noutputs);
#endif
for (i = 0; i < noutputs; i++)
{
e = ee[i];
variable_name = scew_element_contents (scew_element_by_name (e, "variable-name"));
/* Do the outputs include a variable with this name? */
for (j = 0; j < self->outputs->len; j++)
{
output = (RPT_reporting_t *) g_ptr_array_index (self->outputs, j);
if (strcmp (output->name, variable_name) == 0)
break;
}
if (j == self->outputs->len)
g_warning ("no output variable named \"%s\", ignoring", variable_name);
else
{
RPT_reporting_set_frequency (output,
RPT_string_to_frequency (scew_element_contents
(scew_element_by_name
(e, "frequency"))));
#if DEBUG
g_debug ("report \"%s\" %s", variable_name, RPT_frequency_name[output->frequency]);
#endif
}
}
free (ee);
/* Read zone and production type attributes to determine which
* entries should be filled in. */
if (scew_attribute_by_name (params, "zone") != NULL)
zone = naadsm_read_zone_attribute (params, local_data->zones);
else
zone = NULL;
production_type =
naadsm_read_prodtype_attribute (params, "production-type", local_data->production_types);
/* Copy the parameters into the selected zone and production types. */
for (i = 0; i < nprod_types; i++)
{
if (production_type[i])
{
if (has_destruction_cost_params)
{
if (local_data->destruction_cost_params == NULL)
{
local_data->destruction_cost_params =
g_new0 (destruction_cost_data_t *, nprod_types);
g_assert (local_data->destruction_cost_params != NULL);
}
if (local_data->destruction_cost_params[i] == NULL)
{
local_data->destruction_cost_params[i] =
g_new0 (destruction_cost_data_t, 1);
g_assert (local_data->destruction_cost_params[i] != NULL);
}
memcpy (local_data->destruction_cost_params[i],
&destruction_cost_params,
sizeof (destruction_cost_data_t));
}
if (has_vaccination_cost_params)
{
if (local_data->vaccination_cost_params == NULL)
{
local_data->vaccination_cost_params =
g_new0 (vaccination_cost_data_t *, nprod_types);
g_assert (local_data->vaccination_cost_params != NULL);
}
if (local_data->vaccination_cost_params[i] == NULL)
{
local_data->vaccination_cost_params[i] =
g_new0 (vaccination_cost_data_t, 1);
g_assert (local_data->vaccination_cost_params[i] != NULL);
}
memcpy (local_data->vaccination_cost_params[i],
&vaccination_cost_params,
sizeof (vaccination_cost_data_t));
}
if (has_surveillance_cost_param)
{
if (zone)
{
for (j = 0; j < nzones; j++)
{
if (zone[j])
{
if (local_data->surveillance_cost_param == NULL)
{
local_data->surveillance_cost_param =
g_new0 (double *, nzones);
g_assert (local_data->surveillance_cost_param != NULL);
}
if (local_data->surveillance_cost_param[j] == NULL)
{
local_data->surveillance_cost_param[j] =
g_new0 (double, nprod_types);
g_assert (local_data->surveillance_cost_param[j] != NULL);
}
local_data->surveillance_cost_param[j][i] = surveillance_cost_param;
}
}
}
