vscl |
Table of contents
ProcedureVSCL ( Vector scaling, 3 dimensions ) SUBROUTINE VSCL ( S, V1, VOUT ) AbstractMultiply a scalar and a double precision 3-dimensional vector. Required_ReadingNone. KeywordsVECTOR DeclarationsIMPLICIT NONE DOUBLE PRECISION S DOUBLE PRECISION V1 ( 3 ) DOUBLE PRECISION VOUT ( 3 ) Brief_I/OVARIABLE I/O DESCRIPTION -------- --- -------------------------------------------------- S I Scalar to multiply a vector. V1 I Vector to be multiplied. VOUT O Product vector, S * V1. Detailed_InputS is a double precision scalar used to multiply the vector V1. V1 is a double precision 3-dimensional vector, which is to be scaled by S. Detailed_OutputVOUT is a double precision 3-dimensional vector containing the product of the scalar with the vector V1. ParametersNone. ExceptionsError free. FilesNone. ParticularsFor each value of the index I from 1 to 3, this subroutine performs the following multiplication VOUT(I) = S * V1(I) No error checking is performed to guard against numeric overflow or underflow. ExamplesThe numerical results shown for this example may differ across platforms. The results depend on the SPICE kernels used as input, the compiler and supporting libraries, and the machine specific arithmetic implementation. 1) Define a sets of scalar double precision values and use them to scale a given 3-dimensional vector. Example code begins here. PROGRAM VSCL_EX1 IMPLICIT NONE C C Local parameters. C INTEGER SETSIZ PARAMETER ( SETSIZ = 3 ) C C Local variables. C DOUBLE PRECISION S ( SETSIZ ) DOUBLE PRECISION V1 ( 3 ) DOUBLE PRECISION VOUT ( 3 ) INTEGER I C C Define the set of scalars and the input vector. C DATA S / 3.D0, 0.D0, -1.D0 / DATA V1 / 1.D0, 2.D0, -3.D0 / WRITE(*,'(A,3F6.1)') 'Input vector : ', V1 WRITE(*,*) C C Calculate product of each scalar and V1. C DO I=1, SETSIZ CALL VSCL ( S(I), V1, VOUT ) WRITE(*,'(A,F6.1)') 'Scale factor : ', S(I) WRITE(*,'(A,3F6.1)') 'Output vector: ', VOUT WRITE(*,*) END DO END When this program was executed on a Mac/Intel/gfortran/64-bit platform, the output was: Input vector : 1.0 2.0 -3.0 Scale factor : 3.0 Output vector: 3.0 6.0 -9.0 Scale factor : 0.0 Output vector: 0.0 0.0 -0.0 Scale factor : -1.0 Output vector: -1.0 -2.0 3.0 Restrictions1) The user is responsible for insuring that no floating point overflow occurs from multiplying S by any component of V1. No error recovery or reporting scheme is incorporated in this subroutine. Literature_ReferencesNone. Author_and_InstitutionN.J. Bachman (JPL) J. Diaz del Rio (ODC Space) W.M. Owen (JPL) W.L. Taber (JPL) VersionSPICELIB Version 1.1.0, 06-JUL-2021 (JDR) Added IMPLICIT NONE statement. Edited the header to comply with NAIF standard. Added complete code example based on existing example. SPICELIB Version 1.0.2, 23-APR-2010 (NJB) Header correction: assertions that the output can overwrite the input have been removed. SPICELIB Version 1.0.1, 10-MAR-1992 (WLT) Comment section for permuted index source lines was added following the header. SPICELIB Version 1.0.0, 31-JAN-1990 (WMO) |
Fri Dec 31 18:37:06 2021