SHOGUN  4.1.0
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Modules Pages
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 
12 #ifdef HAVE_EIGEN3
14 
15 #if EIGEN_VERSION_AT_LEAST(3,1,0)
16 #include <unsupported/Eigen/MatrixFunctions>
17 #endif // EIGEN_VERSION_AT_LEAST(3,1,0)
18 
19 #include <shogun/lib/SGVector.h>
20 #include <shogun/lib/SGMatrix.h>
26 
27 using namespace Eigen;
28 
29 namespace shogun
30 {
31 
32 CDenseMatrixExactLog::CDenseMatrixExactLog()
33  : COperatorFunction<float64_t>(NULL, NULL, OF_LOG)
34 {
35  SG_GCDEBUG("%s created (%p)\n", this->get_name(), this)
36 }
37 
42 {
43  SG_GCDEBUG("%s created (%p)\n", this->get_name(), this)
44 }
45 
47 {
48  SG_GCDEBUG("%s destroyed (%p)\n", this->get_name(), this)
49 }
50 
51 #if EIGEN_VERSION_AT_LEAST(3,1,0)
53 {
54  SG_DEBUG("Entering...\n");
55 
56  // check for proper downcast
59  REQUIRE(op, "Operator not an instance of DenseMatrixOperator!\n");
61 
62  // compute log(C) using Eigen3
65  Map<MatrixXd> log_mat(log_m.matrix, log_m.num_rows, log_m.num_cols);
66  log_mat=mat.log();
67 
68  // the log(C) is also a linear operator here
69  // reset the operator of this function with log(C)
73 
74  SG_DEBUG("Leaving...\n");
75 }
76 #else
78 {
79  SG_WARNING("Eigen3.1.0 or later required!\n")
80 }
81 #endif // EIGEN_VERSION_AT_LEAST(3,1,0)
82 
84  sample)
85 {
86  SG_DEBUG("Entering...\n");
87 
89  // we don't want the aggregator to be destroyed when the job is unref-ed
90  SG_REF(agg);
92  dynamic_cast<CDenseMatrixOperator<float64_t>*>(m_linear_operator), sample);
93  SG_REF(job);
94  // sanity check
95  REQUIRE(m_computation_engine, "Computation engine is NULL\n");
97  // we can safely unref the job here, computation engine takes it from here
98  SG_UNREF(job);
99 
100  SG_DEBUG("Leaving...\n");
101  return agg;
102 }
103 
104 }
105 #endif // HAVE_EIGEN3
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:378
#define SG_REF(x)
Definition: SGObject.h:51
index_t num_rows
Definition: SGMatrix.h:376
#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:52
#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