ViennaCL - The Vienna Computing Library  1.5.1
toeplitz_matrix_operations.hpp
Go to the documentation of this file.
1 #ifndef VIENNACL_LINALG_TOEPLITZ_MATRIX_OPERATIONS_HPP_
2 #define VIENNACL_LINALG_TOEPLITZ_MATRIX_OPERATIONS_HPP_
3 
4 /* =========================================================================
5  Copyright (c) 2010-2014, Institute for Microelectronics,
6  Institute for Analysis and Scientific Computing,
7  TU Wien.
8  Portions of this software are copyright by UChicago Argonne, LLC.
9 
10  -----------------
11  ViennaCL - The Vienna Computing Library
12  -----------------
13 
14  Project Head: Karl Rupp rupp@iue.tuwien.ac.at
15 
16  (A list of authors and contributors can be found in the PDF manual)
17 
18  License: MIT (X11), see file LICENSE in the base directory
19 ============================================================================= */
20 
25 #include "viennacl/forwards.h"
26 #include "viennacl/ocl/backend.hpp"
27 #include "viennacl/scalar.hpp"
28 #include "viennacl/vector.hpp"
29 #include "viennacl/tools/tools.hpp"
30 #include "viennacl/fft.hpp"
31 
32 namespace viennacl
33 {
34  namespace linalg
35  {
36 
37 
38  // A * x
39 
48  template<class SCALARTYPE, unsigned int ALIGNMENT>
52  {
53  assert(mat.size1() == result.size());
54  assert(mat.size2() == vec.size());
55 
56  viennacl::vector<SCALARTYPE> tmp(vec.size() * 4); tmp.clear();
57  viennacl::vector<SCALARTYPE> tmp2(vec.size() * 4);
58 
59  viennacl::vector<SCALARTYPE> tep(mat.elements().size() * 2);
60  viennacl::detail::fft::real_to_complex(mat.elements(), tep, mat.elements().size());
61 
62 
63 
64  copy(vec, tmp);
65  viennacl::detail::fft::real_to_complex(tmp, tmp2, vec.size() * 2);
66  viennacl::linalg::convolve(tep, tmp2, tmp);
67  viennacl::detail::fft::complex_to_real(tmp, tmp2, vec.size() * 2);
68  copy(tmp2.begin(), tmp2.begin() + static_cast<vcl_ptrdiff_t>(vec.size()), result.begin());
69  }
70 
71  } //namespace linalg
72 
73 
74 
75 } //namespace viennacl
76 
77 
78 #endif
size_type size() const
Returns the length of the vector (cf. std::vector)
Definition: vector.hpp:859
Various little tools used here and there in ViennaCL.
This file provides the forward declarations for the main types used within ViennaCL.
void clear()
Resets all entries to zero. Does not change the size of the vector.
Definition: vector.hpp:885
void copy(std::vector< SCALARTYPE > &cpu_vec, circulant_matrix< SCALARTYPE, ALIGNMENT > &gpu_mat)
Copies a circulant matrix from the std::vector to the OpenCL device (either GPU or multi-core CPU) ...
Definition: circulant_matrix.hpp:150
A vector class representing a linear memory sequence on the GPU. Inspired by boost::numeric::ublas::v...
Definition: forwards.h:208
Implementations of the OpenCL backend, where all contexts are stored in.
iterator begin()
Returns an iterator pointing to the beginning of the vector (STL like)
Definition: vector.hpp:825
vcl_size_t size1() const
Returns the number of rows of the matrix.
Definition: toeplitz_matrix.hpp:93
The vector type with operator-overloads and proxy classes is defined here. Linear algebra operations ...
All routines related to the Fast Fourier Transform. Experimental.
std::ptrdiff_t vcl_ptrdiff_t
Definition: forwards.h:59
A Toeplitz matrix class.
Definition: forwards.h:330
void prod_impl(const matrix_base< NumericT, F > &mat, const vector_base< NumericT > &vec, vector_base< NumericT > &result)
Carries out matrix-vector multiplication.
Definition: matrix_operations.hpp:350
Implementation of the ViennaCL scalar class.
vcl_size_t size2() const
Returns the number of columns of the matrix.
Definition: toeplitz_matrix.hpp:98
viennacl::vector< SCALARTYPE, ALIGNMENT > & elements()
Returns an internal viennacl::vector, which represents a Toeplitz matrix elements.
Definition: toeplitz_matrix.hpp:86