35 #ifndef _BLAZE_MATH_EXPRESSIONS_DVECSVECCROSSEXPR_H_ 36 #define _BLAZE_MATH_EXPRESSIONS_DVECSVECCROSSEXPR_H_ 87 template<
typename VT1
91 :
public CrossExpr< DenseVector< DVecSVecCrossExpr<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 DVecSVecCrossExpr& 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 DVecSVecCrossExpr& 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 DVecSVecCrossExpr& 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 DVecSVecCrossExpr& 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];
493 template<
typename VT1
496 , DisableIf_t< IsZero_v<VT2> >* =
nullptr >
497 inline const DVecSVecCrossExpr<VT1,VT2,TF>
498 dvecsveccross(
const DenseVector<VT1,TF>& lhs,
const SparseVector<VT2,TF>& rhs )
504 return DVecSVecCrossExpr<VT1,VT2,TF>( ~lhs, ~rhs );
523 template<
typename VT1
526 , EnableIf_t< IsZero_v<VT2> >* =
nullptr >
527 inline decltype(
auto)
528 dvecsveccross( const DenseVector<VT1,TF>& lhs, const SparseVector<VT2,TF>& rhs )
536 using ReturnType =
const CrossTrait_t< ResultType_t<VT1>, ResultType_t<VT2> >;
574 template<
typename VT1
577 inline decltype(
auto)
582 if( (~lhs).
size() != 3UL || (~rhs).
size() != 3UL ) {
586 return dvecsveccross( ~lhs, ~rhs );
617 template<
typename VT1
ElementType_t< VT1 > ET1
Element type of the left-hand side dense vector expression.
Definition: DVecSVecCrossExpr.h:102
#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
Expression object for dense vector-sparse vector cross products.The DVecSVecCrossExpr class represent...
Definition: DVecSVecCrossExpr.h:90
Header file for auxiliary alias declarations.
Header file for the UNUSED_PARAMETER function template.
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
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.
Efficient implementation of a compressed matrix.The CompressedMatrix class template is the represent...
Definition: CompressedMatrix.h:224
CompositeType_t< VT1 > CT1
Composite type of the left-hand side dense vector expression.
Definition: DVecSVecCrossExpr.h:100
Base class for all cross product expression templates.The CrossExpr class serves as a tag for all exp...
Definition: CrossExpr.h:66
Constraint on the data type.
Header file for the DenseVector base class.
Header file for the Computation base class.
ReturnType_t< VT1 > RN1
Return type of the left-hand side dense vector expression.
Definition: DVecSVecCrossExpr.h:98
constexpr void UNUSED_PARAMETER(const Args &...)
Suppression of unused parameter warnings.
Definition: Unused.h:81
typename T::ReturnType ReturnType_t
Alias declaration for nested ReturnType type definitions.The ReturnType_t alias declaration provides ...
Definition: Aliases.h:410
static constexpr bool smpAssignable
Compilation switch for the expression template assignment strategy.
Definition: DVecSVecCrossExpr.h:152
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.
decltype((std::declval< RN1 >() *std::declval< RN2 >()) -(std::declval< RN1 >() *std::declval< RN2 >())) ExprReturnType
Expression return type for the subscript operator.
Definition: DVecSVecCrossExpr.h:117
LeftOperand lhs_
Left-hand side dense vector of the cross product expression.
Definition: DVecSVecCrossExpr.h:259
If_t< IsExpression_v< VT2 >, const VT2, const VT2 &> RightOperand
Composite type of the right-hand side sparse vector expression.
Definition: DVecSVecCrossExpr.h:138
Header file for the DisableIf class template.
Header file for the IsTemporary type trait class.
Namespace of the Blaze C++ math library.
Definition: Blaze.h:58
ElementType_t< VT2 > ET2
Element type of the right-hand side sparse vector expression.
Definition: DVecSVecCrossExpr.h:103
Header file for the If class template.
DVecSVecCrossExpr(const VT1 &lhs, const VT2 &rhs) noexcept
Constructor for the DVecSVecCrossExpr class.
Definition: DVecSVecCrossExpr.h:161
#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
TransposeType_t< ResultType > TransposeType
Transpose type for expression template evaluations.
Definition: DVecSVecCrossExpr.h:125
CrossTrait_t< RT1, RT2 > ResultType
Result type for expression template evaluations.
Definition: DVecSVecCrossExpr.h:124
#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
ReturnType operator[](size_t index) const
Subscript operator for the direct access to the vector elements.
Definition: DVecSVecCrossExpr.h:176
static constexpr bool simdEnabled
Compilation switch for the expression template evaluation strategy.
Definition: DVecSVecCrossExpr.h:149
Base class for N-dimensional dense vectors.The DenseVector class is a base class for all arbitrarily ...
Definition: DenseVector.h:76
const If_t< returnExpr, ExprReturnType, ElementType > ReturnType
Return type for expression template evaluations.
Definition: DVecSVecCrossExpr.h:129
#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.
Constraint on the data type.
Header file for the exception macros of the math module.
ElementType_t< ResultType > ElementType
Resulting element type.
Definition: DVecSVecCrossExpr.h:126
Constraint on the data type.
Header file for all forward declarations for expression class templates.
Constraint on the data type.
Header file for the EnableIf class template.
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
ResultType_t< VT1 > RT1
Result type of the left-hand side dense vector expression.
Definition: DVecSVecCrossExpr.h:96
Header file for the CrossExpr base class.
const ResultType CompositeType
Data type for composite expression templates.
Definition: DVecSVecCrossExpr.h:132
ReturnType_t< VT2 > RN2
Return type of the right-hand side sparse vector expression.
Definition: DVecSVecCrossExpr.h:99
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.
bool canAlias(const T *alias) const noexcept
Returns whether the expression can alias with the given address alias.
Definition: DVecSVecCrossExpr.h:240
const StaticVector< ET2, 3UL, TF > RT
Composite type of the right-hand side sparse vector expression.
Definition: DVecSVecCrossExpr.h:144
If_t< IsComputation_v< VT1 >, const StaticVector< ET1, 3UL, TF >, CT1 > LT
Composite type of the left-hand side dense vector expression.
Definition: DVecSVecCrossExpr.h:141
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
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
bool isAliased(const T *alias) const noexcept
Returns whether the expression is aliased with the given address alias.
Definition: DVecSVecCrossExpr.h:252
#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
RightOperand rhs_
Right-hand side sparse vector of the cross product expression.
Definition: DVecSVecCrossExpr.h:260
ResultType_t< VT2 > RT2
Result type of the right-hand side sparse vector expression.
Definition: DVecSVecCrossExpr.h:97
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
Header file for the IsComputation type trait class.
Base class for all compute expression templates.The Computation class serves as a tag for all computa...
Definition: Computation.h:66
LeftOperand leftOperand() const noexcept
Returns the left-hand side dense vector operand.
Definition: DVecSVecCrossExpr.h:218
Base class for sparse vectors.The SparseVector class is a base class for all arbitrarily sized (N-dim...
Definition: Forward.h:138
ReturnType at(size_t index) const
Checked access to the vector elements.
Definition: DVecSVecCrossExpr.h:195
CompositeType_t< VT2 > CT2
Composite type of the right-hand side sparse vector expression.
Definition: DVecSVecCrossExpr.h:101
static constexpr bool returnExpr
Compilation switch for the selection of the subscript operator return type.
Definition: DVecSVecCrossExpr.h:113
RightOperand rightOperand() const noexcept
Returns the right-hand side sparse vector operand.
Definition: DVecSVecCrossExpr.h:228
#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
If_t< IsExpression_v< VT1 >, const VT1, const VT1 &> LeftOperand
Composite type of the left-hand side dense vector expression.
Definition: DVecSVecCrossExpr.h:135
constexpr size_t size() const noexcept
Returns the current size/dimension of the vector.
Definition: DVecSVecCrossExpr.h:208
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