35 #ifndef _BLAZE_MATH_EXPRESSIONS_SVECSVECCROSSEXPR_H_ 36 #define _BLAZE_MATH_EXPRESSIONS_SVECSVECCROSSEXPR_H_ 87 template<
typename VT1
90 class SVecSVecCrossExpr
91 :
public CrossExpr< DenseVector< SVecSVecCrossExpr<VT1,VT2,TF>, TF > >
113 static constexpr
bool returnExpr = ( !IsTemporary_v<RN1> && !IsTemporary_v<RN2> );
116 using ExprReturnType = decltype( ( std::declval<RN1>() * std::declval<RN2>() ) -
117 ( std::declval<RN1>() * std::declval<RN2>() ) );
181 else if( index == 1UL )
199 return (*
this)[index];
208 inline constexpr
size_t size() const noexcept {
239 template<
typename T >
240 inline bool canAlias(
const T* alias )
const noexcept {
241 return (
lhs_.canAlias( alias ) ||
rhs_.canAlias( alias ) );
251 template<
typename T >
252 inline bool isAliased(
const T* alias )
const noexcept {
253 return (
lhs_.isAliased( alias ) ||
rhs_.isAliased( alias ) );
275 template<
typename VT >
286 (~lhs)[0] = x[1UL]*y[2UL] - x[2UL]*y[1UL];
287 (~lhs)[1] = x[2UL]*y[0UL] - x[0UL]*y[2UL];
288 (~lhs)[2] = x[0UL]*y[1UL] - x[1UL]*y[0UL];
305 template<
typename VT >
335 template<
typename VT >
336 friend inline void addAssign( DenseVector<VT,TF>& lhs,
const SVecSVecCrossExpr& rhs )
346 (~lhs)[0] += x[1UL]*y[2UL] - x[2UL]*y[1UL];
347 (~lhs)[1] += x[2UL]*y[0UL] - x[0UL]*y[2UL];
348 (~lhs)[2] += x[0UL]*y[1UL] - x[1UL]*y[0UL];
369 template<
typename VT >
370 friend inline void subAssign( DenseVector<VT,TF>& lhs,
const SVecSVecCrossExpr& rhs )
380 (~lhs)[0] -= x[1UL]*y[2UL] - x[2UL]*y[1UL];
381 (~lhs)[1] -= x[2UL]*y[0UL] - x[0UL]*y[2UL];
382 (~lhs)[2] -= x[0UL]*y[1UL] - x[1UL]*y[0UL];
403 template<
typename VT >
404 friend inline void multAssign( DenseVector<VT,TF>& lhs,
const SVecSVecCrossExpr& rhs )
414 (~lhs)[0] *= x[1UL]*y[2UL] - x[2UL]*y[1UL];
415 (~lhs)[1] *= x[2UL]*y[0UL] - x[0UL]*y[2UL];
416 (~lhs)[2] *= x[0UL]*y[1UL] - x[1UL]*y[0UL];
437 template<
typename VT >
438 friend inline void divAssign( DenseVector<VT,TF>& lhs,
const SVecSVecCrossExpr& rhs )
448 (~lhs)[0] /= x[1UL]*y[2UL] - x[2UL]*y[1UL];
449 (~lhs)[1] /= x[2UL]*y[0UL] - x[0UL]*y[2UL];
450 (~lhs)[2] /= x[0UL]*y[1UL] - x[1UL]*y[0UL];
494 template<
typename VT1
497 , DisableIf_t< IsZero_v<VT1> || IsZero_v<VT2> >* =
nullptr >
498 inline const SVecSVecCrossExpr<VT1,VT2,TF>
499 svecsveccross(
const SparseVector<VT1,TF>& lhs,
const SparseVector<VT2,TF>& rhs )
505 return SVecSVecCrossExpr<VT1,VT2,TF>( ~lhs, ~rhs );
524 template<
typename VT1
527 , EnableIf_t< IsZero_v<VT1> || IsZero_v<VT2> >* =
nullptr >
528 inline decltype(
auto)
529 svecsveccross( const SparseVector<VT1,TF>& lhs, const SparseVector<VT2,TF>& rhs )
537 using ReturnType =
const CrossTrait_t< ResultType_t<VT1>, ResultType_t<VT2> >;
573 template<
typename VT1
576 inline decltype(
auto)
581 if( (~lhs).
size() != 3UL || (~rhs).
size() != 3UL ) {
585 return svecsveccross( ~lhs, ~rhs );
615 template<
typename VT1
#define BLAZE_THROW_INVALID_ARGUMENT(MESSAGE)
Macro for the emission of a std::invalid_argument exception.This macro encapsulates the default way o...
Definition: Exception.h:235
Header file for auxiliary alias declarations.
Header file for the UNUSED_PARAMETER function template.
const StaticVector< ET2, 3UL, TF > RT
Composite type of the right-hand side sparse vector expression.
Definition: SVecSVecCrossExpr.h:144
Header file for basic type definitions.
typename If< Condition, T1, T2 >::Type If_t
Auxiliary alias declaration for the If class template.The If_t alias declaration provides a convenien...
Definition: If.h:109
RightOperand rightOperand() const noexcept
Returns the right-hand side sparse vector operand.
Definition: SVecSVecCrossExpr.h:228
typename T::ResultType ResultType_t
Alias declaration for nested ResultType type definitions.The ResultType_t alias declaration provides ...
Definition: Aliases.h:390
Header file for the serial shim.
bool isAliased(const T *alias) const noexcept
Returns whether the expression is aliased with the given address alias.
Definition: SVecSVecCrossExpr.h:252
Efficient implementation of a compressed matrix.The CompressedMatrix class template is the represent...
Definition: CompressedMatrix.h:224
const StaticVector< ET1, 3UL, TF > LT
Composite type of the left-hand side sparse vector expression.
Definition: SVecSVecCrossExpr.h:141
Constraint on the data type.
Header file for the DenseVector base class.
Header file for the Computation base class.
bool canAlias(const T *alias) const noexcept
Returns whether the expression can alias with the given address alias.
Definition: SVecSVecCrossExpr.h:240
ReturnType_t< VT2 > RN2
Return type of the right-hand side sparse vector expression.
Definition: SVecSVecCrossExpr.h:99
constexpr void UNUSED_PARAMETER(const Args &...)
Suppression of unused parameter warnings.
Definition: Unused.h:81
LeftOperand lhs_
Left-hand side sparse vector of the cross product expression.
Definition: SVecSVecCrossExpr.h:259
typename T::ReturnType ReturnType_t
Alias declaration for nested ReturnType type definitions.The ReturnType_t alias declaration provides ...
Definition: Aliases.h:410
ElementType_t< VT1 > ET1
Element type of the left-hand side sparse vector expression.
Definition: SVecSVecCrossExpr.h:102
typename T::ElementType ElementType_t
Alias declaration for nested ElementType type definitions.The ElementType_t alias declaration provide...
Definition: Aliases.h:170
Efficient implementation of a fixed-sized vector.The StaticVector class template is the representatio...
Definition: Forward.h:61
Constraint on the transpose flag of vector types.
Header file for the DisableIf class template.
TransposeType_t< ResultType > TransposeType
Transpose type for expression template evaluations.
Definition: SVecSVecCrossExpr.h:125
Header file for the IsTemporary type trait class.
Namespace of the Blaze C++ math library.
Definition: Blaze.h:58
Header file for the If class template.
#define BLAZE_CONSTRAINT_MUST_BE_ZERO_TYPE(T)
Constraint on the data type.In case the given data type T is not a zero vector or matrix type...
Definition: Zero.h:61
#define BLAZE_THROW_OUT_OF_RANGE(MESSAGE)
Macro for the emission of a std::out_of_range exception.This macro encapsulates the default way of Bl...
Definition: Exception.h:331
ElementType_t< VT2 > ET2
Element type of the right-hand side sparse vector expression.
Definition: SVecSVecCrossExpr.h:103
Base class for N-dimensional dense vectors.The DenseVector class is a base class for all arbitrarily ...
Definition: DenseVector.h:76
LeftOperand leftOperand() const noexcept
Returns the left-hand side sparse vector operand.
Definition: SVecSVecCrossExpr.h:218
static constexpr bool returnExpr
Compilation switch for the selection of the subscript operator return type.
Definition: SVecSVecCrossExpr.h:113
#define BLAZE_CONSTRAINT_MUST_BE_SPARSE_VECTOR_TYPE(T)
Constraint on the data type.In case the given data type T is not a sparse, N-dimensional vector type...
Definition: SparseVector.h:61
Constraint on the data type.
ReturnType operator[](size_t index) const
Subscript operator for the direct access to the vector elements.
Definition: SVecSVecCrossExpr.h:176
CrossTrait_t< RT1, RT2 > ResultType
Result type for expression template evaluations.
Definition: SVecSVecCrossExpr.h:124
Constraint on the data type.
Header file for the exception macros of the math module.
If_t< IsExpression_v< VT1 >, const VT1, const VT1 &> LeftOperand
Composite type of the left-hand side sparse vector expression.
Definition: SVecSVecCrossExpr.h:135
RightOperand rhs_
Right-hand side sparse vector of the cross product expression.
Definition: SVecSVecCrossExpr.h:260
If_t< IsExpression_v< VT2 >, const VT2, const VT2 &> RightOperand
Composite type of the right-hand side sparse vector expression.
Definition: SVecSVecCrossExpr.h:138
static constexpr bool simdEnabled
Compilation switch for the expression template evaluation strategy.
Definition: SVecSVecCrossExpr.h:149
Constraint on the data type.
Header file for all forward declarations for expression class templates.
Constraint on the data type.
constexpr size_t size() const noexcept
Returns the current size/dimension of the vector.
Definition: SVecSVecCrossExpr.h:208
ResultType_t< VT2 > RT2
Result type of the right-hand side sparse vector expression.
Definition: SVecSVecCrossExpr.h:97
Header file for the EnableIf class template.
const If_t< returnExpr, ExprReturnType, ElementType > ReturnType
Return type for expression template evaluations.
Definition: SVecSVecCrossExpr.h:129
typename CrossTrait< T1, T2 >::Type CrossTrait_t
Auxiliary alias declaration for the CrossTrait class template.The CrossTrait_t alias declaration prov...
Definition: CrossTrait.h:165
Header file for the CrossExpr base class.
Header file for all forward declarations for dense vectors and matrices.
typename T::TransposeType TransposeType_t
Alias declaration for nested TransposeType type definitions.The TransposeType_t alias declaration pro...
Definition: Aliases.h:470
Header file for run time assertion macros.
typename T::CompositeType CompositeType_t
Alias declaration for nested CompositeType type definitions.The CompositeType_t alias declaration pro...
Definition: Aliases.h:90
Header file for the cross product trait.
ReturnType_t< VT1 > RN1
Return type of the left-hand side sparse vector expression.
Definition: SVecSVecCrossExpr.h:98
static constexpr bool smpAssignable
Compilation switch for the expression template assignment strategy.
Definition: SVecSVecCrossExpr.h:152
Header file for the IsZero type trait.
#define BLAZE_FUNCTION_TRACE
Function trace macro.This macro can be used to reliably trace function calls. In case function tracin...
Definition: FunctionTrace.h:94
decltype((std::declval< RN1 >() *std::declval< RN2 >()) -(std::declval< RN1 >() *std::declval< RN2 >())) ExprReturnType
Expression return type for the subscript operator.
Definition: SVecSVecCrossExpr.h:117
constexpr size_t size(const Matrix< MT, SO > &matrix) noexcept
Returns the total number of elements of the matrix.
Definition: Matrix.h:530
decltype(auto) serial(const DenseMatrix< MT, SO > &dm)
Forces the serial evaluation of the given dense matrix expression dm.
Definition: DMatSerialExpr.h:808
#define BLAZE_CONSTRAINT_MUST_NOT_REQUIRE_EVALUATION(T)
Constraint on the data type.In case the given data type T requires an intermediate evaluation within ...
Definition: RequiresEvaluation.h:81
const Type & ReturnType
Return type for expression template evaluations.
Definition: CompressedMatrix.h:3081
#define BLAZE_CONSTRAINT_MUST_BE_DENSE_VECTOR_TYPE(T)
Constraint on the data type.In case the given data type T is not a dense, N-dimensional vector type...
Definition: DenseVector.h:61
const ResultType CompositeType
Data type for composite expression templates.
Definition: SVecSVecCrossExpr.h:132
ElementType_t< ResultType > ElementType
Resulting element type.
Definition: SVecSVecCrossExpr.h:126
ReturnType at(size_t index) const
Checked access to the vector elements.
Definition: SVecSVecCrossExpr.h:195
Header file for the IsComputation type trait class.
Base class for sparse vectors.The SparseVector class is a base class for all arbitrarily sized (N-dim...
Definition: Forward.h:138
SVecSVecCrossExpr(const VT1 &lhs, const VT2 &rhs) noexcept
Constructor for the SVecSVecCrossExpr class.
Definition: SVecSVecCrossExpr.h:161
CompositeType_t< VT1 > CT1
Composite type of the left-hand side sparse vector expression.
Definition: SVecSVecCrossExpr.h:100
ResultType_t< VT1 > RT1
Result type of the left-hand side sparse vector expression.
Definition: SVecSVecCrossExpr.h:96
#define BLAZE_CONSTRAINT_MUST_NOT_BE_ZERO_TYPE(T)
Constraint on the data type.In case the given data type T is a zero vector or matrix type...
Definition: Zero.h:81
#define BLAZE_CONSTRAINT_MUST_BE_VECTOR_WITH_TRANSPOSE_FLAG(T, TF)
Constraint on the data type.In case the given data type T is not a dense or sparse vector type and in...
Definition: TransposeFlag.h:63
#define BLAZE_INTERNAL_ASSERT(expr, msg)
Run time assertion macro for internal checks.In case of an invalid run time expression, the program execution is terminated. The BLAZE_INTERNAL_ASSERT macro can be disabled by setting the BLAZE_USER_ASSERTION flag to zero or by defining NDEBUG during the compilation.
Definition: Assert.h:101
Expression object for sparse vector-sparse vector cross products.The SVecSVecCrossExpr class represen...
Definition: Forward.h:151
CompositeType_t< VT2 > CT2
Composite type of the right-hand side sparse vector expression.
Definition: SVecSVecCrossExpr.h:101
Header file for the IsExpression type trait class.
Header file for the function trace functionality.
decltype(auto) cross(const DenseVector< VT1, TF > &lhs, const DenseVector< VT2, TF > &rhs)
Cross product of two dense vectors ( ).
Definition: DVecDVecCrossExpr.h:544