- L(i,1:i-1))) L(i+1:,i) = (A(i+1:,i) - matmul(
conjg(L(i,1:i-1)), L(i+1:,1:i-1))) / L(i,i) end do
where conjg refers to
complex conjugate of the elements...
-
ZDOTU = X^T * Y, CDOTC,
ZDOTC = X^H * Y
Fortran as dot_product(A,B) or sum(
conjg(A) * B)
Julia as A' * B or
standard library LinearAlgebra as dot(A, B)...
-
matrix AP to A do j=1,n k = j*(j-1)/2 A(1:j,j) = AP(1+k:j+k) A(j,1:j-1) =
conjg(AP(1+k:j-1+k)) end do An
example of
packed storage for
banded matrix: real...
- SIN COS TAN ASIN ACOS ATAN ATAN2 SINH COSH TANH
Complex numbers:
AIMAG CONJG The
following are for characters: LGE LGT LLE LLT
ICHAR CHAR
INDEX The simple...