source: trunk/yat/classifier/DataLookup1D.cc @ 831

Last change on this file since 831 was 831, checked in by Peter, 15 years ago

Refs #185.

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date ID
File size: 3.3 KB
Line 
1// $Id$
2
3/*
4  Copyright (C) 2005 Peter Johansson
5  Copyright (C) 2006 Jari Häkkinen, Markus Ringnér, Peter Johansson
6  Copyright (C) 2007 Peter Johansson
7
8  This file is part of the yat library, http://lev.thep.lu.se/trac/yat
9
10  The yat library is free software; you can redistribute it and/or
11  modify it under the terms of the GNU General Public License as
12  published by the Free Software Foundation; either version 2 of the
13  License, or (at your option) any later version.
14
15  The yat library is distributed in the hope that it will be useful,
16  but WITHOUT ANY WARRANTY; without even the implied warranty of
17  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
18  General Public License for more details.
19
20  You should have received a copy of the GNU General Public License
21  along with this program; if not, write to the Free Software
22  Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA
23  02111-1307, USA.
24*/
25
26#include "DataLookup1D.h"
27#include "DataLookup2D.h"
28#include "MatrixLookup.h"
29
30#include "yat/utility/matrix.h"
31#include "yat/utility/vector.h"
32
33#include <cassert>
34#include <iostream>
35#include <iomanip>
36
37namespace theplu {
38namespace yat {
39namespace classifier {
40
41  DataLookup1D::DataLookup1D(const DataLookup2D& m, const size_t i, 
42                             const bool row_vector)
43    : column_vector_(!row_vector), index_(i), matrix_(&m), owner_(false)
44  {
45    assert( !column_vector_ || i<m.columns());
46    assert( column_vector_ || i<m.rows());
47  }
48 
49
50  DataLookup1D::DataLookup1D(const size_t size, const double value)
51    : column_vector_(false), index_(0), owner_(true)
52  {
53    matrix_ = new MatrixLookup(1,size,value);
54  }
55
56
57  DataLookup1D::DataLookup1D(const DataLookup1D& other)
58    : column_vector_(other.column_vector_), index_(other.index_),
59      matrix_(other.matrix_), owner_(false)
60  {
61  }
62
63
64  DataLookup1D::DataLookup1D(const utility::vector& v, 
65                             const std::vector<size_t>& index)
66    : column_vector_(true), index_(0), owner_(true)
67  {
68    utility::matrix* m = new utility::matrix(1,index.size());
69    for (size_t i=0; i<index.size(); ++i){
70      assert(index[i]<v.size());
71      (*m)(0,i)=v(index[i]);
72    }
73    matrix_ = new MatrixLookup(*m, true);
74  }
75
76
77  DataLookup1D::DataLookup1D(const utility::vector& v)
78    : column_vector_(false), index_(0), owner_(true)
79  {
80    utility::matrix* m = new utility::matrix(1,v.size());
81    for (size_t i=0; i<v.size(); ++i){
82      (*m)(0,i)=v(i);
83    }
84    matrix_ = new MatrixLookup(*m, true);
85  }
86
87
88  DataLookup1D::~DataLookup1D()
89  {
90    if (owner_)
91      delete matrix_;
92  }
93
94
95  size_t DataLookup1D::size(void) const
96  {
97    return column_vector_ ? matrix_->rows() : matrix_->columns();
98  }
99
100
101  double DataLookup1D::operator()(const size_t i) const
102  { 
103    assert(i<size());
104    return column_vector_ ? (*matrix_)(i,index_) : (*matrix_)(index_,i);
105  }
106
107
108  double DataLookup1D::operator[](const size_t i) const
109  { 
110    return this->operator()(i);
111  }
112
113
114  double DataLookup1D::operator*(const DataLookup1D& other) const
115  {
116    assert(other.size()==size());
117    double res=0;
118    for (size_t i = 0; i<size(); i++)
119      res += (*this)(i)*other(i);
120    return res;
121  }
122
123
124  std::ostream& operator<<(std::ostream& os, const DataLookup1D& x)
125  {
126    os.setf(std::ios::dec);
127    os.precision(12);
128   
129    for (size_t i=0; i<x.size(); ++i) {
130      os << x(i);
131      if ((i+1)<x.size())
132        os << os.fill();
133    }
134    return os;
135  }
136
137}}} // of namespace classifier, yat, and theplu
Note: See TracBrowser for help on using the repository browser.