Table of contents
CSPICE_VHAT returns the unit vector along a double precision
3-dimensional vector.
Given:
v1 an arbitrary vector(s).
[3,n] = size(v1); double = class(v1)
the call:
[vout] = cspice_vhat( v1 )
returns:
vout the unit vector(s) in the direction of `v1'.
[3,n] = size(vout); double = class(vout)
If `v1' represents the zero vector, then `vout' will
also be the zero vector.
`vout' returns 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
vector.
Example code begins here.
function vhat_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.
%
fprintf('Scalar case:\n')
for i=1:SETSIZ
[vout] = cspice_vhat( 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( ' \n' )
end
%
% Repeat the operation with one single call to cspice_vhat.
%
[vout] = cspice_vhat( 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( ' \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
Vector : 0.00000010 0.00000020 0.00000020
Unit vector: 0.33333333 0.66666667 0.66666667
Vectorized case:
Vector : 5.00000000 12.00000000 0.00000000
Unit vector: 0.38461538 0.92307692 0.00000000
Vector : 0.00000010 0.00000020 0.00000020
Unit vector: 0.33333333 0.66666667 0.66666667
cspice_vhat determines the magnitude of `v1' and then divides each
component of `v1' by the magnitude. This process is highly stable
over the whole range of 3-dimensional vectors.
1) If `v1' represents the zero vector, then `vout' will also be the
zero vector.
2) If the input argument `v1' is undefined, an error is signaled
by the Matlab error handling system.
3) 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.
1) There is no known case whereby floating point overflow may
occur. Thus, no error recovery or reporting scheme is
incorporated into this routine.
MICE.REQ
None.
J. Diaz del Rio (ODC Space)
E.D. Wright (JPL)
-Mice Version 1.1.0, 10-AUG-2021 (EDW) (JDR)
Updated 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, 18-DEC-2014 (EDW)
Edited -I/O section to conform to NAIF standard for Mice
documentation.
-Mice Version 1.0.1, 30-DEC-2008 (EDW)
Corrected misspellings.
-Mice Version 1.0.0, 25-APR-2006 (EDW)
unitize a 3-dimensional vector
|