ObjStore::Lib::PDL - Persistent PDL-compatible matrices

ObjStore-Lib-PDL documentation Contained in the ObjStore-Lib-PDL distribution.

Index

Code Index:

NAME

Top

ObjStore::Lib::PDL - Persistent PDL-compatible matrices

SYNOPSIS

Top

    use PDL::Lite;
    use ObjStore::Lib::PDL;

    begin 'update', sub {
      my $pdl = ObjStore::Lib::PDL->new($near,
			  { Datatype => PDL::float(), Dims => [3,3] });

      $pdl->slice(":,4")->clump(2) *= 2;  #or whatever
    };
    die if $@;

DESCRIPTION

Top

The main thing of interest is that dimensions are arranged in memory such that the first dimension is the most packed. For example, in a PDL of dimensions [2,3] the layout is as follows:

  [0,1]
  [2,3]
  [4,5]

Be aware that this memory layout convention is dependent on the implementation of PDL. It is, however, very unlikely to change.

SEE ALSO

Top

PDL

AUTHOR

Top

Copyright © 1999 Joshua Nathaniel Pritikin. All rights reserved.

This package is free software and is provided "as is" without express or implied warranty. It may be used, redistributed and/or modified under the terms of the Perl Artistic License (see http://www.perl.com/perl/misc/Artistic.html)

ObjStore-Lib-PDL documentation Contained in the ObjStore-Lib-PDL distribution. 
use strict;
package ObjStore::Lib::PDL;
use Carp;
use ObjStore;
use base ('ObjStore::UNIVERSAL','PDL','DynaLoader');
use vars qw($VERSION @OVERLOAD);
$VERSION = '0.92';
BEGIN {
    require PDL::Lite;
    my @ops;
    for (values %overload::ops) { push @ops, split /\s/ }
    for (@ops) {
	my $meth;
	if ($_ eq '""') {
	    $meth = overload::Method('ObjStore::UNIVERSAL', $_);
	} else {
	    $meth = overload::Method('PDL', $_);
	}
	if ($meth) {
	    push @OVERLOAD, $_, $meth;
	} else {
	    #warn "PDL does not overload '$_'\n";
	}
    }
}
use overload @OVERLOAD;

__PACKAGE__->bootstrap($VERSION);
$ObjStore::SCHEMA{'ObjStore::Lib::PDL'}->
    load($ObjStore::Config::SCHEMA_DBDIR."/Lib-PDL-01.adb");

sub new {
    my ($this, $near, $how) = @_;
    $near = $near->segment_of if ref $near;
    my $class = ref $this || $this;
    my $o = _allocate($class, $near);
    if ($how) {
	if (exists $how->{Datatype}) {
	    my $dt = $how->{Datatype};
	    $dt = $dt->[0] if ref $dt;
	    $o->set_datatype($dt);
	}
	$o->setdims($how->{Dims})
	    if exists $how->{Dims};
    }
    $o;
}

#*initialize = \&PDL::initialize;

sub initialize {
#    my ($class) = @_;
    PDL->initialize();
}

# ObjStore::UNIVERSAL::isa is naughty!
*isa = \&UNIVERSAL::isa;

$ObjStore::STARGATE{PDL} = sub {
    my ($class, $sv, $seg) = @_;
    my $pt = $class eq 'PDL' ? 'ObjStore::Lib::PDL' : $class;
    my @dims = $sv->dims;
    my $dt = $sv->get_datatype;
    my $spec = { Datatype => $dt, Dims => \@dims };
    my $o = $pt->new($seg, $spec);
    $o .= $sv;
    $o;
};

END { _PurgeFreelist() }

1;