ViennaCL - The Vienna Computing Library  1.5.2
execute_scalar_dispatcher.hpp
Go to the documentation of this file.
1 #ifndef VIENNACL_SCHEDULER_EXECUTE_SCALAR_DISPATCHER_HPP
2 #define VIENNACL_SCHEDULER_EXECUTE_SCALAR_DISPATCHER_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 
21 
26 #include <assert.h>
27 
28 #include "viennacl/forwards.h"
32 
33 namespace viennacl
34 {
35  namespace scheduler
36  {
37  namespace detail
38  {
40  template <typename ScalarType1>
41  void as(lhs_rhs_element & s1,
42  lhs_rhs_element const & s2, ScalarType1 const & alpha, vcl_size_t len_alpha, bool reciprocal_alpha, bool flip_sign_alpha)
43  {
46  && bool("Arguments are not vector types!"));
47 
48  switch (s1.numeric_type)
49  {
50  case FLOAT_TYPE:
51  assert(s2.numeric_type == FLOAT_TYPE && bool("Vectors do not have the same scalar type"));
53  *s2.vector_float, convert_to_float(alpha), len_alpha, reciprocal_alpha, flip_sign_alpha);
54  break;
55  case DOUBLE_TYPE:
56  assert(s2.numeric_type == DOUBLE_TYPE && bool("Vectors do not have the same scalar type"));
58  *s2.vector_double, convert_to_double(alpha), len_alpha, reciprocal_alpha, flip_sign_alpha);
59  break;
60  default:
61  throw statement_not_supported_exception("Invalid arguments in scheduler when calling av()");
62  }
63  }
64 
66  template <typename ScalarType1, typename ScalarType2>
67  void asbs(lhs_rhs_element & s1,
68  lhs_rhs_element const & s2, ScalarType1 const & alpha, vcl_size_t len_alpha, bool reciprocal_alpha, bool flip_sign_alpha,
69  lhs_rhs_element const & s3, ScalarType2 const & beta, vcl_size_t len_beta, bool reciprocal_beta, bool flip_sign_beta)
70  {
74  && bool("Arguments are not vector types!"));
75 
76  switch (s1.numeric_type)
77  {
78  case FLOAT_TYPE:
79  assert(s2.numeric_type == FLOAT_TYPE && s3.numeric_type == FLOAT_TYPE && bool("Vectors do not have the same scalar type"));
81  *s2.vector_float, convert_to_float(alpha), len_alpha, reciprocal_alpha, flip_sign_alpha,
82  *s3.vector_float, convert_to_float(beta), len_beta, reciprocal_beta, flip_sign_beta);
83  break;
84  case DOUBLE_TYPE:
85  assert(s2.numeric_type == DOUBLE_TYPE && s3.numeric_type == DOUBLE_TYPE && bool("Vectors do not have the same scalar type"));
87  *s2.vector_double, convert_to_double(alpha), len_alpha, reciprocal_alpha, flip_sign_alpha,
88  *s3.vector_double, convert_to_double(beta), len_beta, reciprocal_beta, flip_sign_beta);
89  break;
90  default:
91  throw statement_not_supported_exception("Invalid arguments in scheduler when calling avbv()");
92  }
93  }
94 
96  template <typename ScalarType1, typename ScalarType2>
98  lhs_rhs_element const & s2, ScalarType1 const & alpha, vcl_size_t len_alpha, bool reciprocal_alpha, bool flip_sign_alpha,
99  lhs_rhs_element const & s3, ScalarType2 const & beta, vcl_size_t len_beta, bool reciprocal_beta, bool flip_sign_beta)
100  {
104  && bool("Arguments are not vector types!"));
105 
106  switch (s1.numeric_type)
107  {
108  case FLOAT_TYPE:
109  assert(s2.numeric_type == FLOAT_TYPE && s3.numeric_type == FLOAT_TYPE && bool("Vectors do not have the same scalar type"));
111  *s2.vector_float, convert_to_float(alpha), len_alpha, reciprocal_alpha, flip_sign_alpha,
112  *s3.vector_float, convert_to_float(beta), len_beta, reciprocal_beta, flip_sign_beta);
113  break;
114  case DOUBLE_TYPE:
115  assert(s2.numeric_type == DOUBLE_TYPE && s3.numeric_type == DOUBLE_TYPE && bool("Vectors do not have the same scalar type"));
117  *s2.vector_double, convert_to_double(alpha), len_alpha, reciprocal_alpha, flip_sign_alpha,
118  *s3.vector_double, convert_to_double(beta), len_beta, reciprocal_beta, flip_sign_beta);
119  break;
120  default:
121  throw statement_not_supported_exception("Invalid arguments in scheduler when calling avbv_v()");
122  }
123  }
124 
125 
126  } // namespace detail
127  } // namespace scheduler
128 } // namespace viennacl
129 
130 #endif
131 
statement_node_subtype subtype
Definition: forwards.h:270
std::size_t vcl_size_t
Definition: forwards.h:58
void asbs(lhs_rhs_element &s1, lhs_rhs_element const &s2, ScalarType1 const &alpha, vcl_size_t len_alpha, bool reciprocal_alpha, bool flip_sign_alpha, lhs_rhs_element const &s3, ScalarType2 const &beta, vcl_size_t len_beta, bool reciprocal_beta, bool flip_sign_beta)
Wrapper for viennacl::linalg::avbv(), taking care of the argument unwrapping.
Definition: execute_scalar_dispatcher.hpp:67
Implementations of scalar operations.
Definition: forwards.h:217
Definition: forwards.h:185
void avbv_v(vector_base< T > &vec1, vector_base< T > const &vec2, ScalarType1 const &alpha, vcl_size_t len_alpha, bool reciprocal_alpha, bool flip_sign_alpha, vector_base< T > const &vec3, ScalarType2 const &beta, vcl_size_t len_beta, bool reciprocal_beta, bool flip_sign_beta)
Definition: vector_operations.hpp:115
double convert_to_double(float d)
Definition: execute_util.hpp:88
void as(lhs_rhs_element &s1, lhs_rhs_element const &s2, ScalarType1 const &alpha, vcl_size_t len_alpha, bool reciprocal_alpha, bool flip_sign_alpha)
Wrapper for viennacl::linalg::av(), taking care of the argument unwrapping.
Definition: execute_scalar_dispatcher.hpp:41
This file provides the forward declarations for the main types used within ViennaCL.
A class representing the 'data' for the LHS or RHS operand of the respective node.
Definition: forwards.h:267
void asbs_s(lhs_rhs_element &s1, lhs_rhs_element const &s2, ScalarType1 const &alpha, vcl_size_t len_alpha, bool reciprocal_alpha, bool flip_sign_alpha, lhs_rhs_element const &s3, ScalarType2 const &beta, vcl_size_t len_beta, bool reciprocal_beta, bool flip_sign_beta)
Wrapper for viennacl::linalg::avbv_v(), taking care of the argument unwrapping.
Definition: execute_scalar_dispatcher.hpp:97
Definition: forwards.h:170
viennacl::vector_base< float > * vector_float
Definition: forwards.h:315
statement_node_numeric_type numeric_type
Definition: forwards.h:271
viennacl::vector_base< double > * vector_double
Definition: forwards.h:316
Definition: forwards.h:184
Provides the datastructures for dealing with a single statement such as 'x = y + z;'.
void avbv(vector_base< T > &vec1, vector_base< T > const &vec2, ScalarType1 const &alpha, vcl_size_t len_alpha, bool reciprocal_alpha, bool flip_sign_alpha, vector_base< T > const &vec3, ScalarType2 const &beta, vcl_size_t len_beta, bool reciprocal_beta, bool flip_sign_beta)
Definition: vector_operations.hpp:78
void av(vector_base< T > &vec1, vector_base< T > const &vec2, ScalarType1 const &alpha, vcl_size_t len_alpha, bool reciprocal_alpha, bool flip_sign_alpha)
Definition: vector_operations.hpp:49
float convert_to_float(float f)
Definition: execute_util.hpp:75
statement_node_type_family type_family
Definition: forwards.h:269
Definition: forwards.h:216
Provides various utilities for implementing the execution of statements.
Exception for the case the scheduler is unable to deal with the operation.
Definition: forwards.h:36