SHOGUN
3.2.1
首页
相关页面
模块
类
文件
文件列表
文件成员
全部
类
命名空间
文件
函数
变量
类型定义
枚举
枚举值
友元
宏定义
组
页
src
shogun
kernel
HistogramIntersectionKernel.cpp
浏览该文件的文档.
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) 2010 Koen van de Sande
8
* Copyright (C) 2010 Koen van de Sande / University of Amsterdam
9
*/
10
11
#include <
shogun/lib/common.h
>
12
#include <
shogun/kernel/HistogramIntersectionKernel.h
>
13
#include <
shogun/features/Features.h
>
14
#include <
shogun/io/SGIO.h
>
15
16
using namespace
shogun;
17
18
CHistogramIntersectionKernel::CHistogramIntersectionKernel
()
19
:
CDotKernel
(0), m_beta(1.0)
20
{
21
register_params
();
22
}
23
24
CHistogramIntersectionKernel::CHistogramIntersectionKernel
(int32_t size)
25
:
CDotKernel
(size), m_beta(1.0)
26
{
27
register_params
();
28
}
29
30
CHistogramIntersectionKernel::CHistogramIntersectionKernel
(
31
CDenseFeatures<float64_t>
* l,
CDenseFeatures<float64_t>
* r,
32
float64_t
beta, int32_t size)
33
:
CDotKernel
(size), m_beta(beta)
34
{
35
init
(l,r);
36
register_params
();
37
}
38
39
CHistogramIntersectionKernel::~CHistogramIntersectionKernel
()
40
{
41
cleanup
();
42
}
43
44
bool
CHistogramIntersectionKernel::init(
CFeatures
* l,
CFeatures
* r)
45
{
46
bool
result=CDotKernel::init(l,r);
47
init_normalizer
();
48
return
result;
49
}
50
51
float64_t
CHistogramIntersectionKernel::compute
(int32_t idx_a, int32_t idx_b)
52
{
53
int32_t alen, blen;
54
bool
afree, bfree;
55
56
float64_t
* avec=
57
((
CDenseFeatures<float64_t>
*)
lhs
)->get_feature_vector(idx_a, alen, afree);
58
float64_t
* bvec=
59
((
CDenseFeatures<float64_t>
*)
rhs
)->get_feature_vector(idx_b, blen, bfree);
60
ASSERT
(alen==blen)
61
62
float64_t
result=0;
63
64
// checking if beta is default or not
65
if
(
m_beta
== 1.0)
66
{
67
// compute standard histogram intersection kernel
68
for
(int32_t i=0; i<alen; i++)
69
result += (avec[i] < bvec[i]) ? avec[i] : bvec[i];
70
}
71
else
72
{
73
//compute generalized histogram intersection kernel
74
for
(int32_t i=0; i<alen; i++)
75
result +=
CMath::min
(
CMath::pow
(avec[i],
m_beta
),
CMath::pow
(bvec[i],
m_beta
));
76
}
77
((
CDenseFeatures<float64_t>
*)
lhs
)->free_feature_vector(avec, idx_a, afree);
78
((
CDenseFeatures<float64_t>
*)
rhs
)->free_feature_vector(bvec, idx_b, bfree);
79
80
return
result;
81
}
82
83
void
CHistogramIntersectionKernel::register_params
()
84
{
85
SG_ADD
(&
m_beta
,
"beta"
,
"the beta parameter of the kernel"
,
MS_AVAILABLE
);
86
}
SHOGUN
机器学习工具包 - 项目文档