Table of contents
CSPICE_UNORM normalizes a double precision 3-vector and returns its
magnitude.
Given:
v1 arbitrary 3-vector(s), including the zero vector.
[3,n] = size(v1); double = class(v1)
the call:
[vout, vmag] = cspice_unorm( v1 )
returns:
vout the unit vector(s) in the direction of `v1'.
[3,n] = size(vout); double = class(vout)
If `v1' is the zero vector, then `vout' will also be the
zero vector.
vmag the magnitude(s) of `v1'.
[1,n] = size(vmag); double = class(vmag)
`vout' and `vmag' return with the same vectorization measure,
N, as `v1'.
None.
Any numerical results shown for this example may differ between
platforms as the results depend on the SPICE kernels used as input
and the machine specific arithmetic implementation.
1) Define a set of vectors and compute their corresponding
unit vectors and magnitudes.
Example code begins here.
function unorm_ex1()
%
% Local parameters.
%
SETSIZ = 2;
%
% Define the vector set.
%
seta = [ [5.0, 12.0, 0.0]', [1.e-7, 2.e-7, 2.e-7]' ];
%
% Calculate the unit vectors and magnitudes.
%
fprintf('Scalar case:\n')
for i=1:SETSIZ
[vout, vmag] = cspice_unorm( seta(:,i) );
fprintf( 'Vector : %12.8f %12.8f %12.8f\n', ...
seta(1,i), seta(2,i), seta(3,i) )
fprintf( 'Unit vector: %12.8f %12.8f %12.8f\n', ...
vout(1), vout(2), vout(3) )
fprintf( 'Magnitude : %12.8f\n', vmag )
fprintf( ' \n' )
end
%
% Repeat the operation with one single call to cspice_vhat.
%
[vout, vmag] = cspice_unorm( seta );
fprintf('Vectorized case:\n')
for i=1:SETSIZ
fprintf( 'Vector : %12.8f %12.8f %12.8f\n', ...
seta(1,i), seta(2,i), seta(3,i) )
fprintf( 'Unit vector: %12.8f %12.8f %12.8f\n', ...
vout(1,i), vout(2,i), vout(3,i) )
fprintf( 'Magnitude : %12.8f\n', vmag(i) )
fprintf( ' \n' )
end
When this program was executed on a Mac/Intel/Octave6.x/64-bit
platform, the output was:
Scalar case:
Vector : 5.00000000 12.00000000 0.00000000
Unit vector: 0.38461538 0.92307692 0.00000000
Magnitude : 13.00000000
Vector : 0.00000010 0.00000020 0.00000020
Unit vector: 0.33333333 0.66666667 0.66666667
Magnitude : 0.00000030
Vectorized case:
Vector : 5.00000000 12.00000000 0.00000000
Unit vector: 0.38461538 0.92307692 0.00000000
Magnitude : 13.00000000
Vector : 0.00000010 0.00000020 0.00000020
Unit vector: 0.33333333 0.66666667 0.66666667
Magnitude : 0.00000030
cspice_unorm references a function called cspice_vnorm (which itself is
numerically stable) to calculate the norm of the input vector `v1'.
If the norm is equal to zero, then each component of the output
vector `vout' is set to zero. Otherwise, `vout' is calculated by
dividing `v1' by the norm.
1) If the input argument `v1' is undefined, an error is signaled
by the Matlab error handling system.
2) If the input argument `v1' is not of the expected type, or it
does not have the expected dimensions and size, an error is
signaled by the Mice interface.
None.
None.
MICE.REQ
None.
J. Diaz del Rio (ODC Space)
E.D. Wright (JPL)
-Mice Version 1.1.0, 27-AUG-2021 (EDW) (JDR)
Edited the header to comply with NAIF standard. Added
complete code example to -Examples section.
Added -Parameters, -Exceptions, -Files, -Restrictions,
-Literature_References and -Author_and_Institution sections, and
completed -Particulars section.
Eliminated use of "lasterror" in rethrow.
Removed reference to the function's corresponding CSPICE header from
-Required_Reading section.
-Mice Version 1.0.2, 09-NOV-2012 (EDW)
Edited -I/O section to conform to NAIF standard for Mice
documentation.
-Mice Version 1.0.1, 21-APR-2010 (EDW)
Corrected description of call example to list all output arguments,
from:
vout = cspice_unorm(v1)
to
[vout, vmag] = cspice_unorm(v1)
-Mice Version 1.0.0, 25-APR-2006 (EDW)
3-dimensional unit vector and norm
|