cond_ellipsoid.c

Go to the documentation of this file.

/*
 * This program is free software: you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation, either version 3 of the License, or
 * (at your option) any later version.
 *
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with this program.  If not, see <http://www.gnu.org/licenses/>.
 */
 
#include "cond_ellipsoid.h"
#include "ea.h"
#include "sam.h"
#include "utils.h"
 
#define N_MU 2 
 
static const int MU_TYPE[N_MU] = { SAM_LOG_NORMAL, SAM_LOG_NORMAL };
 
static double
cond_ellipsoid_dist(const struct XCSF *xcsf, const struct Cl *c,
                    const double *x)
{
    const struct CondEllipsoid *cond = c->cond;
    double dist = 0;
    for (int i = 0; i < xcsf->x_dim; ++i) {
        const double d = (x[i] - cond->center[i]) / cond->spread[i];
        dist += d * d;
    }
    return dist;
}
 
void
cond_ellipsoid_init(const struct XCSF *xcsf, struct Cl *c)
{
    struct CondEllipsoid *new = malloc(sizeof(struct CondEllipsoid));
    new->center = malloc(sizeof(double) * xcsf->x_dim);
    new->spread = malloc(sizeof(double) * xcsf->x_dim);
    new->mu = malloc(sizeof(double) * N_MU);
    const double spread_max = fabs(xcsf->cond->max - xcsf->cond->min);
    for (int i = 0; i < xcsf->x_dim; ++i) {
        new->center[i] = rand_uniform(xcsf->cond->min, xcsf->cond->max);
        new->spread[i] = rand_uniform(xcsf->cond->spread_min, spread_max);
    }
    if (xcsf->cond->sam) { // self-adaptive mutation
        sam_init(new->mu, N_MU, MU_TYPE);
    } else { // constant mutation rate
        new->mu[0] = xcsf->cond->p_mu;
        new->mu[1] = xcsf->cond->mu;
    }
    c->cond = new;
}
 
void
cond_ellipsoid_free(const struct XCSF *xcsf, const struct Cl *c)
{
    (void) xcsf;
    const struct CondEllipsoid *cond = c->cond;
    free(cond->center);
    free(cond->spread);
    free(cond->mu);
    free(c->cond);
}
 
void
cond_ellipsoid_copy(const struct XCSF *xcsf, struct Cl *dest,
                    const struct Cl *src)
{
    struct CondEllipsoid *new = malloc(sizeof(struct CondEllipsoid));
    const struct CondEllipsoid *src_cond = src->cond;
    new->center = malloc(sizeof(double) * xcsf->x_dim);
    new->spread = malloc(sizeof(double) * xcsf->x_dim);
    new->mu = malloc(sizeof(double) * N_MU);
    memcpy(new->center, src_cond->center, sizeof(double) * xcsf->x_dim);
    memcpy(new->spread, src_cond->spread, sizeof(double) * xcsf->x_dim);
    memcpy(new->mu, src_cond->mu, sizeof(double) * N_MU);
    dest->cond = new;
}
 
void
cond_ellipsoid_cover(const struct XCSF *xcsf, const struct Cl *c,
                     const double *x)
{
    const struct CondEllipsoid *cond = c->cond;
    const double spread_max = fabs(xcsf->cond->max - xcsf->cond->min);
    for (int i = 0; i < xcsf->x_dim; ++i) {
        cond->center[i] = x[i];
        cond->spread[i] = rand_uniform(xcsf->cond->spread_min, spread_max);
    }
}
 
void
cond_ellipsoid_update(const struct XCSF *xcsf, const struct Cl *c,
                      const double *x, const double *y)
{
    (void) y;
    if (xcsf->cond->eta > 0) {
        const struct CondEllipsoid *cond = c->cond;
        for (int i = 0; i < xcsf->x_dim; ++i) {
            cond->center[i] += xcsf->cond->eta * (x[i] - cond->center[i]);
        }
    }
}
 
bool
cond_ellipsoid_match(const struct XCSF *xcsf, const struct Cl *c,
                     const double *x)
{
    return (cond_ellipsoid_dist(xcsf, c, x) < 1);
}
 
bool
cond_ellipsoid_crossover(const struct XCSF *xcsf, const struct Cl *c1,
                         const struct Cl *c2)
{
    const struct CondEllipsoid *cond1 = c1->cond;
    const struct CondEllipsoid *cond2 = c2->cond;
    bool changed = false;
    if (rand_uniform(0, 1) < xcsf->ea->p_crossover) {
        for (int i = 0; i < xcsf->x_dim; ++i) {
            if (rand_uniform(0, 1) < 0.5) {
                const double tmp = cond1->center[i];
                cond1->center[i] = cond2->center[i];
                cond2->center[i] = tmp;
                changed = true;
            }
            if (rand_uniform(0, 1) < 0.5) {
                const double tmp = cond1->spread[i];
                cond1->spread[i] = cond2->spread[i];
                cond2->spread[i] = tmp;
                changed = true;
            }
        }
    }
    return changed;
}
 
