magnitude | R Documentation |
Following Perwass, the magnitude of a multivector is defined as
\mjdeqn\left|\left|A\right|\right| = \sqrtA\ast Aomitted; see PDF
Where \mjeqnA\ast A? denotes the Euclidean scalar product
eucprod()
. Recall that the Euclidean scalar product is never
negative (the function body is sqrt(abs(eucprod(z)))
; the
abs()
is needed to avoid numerical roundoff errors in
eucprod()
giving a negative value).
## S3 method for class 'clifford' Mod(z)
z |
Clifford objects |
If you want the square,
\mjeqn\left|\left|A\right|\right|^2||A||^2 and not
\mjeqn\left|\left|A\right|\right|||A||, it is faster and more accurate
to use eucprod(A)
, because this avoids a needless square root.
There is a nice example of scalar product at rcliff.Rd
.
Robin K. S. Hankin
Ops.clifford
,
Conj
,
rcliff
Mod(rcliff()) # Perwass, p68, asserts that if A is a k-blade, then (in his notation) # AA == A*A. # In package idiom, A*A == A %star% A: A <- rcliff() Mod(A*A - A %star% A) # meh A <- rblade() Mod(A*A - A %star% A) # should be small
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.