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)
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);
209 "Invalid block begin index (%d, %d)!\n",
210 block_begin_row, block_begin_col)
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");
248 "Invalid block begin index (%d, %d)!\n", block_begin, block_begin)
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);
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);
283 "Invalid block begin index (%d, %d)!\n", block_begin, block_begin)
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);
299 block_begin, block_begin, block_size, block_size)), no_diag);
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);
329 "Invalid block begin index (%d, %d)!\n",
330 block_begin_row, block_begin_col)
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 EFeatureType get_feature_type() const =0
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 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)
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()
virtual EFeatureClass get_feature_class() const =0
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 int32_t get_num_vectors() const =0
Block< Matrix > block(Matrix matrix, index_t row_begin, index_t col_begin, index_t row_size, index_t col_size)