ViennaCL - The Vienna Computing Library  1.5.2
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
Various little tools used here and there in ViennaCL.
viennacl::vector< SCALARTYPE, ALIGNMENT > & elements()
Returns an internal viennacl::vector, which represents a Toeplitz matrix elements.
Definition: toeplitz_matrix.hpp:86
This file provides the forward declarations for the main types used within ViennaCL.
vcl_size_t size2() const
Returns the number of columns of the matrix.
Definition: toeplitz_matrix.hpp:98
Main namespace in ViennaCL. Holds all the basic types such as vector, matrix, etc. and defines operations upon them.
Definition: cpu_ram.hpp:29
A Toeplitz matrix class.
Definition: forwards.h:330
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
iterator begin()
Returns an iterator pointing to the beginning of the vector (STL like)
Definition: vector.hpp:803
size_type size() const
Returns the length of the vector (cf. std::vector)
Definition: vector.hpp:837
Implementations of the OpenCL backend, where all contexts are stored in.
A vector class representing a linear memory sequence on the GPU. Inspired by boost::numeric::ublas::v...
Definition: forwards.h:208
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
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.
void clear()
Resets all entries to zero. Does not change the size of the vector.
Definition: vector.hpp:863
vcl_size_t size1() const
Returns the number of rows of the matrix.
Definition: toeplitz_matrix.hpp:93