Matrix(3pm) User Contributed Perl Documentation Matrix(3pm)

# NAME¶

PDL::Matrix -- a convenience matrix class for column-major access

# VERSION¶

This document refers to version PDL::Matrix 0.5 of PDL::Matrix

# SYNOPSIS¶

```  use PDL::Matrix;
\$m = mpdl [[1,2,3],[4,5,6]];
\$m = PDL::Matrix->pdl([[1,2,3],[4,5,6]]);
\$m = msequence(4,3);
@dimsa = \$x->mdims; # 'dims' is not overloaded
\$v = vpdl [0,1,2,3]
\$v = vzeroes(4);
```

# DESCRIPTION¶

## Overview¶

This package tries to help people who want to use PDL for 2D matrix computation with lots of indexing involved. It provides a PDL subclass so one- and two-dimensional ndarrays that are used as vectors resp and matrices can be typed in using traditional matrix convention.

If you want to know more about matrix operation support in PDL, you want to read PDL::MatrixOps or PDL::Slatec.

The original pdl class refers to the first index as the first row, the second index as the first column of a matrix. Consider

```  print \$B = sequence(3,2)
[
[0 1 2]
[3 4 5]
]
```

which gives a 2x3 matrix in terms of the matrix convention, but the constructor used (3,2). This might get more confusing when using slices like sequence(3,2)->slice("1:2,(0)") : with traditional matrix convention one would expect [2 4] instead of [1 2].

This subclass PDL::Matrix overloads the constructors and indexing functions of pdls so that they are compatible with the usual matrix convention, where the first dimension refers to the row of a matrix. So now, the above example would be written as

```  print \$B = PDL::Matrix->sequence(3,2) # or \$B = msequence(3,2)
[
[0 1]
[2 3]
[4 5]
]
```

Routines like eigens or inv can be used without any changes.

Furthermore one can construct and use vectors as n x 1 matrices without mentioning the second index '1'.

## Implementation¶

"PDL::Matrix" works by overloading a number of PDL constructors and methods such that first and second args (corresponding to first and second dims of corresponding matrices) are effectively swapped. It is not yet clear if PDL::Matrix achieves a consistent column-major look-and-feel in this way.

As of version 0.5 (rewrite by CED) the matrices are stored in the usual way, just constructed and stringified differently. That way indexing and everything else works the way you think it should.

 2022-10-26 perl v5.36.0