SHOGUN  5.0.0
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Modules
DenseMatrixExactLog.cpp
Go to the documentation of this file.
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) 2013 Soumyajit De
8  */
9 
10 #include <shogun/lib/common.h>
11 
13 
14 #if EIGEN_VERSION_AT_LEAST(3,1,0)
15 #include <unsupported/Eigen/MatrixFunctions>
16 #endif // EIGEN_VERSION_AT_LEAST(3,1,0)
17 
18 #include <shogun/lib/SGVector.h>
19 #include <shogun/lib/SGMatrix.h>
25 
26 using namespace Eigen;
27 
28 namespace shogun
29 {
30 
31 CDenseMatrixExactLog::CDenseMatrixExactLog()
32  : COperatorFunction<float64_t>(NULL, NULL, OF_LOG)
33 {
34  SG_GCDEBUG("%s created (%p)\n", this->get_name(), this)
35 }
36 
41 {
42  SG_GCDEBUG("%s created (%p)\n", this->get_name(), this)
43 }
44 
46 {
47  SG_GCDEBUG("%s destroyed (%p)\n", this->get_name(), this)
48 }
49 
50 #if EIGEN_VERSION_AT_LEAST(3,1,0)
52 {
53  SG_DEBUG("Entering...\n");
54 
55  // check for proper downcast
58  REQUIRE(op, "Operator not an instance of DenseMatrixOperator!\n");
60 
61  // compute log(C) using Eigen3
64  Map<MatrixXd> log_mat(log_m.matrix, log_m.num_rows, log_m.num_cols);
65 #if EIGEN_WITH_LOG_BUG_1229
66  MatrixXd tmp = mat;
67  log_mat=tmp.log();
68 #else
69  log_mat=mat.log();
70 #endif
71 
72  // the log(C) is also a linear operator here
73  // reset the operator of this function with log(C)
77 
78  SG_DEBUG("Leaving...\n");
79 }
80 #else
82 {
83  SG_WARNING("Eigen3.1.0 or later required!\n")
84 }
85 #endif // EIGEN_VERSION_AT_LEAST(3,1,0)
86 
88  sample)
89 {
90  SG_DEBUG("Entering...\n");
91 
93  // we don't want the aggregator to be destroyed when the job is unref-ed
94  SG_REF(agg);
96  dynamic_cast<CDenseMatrixOperator<float64_t>*>(m_linear_operator), sample);
97  SG_REF(job);
98  // sanity check
99  REQUIRE(m_computation_engine, "Computation engine is NULL\n");
101  // we can safely unref the job here, computation engine takes it from here
102  SG_UNREF(job);
103 
104  SG_DEBUG("Leaving...\n");
105  return agg;
106 }
107 
108 }
SGMatrix< T > get_matrix_operator() const
CIndependentComputationEngine * m_computation_engine
Abstract template base class for computing for a linear operator C and a vector s. submit_jobs method creates a bunch of jobs needed to solve for this particular and attaches one unique job aggregator to each of them, then submits them all to the computation engine.
Class that represents the job of applying the log of a CDenseMatrixOperator on a real vector...
Definition: SGMatrix.h:20
#define REQUIRE(x,...)
Definition: SGIO.h:206
index_t num_cols
Definition: SGMatrix.h:376
#define SG_REF(x)
Definition: SGObject.h:54
index_t num_rows
Definition: SGMatrix.h:374
#define SG_GCDEBUG(...)
Definition: SGIO.h:102
Template class that aggregates scalar job results in each submit_result call, finalize then transform...
double float64_t
Definition: common.h:50
Abstract base class that provides an interface for computing an aggeregation of the job results of in...
virtual const char * get_name() const
#define SG_UNREF(x)
Definition: SGObject.h:55
#define SG_DEBUG(...)
Definition: SGIO.h:107
all of classes and functions are contained in the shogun namespace
Definition: class_list.h:18
CLinearOperator< float64_t > * m_linear_operator
Abstract template base class that represents a linear operator, e.g. a matrix.
Abstract base class for solving multiple independent instances of CIndependentJob. It has one method, submit_job, which may add the job to an internal queue and might block if there is yet not space in the queue. After jobs are submitted, it might not yet be ready. wait_for_all waits until all jobs are completed, which *must be* called to guarantee that all jobs are finished.
#define SG_WARNING(...)
Definition: SGIO.h:128
virtual void submit_job(CIndependentJob *job)=0
virtual CJobResultAggregator * submit_jobs(SGVector< float64_t > sample)

SHOGUN Machine Learning Toolbox - Documentation