Foro - Perl en Español

por **LAS_GC** » 2011-07-22 06:51 @327

Hola a todos,

Os escribo porque no tengo mucha experiencia en Perl y necesito optimizar una función, por ejemplo, f(x1,x2)=x1**2+x2**2 a partir de (4,5). He encontrado la librería SIMPLEX, pero no consigo que me funcione y creo que la culpa es del vector $init.

¿Alguien puede ayudarme?

Información sobre PDL::Opt::Simplex
http://search.cpan.org/~chm/PDL-2.4.9_0 ... Simplex.pm

Esto es lo que he escrito:

Sintáxis: [ Descargar ] [ Ocultar ]

Using perl Syntax Highlighting

#!/usr/bin/perl
 
use PDL;
use PDL::Opt::Simplex;
 
$init = ( 4, 5 );
 
$initsize = 2;
$minsize  = 1e-6;
$maxiter  = 500;
 
sub evaluate_func_at {
    my ($xv) = @_;
    my $x1   = $xv->slice("(0)");
    my $x2   = $xv->slice("(1)");
    return $x1**2 + $x2**2;
}
 
( $optimum, $ssize, $optval ) = simplex(
    $init, $initsize, $minsize, $maxiter,
    sub { evaluate_func_at( $_[0] ) },
    sub { display_simplex( $_[0] ) }
);
Coloreado en 0.004 segundos,  usando GeSHi 1.0.8.4

Error: Can't call method "getdim" on an undefined value at /usr/lib64/perl5/PDL/Opt/Simplex.pm line 116.

Muchas gracias de antemano
Verónica

Y vemos lo que pasa: el $init que le has pasado a simplex() debería ser un objeto pdl, pero no lo es. Has creado el $init como una referencia a un array de dos elementos, pero no como objeto pdl.

Prueba a cambiar la línea 6 por esta:

$init = pdl( 4, 5 );

Tienes más información en la sección 'pdl' de la página de manual de PDL::Core.

A propósito, me parece que tienes una versión antigua de PDL. Al menos, no coincide con la versión v2.4.9_003 que está en CPAN.

por **LAS_GC** » 2011-07-22 07:31 @354

Hola de nuevo y muchas gracias, explorer, por responder tan rápido.

He probado lo que me has dicho y me sale otro error:
Undefined subroutine &main::display_simplex called at ./ejemplo.pl line 21.

Este error es lógico porque no tengo esta subrutina añadida sub {display_simplex($_[0])) y no tengo ni idea de cómo hacerla

En la ayuda no he visto nada sobre esto. Igual es tan básico que ni se molestaron.

Yo he utilizado en C++ esta librería y normalmente el vector de entrada inicial son N+1 valores para cada variable. Pero aquí parece ser que solo es uno. Me refiero por ejemplo a x1=(4,1,3) y x2=(5,2,6).
$init=[(4,5) (1,2) (3,6)]. ¿A lo mejor el error viene de ahí?

Saludos
Verónica

por **explorer** » 2011-07-22 11:32 @522

Según el código del módulo, esa subrutina es opcional, así que la puedes quitar.

Sintáxis: [ Descargar ] [ Ocultar ]

Using perl Syntax Highlighting

( $optimum, $ssize, $optval ) = simplex(
    $init, $initsize, $minsize, $maxiter,
    sub { evaluate_func_at( $_[0] ) },
);
Coloreado en 0.001 segundos,  usando GeSHi 1.0.8.4

Y sí, la documentación de PDL es muy escasa... casi la única manera de obtener asistencia -casi- inmediata es uniéndose a su lista de correo.

En cuanto a las dimensiones de $init, según el código (te aconsejo que, al menos, lo mires), si el cardinal de la segunda dimensión de $init es 1, entonces lo resuelve de la forma clásica (construir un tetraedro, bucle por los valores, y reduciendo el tamaño del tetraedro hasta encontrar el punto de equilibrio). Si el cardinal de la segunda dimensión es igual al cardinal de la primera dimensión más uno, entonces aplica otro método más genérico (para más dimensiones). Si no se cumple la condición de $nd2 == $nd +1, entonces sale sin calcular nada.

Te aconsejo que te unas a la lista de correo de PDL. Es en inglés, pero hay gente muy buena, y suele dar respuestas rápidas.

por **LAS_GC** » 2011-07-25 07:22 @348

Hola de nuevo,
Con tu ayuda ya lo he resuelto. He quitado esa subrutina y lo que hago es un simple print() de x1, x2 y la función, que al final es lo que me interesa ver.

En cuanto a $initsize, fue un poco más complicado. Con este método se dibuja inicialmente un triángulo. Para cada uno de los tres vértices se calcula la función y el punto de mayor valor se cambia (a través de expansiones). Con $init le damos el centro de este triángulo y se utiliza $initsize para encontrar los tres vértices.

Muchas gracias de nuevo. Gracias por tus consejos.
Saludos

LAS_GC

--------------------------------------------------

El script final es:

Sintáxis: [ Descargar ] [ Ocultar ]

Using perl Syntax Highlighting

#!/usr/bin/perl
 
use PDL;
use PDL::Opt::Simplex;
 
my $Bext   = -1e-3;
my $theta0 = -0.20;
my $DB     = 0.9e-3;
my $Dt     = 0.02;
 
my $init = pdl( $Bext, $theta0 );
 
my $initsize = pdl( $DB, $Dt );
my $minsize  = 1e-6;
my $maxiter  = 500;
 
sub evaluate_func_at {
    my ($xv) = @_;
 
    my $x1 = $xv->slice("(0)");
    print "x1 = $x1 \n";
    my $x2 = $xv->slice("(1)");
    print "x2 = $x2 \n";
 
    $rfx = $rx1**2 + $rx2**2;
    print "F(x1,x2) = $rfx \n";
 
    return $rfx;
 
}
 
( $optimum, $ssize, $optval )
    = simplex( $init, $initsize, $minsize, $maxiter,
    sub { evaluate_func_at( $_[0] ) } );
 
print "OPTIMUM = $optimum \n";
print "SSIZE   = $ssize\n";
print "OPTVAL  = $optval\n";
Coloreado en 0.002 segundos,  usando GeSHi 1.0.8.4

Foro - Perl en Español

Error getdim en PDL::Opt::Simplex

Error getdim en PDL::Opt::Simplex

Publicidad

Re: Error getdim en PDL::Opt::Simplex

Re: Error getdim en PDL::Opt::Simplex

Re: Error getdim en PDL::Opt::Simplex

Re: Error getdim en PDL::Opt::Simplex

¿Quién está conectado?