RM Procedure

Reads data into a one- or two-dimensional matrix.

RM, a [, rows, columns]

rows (optional) Number of rows in the matrix.

columns (optional) Number of columns in the matrix.

a Named variable into which the data is stored.

Complex If present and nonzero, creates a single-precision complex matrix.

Dcomplex If present and nonzero, creates a double-precision complex matrix.

Double If present and nonzero, creates a double-precision matrix.

Procedure RM is used to read data into matrices according to the PV-WAVE matrix-storage scheme. If rows and columns are not specified, RM attempts to use the current definition of a to determine the number of rows and columns of a. If rows and columns are not specified and a is undefined or a scalar, an error is issued.

Upon invoking RM, the user is prompted with the row number for which input is expected. The prompt for the next row to be filled does not appear until the current row is filled. If the amount of data input for a particular row is larger than the defined number of columns of a, then the extra trailing input is ignored, and the prompt for the next row is given.

The matrix-printing procedures PM or PMF must be used to correctly print a matrix read in with RM.

This example reads a 2-by-3 matrix and prints the results using the matrix-printing procedure PM.

RM, a, 2, 3

row 0: 11 22 33
row 1: 40 50 60

; Read a 2-by-3 matrix.

PM, a

11.0000           22.0000           33.0000
40.0000           50.0000           60.0000

; Output the matrix.

In this example, a complex matrix is read. Notice that the elements input as integers are promoted to type complex with the value of the imaginary part set to zero.

RM, a, 3, 2, /Complex

row 0: (1, 0)         (1, 1)
row 1:   1                 -1
row 2: (10.0, 0)   (0, -10)

; Read the matrix; note that keyword Complex is set.

PM, a

(  1.00000, 0.00000)(  1.00000,  1.00000)
(  1.00000, 0.00000)( -1.00000,  0.00000)
( 10.0000,  0.00000)(  0.00000, -10.0000)

; Print the result.

In this example, a 3-by-3 matrix and a 3-by-1 matrix are read. These matrices are then used in a call to the PV-WAVE:IMSL Mathematics Toolkit function LUSOL.

RM, a, 3, 3

row 0: 1 3 3
row 1: 1 3 4
row 2: 1 4 3

; Read the coefficient matrix.

RM, b, 3, 1

row 0:  1
row 1:  4
row 2: -1

; Read the right-hand side.

x = LUSOL(b, a)

; Call LUSOL to compute the solution.

PM, x

-2.00000
-2.00000
  3.00000

; Output the results.

PM, a#x - b

0.00000
0.00000
0.00000

PM, PMF, RMF

See Matrices in Chapter 5 of the PV-WAVE Programmer's Guide for more information.