1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
|
/**
* @file msvmmaj_init.c
* @author Gertjan van den Burg
* @date January 7, 2014
* @brief Functions for initializing model and data structures
*
* @details
* This file contains functions for initializing a MajModel instance
* and a MajData instance. In addition, default values for these
* structures are defined here (and only here). Functions for allocating
* memory for the model structure and freeing of the model and data structures
* are also included.
*
*/
#include <math.h>
#include "msvmmaj.h"
#include "msvmmaj_init.h"
/**
* @brief Initialize a MajModel structure
*
* @details
* A MajModel structure is initialized and the default value for the
* parameters are set. A pointer to the initialized model is returned.
*
* @returns initialized MajModel
*/
struct MajModel *msvmmaj_init_model()
{
struct MajModel *model = Malloc(struct MajModel, 1);
// set default values
model->p = 1.0;
model->lambda = pow(2, -8.0);
model->epsilon = 1e-6;
model->kappa = 0.0;
model->weight_idx = 1;
model->kerneltype = K_LINEAR;
return model;
}
/**
* @brief Initialize a MajData structure
*
* @details
* A MajData structure is initialized and default values are set.
* A pointer to the initialized data is returned.
*
* @returns initialized MajData
*
*/
struct MajData *msvmmaj_init_data()
{
struct MajData *data = Malloc(struct MajData, 1);
// set default values
data->kerneltype = K_LINEAR;
return data;
}
/**
* @brief Allocate memory for a MajModel
*
* @details
* This function can be used to allocate the memory needed for a MajModel. All
* arrays in the model are specified and initialized to 0.
*
* @param[in] model MajModel to allocate
*
*/
void msvmmaj_allocate_model(struct MajModel *model)
{
long n = model->n;
long m = model->m;
long K = model->K;
model->W = Calloc(double, m*(K-1));
if (model->W == NULL) {
fprintf(stderr, "Failed to allocate memory for W.\n");
exit(1);
}
model->t = Calloc(double, K-1);
if (model->t == NULL) {
fprintf(stderr, "Failed to allocate memory for t.\n");
exit(1);
}
model->V = Calloc(double, (m+1)*(K-1));
if (model->V == NULL) {
fprintf(stderr, "Failed to allocate memory for V.\n");
exit(1);
}
model->Vbar = Calloc(double, (m+1)*(K-1));
if (model->Vbar == NULL) {
fprintf(stderr, "Failed to allocate memory for Vbar.\n");
exit(1);
}
model->U = Calloc(double, K*(K-1));
if (model->U == NULL) {
fprintf(stderr, "Failed to allocate memory for U.\n");
exit(1);
}
model->UU = Calloc(double, n*K*(K-1));
if (model->UU == NULL) {
fprintf(stderr, "Failed to allocate memory for UU.\n");
exit(1);
}
model->Q = Calloc(double, n*K);
if (model->Q == NULL) {
fprintf(stderr, "Failed to allocate memory for Q.\n");
exit(1);
}
model->H = Calloc(double, n*K);
if (model->H == NULL) {
fprintf(stderr, "Failed to allocate memory for H.\n");
exit(1);
}
model->R = Calloc(double, n*K);
if (model->R == NULL) {
fprintf(stderr, "Failed to allocate memory for R.\n");
exit(1);
}
model->rho = Calloc(double, n);
if (model->rho == NULL) {
fprintf(stderr, "Failed to allocate memory for rho.\n");
exit(1);
}
}
/**
* @brief Free allocated MajModel struct
*
* @details
* Simply free a previously allocated MajModel by freeing all its component
* arrays. Note that the model struct itself is also freed here.
*
* @param[in] model MajModel to free
*
*/
void msvmmaj_free_model(struct MajModel *model)
{
free(model->W);
free(model->t);
free(model->V);
free(model->Vbar);
free(model->U);
free(model->UU);
free(model->Q);
free(model->H);
free(model->rho);
free(model->R);
free(model);
}
/**
* @brief Free allocated MajData struct
*
* @details
* Simply free a previously allocated MajData struct by freeing all its
* components. Note that the data struct itself is also freed here.
*
* @param[in] data MajData struct to free
*
*/
void msvmmaj_free_data(struct MajData *data)
{
free(data->Z);
free(data->y);
free(data);
}
|
