gensvm_simplex.c

Go to the documentation of this file.

#include "gensvm_simplex.h"

void gensvm_simplex(struct GenModel *model)
{
    long i, j, K = model->K;

    for (i=0; i<K; i++) {
        for (j=0; j<K-1; j++) {
            if (i <= j) {
                matrix_set(model->U, K-1, i, j,
                            -1.0/sqrt(2.0*(j+1)*(j+2)));
            } else if (i == j+1) {
                matrix_set(model->U, K-1, i, j,
                            sqrt((j+1)/(2.0*(j+2))));
            } else {
                matrix_set(model->U, K-1, i, j, 0.0);
            }
        }
    }
}

void gensvm_simplex_diff(struct GenModel *model)
{
    long i, j, l, K = model->K;
    double value;

    // UU is a 2D block matrix, where block i has the differences:
    // U(i, :) - U(j, :) for all j
    for (i=0; i<K; i++) {
        for (j=0; j<K; j++) {
            for (l=0; l<K-1; l++) {
                value = matrix_get(model->U, K-1, i, l);
                value -= matrix_get(model->U, K-1, j, l);
                matrix_set(model->UU, K-1, i*K+j, l, value);
            }
        }
    }
}