Museum

Home

Lab Overview

Retrotechnology Articles

⇒ Online Manual

Media Vault

Software Library

Restoration Projects

Artifacts Sought

VEC_$DREC3C                       Domain/OS                        VEC_$DREC3C


NAME
     vec_$drec3c - calculate third-order linear recurrence of single-precision
     vector with constant coefficients

SYNOPSIS (C)
     #include <apollo/base.h>
     #include <apollo/vec.h>

     void vec_$drec3c(
          double *start_vec,
          long int &count,
          double &mult_constant,
          double &mult_constant_2,
          double &mult_constant_3,
          double *result_vec)

SYNOPSIS (Pascal)
     %include '/sys/ins/base.ins.pas';
     %include '/sys/ins/vec.ins.pas';

     procedure vec_$drec3c(
          in start_vec: univ vec_$double_vector;
          in count: integer32;
          in mult_constant: double;
          in mult_constant_2: double;
          in mult_constant_3: double;
          var result_vec: univ vec_$double_vector);

SYNOPSIS (FORTRAN)
     %include '/sys/ins/base.ins.ftn'
     %include '/sys/ins/vec.ins.ftn'

           parameter (nvec = 10)

           real*8 start_vec(nvec), result_vec(nvec)
           real*8 mult_constant, mult_constant_2, mult_constant_3
           integer*4 count

           call vec_$drec3c(start_vec, count, mult_constant,
          & mult_constant_2, mult_constant_3, result_vec)

DESCRIPTION
     Vec_$drec3c calculates a third-order linear recurrence with constant
     coefficients based on start_vec.  It differs from vec_$rec3c in that the
     vectors being handled are double-precision floating point.

     This routine reads the first three entries in the array result_vec; it
     then writes count entries into result_vec based on the following formula:

     For each integer I such that 1 <= I <= count,
     result_vec(I+3) = start_vec(I) + (mult_constant x result_vec(I+2))
                       + (mult_constant_2 x result_vec(I+1))
                       + (mult_constant_3 x result_vec(I))

     Notice that result_vec is used both for input and output, and that it
     must be large enough to hold count + 3 entries.

     start_vec
          An input vector.

     mult_constant
          A scalar multiplier.

     mult_constant_2
          A scalar multiplier.

     mult_constant_3
          A scalar multiplier.

     count
          The number of elements to be operated on.

     result_vec
          The vector containing two input values and the result of the
          recurrence calculation.

NOTES
     When vec_$drec3c is used to operate on matrixes in C and Pascal,
     start_vec and result_vec are row vectors; in FORTRAN, they are column
     vectors.

     As in all the vec_$ calls,  the result array must not overlap any of the
     input arrays; the result array may be identical to an input, but must not
     contain any smaller subset of an input.   Because of pipelining, using
     overlapping input and output arrays may cause incorrect results.


SEE ALSO
     vec_$rec3, vec_$rec3c, vec_$rec3c_i, vec_$drec3c_i, vec_$irec3c,
     vec_$irec3c_i, vec_$irec3c16, vec_$irec3c16_i.

Typewritten Software • bear@typewritten.org • Edmonds, WA 98026