NAME

Math::Matrix::MaybeGSL - Uniform use of Math::MatrixReal and Math::GSL::Matrix.

VERSION

version 0.008

SYNOPSIS

use Math::Matrix::MaybeGSL;

my $matrix = Matrix->new(3, 4);

# puts first position of matrix with value 10
$matrix->assign(1, 1, 10);

# gets last position of matrix (should hold 0)
my $l = $matrix->element(3, 4);

DESCRIPTION

This module interfaces with Math::GSL::Matrix or, if that is not available, Math::MatrixReal. The idea behind this module is to allow the development of tools that use matrices that will work in pure Perl (using Math::MatrixReal) or with extra efficiency using Math::GSL::Matrix.

Given the two modules API is quite distinct, this module defines its own API, stealing method names from both these modules.

METHODS

Matrix

This is a simple function that returns this package name: Math::Matrix::MaybeGSL. It allows a simple interface as shown below for the constructors.

isGSL

Returns a true value is running over Math::GSL backend.

if (Matrix->isGSL) { ... }

new

Construct a new matrix object. Receives as arguments the number of rows and columns of the matrix being created.

my $matrix = Matrix->new(20, 30);

Yes, although the module name is Math::Matrix::MaybeGSL, the Matrix subroutine will make it easier to use (shorter name).

new_from_cols

Receives a nested list with the matrix elements, one column at a time:

my $matrix = Matrix->new_from_cols( [[1, 2], [3, 4]]);

returns  [ 1 3 ]
         [ 2 4 ]

new_from_rows

Receives a nested list with the matrix elements, one row at a time:

my $matrix = Matrix->new_from_rows( [[1, 2], [3, 4]]);

returns  [ 1 2 ]
         [ 3 4 ]

dim

Returns a list (a pair) with the number of lines and columns of the matrix.

my ($rows, $columns) = $matrix->dim();

assign

Sets a value in a specific position. Note that indexes start at 1 unlike Perl and some other programming languages.

# sets the first element of the matrix to 10.
$matrix->assign(1, 1, 10);

element

Retrieves a value from a specific position of the matrix. Note that indexes start at 1 unlike Perl and some other programming languages.

# retrieves the second element of the first row
my $val = $matrix->element(1, 2);

each

Apply a specific function to every element of the matrix, returning a new one.

# square all elements
$squared_matrix = $matrix->each( { shift ** 2 } );

hconcat

Concatenates two matrices horizontally. Note they must have the same number of rows.

$C = $a->hconcat($b);

if A = [ 1 2 ]  and B = [ 5 6 ]  then C = [ 1 2 5 6 ]
       [ 3 4 ]          [ 7 8 ]           [ 3 4 7 8 ]

vconcat

Concatenates two matrices horizontally. Note they must have the same number of rows.

$C = $a->vconcat($b);

if A = [ 1 2 ]  and B = [ 5 6 ]  then C = [ 1 2 ]
       [ 3 4 ]          [ 7 8 ]           [ 3 4 ]
                                          [ 5 6 ]
                                          [ 7 8 ]

max

Returns the maximum value of the matrix. In scalar context the position is also returned. For vectors (matrices whose number of rows or columns is 1) only a position value is returned.

$max = $matrix->max();
($max, $row, $col) = $matrix->max();

min

Returns the minimum value of the matrix. In scalar context the position is also returned. For vectors (matrices whose number of rows or columns is 1) only a position value is returned.

$min = $matrix->min();
($min, $row, $col) = $matrix->min();

det

Returns the determinant of the matrix, without going through the rigamarole of computing a LR decomposition.

as_list

Get the contents of a matrix instance as a Perl list.

write

Given a matrix and a filename, writes that matrix to the file. Note that if the file exists it will be overwritten. Also, files written by Math::GSL will not be compatible with files written by Math::MatrixReal nor vice-versa.

$matrix->write("my_matrix.dat");

read

Reads a matrix written by the write method. Note that it will only read matrices written by the same back-end that is being used for reading.

my $matrix = Matrix->load("my_matrix.dat");

row

Returns the selected row in a matrix as a new matrix object. Note that indexes start at 1 unlike Perl and some other programming languages.

my $row = $matrix->row(1);

find_zeros

Given a matrix, returns a nested list of indices corresponding to zero values in the given matrix. Note that indexes start at 1 unlike Perl and some other programming languages.

my @indices = $matrix->find_zeros();

transpose

Returns transposed matrix.

OVERLOAD

For now only matrix multiplication and addition are overloaded, in the usual operators, * and +, correspondingly. Take attention that these operations only work if the matrix dimensions are compatible.

$m = $a * $b;
$n = $a + $b;

BUGS

At this initial stage of this module, only the methods that I am really needing for my depending applications are implemented. Therefore, it might not include the method that you were looking for. Nevertheless, send me an e-mail (or open an issue on GitHub) and I'll be happy to include it (given the two modules support it).

SEE ALSO

Check Math::MatrixReal and Math::GSL::Matrix documentation.

CONTRIBUTORS

  • Andrius Merkys <merkys@cpan.org>

  • Ivan Baidakou

  • Gabor Szabo

AUTHOR

Alberto Simões <ambs@cpan.org>

COPYRIGHT AND LICENSE

This software is copyright (c) 2014-2023 by Alberto Simões.

This is free software; you can redistribute it and/or modify it under the same terms as the Perl 5 programming language system itself.