#include <deal.II/lac/cuda_vector.h>

Inheritance diagram for LinearAlgebra::CUDAWrappers::Vector< Number >:

Public Member Functions
	Vector ()

	Vector (const Vector< Number > &V)

	Vector (const size_type n)

	~Vector ()

void	reinit (const size_type n, const bool omit_zeroing_entries=false)

virtual void	reinit (const VectorSpaceVector< Number > &V, const bool omit_zeroing_entries=false) override

virtual void	import (const ReadWriteVector< Number > &V, VectorOperation::values operation, std::shared_ptr< const CommunicationPatternBase > communication_pattern=std::shared_ptr< const CommunicationPatternBase >()) override

virtual Vector< Number > &	operator= (const Number s) override

virtual Vector< Number > &	operator*= (const Number factor) override

virtual Vector< Number > &	operator/= (const Number factor) override

virtual Vector< Number > &	operator+= (const VectorSpaceVector< Number > &V) override

virtual Vector< Number > &	operator-= (const VectorSpaceVector< Number > &V) override

virtual Number	operator* (const VectorSpaceVector< Number > &V) const override

virtual void	add (const Number a) override

virtual void	add (const Number a, const VectorSpaceVector< Number > &V) override

virtual void	add (const Number a, const VectorSpaceVector< Number > &V, const Number b, const VectorSpaceVector< Number > &W) override

virtual void	sadd (const Number s, const Number a, const VectorSpaceVector< Number > &V) override

virtual void	scale (const VectorSpaceVector< Number > &scaling_factors) override

virtual void	equ (const Number a, const VectorSpaceVector< Number > &V) override

virtual bool	all_zero () const override

virtual value_type	mean_value () const override

virtual real_type	l1_norm () const override

virtual real_type	l2_norm () const override

virtual real_type	linfty_norm () const override

virtual Number	add_and_dot (const Number a, const VectorSpaceVector< Number > &V, const VectorSpaceVector< Number > &W) override

Number *	get_values () const

virtual size_type	size () const override

virtual ::IndexSet	locally_owned_elements () const override

virtual void	print (std::ostream &out, const unsigned int precision=2, const bool scientific=true, const bool across=true) const override

virtual std::size_t	memory_consumption () const override

Public Member Functions inherited from LinearAlgebra::VectorSpaceVector< Number >
virtual void	compress (VectorOperation::values)

virtual	~VectorSpaceVector ()=default

Static Public Member Functions
static::ExceptionBase &	ExcVectorTypeNotCompatible ()

Private Attributes
Number *	val

size_type	n_elements

Detailed Description

template<typename Number>
class LinearAlgebra::CUDAWrappers::Vector< Number >

This class implements a vector using CUDA for use on Nvidia GPUs. This class is derived from the LinearAlgebra::VectorSpaceVector class.

Note: Only float and double are supported.

See also: CUDAWrappers

Author: Karl Ljungkvist, Bruno Turcksin, 2016

Definition at line 53 of file cuda_vector.h.

Constructor & Destructor Documentation

template<typename Number>

LinearAlgebra::CUDAWrappers::Vector< Number >::Vector ( )

Constructor. Create a vector of dimension zero.

template<typename Number>

LinearAlgebra::CUDAWrappers::Vector< Number >::Vector ( const Vector< Number > & V )

Copy constructor.

template<typename Number>

LinearAlgebra::CUDAWrappers::Vector< Number >::Vector ( const size_type n )

explicit

Constructor. Set dimension to n and initialize all elements with zero.

The constructor is made explicit to avoid accident like this: v=0;. Presumably, the user wants to set every elements of the vector to zero, but instead, what happens is this call: v=Vector<Number>(0);, i.e. the vector is replaced by one of length zero.

template<typename Number>

LinearAlgebra::CUDAWrappers::Vector< Number >::~Vector ( )

Destructor.

Member Function Documentation

template<typename Number>

void LinearAlgebra::CUDAWrappers::Vector< Number >::reinit	(	const size_type	n,
		const bool	omit_zeroing_entries = `false`
	)

