35 #ifndef _BLAZE_MATH_EXPRESSIONS_SVECDVECCROSSEXPR_H_
36 #define _BLAZE_MATH_EXPRESSIONS_SVECDVECCROSSEXPR_H_
83 template<
typename VT1
85 class SVecDVecCrossExpr :
public DenseVector< SVecDVecCrossExpr<VT1,VT2>, false >
143 enum { vectorizable = 0 };
146 enum { smpAssignable = 0 };
175 else if( index == 1UL )
176 return lhs_[2UL] * rhs_[0UL] -
lhs_[0UL] * rhs_[2UL];
178 return lhs_[0UL] * rhs_[1UL] -
lhs_[1UL] * rhs_[0UL];
189 inline ReturnType
at(
size_t index )
const {
193 return (*
this)[index];
233 template<
typename T >
235 return (
lhs_.canAlias( alias ) ||
rhs_.canAlias( alias ) );
245 template<
typename T >
247 return (
lhs_.isAliased( alias ) ||
rhs_.isAliased( alias ) );
269 template<
typename VT >
277 LT x(
serial( rhs.lhs_ ) );
278 RT y(
serial( rhs.rhs_ ) );
280 (~lhs)[0] = x[1UL]*y[2UL] - x[2UL]*y[1UL];
281 (~lhs)[1] = x[2UL]*y[0UL] - x[0UL]*y[2UL];
282 (~lhs)[2] = x[0UL]*y[1UL] - x[1UL]*y[0UL];
299 template<
typename VT >
311 const ResultType tmp(
serial( rhs ) );
329 template<
typename VT >
330 friend inline void addAssign( DenseVector<VT,false>& lhs,
const SVecDVecCrossExpr& rhs )
337 LT x(
serial( rhs.lhs_ ) );
338 RT y(
serial( rhs.rhs_ ) );
340 (~lhs)[0] += x[1UL]*y[2UL] - x[2UL]*y[1UL];
341 (~lhs)[1] += x[2UL]*y[0UL] - x[0UL]*y[2UL];
342 (~lhs)[2] += x[0UL]*y[1UL] - x[1UL]*y[0UL];
363 template<
typename VT >
364 friend inline void subAssign( DenseVector<VT,false>& lhs,
const SVecDVecCrossExpr& rhs )
371 LT x(
serial( rhs.lhs_ ) );
372 RT y(
serial( rhs.rhs_ ) );
374 (~lhs)[0] -= x[1UL]*y[2UL] - x[2UL]*y[1UL];
375 (~lhs)[1] -= x[2UL]*y[0UL] - x[0UL]*y[2UL];
376 (~lhs)[2] -= x[0UL]*y[1UL] - x[1UL]*y[0UL];
397 template<
typename VT >
398 friend inline void multAssign( DenseVector<VT,false>& lhs,
const SVecDVecCrossExpr& rhs )
405 LT x(
serial( rhs.lhs_ ) );
406 RT y(
serial( rhs.rhs_ ) );
408 (~lhs)[0] *= x[1UL]*y[2UL] - x[2UL]*y[1UL];
409 (~lhs)[1] *= x[2UL]*y[0UL] - x[0UL]*y[2UL];
410 (~lhs)[2] *= x[0UL]*y[1UL] - x[1UL]*y[0UL];
466 template<
typename T1
468 inline const SVecDVecCrossExpr<T1,T2>
473 if( (~lhs).
size() != 3UL || (~rhs).
size() != 3UL ) {
492 template<
typename VT1,
typename VT2 >
493 struct Size< SVecDVecCrossExpr<VT1,VT2> > :
public SizeT<3UL>
SelectType< IsExpression< VT1 >::value, const VT1, const VT1 & >::Type LeftOperand
Composite type of the left-hand side sparse vector expression.
Definition: SVecDVecCrossExpr.h:129
#define BLAZE_THROW_INVALID_ARGUMENT(MESSAGE)
Macro for the emission of a std::invalid_argument exceptionThis macro encapsulates the default way of...
Definition: Exception.h:187
CrossTrait< RT1, RT2 >::Type ResultType
Result type for expression template evaluations.
Definition: SVecDVecCrossExpr.h:118
Compile time check whether the given type is a temporary vector or matrix type.This type trait class ...
Definition: IsTemporary.h:87
const StaticVector< ET1, 3UL, false > LT
Composite type of the left-hand side sparse vector expression.
Definition: SVecDVecCrossExpr.h:135
Header file for basic type definitions.
BLAZE_ALWAYS_INLINE size_t size(const Vector< VT, TF > &vector)
Returns the current size/dimension of the vector.
Definition: Vector.h:252
LeftOperand lhs_
Left-hand side sparse vector of the cross product expression.
Definition: SVecDVecCrossExpr.h:253
SVecDVecCrossExpr< VT1, VT2 > This
Type of this SVecDVecCrossExpr instance.
Definition: SVecDVecCrossExpr.h:117
ResultType::TransposeType TransposeType
Transpose type for expression template evaluations.
Definition: SVecDVecCrossExpr.h:119
#define BLAZE_CONSTRAINT_MUST_BE_COLUMN_VECTOR_TYPE(T)
Constraint on the data type.In case the given data type T is not a column dense or sparse vector type...
Definition: ColumnVector.h:79
Efficient implementation of a compressed matrix.The CompressedMatrix class template is the represent...
Definition: CompressedMatrix.h:207
size_t size() const
Returns the current size/dimension of the vector.
Definition: SVecDVecCrossExpr.h:202
Expression object for sparse vector-dense vector cross products.The SVecDVecCrossExpr class represent...
Definition: Forward.h:121
const This & CompositeType
Data type for composite expression templates.
Definition: CompressedMatrix.h:2588
Type ElementType
Type of the sparse matrix elements.
Definition: CompressedMatrix.h:259
Header file for the DenseVector base class.
const DMatSerialExpr< MT, SO > serial(const DenseMatrix< MT, SO > &dm)
Forces the serial evaluation of the given dense matrix expression dm.
Definition: DMatSerialExpr.h:721
Header file for the Computation base class.
Header file for the SizeT class template.
VT1::ReturnType RN1
Return type of the left-hand side sparse vector expression.
Definition: SVecDVecCrossExpr.h:93
Constraint on the data type.
Base template for the CrossTrait class.
Definition: CrossTrait.h:110
Efficient implementation of a fixed-sized vector.The StaticVector class template is the representatio...
Definition: Forward.h:61
VT2::ElementType ET2
Element type of the right-hand side dense vector expression.
Definition: SVecDVecCrossExpr.h:98
ReturnType operator[](size_t index) const
Subscript operator for the direct access to the vector elements.
Definition: SVecDVecCrossExpr.h:170
Constraint on the transpose flag of vector types.
Header file for the MultExprTrait class template.
LeftOperand leftOperand() const
Returns the left-hand side sparse vector operand.
Definition: SVecDVecCrossExpr.h:212
Compile time type selection.The SelectType class template selects one of the two given types T1 and T...
Definition: SelectType.h:59
Header file for the IsTemporary type trait class.
Namespace of the Blaze C++ math library.
Definition: Blaze.h:57
VT2::ReturnType RN2
Return type of the right-hand side dense vector expression.
Definition: SVecDVecCrossExpr.h:94
RightOperand rightOperand() const
Returns the right-hand side dense vector operand.
Definition: SVecDVecCrossExpr.h:222
#define BLAZE_THROW_OUT_OF_RANGE(MESSAGE)
Macro for the emission of a std::out_of_range exceptionThis macro encapsulates the default way of Bla...
Definition: Exception.h:331
SelectType< IsExpression< VT2 >::value, const VT2, const VT2 & >::Type RightOperand
Composite type of the right-hand side dense vector expression.
Definition: SVecDVecCrossExpr.h:132
RightOperand rhs_
Right-hand side dense vector of the cross product expression.
Definition: SVecDVecCrossExpr.h:254
const DVecDVecCrossExpr< T1, T2 > operator%(const DenseVector< T1, false > &lhs, const DenseVector< T2, false > &rhs)
Operator for the cross product of two dense vectors ( ).
Definition: DVecDVecCrossExpr.h:466
Base class for N-dimensional dense vectors.The DenseVector class is a base class for all arbitrarily ...
Definition: DenseVector.h:70
#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:79
SelectType< IsComputation< VT2 >::value, const StaticVector< ET2, 3UL, false >, CT2 >::Type RT
Composite type of the right-hand side dense vector expression.
Definition: SVecDVecCrossExpr.h:138
#define BLAZE_CONSTRAINT_MUST_BE_REFERENCE_TYPE(T)
Constraint on the data type.In case the given data type T is not a reference type, a compilation error is created.
Definition: Reference.h:78
VT1::ElementType ET1
Element type of the left-hand side sparse vector expression.
Definition: SVecDVecCrossExpr.h:97
bool isAliased(const T *alias) const
Returns whether the expression is aliased with the given address alias.
Definition: SVecDVecCrossExpr.h:246
Constraint on the data type.
Evaluation of the return type of a subtraction expression.Via this type trait it is possible to evalu...
Definition: SubExprTrait.h:104
VT1::ResultType RT1
Result type of the left-hand side sparse vector expression.
Definition: SVecDVecCrossExpr.h:91
Type ElementType
Type of the sparse matrix elements.
Definition: CompressedMatrix.h:2586
Header file for the SelectType class template.
Header file for all forward declarations for expression class templates.
Constraint on the data type.
const ResultType CompositeType
Data type for composite expression templates.
Definition: SVecDVecCrossExpr.h:126
Header file for the serial shim.
Header file for the CrossExpr base class.
VT1::CompositeType CT1
Composite type of the left-hand side sparse vector expression.
Definition: SVecDVecCrossExpr.h:95
Header file for all forward declarations for dense vectors and matrices.
const Type & ReturnType
Return type for expression template evaluations.
Definition: CompressedMatrix.h:2587
SubExprTrait< typename MultExprTrait< RN1, RN2 >::Type, typename MultExprTrait< RN1, RN2 >::Type >::Type ExprReturnType
Expression return type for the subscript operator.
Definition: SVecDVecCrossExpr.h:112
Header file for run time assertion macros.
ReturnType at(size_t index) const
Checked access to the vector elements.
Definition: SVecDVecCrossExpr.h:189
Header file for the cross product trait.
VT2::ResultType RT2
Result type of the right-hand side dense vector expression.
Definition: SVecDVecCrossExpr.h:92
VT2::CompositeType CT2
Composite type of the right-hand side dense vector expression.
Definition: SVecDVecCrossExpr.h:96
#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:79
ResultType::ElementType ElementType
Resulting element type.
Definition: SVecDVecCrossExpr.h:120
Header file for the IsComputation type trait class.
bool canAlias(const T *alias) const
Returns whether the expression can alias with the given address alias.
Definition: SVecDVecCrossExpr.h:234
Base class for sparse vectors.The SparseVector class is a base class for all arbitrarily sized (N-dim...
Definition: Forward.h:118
#define BLAZE_FUNCTION_TRACE
Function trace macro.This macro can be used to reliably trace function calls. In case function tracin...
Definition: FunctionTrace.h:157
This ResultType
Result type for expression template evaluations.
Definition: CompressedMatrix.h:2583
Header file for exception macros.
Header file for the SubExprTrait class template.
Header file for the Size type trait.
Evaluation of the resulting expression type of a multiplication.Via this type trait it is possible to...
Definition: MultExprTrait.h:137
SVecDVecCrossExpr(const VT1 &lhs, const VT2 &rhs)
Constructor for the SVecDVecCrossExpr class.
Definition: SVecDVecCrossExpr.h:155
#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
const SelectType< returnExpr, ExprReturnType, ElementType >::Type ReturnType
Return type for expression template evaluations.
Definition: SVecDVecCrossExpr.h:123
Header file for the IsExpression type trait class.
Header file for the FunctionTrace class.