BLAS 1 Feature completeness

[crtrott]

This is the issue to track progress for BLAS 1 feature completeness.

Function	Operation Performed	BLAS Name	KK impl done	TPL hooks	Supports MV
abs	Y(i) = abs(X(i))	--	X		X
axpby	Y(i) = a*X(i) + b*Y(i)	--	X		X
axpy	Y(i) += a*X(i)	axpy	X	X	X
dot	Sum( X(i)*Y(i) )	dot	X	X	X
iamax	MaxIndex(abs(X(i)))	iamax	X
mult	Z(i)=a*A(i)*X(i)+c*Z(i)	--	X		X
nrm1	Sum( abs(X(i))	asum	X	X	X
nrm2	sqrt(Sum(abs(X(i))*abs(X(i))))	nrm2	X	X	X
nrm2_squared	Sum ( abs(X(i))*abs(X(i)) )	--	X		X
nrminf	Max(abs(X(i)))	--	X		X
reciprocal	Y(i) = 1/X(i)	--	X		X
rot		rot	X	X
rotm		rotm	X	X
rotg		rotg	X	X
rotmg		rotmg	X	X
scal	Y(i) = a*X(i)	scal*	X	X	X
sum	Sum( X(i) )	--	X		X
update	Z(i)=a*X(i)+b*Y(i)+c*Z(i)	--	X		X

[crtrott]

There are some inconsistencies in our API. Blas generally has the output array as the last one for example: axpy(a,X,Y) : Y += aX. We do that for some like axpby(a,X,b,Y) and update(a,X,b,Y,c,Z) but not mult(c,Y,a,A,X) : Y = cY + bAX or the multivector dot(R,X,Y) where R is the result vector.

[crtrott]

rot, rotm, rotg, rotmg nobody uses them and we do that on 2x2 so nobody cares.

[crtrott]

Siva does iamax

[kliegeois]

@lucbv can we close this one as it is completed?

[lucbv]

Sure