Reinit functionality. The flag omit_zeroing_entries determines whether the vector should be filled with zero (false) or left untouched (true).

template<typename Number>

virtual void LinearAlgebra::CUDAWrappers::Vector< Number >::reinit	(	const VectorSpaceVector< Number > &	V,
		const bool	omit_zeroing_entries = `false`
	)

overridevirtual

Change the dimension to that of the vector V. The elements of V are not copied.

Implements LinearAlgebra::VectorSpaceVector< Number >.

template<typename Number>

virtual void LinearAlgebra::CUDAWrappers::Vector< Number >::import	(	const ReadWriteVector< Number > &	V,
		VectorOperation::values	operation,
		std::shared_ptr< const CommunicationPatternBase >	communication_pattern = `std::shared_ptr< const CommunicationPatternBase >()`
	)

overridevirtual

Import all the element from the input vector V. VectorOperation::values operation is used to decide if the elements int V should be added to the current vector or replace the current elements. The last parameter is not used. It is only used for distributed vectors. This is the function that should be used to copy a vector to the GPU.

Implements LinearAlgebra::VectorSpaceVector< Number >.

template<typename Number>

virtual Vector<Number>& LinearAlgebra::CUDAWrappers::Vector< Number >::operator= ( const Number s )

overridevirtual

Sets all elements of the vector to the scalar s. This operation is only allowed if s is equal to zero.

Implements LinearAlgebra::VectorSpaceVector< Number >.

template<typename Number>

virtual Vector<Number>& LinearAlgebra::CUDAWrappers::Vector< Number >::operator*= ( const Number factor )

overridevirtual

Multiply the entive vector by a fixed factor.

Implements LinearAlgebra::VectorSpaceVector< Number >.

template<typename Number>

virtual Vector<Number>& LinearAlgebra::CUDAWrappers::Vector< Number >::operator/= ( const Number factor )

overridevirtual

Divide the entire vector by a fixed factor.

Implements LinearAlgebra::VectorSpaceVector< Number >.

template<typename Number>

virtual Vector<Number>& LinearAlgebra::CUDAWrappers::Vector< Number >::operator+= ( const VectorSpaceVector< Number > & V )

overridevirtual

Add the vector V to the present one.

Implements LinearAlgebra::VectorSpaceVector< Number >.

template<typename Number>

virtual Vector<Number>& LinearAlgebra::CUDAWrappers::Vector< Number >::operator-= ( const VectorSpaceVector< Number > & V )

overridevirtual

Subtract the vector V from the present one.

Implements LinearAlgebra::VectorSpaceVector< Number >.

template<typename Number>

virtual Number LinearAlgebra::CUDAWrappers::Vector< Number >::operator* ( const VectorSpaceVector< Number > & V ) const

overridevirtual

Return the scalar product of two vectors.

Implements LinearAlgebra::VectorSpaceVector< Number >.

template<typename Number>

virtual void LinearAlgebra::CUDAWrappers::Vector< Number >::add ( const Number a )

overridevirtual

Add to all components. Note that a is a scalar not a vector.

Implements LinearAlgebra::VectorSpaceVector< Number >.

template<typename Number>

virtual void LinearAlgebra::CUDAWrappers::Vector< Number >::add	(	const Number	a,
		const VectorSpaceVector< Number > &	V
	)

overridevirtual

Simple addition of a multiple of a vector, i.e. *this += a*V.

Implements LinearAlgebra::VectorSpaceVector< Number >.

template<typename Number>

virtual void LinearAlgebra::CUDAWrappers::Vector< Number >::add	(	const Number	a,
		const VectorSpaceVector< Number > &	V,
		const Number	b,
		const VectorSpaceVector< Number > &	W
	)

overridevirtual

Multiple addition of scaled vectors, i.e. *this += a*V.

Implements LinearAlgebra::VectorSpaceVector< Number >.

template<typename Number>

virtual void LinearAlgebra::CUDAWrappers::Vector< Number >::sadd	(	const Number	s,
		const Number	a,
		const VectorSpaceVector< Number > &	V
	)

overridevirtual

Scaling and simple addition of a multiple of a vector, i.e. *this = s*(*this)+a*V

