SHOGUN
3.2.1
首页
相关页面
模块
类
文件
文件列表
文件成员
全部
类
命名空间
文件
函数
变量
类型定义
枚举
枚举值
友元
宏定义
组
页
src
shogun
kernel
normalizer
SqrtDiagKernelNormalizer.h
浏览该文件的文档.
1
/*
2
* This program is free software; you can redistribute it and/or modify
3
* it under the terms of the GNU General Public License as published by
4
* the Free Software Foundation; either version 3 of the License, or
5
* (at your option) any later version.
6
*
7
* Written (W) 2009 Soeren Sonnenburg
8
* Copyright (C) 2009 Fraunhofer Institute FIRST and Max-Planck-Society
9
*/
10
11
#ifndef _SQRTDIAGKERNELNORMALIZER_H___
12
#define _SQRTDIAGKERNELNORMALIZER_H___
13
14
#include <
shogun/kernel/normalizer/KernelNormalizer.h
>
15
#include <
shogun/kernel/string/CommWordStringKernel.h
>
16
17
namespace
shogun
18
{
29
class
CSqrtDiagKernelNormalizer
:
public
CKernelNormalizer
30
{
31
public
:
36
CSqrtDiagKernelNormalizer
(
bool
use_opt_diag=
false
):
CKernelNormalizer
(),
37
sqrtdiag_lhs
(NULL),
num_sqrtdiag_lhs
(0),
38
sqrtdiag_rhs
(NULL),
num_sqrtdiag_rhs
(0),
39
use_optimized_diagonal_computation
(use_opt_diag)
40
{
41
m_parameters
->
add_vector
(&
sqrtdiag_lhs
, &
num_sqrtdiag_lhs
,
"sqrtdiag_lhs"
,
42
"sqrt(K(x,x)) for left hand side examples."
);
43
m_parameters
->
add_vector
(&
sqrtdiag_rhs
, &
num_sqrtdiag_rhs
,
"sqrtdiag_rhs"
,
44
"sqrt(K(x,x)) for right hand side examples."
);
45
SG_ADD
(&
use_optimized_diagonal_computation
,
46
"use_optimized_diagonal_computation"
,
47
"flat if optimized diagonal computation is used"
,
MS_NOT_AVAILABLE
);
48
}
49
51
virtual
~CSqrtDiagKernelNormalizer
()
52
{
53
SG_FREE(
sqrtdiag_lhs
);
54
SG_FREE(
sqrtdiag_rhs
);
55
}
56
59
virtual
bool
init
(
CKernel
* k)
60
{
61
ASSERT
(k)
62
num_sqrtdiag_lhs
=k->
get_num_vec_lhs
();
63
num_sqrtdiag_rhs
=k->
get_num_vec_rhs
();
64
ASSERT
(
num_sqrtdiag_lhs
>0)
65
ASSERT
(
num_sqrtdiag_rhs
>0)
66
67
CFeatures
* old_lhs=k->
lhs
;
68
CFeatures
* old_rhs=k->
rhs
;
69
70
k->
lhs
=old_lhs;
71
k->
rhs
=old_lhs;
72
bool
r1=
alloc_and_compute_diag
(k,
sqrtdiag_lhs
,
num_sqrtdiag_lhs
);
73
74
k->
lhs
=old_rhs;
75
k->
rhs
=old_rhs;
76
bool
r2=
alloc_and_compute_diag
(k,
sqrtdiag_rhs
,
num_sqrtdiag_rhs
);
77
78
k->
lhs
=old_lhs;
79
k->
rhs
=old_rhs;
80
81
return
r1 && r2;
82
}
83
89
virtual
float64_t
normalize
(
90
float64_t
value, int32_t idx_lhs, int32_t idx_rhs)
91
{
92
float64_t
sqrt_both=
sqrtdiag_lhs
[idx_lhs]*
sqrtdiag_rhs
[idx_rhs];
93
return
value/sqrt_both;
94
}
95
100
virtual
float64_t
normalize_lhs
(
float64_t
value, int32_t idx_lhs)
101
{
102
return
value/
sqrtdiag_lhs
[idx_lhs];
103
}
104
109
virtual
float64_t
normalize_rhs
(
float64_t
value, int32_t idx_rhs)
110
{
111
return
value/
sqrtdiag_rhs
[idx_rhs];
112
}
113
114
public
:
119
bool
alloc_and_compute_diag
(
CKernel
* k,
float64_t
* &v, int32_t num)
120
{
121
SG_FREE(v);
122
v=SG_MALLOC(
float64_t
, num);
123
124
for
(int32_t i=0; i<num; i++)
125
{
126
if
(k->
get_kernel_type
() ==
K_COMMWORDSTRING
)
127
{
128
if
(
use_optimized_diagonal_computation
)
129
v[i]=sqrt(((
CCommWordStringKernel
*) k)->compute_diag(i));
130
else
131
v[i]=sqrt(((
CCommWordStringKernel
*) k)->compute_helper(i,i,
true
));
132
}
133
else
134
v[i]=sqrt(k->
compute
(i,i));
135
136
if
(v[i]==0.0)
137
v[i]=1e-16;
/* avoid divide by zero exception */
138
}
139
140
return
(v!=NULL);
141
}
142
144
virtual
const
char
*
get_name
()
const
{
return
"SqrtDiagKernelNormalizer"
; }
145
146
protected
:
148
float64_t
*
sqrtdiag_lhs
;
149
151
int32_t
num_sqrtdiag_lhs
;
152
154
float64_t
*
sqrtdiag_rhs
;
155
157
int32_t
num_sqrtdiag_rhs
;
158
160
bool
use_optimized_diagonal_computation
;
161
};
162
}
163
#endif
SHOGUN
机器学习工具包 - 项目文档