21 #ifdef HAVE_LINALG_LIB
24 using namespace linalg;
25 #endif // HAVE_LINALG_LIB
27 void CCustomKernel::init()
39 "Subset stack of columns", MS_NOT_AVAILABLE);
40 SG_ADD(&m_free_km, "free_km", "Whether kernel matrix should be freed in "
41 "destructor", MS_NOT_AVAILABLE);
42 SG_ADD(&m_is_symmetric, "is_symmetric", "Whether kernel matrix is symmetric",
44 SG_ADD(&kmatrix, "kmatrix", "Kernel matrix.", MS_NOT_AVAILABLE);
45 SG_ADD(&upper_diagonal, "upper_diagonal", "Upper diagonal", MS_NOT_AVAILABLE);
49 :
CKernel(10), kmatrix(), upper_diagonal(false)
77 :
CKernel(10), upper_diagonal(false)
86 :
CKernel(10), upper_diagonal(false)
119 REQUIRE(l,
"CFeatures l should not be NULL\n")
120 REQUIRE(r,
"CFeatures r should not be NULL\n")
124 "Different FeatureClass: l is %d, r is %d\n",
127 "Different FeatureType: l is %d, r is %d\n",
163 #ifdef HAVE_LINALG_LIB
166 index_t block_size,
bool no_diag)
172 SG_INFO(
"Row/col subsets initialized! Falling back to "
173 "CKernel::sum_symmetric_block (slower)!\n");
180 "Invalid
block begin index (%d, %d)!\n", block_begin, block_begin)
182 "Invalid
block size (%d) at starting index (%d, %d)! "
183 "Please use smaller blocks!", block_size, block_begin, block_begin)
184 REQUIRE(block_size>=1, "Invalid
block size (%d)!\n", block_size)
188 return sum_symmetric<Backend::EIGEN3>(
block(
kmatrix, block_begin,
189 block_begin, block_size, block_size), no_diag);
194 index_t block_size_col,
bool no_diag)
200 SG_INFO(
"Row/col subsets initialized! Falling back to "
201 "CKernel::sum_block (slower)!\n");
203 block_size_row, block_size_col, no_diag);
206 REQUIRE(kmatrix.matrix,
"The kernel matrix is not initialized!\n")
207 REQUIRE(block_begin_row>=0 && block_begin_row<kmatrix.num_rows &&
208 block_begin_col>=0 && block_begin_col<kmatrix.num_cols,
209 "Invalid
block begin index (%d, %d)!\n",
210 block_begin_row, block_begin_col)
211 REQUIRE(block_begin_row+block_size_row<=kmatrix.num_rows &&
212 block_begin_col+block_size_col<=kmatrix.num_cols,
213 "Invalid
block size (%d, %d) at starting index (%d, %d)! "
214 "Please use smaller blocks!", block_size_row, block_size_col,
215 block_begin_row, block_begin_col)
216 REQUIRE(block_size_row>=1 && block_size_col>=1,
217 "Invalid
block size (%d, %d)!\n", block_size_row, block_size_col)
220 if (no_diag && block_size_row!=block_size_col)
222 SG_WARNING(
"Not removing the main diagonal since block is not square!\n");
228 return sum<Backend::EIGEN3>(
block(kmatrix, block_begin_row, block_begin_col,
229 block_size_row, block_size_col), no_diag);
233 block_begin,
index_t block_size,
bool no_diag)
239 SG_INFO(
"Row/col subsets initialized! Falling back to "
240 "CKernel::row_wise_sum_symmetric_block (slower)!\n");
245 REQUIRE(kmatrix.matrix,
"The kernel matrix is not initialized!\n")
246 REQUIRE(m_is_symmetric, "The
kernel matrix is not symmetric!\n")
247 REQUIRE(block_begin>=0 && block_begin<kmatrix.num_cols,
248 "Invalid
block begin index (%d, %d)!\n", block_begin, block_begin)
249 REQUIRE(block_begin+block_size<=kmatrix.num_cols,
250 "Invalid
block size (%d) at starting index (%d, %d)! "
251 "Please use smaller blocks!", block_size, block_begin, block_begin)
252 REQUIRE(block_size>=1, "Invalid
block size (%d)!\n", block_size)
255 block_begin, block_size, block_size), no_diag);
259 for (
index_t i=0; i<s.vlen; ++i)
274 SG_INFO(
"Row/col subsets initialized! Falling back to "
275 "CKernel::row_wise_sum_squared_sum_symmetric_block (slower)!\n");
277 block_size, no_diag);
280 REQUIRE(kmatrix.matrix,
"The kernel matrix is not initialized!\n")
281 REQUIRE(m_is_symmetric, "The
kernel matrix is not symmetric!\n")
282 REQUIRE(block_begin>=0 && block_begin<kmatrix.num_cols,
283 "Invalid
block begin index (%d, %d)!\n", block_begin, block_begin)
284 REQUIRE(block_begin+block_size<=kmatrix.num_cols,
285 "Invalid
block size (%d) at starting index (%d, %d)! "
286 "Please use smaller blocks!", block_size, block_begin, block_begin)
287 REQUIRE(block_size>=1, "Invalid
block size (%d)!\n", block_size)
295 block_begin, block_begin, block_size, block_size), no_diag);
298 elementwise_square<Backend::EIGEN3>(
block(kmatrix,
299 block_begin, block_begin, block_size, block_size)), no_diag);
301 for (
index_t i=0; i<sum.vlen; ++i)
302 row_sum(i, 0)=sum[i];
304 for (
index_t i=0; i<sq_sum.vlen; ++i)
305 row_sum(i, 1)=sq_sum[i];
314 index_t block_size_col,
bool no_diag)
320 SG_INFO(
"Row/col subsets initialized! Falling back to "
321 "CKernel::row_col_wise_sum_block (slower)!\n");
323 block_size_row, block_size_col, no_diag);
326 REQUIRE(kmatrix.matrix,
"The kernel matrix is not initialized!\n")
327 REQUIRE(block_begin_row>=0 && block_begin_row<kmatrix.num_rows &&
328 block_begin_col>=0 && block_begin_col<kmatrix.num_cols,
329 "Invalid
block begin index (%d, %d)!\n",
330 block_begin_row, block_begin_col)
331 REQUIRE(block_begin_row+block_size_row<=kmatrix.num_rows &&
332 block_begin_col+block_size_col<=kmatrix.num_cols,
333 "Invalid
block size (%d, %d) at starting index (%d, %d)! "
334 "Please use smaller blocks!", block_size_row, block_size_col,
335 block_begin_row, block_begin_col)
336 REQUIRE(block_size_row>=1 && block_size_col>=1,
337 "Invalid
block size (%d, %d)!\n", block_size_row, block_size_col)
340 if (no_diag && block_size_row!=block_size_col)
342 SG_WARNING(
"Not removing the main diagonal since block is not square!\n");
352 block_begin_row, block_begin_col, block_size_row,
353 block_size_col), no_diag);
356 block_begin_row, block_begin_col, block_size_row,
357 block_size_col), no_diag);
362 for (
index_t i=0; i<colwise.vlen; ++i)
363 sum[i+rowwise.
vlen]=colwise[i];
369 #endif // HAVE_EIGEN3
370 #endif // HAVE_LINALG_LIB
virtual void add_row_subset(SGVector< index_t > subset)
SGMatrix< float32_t > kmatrix
int32_t num_rhs
number of feature vectors on right hand side
The Custom Kernel allows for custom user provided kernel matrices.
The class IndexFeatures implements features that contain the index of the features. This features used in the CCustomKernel::init to make the subset of the kernel matrix. Initial CIndexFeature of row_idx and col_idx, pass them to the CCustomKernel::init(row_idx, col_idx), then use CCustomKernel::get_kernel_matrix() will get the sub kernel matrix specified by the row_idx and col_idx.
SGMatrix< float32_t > get_float32_kernel_matrix()
virtual float64_t sum_block(index_t block_begin_row, index_t block_begin_col, index_t block_size_row, index_t block_size_col, bool no_diag=false)
bool get_lhs_equals_rhs()
CSubsetStack * m_col_subset_stack
virtual int32_t get_num_vectors() const =0
float64_t kernel(int32_t idx_a, int32_t idx_b)
virtual void remove_all_row_subsets()
SGMatrix< float64_t > get_kernel_matrix()
class to add subset support to another class. A CSubsetStackStack instance should be added and wrappe...
virtual SGMatrix< float64_t > row_wise_sum_squared_sum_symmetric_block(index_t block_begin, index_t block_size, bool no_diag=true)
virtual void remove_col_subset()
SGVector< index_t > get_feature_index()
virtual void remove_row_subset()
virtual void add_col_subset(SGVector< index_t > subset)
virtual void add_subset(SGVector< index_t > subset)
Class SGObject is the base class of all shogun objects.
virtual void row_subset_changed_post()
virtual void remove_all_subsets()
virtual SGVector< float64_t > row_col_wise_sum_block(index_t block_begin_row, index_t block_begin_col, index_t block_size_row, index_t block_size_col, bool no_diag=false)
virtual void add_subset_in_place(SGVector< index_t > subset)
virtual float64_t sum_symmetric_block(index_t block_begin, index_t block_size, bool no_diag=true)
bool set_full_kernel_matrix_from_full(SGMatrix< float32_t > full_kernel_matrix, bool check_symmetry=false)
virtual SGVector< float64_t > row_wise_sum_symmetric_block(index_t block_begin, index_t block_size, bool no_diag=true)
virtual EFeatureClass get_feature_class() const =0
int32_t num_lhs
number of feature vectors on left hand side
virtual bool dummy_init(int32_t rows, int32_t cols)
virtual void remove_all_col_subsets()
virtual bool init_normalizer()
CFeatures * rhs
feature vectors to occur on right hand side
The class DummyFeatures implements features that only know the number of feature objects (but don't a...
all of classes and functions are contained in the shogun namespace
virtual EKernelType get_kernel_type()=0
CFeatures * lhs
feature vectors to occur on left hand side
The class Features is the base class of all feature objects.
virtual void col_subset_changed_post()
virtual bool has_subsets() const
virtual void remove_subset()
virtual void add_col_subset_in_place(SGVector< index_t > subset)
CSubsetStack * m_row_subset_stack
virtual void add_row_subset_in_place(SGVector< index_t > subset)
virtual EFeatureType get_feature_type() const =0
Block< Matrix > block(Matrix matrix, index_t row_begin, index_t col_begin, index_t row_size, index_t col_size)