Ejemplo: El Producto de Matrices

El siguiente ejemplo muestra una subrutina matrixProd que realiza el producto de matrices. Veamos un ejemplo de llamada:

lhp@nereida:~/Lperl/src$ sed -ne '58,64p' matrixproduct.pl | cat -n
     1  my $A = [[1,2,3],[2,4,6],[3,6,9]];
     2  my $B = [[1,2],[2,4],[3,6]];
     3  my $C = matrixProd($A,$B);
     4
     5  print "Matriz A \n"; printMat($A);
     6  print "\nMatriz B \n"; printMat($B);
     7  print "\nMatriz C \n"; printMat($C);

La ejecución del programa produce la siguiente salida:

lhp@nereida:~/Lperl/src$ matrixproduct.pl
Matriz A
1       2       3
2       4       6
3       6       9

Matriz B
1       2
2       4
3       6

Matriz C
14      28
28      56
42      84

Producto de Matrices Usando Referencias

Sigue el código e la función matrixProd:

lhp@europa:~/projects/perl/src/perltesting$  sed -ne '1,26p' matrixproduct.pl | cat -n
 1  #!/usr/bin/perl -w
 2  use strict;
 3  use List::Util qw(sum);
 4  use List::MoreUtils qw(any);
 5  use Scalar::Util qw{reftype};
 6
 7  sub matrixProd {
 8    my ($A, $B) = @_;
 9    die "Error. Se esperaban dos matrices" if any { !defined($_) or reftype($_) ne 'ARRAY' } ($A, $B, @$A, @$B);
10
11    my $nrA = @$A;
12    my $nrB = @$B;
13    my $ncB = @{$B->[0]};
14    my @C;
15
16    die ("Las matrices no son multiplicables\n") if any { @$_ != $nrB } @$A;
17
18    for(my $i = 0; $i < $nrA; $i++)  {
19      for(my $j = 0; $j < $ncB; $j++) {
20        my $k = 0;
21        $C[$i][$j] = sum(map { $A->[$i][$k++] * $_ } map { $_->[$k] } @$B);
22      }
23    }
24
25    return \@C;
26  }

El Módulo PDL

Todo queda mucho mas simple usando el módulo PDL (PDL significa Perl Data Language) el cual da soporte al cálculo científico en Perl:

1	#!/usr/bin/perl -w
2	use PDL;
3	
4	$a = pdl [[1,2,3],[2,4,6],[3,6,9]];
5	$b = pdl [[1,2],[2,4],[3,6]];
6	$c = $a x $b; # x esta sobrecargado
7	
8	print "a = $a\nb = $b\nc = $c\n";

PDL es un módulo diseñado para el cálculo y visualización de datos científicos.

$b->[1] } 22 map { [ $_, -M $_ ] } @_ 23 } 24 25 sub latest { 26 (sort_by_time( @_ ) )[0]; 27 } 28 29 sub dirs { 30 grep { -d $_ } @_; 31 } La principal novedad en este código se refiere al control (línea 11) de que comandos están disponibles. La función can puede ser llamada de la forma:

paquete->can('nombre_de_subrutina')

devuelve cierto si existe una subrutina con ese paquete y falso en caso contrario. De hecho devuelve una referencia (dura/no simbólica) a la subrutina si esta existe.

Recuerde

• En un programa Perl pueden existir varias tablas de símbolos. Una declaración package cambia el espacio de nombres hasta que encontremos una nueva declaración package, o hasta el final del bloque actual.
• De hecho, las variables en Perl se clasifican como package variables y lexical variables.
• Estas últimas, las léxicas, son las declaradas con my.
• Las variables package pertenecen, naturalmente, a un package (normalmente el actual).
• El package inicial es el package main.
• Cuando sea necesario hacer explícito a que package pertenece la variable, puede hacerse prefijando su nombre con el del package, separado por ::.