diff options
author | Paul Duncan <pabs@pablotron.org> | 2019-02-03 18:36:52 -0500 |
---|---|---|
committer | Paul Duncan <pabs@pablotron.org> | 2019-02-03 18:36:52 -0500 |
commit | 33a722132491ebdd31722f0cada0f81f6b082282 (patch) | |
tree | 389ae93db5d0a1ff8085ea5f38026d5e060deea1 /km-set.c | |
parent | aa74bd04f66217ff4d617924630b13d721578159 (diff) | |
download | kmeans-33a722132491ebdd31722f0cada0f81f6b082282.tar.bz2 kmeans-33a722132491ebdd31722f0cada0f81f6b082282.zip |
cluster init refactoring, fix best sort, add km_score()
Diffstat (limited to 'km-set.c')
-rw-r--r-- | km-set.c | 83 |
1 files changed, 0 insertions, 83 deletions
@@ -265,86 +265,3 @@ km_set_get_row( const size_t num_floats = set->shape.num_floats; return set->floats + i * num_floats; } - -// init a set with num_clusters clusters of shape num_floats by picking -// random cluster centers -bool -km_set_init_rand_clusters( - km_set_t * const cs, - const size_t num_floats, - const size_t num_clusters, - km_rand_t * const rs -) { - // init cluster shape - const km_shape_t shape = { - .num_floats = num_floats, - .num_ints = 1, - }; - - // generate random cluster centers - float floats[num_floats * num_clusters]; - if (!km_rand_fill(rs, num_floats * num_clusters, floats)) { - // return failure - return false; - } - - // FIXME: should probably be heap-allocated - int ints[num_clusters]; - memset(ints, 0, sizeof(ints)); - - // init cluster set - if (!km_set_init(cs, &shape, num_clusters)) { - // return failure - return false; - } - - // add data, return result - return km_set_push(cs, num_clusters, floats, ints); -} - -// init a set with num_clusters clusters of shape num_floats by picking -// random initial points from the set -bool -km_set_init_rand_points( - km_set_t * const cs, - const km_set_t * const set, - const size_t num_clusters, - km_rand_t * const rs -) { - const size_t num_floats = set->shape.num_floats, - stride = sizeof(float) * num_floats; - - // init cluster shape - const km_shape_t shape = { - .num_floats = num_floats, - .num_ints = 1, - }; - - // get random row offsets - size_t rows[num_clusters]; - if (!km_rand_fill_sizes(rs, num_clusters, rows)) { - // return failure - return false; - } - - // generate random cluster centers - float floats[num_floats * num_clusters]; - for (size_t i = 0; i < num_clusters; i++) { - const size_t row_num = rows[i] % set->num_rows; - const float * const row_floats = km_set_get_row(set, row_num); - memcpy(floats + i * num_floats, row_floats, stride); - } - - // FIXME: should probably be heap-allocated - int ints[num_clusters]; - memset(ints, 0, sizeof(ints)); - - // init cluster set - if (!km_set_init(cs, &shape, num_clusters)) { - // return failure - return false; - } - - // add data, return result - return km_set_push(cs, num_clusters, floats, ints); -} |