FMSMM

PURPOSE: Matrix Multiply

[C] = [D] - [A]T{B}, IACCUM = -1, TRANSA='T'
[C] = [D] - [A] {B}, IACCUM = -1, TRANSA='N'


[C] =       [A]T{B}, IACCUM =  0, TRANSA='T'
[C] =       [A] {B}, IACCUM =  0, TRANSA='N'


[C] = [D] + [A]T{B}, IACCUM = +1, TRANSA='T'
[C] = [D] + [A] {B}, IACCUM = +1, TRANSA='N'

SYNOPSIS

INPUT PARAMETERS:

• LUC(25) = Integer array.
  Vector or matrix file attributes for the output array [C].

• LUD(25) = Integer array.
  Vector or matrix file attributes for the input array [D].
  This array is used to initialize the output array [C] when IACCUM is -1 or 1.
  This parameter must be the same type (matrix, vector) and size as LUC.
  This parameter may be the same as LUC.

• IACCUM = Integer.
  Determines if the matrix product of [A] and [B] is subtracted, stored into or added to [C].
  This value overrides the FMS Parameter IACCUM.
  Possible values are -1, 0 and 1.

• TRANSA = Character*1.
  Determines if the Vector file {A} is transposed.
  Possible values are 'N' and 'T'.
  NOTE: TRANSA is always 'N' when [A] is a matrix or submatrix file.

• LUA(25) = Integer array.
  Vector, matrix or submatrix file attributes for the input array [A].

• LUB(25) = Integer array.
  Vector file attributes for the input array {B}.

  OUTPUT PARAMETERS:

  None.

  FMS PARAMETERS:

  The following FMS Parameters are especially important to this routine:

  Parameter	Description
  MMROW	Starting row in [C] to store the product of [A] and [B]. Default value = 0
  MMCOL	Starting column in [C] to store the product of [A] and [B]. Default value = 0
  MMKA	Number of [0] columns to add to the beginning of [A] Default value = 0
  IPRMV	Matrix-vectors multiply print code
  IPRVV	Vectors-vectors multiply print code

  DESCRIPTION:

  This subroutine performs one of the following multiply operations:

  CASE	LUC and LUD	IACCUM	TRANSA	LUA	LUB	Compute
  1	Vector File	-1 0 1	'N'	Matrix File	Vector File	{C} = {D} - [A]{B} {C} = [A]{B} {C} = {D} + [A]{B}
  2	Vector File	-1 0 1	'N'	Submatrix File	Vector File	{C} = {D} - [Si]{B} {C} = [Si]{B} {C} = {D} + [Si]{B}
  3	Vector File	-1 0 1	'N'	Vector File	Vector File	{C} = {D} - {A}{B} {C} = {A}{B} {C} = {D} + {A}{B}
  4	Vector File	-1 0 1	'T'	Vector File	Vector File	{C} = {D} - {A}T{B} {C} = {A}T{B} {C} = {D} + {A}T{B}
  5	Matrix File	-1 0 1	'T'	Vector File	Vector File	[C] = [D] - {A}T{B} [C] = {A}T{B} [C] = [D] + {A}T{B}

  The calculation proceeds in two steps. First, the output array [C] is initialized to the input array [D] when IACCUM is -1 or 1. When IACCUM=0, [C] is initialized to [0].

  The matrix product [A][B] or [A]T[B] is then subtracted from [C], stored into [C] or added to [C], according to the value of IACCUM (-1, 0, or 1). The matrix [A] may be transposed (TRANSA='T') or not (TRANSA='N').

  Normally the product of [A] and [B] is stored in [C] starting at the first row and column of matrix [C]. However, the FMS Parameters MMROW and MMCOL may be used to shift the storage to start at C(1+MMROW,1+MMCOL). The default values of MMROW and MMCOL are 0. This option only applies when [C] is a matrix file in block or slab format.

  If any part of the product of [A] and [B], including any shifting by MMROW and MMCOL, falls outside of the matrix [C], those values are not computed. If the product of [A] and [B] is smaller than [C], only the appropriate terms in [C] are modified.

  The FMS Parameter MMKA may be used to shift the matrix [A] by adding MMKA columns of zero's. The shifted matrix [A] then becomes [0|A]. The default value of MMKA is 0.

  EXAMPLE_15 illustrates how to use this subroutine and the FMS Parameters MMROW, MMCOL and MMKA.

  ---

  
  Copyright © Multipath Corporation