--- a/LibDBlas.def
+++ b/LibDBlas.def
@@ -40,13 +40,14 @@
(* *)
(* zswap : swap two vectors *)
(* zcopy : copy vectors to vector *)
+ (* zdotu : dot product of two vectors *)
(* zdotc : dot product of two vectors *)
(* dznrm2 : euklidian norm of a vector *)
(* zscal : scale vector by a constant *)
(* zaxpy : constant times a vector plus a vector *)
(* zdrot : applies a plane rotation *)
(* *)
- (* zgemv : matrix matrix operations (complex) *)
+ (* zgemv : matrix vector operations (complex) *)
(* zgemm : matrix matrix operations (complex) *)
(* *)
(*------------------------------------------------------------------------*)
@@ -405,7 +406,7 @@
(* *)
(* A := Alpha*(x;y) + A, *)
(* *)
- (* where Alpha is a scalar, x is an m element vector, y is *)
+ (* where Alpha is a scalar, X is an m element vector, Y is *)
(* an n element vector and A is an m by n matrix. *)
(* (x;y) denotes a dyadic product of the vectors x and y. *)
(* *)
@@ -464,6 +465,17 @@
(* MRi, Modula-2 10.04.16 | 09.09.18 (complex version) *)
(*----------------------------------------------------------------*)
+PROCEDURE zdotu( N : INTEGER;
+ VAR X : (* ARRAY OF *) CFLOAT;
+ IncX : INTEGER;
+ VAR Y : (* ARRAY OF *) CFLOAT;
+ IncY : INTEGER) : CFLOAT;
+
+ (*----------------------------------------------------------------*)
+ (* Forms the dot product of two vectors X*Y *)
+ (* Uses unrolled loops for increments equal to one. *)
+ (*----------------------------------------------------------------*)
+
PROCEDURE zdotc( N : INTEGER;
VAR X : (* ARRAY OF *) CFLOAT;
IncX : INTEGER;
@@ -471,10 +483,9 @@
IncY : INTEGER) : CFLOAT;
(*----------------------------------------------------------------*)
- (* Forms the dot product of two vectors. Uses unrolled loops for *)
- (* increments equal to one. *)
- (*----------------------------------------------------------------*)
-
+ (* Forms the dot product of two vectors conj(X)*Y *)
+ (* Uses unrolled loops for increments equal to one. *)
+ (*----------------------------------------------------------------*)
PROCEDURE dznrm2( N : INTEGER;
VAR X : (* ARRAY OF *) CFLOAT;
@@ -482,7 +493,7 @@
(*----------------------------------------------------------------*)
(* dznrm2 returns the euclidean norm of a vector so that *)
- (* dznrm2 := sqrt( X**H*X ) *)
+ (* dznrm2 := sqrt(conj(X')*X) *)
(*----------------------------------------------------------------*)
PROCEDURE zscal( n : INTEGER;