bool
cond_ellipsoid_mutate(const struct XCSF *xcsf, const struct Cl *c)
{
    bool changed = false;
    const struct CondEllipsoid *cond = c->cond;
    double *center = cond->center;
    double *spread = cond->spread;
 
    if (xcsf->cond->sam) {
        sam_adapt(cond->mu, N_MU, MU_TYPE);
    }
 
    for (int i = 0; i < xcsf->x_dim; ++i) {
        if (rand_uniform(0, 1) < cond->mu[0]) {
            double orig = center[i];
            center[i] += rand_normal(0, cond->mu[1]);
            center[i] = clamp(center[i], xcsf->cond->min, xcsf->cond->max);
            if (orig != center[i]) {
                changed = true;
            }
 
            orig = spread[i];
            spread[i] += rand_normal(0, cond->mu[1]);
            spread[i] = fmax(DBL_EPSILON, spread[i]);
            if (orig != spread[i]) {
                changed = true;
            }
        }
    }
    return changed;
}
 
bool
cond_ellipsoid_general(const struct XCSF *xcsf, const struct Cl *c1,
                       const struct Cl *c2)
{
    const struct CondEllipsoid *cond1 = c1->cond;
    const struct CondEllipsoid *cond2 = c2->cond;
    double *temp = malloc(sizeof(double) * xcsf->x_dim);
    memcpy(temp, cond2->center, sizeof(double) * xcsf->x_dim);
    for (int i = 0; i < xcsf->x_dim; ++i) {
        if (cond1->center[i] != cond2->center[i] ||
            cond1->spread[i] != cond2->spread[i]) {
            temp[i] += cond2->spread[i];
            if (cond_ellipsoid_dist(xcsf, c1, temp) > 1) {
                free(temp);
                return false;
            }
            temp[i] -= 2 * cond2->spread[i];
            if (cond_ellipsoid_dist(xcsf, c1, temp) > 1) {
                free(temp);
                return false;
            }
            temp[i] = cond2->center[i];
        }
    }
    free(temp);
    return true;
}
 
void
cond_ellipsoid_print(const struct XCSF *xcsf, const struct Cl *c)
{
    char *json_str = cond_ellipsoid_json_export(xcsf, c);
    printf("%s\n", json_str);
    free(json_str);
}
 
double
cond_ellipsoid_size(const struct XCSF *xcsf, const struct Cl *c)
{
    (void) c;
    return xcsf->x_dim;
}
 
size_t
cond_ellipsoid_save(const struct XCSF *xcsf, const struct Cl *c, FILE *fp)
{
    size_t s = 0;
    const struct CondEllipsoid *cond = c->cond;
    s += fwrite(cond->center, sizeof(double), xcsf->x_dim, fp);
    s += fwrite(cond->spread, sizeof(double), xcsf->x_dim, fp);
    s += fwrite(cond->mu, sizeof(double), N_MU, fp);
    return s;
}
 
size_t
cond_ellipsoid_load(const struct XCSF *xcsf, struct Cl *c, FILE *fp)
{
    size_t s = 0;
    struct CondEllipsoid *new = malloc(sizeof(struct CondEllipsoid));
    new->center = malloc(sizeof(double) * xcsf->x_dim);
    new->spread = malloc(sizeof(double) * xcsf->x_dim);
    new->mu = malloc(sizeof(double) * N_MU);
    s += fread(new->center, sizeof(double), xcsf->x_dim, fp);
    s += fread(new->spread, sizeof(double), xcsf->x_dim, fp);
    s += fread(new->mu, sizeof(double), N_MU, fp);
    c->cond = new;
    return s;
}
 
char *
cond_ellipsoid_json_export(const struct XCSF *xcsf, const struct Cl *c)
{
    const struct CondEllipsoid *cond = c->cond;
    cJSON *json = cJSON_CreateObject();
    cJSON_AddStringToObject(json, "type", "hyperellipsoid");
    cJSON *center = cJSON_CreateDoubleArray(cond->center, xcsf->x_dim);
    cJSON *spread = cJSON_CreateDoubleArray(cond->spread, xcsf->x_dim);
    cJSON *mutation = cJSON_CreateDoubleArray(cond->mu, N_MU);
    cJSON_AddItemToObject(json, "center", center);
    cJSON_AddItemToObject(json, "spread", spread);
    cJSON_AddItemToObject(json, "mutation", mutation);
    char *string = cJSON_Print(json);
    cJSON_Delete(json);
    return string;
}
 
void
cond_ellipsoid_json_import(const struct XCSF *xcsf, struct Cl *c,
                           const cJSON *json)
{
    struct CondEllipsoid *cond = c->cond;
    const cJSON *item = cJSON_GetObjectItem(json, "center");
    if (item != NULL && cJSON_IsArray(item)) {
        if (cJSON_GetArraySize(item) == xcsf->x_dim) {
            for (int i = 0; i < xcsf->x_dim; ++i) {
                const cJSON *item_i = cJSON_GetArrayItem(item, i);
                cond->center[i] = item_i->valuedouble;
            }
        } else {
            printf("Import error: center length mismatch\n");
            exit(EXIT_FAILURE);
        }
    }
    item = cJSON_GetObjectItem(json, "spread");
    if (item != NULL && cJSON_IsArray(item)) {
        if (cJSON_GetArraySize(item) == xcsf->x_dim) {
            for (int i = 0; i < xcsf->x_dim; ++i) {
                const cJSON *item_i = cJSON_GetArrayItem(item, i);
                cond->spread[i] = item_i->valuedouble;
            }
        } else {
            printf("Import error: spread length mismatch\n");
            exit(EXIT_FAILURE);
        }
    }
    sam_json_import(cond->mu, N_MU, json);
}