GENERIC INTERFACE MatrixFast(R, V, MS);

Arithmetic for Modula-3, see doc for details

TYPE T = MS.T;

PROCEDURE Add (x, y: T; ): T;    (* x + y *)
PROCEDURE Sub (x, y: T; ): T;    (* x - y *)
PROCEDURE IsZero (x: T; ): BOOLEAN;
PROCEDURE Equal (x, y: T; ): BOOLEAN; (* return v1=v2 *)

PROCEDURE Scale (x: T; y: R.T; ): T; (* x:=x*factor *)
PROCEDURE Mul (x, y: T; ): T;    (* x * y *)
PROCEDURE MulV (x: T; y: V.T; ): V.T; (* A * b *)
PROCEDURE MulTV (x: T; y: V.T; ): V.T; (* A^T * b or b^T*A *)
PROCEDURE MulMMA (x: T; ): T;    (* x*x^* *)
PROCEDURE MulMAM (x: T; ): T;    (* x^**x *)

CONST Adjoint = MS.Transpose;

PROCEDURE Trace (x: T; ): R.T;   (* sum of the diagonal elements *)
PROCEDURE Determinant (x: T; ): R.T;

END MatrixFast.