VDmath Namespace Reference

Functions
template<class T >
T	Ddot2 (Nektar::LibUtilities::CommSharedPtr &pComm, int n, const T *w, const T *x, const int *y)
	dot2 (vector times vector times vector): z = w*x*y More...
template<class T >
T	Ddot2 (Nektar::LibUtilities::CommSharedPtr &pComm, int n, const T *w, const int incw, const T *x, const int incx, const int *y, const int incy)
	dot2 (vector times vector times vector): z = w*x*y More...
template<class T >
T	Ddot2 (LibUtilities::CommSharedPtr &pComm, int n, const Array< OneD, const T > &w, const Array< OneD, const T > &x, const Array< OneD, const int > &y)
template<class T >
T	Ddot2 (LibUtilities::CommSharedPtr &pComm, int n, const Array< OneD, const T > &w, const int incw, const Array< OneD, const T > &x, const int incx, const Array< OneD, const int > &y, const int incy)

Function Documentation

Ddot2() [1/4]

template<class T >

T VDmath::Ddot2	(	LibUtilities::CommSharedPtr &	pComm,
		int	n,
		const Array< OneD, const T > &	w,
		const Array< OneD, const T > &	x,
		const Array< OneD, const int > &	y
	)

Definition at line 48 of file VDmathArray.hpp.

{
    ASSERTL1(n <= w.size() + w.GetOffset(), "Array out of bounds");
    ASSERTL1(n <= x.size() + x.GetOffset(), "Array out of bounds");
    ASSERTL1(n <= y.size() + y.GetOffset(), "Array out of bounds");
 
    return Ddot2(pComm, n, &w[0], &x[0], &y[0]);
}

References ASSERTL1, and Ddot2().

Ddot2() [2/4]

template<class T >

T VDmath::Ddot2	(	LibUtilities::CommSharedPtr &	pComm,
		int	n,
		const Array< OneD, const T > &	w,
		const int	incw,
		const Array< OneD, const T > &	x,
		const int	incx,
		const Array< OneD, const int > &	y,
		const int	incy
	)

Definition at line 61 of file VDmathArray.hpp.

{
    ASSERTL1(n * incw <= w.size() + w.GetOffset(), "Array out of bounds");
    ASSERTL1(n * incx <= x.size() + x.GetOffset(), "Array out of bounds");
    ASSERTL1(n * incy <= y.size() + y.GetOffset(), "Array out of bounds");
 
    return Ddot2(pComm, n, &w[0], incw, &x[0], incx, &y[0], incy);
}

References ASSERTL1, and Ddot2().

Ddot2() [3/4]

template<class T >

T VDmath::Ddot2	(	Nektar::LibUtilities::CommSharedPtr &	pComm,
		int	n,
		const T *	w,
		const int	incw,
		const T *	x,
		const int	incx,
		const int *	y,
		const int	incy
	)

dot2 (vector times vector times vector): z = w*x*y

Definition at line 63 of file VDmath.hpp.

{
    T sum = 0;
 
    while (n--)
    {
        sum += (*y == 1 ? (*w) * (*x) : 0.0);
        w += incw;
        x += incx;
        y += incy;
    }
    pComm->AllReduce(sum, Nektar::LibUtilities::ReduceSum);
    return sum;
}

References Nektar::LibUtilities::ReduceSum.

Ddot2() [4/4]

template<class T >

T VDmath::Ddot2	(	Nektar::LibUtilities::CommSharedPtr &	pComm,
		int	n,
		const T *	w,
		const T *	x,
		const int *	y
	)

dot2 (vector times vector times vector): z = w*x*y

Definition at line 45 of file VDmath.hpp.

{
    T sum = 0;
 
    while (n--)
    {
        sum += (*y == 1 ? (*w) * (*x) : 0);
        ++w;
        ++x;
        ++y;
    }
    pComm->AllReduce(sum, Nektar::LibUtilities::ReduceSum);
    return sum;
}

References Nektar::LibUtilities::ReduceSum.

Referenced by Ddot2().