All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
circulant_matrix_operations.hpp
Go to the documentation of this file.
#ifndef VIENNACL_LINALG_CIRCULANT_MATRIX_OPERATIONS_HPP_
#define VIENNACL_LINALG_CIRCULANT_MATRIX_OPERATIONS_HPP_
/* =========================================================================
Copyright (c) 2010-2016, Institute for Microelectronics,
Institute for Analysis and Scientific Computing,
TU Wien.
Portions of this software are copyright by UChicago Argonne, LLC.
-----------------
ViennaCL - The Vienna Computing Library
-----------------
Project Head: Karl Rupp rupp@iue.tuwien.ac.at
(A list of authors and contributors can be found in the manual)
License: MIT (X11), see file LICENSE in the base directory
============================================================================= */
#include "viennacl/forwards.h"
#include "viennacl/ocl/backend.hpp"
#include "viennacl/scalar.hpp"
#include "viennacl/vector.hpp"
#include "viennacl/tools/tools.hpp"
#include "viennacl/fft.hpp"
//#include "viennacl/linalg/kernels/coordinate_matrix_kernels.h"
namespace viennacl
{
namespace linalg
{
// A * x
template<typename NumericT, unsigned int AlignmentV>
void prod_impl(viennacl::circulant_matrix<NumericT, AlignmentV> const & mat,
viennacl::vector_base<NumericT> const & vec,
viennacl::vector_base<NumericT> & result)
{
assert(mat.size1() == result.size() && bool("Dimension mismatch"));
assert(mat.size2() == vec.size() && bool("Dimension mismatch"));
//result.clear();
//std::cout << "prod(circulant_matrix" << ALIGNMENT << ", vector) called with internal_nnz=" << mat.internal_nnz() << std::endl;
viennacl::vector<NumericT> circ(mat.elements().size() * 2);
viennacl::linalg::real_to_complex(mat.elements(), circ, mat.elements().size());
viennacl::vector<NumericT> tmp(vec.size() * 2);
viennacl::vector<NumericT> tmp2(vec.size() * 2);
viennacl::linalg::real_to_complex(vec, tmp, vec.size());
viennacl::linalg::convolve(circ, tmp, tmp2);
viennacl::linalg::complex_to_real(tmp2, result, vec.size());
}
} //namespace linalg
} //namespace viennacl
#endif
viennacl::circulant_matrix::elements
viennacl::vector< NumericT, AlignmentV > & elements()
Returns an internal viennacl::vector, which represents a circulant matrix elements.
Definition: circulant_matrix.hpp:86
tools.hpp
Various little tools used here and there in ViennaCL.
viennacl::linalg::real_to_complex
void real_to_complex(viennacl::vector_base< NumericT > const &in, viennacl::vector_base< NumericT > &out, vcl_size_t size)
Create complex vector from real vector (even elements(2*k) = real part, odd elements(2*k+1) = imagina...
Definition: fft_operations.hpp:390
viennacl::linalg::complex_to_real
void complex_to_real(viennacl::vector_base< NumericT > const &in, viennacl::vector_base< NumericT > &out, vcl_size_t size)
Create real vector from complex vector (even elements(2*k) = real part, odd elements(2*k+1) = imagina...
Definition: fft_operations.hpp:421
forwards.h
This file provides the forward declarations for the main types used within ViennaCL.
viennacl::circulant_matrix::size1
vcl_size_t size1() const
Returns the number of rows of the matrix.
Definition: circulant_matrix.hpp:92
viennacl::vector_base< NumericT >
viennacl::vector< NumericT >
backend.hpp
Implementations of the OpenCL backend, where all contexts are stored in.
vector.hpp
The vector type with operator-overloads and proxy classes is defined here. Linear algebra operations ...
viennacl::vector_base::size
size_type size() const
Returns the length of the vector (cf. std::vector)
Definition: vector_def.hpp:118
fft.hpp
All routines related to the Fast Fourier Transform. Experimental.
viennacl::circulant_matrix
A Circulant matrix class.
Definition: circulant_matrix.hpp:41
convolve
ScalarType convolve(std::vector< ScalarType > &in1, std::vector< ScalarType > &in2, unsigned int, unsigned int, unsigned int)
Definition: fft_1d.cpp:944
viennacl::circulant_matrix::size2
vcl_size_t size2() const
Returns the number of columns of the matrix.
Definition: circulant_matrix.hpp:97
viennacl::linalg::prod_impl
void prod_impl(const matrix_base< NumericT > &mat, const vector_base< NumericT > &vec, vector_base< NumericT > &result)
Carries out matrix-vector multiplication.
Definition: matrix_operations.hpp:438
scalar.hpp
Implementation of the ViennaCL scalar class.