Implements LinearAlgebra::VectorSpaceVector< Number >.

template<typename Number>

virtual void LinearAlgebra::CUDAWrappers::Vector< Number >::scale ( const VectorSpaceVector< Number > & scaling_factors )

overridevirtual

Scale each element of this vector by the corresponding element in the argument. This function is mostly meant to simulate multiplication (and immediate re-assignment) by a diagonal scaling matrix.

Implements LinearAlgebra::VectorSpaceVector< Number >.

template<typename Number>

virtual void LinearAlgebra::CUDAWrappers::Vector< Number >::equ	(	const Number	a,
		const VectorSpaceVector< Number > &	V
	)

overridevirtual

Assignment *this = a*V.

Implements LinearAlgebra::VectorSpaceVector< Number >.

template<typename Number>

virtual bool LinearAlgebra::CUDAWrappers::Vector< Number >::all_zero ( ) const

overridevirtual

Return whether the vector contains only elements with value zero.

Implements LinearAlgebra::VectorSpaceVector< Number >.

template<typename Number>

virtual value_type LinearAlgebra::CUDAWrappers::Vector< Number >::mean_value ( ) const

overridevirtual

Return the mean value of all the entries of this vector.

Implements LinearAlgebra::VectorSpaceVector< Number >.

template<typename Number>

virtual real_type LinearAlgebra::CUDAWrappers::Vector< Number >::l1_norm ( ) const

overridevirtual

Return the l₁ norm of the vector (i.e., the sum of the absolute values of all entries among all processors).

Implements LinearAlgebra::VectorSpaceVector< Number >.

template<typename Number>

virtual real_type LinearAlgebra::CUDAWrappers::Vector< Number >::l2_norm ( ) const

overridevirtual

Return the l₂ norm of the vector (i.e., the square root of the sum of the square of all entries among all processors).

Implements LinearAlgebra::VectorSpaceVector< Number >.

template<typename Number>

virtual real_type LinearAlgebra::CUDAWrappers::Vector< Number >::linfty_norm ( ) const

overridevirtual

Return the maximum norm of the vector (i.e., the maximum absolute value among all entries and among all processors).

Implements LinearAlgebra::VectorSpaceVector< Number >.

template<typename Number>

virtual Number LinearAlgebra::CUDAWrappers::Vector< Number >::add_and_dot	(	const Number	a,
		const VectorSpaceVector< Number > &	V,
		const VectorSpaceVector< Number > &	W
	)

overridevirtual

Perform a combined operation of a vector addition and a subsequent inner product, returning the value of the inner product. In other words, the result of this function is the same as if the user called

this->add(a, V);

return_value = *this * W;

The reason this function exists is that this operation involves less memory transfer than calling the two functions separately. This method only needs to load three vectors, this, V, W, whereas calling separate methods means to load the calling vector this twice. Since most vector operations are memory transfer limited, this reduces the time by 25% (or 50% if W equals this).

For complex-valued vectors, the scalar product in the second step is implemented as \(\left<v,w\right>=\sum_i v_i \bar{w_i}\).

Implements LinearAlgebra::VectorSpaceVector< Number >.

template<typename Number >

Number * Vector< Number >::get_values ( ) const

inline

Return the pointer to the underlying array.

Definition at line 314 of file cuda_vector.h.

template<typename Number >

Vector< Number >::size_type Vector< Number >::size ( ) const

inlineoverridevirtual

Return the size of the vector.

Implements LinearAlgebra::VectorSpaceVector< Number >.

Definition at line 323 of file cuda_vector.h.

template<typename Number >

IndexSet Vector< Number >::locally_owned_elements ( ) const

inlineoverridevirtual

Return an index set that describe which elements of this vector are owned by the current processor, i.e. [0, size).

Implements LinearAlgebra::VectorSpaceVector< Number >.

Definition at line 331 of file cuda_vector.h.

template<typename Number>

virtual void LinearAlgebra::CUDAWrappers::Vector< Number >::print	(	std::ostream &	out,
		const unsigned int	precision = `2`,
		const bool	scientific = `true`,
		const bool	across = `true`
	)		const