Skip to contents

Magnitude of a clifford object

magnitude.Rd

Following Perwass, the magnitude of a multivector is defined as

||A||=A∗A

Where A∗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).

Usage

# S3 method for class 'clifford'
Mod(z)

Arguments

z

Clifford objects

Author

Robin K. S. Hankin

Note

If you want the square, ||A||2 and not ||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.

See also

Ops.clifford, Conj, rcliff

Examples



Mod(rcliff())
#> [1] 14.38749


# 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
#> [1] 82.19489

A <- rblade()
Mod(A*A - A %star% A)  # should be small
#> [1] 0