http://www.perlmonks.org?node_id=1197401


in reply to Re^15: High Performance Game of Life (updated)
in thread High Performance Game of Life

Update 1: Added results for the original version and initial improvements.
Update 2: Added results from tbench1-infinite.pl using Game::Life::Infinite::Board.
Update 3: Added results for C++, Organism.h and tbench1.cpp.

In this thread, we've learned that pack i2 is faster than pack ii. Furthermor, functions having inner loops benefit greatly by inlining critical paths. In essence this is what we've done. The pack i2 solution fully optimized is found here, by tybalt89. The mapping of two numbers into one can be found here. Similarly, a shorter implementation by tybalt89 is found here. The fastest time was noted for each run with nothing else running in the background.

The maximum key length for $c, obtained separately, is 7, 12, and 18 for pack i2, the mapping of two numbers, and the shorter solution by tybalt89, respectively. The fix for the latter solution, when running cperl on 64-bit hardware, is using 30 bits instead of 32, described here.

The x and y tmp files were made using eyepopslikeamosquito's createblinker.pl script, found at the top of this thread.

$ perl createblinker.pl 500000 -900000 100 >x.tmp 2>y.tmp

[download]

Benchmark results were captured on a 2.6 GHz laptop (i7 Haswell). I apologize for not having something older to benchmark on. If anything, the OP's machine is faster than mine. ;-)

Results:

c++              138 MB    4 secs

               mem size   v5.26.0   v5.24.2   v5.22.4

infinite       7,548 MB  122 secs  119 secs  133 secs
original         704 MB  167 secs  168 secs  180 secs
improvements     744 MB   57 secs   57 secs   63 secs
optimized i2   1,543 MB   38 secs   39 secs   40 secs
2 nums into 1  1,510 MB   39 secs   40 secs   42 secs
shorter impl   1,661 MB   40 secs   42 secs   44 secs

[download]

cperl:

$ /opt/cperl-5.24.3c/bin/cperl -I. tbench1.pl x.tmp 2
cell count at start = 1500000
run benchmark for 2 ticks
cell count at end = 1500000
time taken: 116 secs   - infinite board
time taken: 163 secs   - original
time taken:  54 secs   - improvements
time taken:  36 secs   - optimized pack i2
time taken:  37 secs   - two numbers into one
time taken:  37 secs   - shorter implementation

[download]

Pack returns unreadable data, but is fast. Readable keys may be preferred for storing into a DB. Stringification "$x:$y" is one way. Unfortunately, that requires split to extract the values and a text field versus numeric if storing into a DB. Bit manipulation is another way.

Running Game::Life::Infinite::Board consumes lots of memory. If possible, check for 8 ~ 10 gigabytes of available memory to minimize the OS from swapping. It also takes ~ 10 seconds during global cleanup while exiting.

eyepopslikeamosquito, infinite isn't running ~ 2x slower as reported here for 1.5 million cells. My laptop has 1600 MHz RAM and verified available memory before running.

Regards, Mario

Replies are listed 'Best First'.
Re^17: High Performance Game of Life (updated - results)
by tybalt89 (Monsignor) on Aug 15, 2017 at 22:06 UTC

    New entry.

    Hey, it passes all the tests :)

    hehehe 

    package Organism;

# based on http://perlmonks.org/?node_id=1197284

use strict;
use warnings;

sub count
  {
  return shift->{config}[0] =~ tr/1//;
  }

# Input a list of [ x, y ] coords
sub insert_cells
  {
  my $extra = 3;
  my $self = shift;
  my $xl = my $xh = $_[0][0]; # find cell limits
  my $yl = my $yh = $_[0][1];
  for (@_)
    {
    my ($x, $y) = @$_;
    $xl > $x and $xl = $x;
    $xh < $x and $xh = $x;
    $yl > $y and $yl = $y;
    $yh < $y and $yh = $y;
    }
  my $xoffset = $xl - $extra; # get sizes and insert live cells
  my $w = $xh - $xl + 2 * $extra;
  my $yoffset = $yl - $extra;
  my $h = $yh - $yl + 2 * $extra;
  my $grid = '0' x $w x $h;
  for (@_)
    {
    my ($x, $y) = @$_;
    substr $grid, $x - $xoffset + ($y - $yoffset) * $w, 1, '1';
    }
  $self->{config} = [ $grid, $w, $h, $xoffset, $yoffset ];
  }

# Return sorted list of cells in the Organism.
# Used for verification and testing the state of the organism.
sub get_live_cells
  {
  my $self = shift;
  my ( $grid, $w, $h, $xoffset, $yoffset ) = @{ $self->{config} };
  my @cells;
  push @cells, [ $-[0] % $w + $xoffset, int( $-[0] / $w ) + $yoffset ]
    while $grid =~ /1/g;
  sort { $a->[0] <=> $b->[0] || $a->[1] <=> $b->[1] } @cells;
  }

sub tick
  {
  my $self = shift;
  my ( $grid, $w, $h ) = @{ $self->{config} };
  my $all = '0' x ($w + 1) . $grid;
  my $sum = $all =~ tr/1/2/r;
  ( $sum |= substr $all, $_ ) =~ tr/1357/2468/ for
    1, 2, $w, $w + 2, $w * 2, $w * 2 + 1, $w * 2 + 2; # other 7 neighb
+ors
  $grid = substr $grid | $sum, 0, $w * $h;
  $self->{config}[0] = $grid =~ tr/1-9/000011100/r; # dead or alive
  }

sub new {
   my $class = shift;
   my %init_self = ( );
   bless \%init_self, $class;
}

1;

    [download]
[reply]
[d/l]

      Hi tybalt89,

      Wow! Your new entry runs faster than C++. Also, memory consumption is less than 500 MB ;-)

      $ perl createblinker.pl 500000 -900000 100 >x.tmp 2>y.tmp
$ g++ -o tbench1 -std=c++11 -Wall -O3 tbench1.cpp

$ time ./tbench1 x.tmp 2
cell count at start = 1500000
run benchmark for 2 ticks
cell count at end = 1500000
time taken 4 secs

real  0m5.240s   mem 139 MB
user  0m5.149s
sys   0m0.085s

$ time /opt/perl-5.26.0/bin/perl -I. tbench1.pl x.tmp 2
cell count at start = 1500000
run benchmark for 2 ticks
cell count at end = 1500000
time taken: 1 secs

real  0m3.482s   mem 492 MB
user  0m3.242s
sys   0m0.233s

      [download]

      Micro-optimization may be a subjective matter. At this level, one may want to for 2%.

      I've replaced 3 multiplications ( $w * 2 ) with ( $w << 1 ).

      ( $sum |= substr $all, $_ ) =~ tr/1357/2468/ for
  1, 2, $w, $w + 2, ($w << 1), ($w << 1) + 1, ($w << 1) + 2; # other 7
+ neighbors

      [download]
      $ time /opt/perl-5.26.0/bin/perl -I. tbench1.pl x.tmp 2
cell count at start = 1500000
run benchmark for 2 ticks
cell count at end = 1500000
time taken: 1 secs

real  0m3.420s   mem 492 MB
user  0m3.203s
sys   0m0.205s

      [download]

      Regards, Mario

[reply]
[d/l]
[select]
[reply]

        Thanks for the new test.

        Here's my "new entry" fixed to pass it. 

        package Organism;

# based on http://perlmonks.org/?node_id=1197284

use strict;
use warnings;

sub count
  {
  return shift->{config}[0] =~ tr/1//;
  }

# Input a list of [ x, y ] coords
sub insert_cells
  {
  my $extra = 3;
  my $self = shift;
  my $xl = my $xh = $_[0][0]; # find cell limits
  my $yl = my $yh = $_[0][1];
  for (@_)
    {
    my ($x, $y) = @$_;
    $xl > $x and $xl = $x;
    $xh < $x and $xh = $x;
    $yl > $y and $yl = $y;
    $yh < $y and $yh = $y;
    }
  my $xoffset = $xl - $extra; # get sizes and insert live cells
  my $w = $xh - $xl + 2 * $extra;
  my $yoffset = $yl - $extra;
  my $h = $yh - $yl + 2 * $extra;
  my $grid = '0' x $w x $h;
  for (@_)
    {
    my ($x, $y) = @$_;
    substr $grid, $x - $xoffset + ($y - $yoffset) * $w, 1, '1';
    }
  $self->{config} = [ $grid, $w, $h, $xoffset, $yoffset ];
  }

# Return sorted list of cells in the Organism.
# Used for verification and testing the state of the organism.
sub get_live_cells
  {
  my $self = shift;
  my ( $grid, $w, $h, $xoffset, $yoffset ) = @{ $self->{config} };
  my @cells;
  push @cells, [ $-[0] % $w + $xoffset, int( $-[0] / $w ) + $yoffset ]
    while $grid =~ /1/g;
  sort { $a->[0] <=> $b->[0] || $a->[1] <=> $b->[1] } @cells;
  }

sub tick
  {
  my $self = shift;
  my ( $grid, $w, $h, $xoffset, $yoffset ) = @{ $self->{config} };
  # expand
  $grid = join '00', '0' x ($w + 1), unpack("(a$w)*", $grid), '0' x ($
+w + 1);
  $w += 2;
  $h += 2;
  $xoffset--;
  $yoffset--;
  # now get new generation
  my $all = '0' x ($w + 1) . $grid;
  my $sum = $all =~ tr/1/2/r;
  ( $sum |= substr $all, $_ ) =~ tr/1357/2468/ for
    1, 2, $w, $w + 2, $w * 2, $w * 2 + 1, $w * 2 + 2; # other 7 neighb
+ors
  $grid = substr $grid | $sum, 0, $w * $h;
  $self->{config} = [ $grid =~ tr/1-9/000011100/r, $w, $h, $xoffset, $
+yoffset ];
  }

sub new {
   my $class = shift;
   my %init_self = ( );
   bless \%init_self, $class;
}

1;

        [download]

        Tweak a couple of lines and add only 5 new ones and it's fixed. I really love perl!

[reply]
[d/l]
Re^17: High Performance Game of Life (updated - results)
by eyepopslikeamosquito (Archbishop) on Aug 15, 2017 at 22:38 UTC

    Re having an old PC - and rubbing salt into tybalt89's wounds - please note that my (Haswell 4770K-powered) home PC is unchanged from when I mentioned it over three years ago in The 10**21 Problem (Part 2) and so is getting rather old and creaky now. Definitely due for an upgrade. :)

    Though my ancient PC does have 32 GB of memory (and so there is no swapping), given my very slow benchmark results with Game::Life::Infinite::Board, I'm suspicious that my doddering, ancient memory is starting to fail, due to being so old. This may be contributing to the significantly different results I am seeing.

    When running: 

    perl tbench1.pl x.tmp 2

    [download]
    with my originally shortened Organism.pm: 
    package Organism;

use strict;
# use warnings;

sub count {
   return scalar keys %{ shift->{Cells} };
}

# Input a list of [ x, y ] coords
sub insert_cells {
   my $cells = shift->{Cells};
   for my $r (@_) { $cells->{ pack 'i2', @{$r} } = undef }
}

# Return sorted list of cells in the Organism.
# Used for verification and testing the state of the organism.
sub get_live_cells {
   sort { $a->[0] <=> $b->[0] || $a->[1] <=> $b->[1] }
      map { [ unpack 'i2', $_ ] }
      keys %{ shift->{Cells} };
}

sub tick {
   my $self  = shift;
   my $cells = $self->{Cells};
   my ( $k1, $k2, $k3, $k4, $k5, $k6, $k7, $k8,
        $x0, $x1, $x2, $y0, $y1, $y2, %new_cells );

   for my $c (keys %{ $cells }) {
      # Get the (up to 8) dead cells surrounding the cell
      ( $x0, $y0 ) = unpack 'i2', $c;
      ( $x1, $x2, $y1, $y2 ) = ( $x0 - 1, $x0 + 1, $y0 - 1, $y0 + 1 );
      my @zcells = (
         ($k1 = pack 'i2', $x1, $y1) x !(exists $cells->{$k1}),
         ($k2 = pack 'i2', $x1, $y0) x !(exists $cells->{$k2}),
         ($k3 = pack 'i2', $x1, $y2) x !(exists $cells->{$k3}),
         ($k4 = pack 'i2', $x0, $y1) x !(exists $cells->{$k4}),
         ($k5 = pack 'i2', $x0, $y2) x !(exists $cells->{$k5}),
         ($k6 = pack 'i2', $x2, $y1) x !(exists $cells->{$k6}),
         ($k7 = pack 'i2', $x2, $y0) x !(exists $cells->{$k7}),
         ($k8 = pack 'i2', $x2, $y2) x !(exists $cells->{$k8}) );

      # Check the live cell
      # Note: next line equivalent to nlive == 2 || nlive == 3
      @zcells == 5 || @zcells == 6 and $new_cells{$c} = undef;

      # Check the dead cells
      for my $z (@zcells) {
         ( $x0, $y0 ) = unpack 'i2', $z;
         ( $x1, $x2, $y1, $y2 ) = ( $x0 - 1, $x0 + 1, $y0 - 1, $y0 + 1
+ );
         # Get num live
         (   ( exists $cells->{ pack 'i2', $x1, $y1 } )
           + ( exists $cells->{ pack 'i2', $x1, $y0 } )
           + ( exists $cells->{ pack 'i2', $x1, $y2 } )
           + ( exists $cells->{ pack 'i2', $x0, $y1 } )
           + ( exists $cells->{ pack 'i2', $x0, $y2 } )
           + ( exists $cells->{ pack 'i2', $x2, $y1 } )
           + ( exists $cells->{ pack 'i2', $x2, $y0 } )
           + ( exists $cells->{ pack 'i2', $x2, $y2 } )
         ) == 3 and $new_cells{$z} = undef;
      }
   }

   $self->{Cells} = \%new_cells;
}

sub new {
   my $class = shift;
   my %init_self = ( Cells => {} );
   bless \%init_self, $class;
}

1;

    [download]
    I get 62 seconds, while after making the excellent improvements suggested by tybalt89 I get 80 seconds. 
    package Organism;

use strict;
# use warnings;

sub count {
   return scalar keys %{ shift->{Cells} };
}

# Input a list of [ x, y ] coords
sub insert_cells {
   my $cells = shift->{Cells};
   for my $r (@_) { $cells->{ pack 'i2', @{$r} } = undef }
}

# Return sorted list of cells in the Organism.
# Used for verification and testing the state of the organism.
sub get_live_cells {
   sort { $a->[0] <=> $b->[0] || $a->[1] <=> $b->[1] }
      map { [ unpack 'i2', $_ ] }
      keys %{ shift->{Cells} };
}

sub tick {
   my $self  = shift;
   my $cells = $self->{Cells};
   my ( $k1, $k2, $k3, $k4, $k5, $k6, $k7, $k8,
        $x0, $x1, $x2, $y0, $y1, $y2, %new_cells, %dead_cells );

   for my $c (keys %{ $cells }) {
      # Get the (up to 8) dead cells surrounding the cell
      ( $x0, $y0 ) = unpack 'i2', $c;
      ( $x1, $x2, $y1, $y2 ) = ( $x0 - 1, $x0 + 1, $y0 - 1, $y0 + 1 );
      $dead_cells{$_}++ for my @zcells = (
         ($k1 = pack 'i2', $x1, $y1) x !(exists $cells->{$k1}),
         ($k2 = pack 'i2', $x1, $y0) x !(exists $cells->{$k2}),
         ($k3 = pack 'i2', $x1, $y2) x !(exists $cells->{$k3}),
         ($k4 = pack 'i2', $x0, $y1) x !(exists $cells->{$k4}),
         ($k5 = pack 'i2', $x0, $y2) x !(exists $cells->{$k5}),
         ($k6 = pack 'i2', $x2, $y1) x !(exists $cells->{$k6}),
         ($k7 = pack 'i2', $x2, $y0) x !(exists $cells->{$k7}),
         ($k8 = pack 'i2', $x2, $y2) x !(exists $cells->{$k8}) );

      # Check the live cell
      # Note: next line equivalent to nlive == 2 || nlive == 3
      @zcells == 5 || @zcells == 6 and $new_cells{$c} = undef;
   }
   $dead_cells{$_} == 3 and $new_cells{$_} = undef for keys %dead_cell
+s;
   $self->{Cells} = \%new_cells;
}

sub new {
   my $class = shift;
   my %init_self = ( Cells => {} );
   bless \%init_self, $class;
}

1;

    [download]

    Update: These tests were run with perl v5.24.0: 

    This is perl 5, version 24, subversion 0 (v5.24.0) built for MSWin32-x
+64-multi-thread

    [download]
    Need to upgrade that ancient Perl. :)

    Update: After upgrading to perl v5.26.0, both are faster: (62 secs, 80 secs) improved to (53 secs, 71 secs).

[reply]
[d/l]
[select]

      Hi eyepopslikeamosquito,

      Oh I wished I had something older to run on. Unfortunately, my laptop is the slowest machine I have. The Haswell/Crystalwell chip has 128MB of eDRAM. It's purpose is for the GPU. But, the 4 cores has access to it as well when the GPU is idle, which is the case while benchmarking Perl. The late 2013 i7-4960HQ CPU runs at 2.6 GHz with turbo-boost at 3.8 GHz. I've forgotten that it could run that high on one core with nothing else running.

      I have some news to share. It's possible for tybalt89's excellent optimization to run faster. The extra optimization shaves 6 ~ 7 seconds. Memory consumption reduced as well. I applied the same optimization to the two bit implementations, mine and tybalt89's shorter implementation. For the latter, also applied a change mentioned earlier to not fail with cperl.

      # use 30-bits on 64-bit hw for cperl compatibility
my $half = ( ( log(~0 + 1) / log(2) ) >= 64 ) ? 30 : 16;

      [download]

      Before:

      $ perl createblinker.pl 500000 -900000 100 >x.tmp 2>y.tmp
$ perl tbench1.pl x.tmp 2

               mem size   v5.26.0   v5.24.2   v5.22.4   cperl
infinite       7,548 MB  122 secs  119 secs  133 secs  116 secs
original         704 MB  167 secs  168 secs  180 secs  163 secs
improvements     744 MB   57 secs   57 secs   63 secs   54 secs
optimized i2   1,543 MB   38 secs   39 secs   40 secs   36 secs
2 nums into 1  1,510 MB   39 secs   40 secs   42 secs   37 secs
shorter impl   1,661 MB   40 secs   42 secs   44 secs   37 secs

      [download]

      After:

                     mem size   v5.26.0   v5.24.2   v5.22.4   cperl
optimized i2   1,326 MB   31 secs   32 secs   33 secs   30 secs
2 nums into 1  1,292 MB   33 secs   34 secs   35 secs   31 secs
shorter impl   1,443 MB   35 secs   36 secs   37 secs   31 secs

      [download]

      The optimization was made to "Check dead cells".

      package Organism;

use strict;
# use warnings;

sub count {
   return scalar keys %{ shift->{Cells} };
}

# Input a list of [ x, y ] coords
sub insert_cells {
   my $cells = shift->{Cells};
   for my $r (@_) { $cells->{ pack 'i2', @{$r} } = undef }
}

# Return sorted list of cells in the Organism.
# Used for verification and testing the state of the organism.
sub get_live_cells {
   sort { $a->[0] <=> $b->[0] || $a->[1] <=> $b->[1] }
      map { [ unpack 'i2', $_ ] }
      keys %{ shift->{Cells} };
}

sub tick {
   my $self  = shift;
   my $cells = $self->{Cells};
   my ( $k1, $k2, $k3, $k4, $k5, $k6, $k7, $k8,
        $x0, $x1, $x2, $y0, $y1, $y2, %new_cells, %dead_cells );

   for my $c (keys %{ $cells }) {
      # Get the (up to 8) dead cells surrounding the cell
      ( $x0, $y0 ) = unpack 'i2', $c;
      ( $x1, $x2, $y1, $y2 ) = ( $x0 - 1, $x0 + 1, $y0 - 1, $y0 + 1 );

      my @zcells = (
         ($k1 = pack 'i2', $x1, $y1) x !(exists $cells->{$k1}),
         ($k2 = pack 'i2', $x1, $y0) x !(exists $cells->{$k2}),
         ($k3 = pack 'i2', $x1, $y2) x !(exists $cells->{$k3}),
         ($k4 = pack 'i2', $x0, $y1) x !(exists $cells->{$k4}),
         ($k5 = pack 'i2', $x0, $y2) x !(exists $cells->{$k5}),
         ($k6 = pack 'i2', $x2, $y1) x !(exists $cells->{$k6}),
         ($k7 = pack 'i2', $x2, $y0) x !(exists $cells->{$k7}),
         ($k8 = pack 'i2', $x2, $y2) x !(exists $cells->{$k8}) );

      # Check the live cell
      # Note: next line equivalent to nlive == 2 || nlive == 3
      @zcells == 5 || @zcells == 6 and $new_cells{$c} = undef;

      # Check the dead cells
      for my $z ( @zcells ) {
         $new_cells{$z} = undef if ++$dead_cells{$z} == 3;
      }
   }

   $self->{Cells} = \%new_cells;
}

sub new {
   my $class = shift;
   my %init_self = ( Cells => {} );
   bless \%init_self, $class;
}

1;

      [download]

      Regards, Mario

[reply]
[d/l]
[select]

        I have some news to share. It's possible for tybalt89's excellent optimization to run faster. ... The optimization was made to "Check dead cells".

        I don't believe your optimization of tybalt89's original code is correct. At least, it does not pass my more rigorous tgol3.t test program, shown below.

        Please test your solutions for correctness against this more thorough test (in addition to the original simple tgol.t blinker test). 

        # tgol3.t - Simple lidka test of Conway Game of Life Organism class
use strict;
use warnings;
use Organism;
use Test::More;
my $nticks = 100;
my $ntests = 3 + ( $nticks + 1) * 2;
plan tests => $ntests;

sub test_one {
   my $org = shift;           # Organism handle
   my $desc = shift;          # Test description
   my $nexpected = shift;     # Expected cell count
   # my $expected = shift;    # Array ref of (sorted) expected cells
   my @cells = $org->get_live_cells();
   my $ncells = $org->count();
   cmp_ok( $ncells, '==', $nexpected, "$desc cell count ($ncells)" );
   cmp_ok( scalar(@cells), '==', $nexpected, "$desc cell array count" 
+);
   # is_deeply( \@cells, $expected, "$desc cell array" );
}

# Test first 100 ticks from famous lidka methuselah
# See: http://conwaylife.com/wiki/Lidka
my @lidka_ticks = (
   13,    # [ 0] initial count

   15,    # [ 1]
   15,    # [ 2]
   19,    # [ 3]
   19,    # [ 4]
   23,    # [ 5]
   23,    # [ 6]
   32,    # [ 7]
   29,    # [ 8]
   47,    # [ 9]
   27,    # [10]

   32,    # [11]
   36,    # [12]
   42,    # [13]
   48,    # [14]
   48,    # [15]
   46,    # [16]
   60,    # [17]
   54,    # [18]
   56,    # [19]
   64,    # [20]

   86,    # [21]
   64,    # [22]
   74,    # [23]
   70,    # [24]
   68,    # [25]
   52,    # [26]
   58,    # [27]
   50,    # [28]
   44,    # [29]
   50,    # [30]

   54,    # [31]
   80,    # [32]
   50,    # [33]
   54,    # [34]
   56,    # [35]
   54,    # [36]
   62,    # [37]
   50,    # [38]
   58,    # [39]
   56,    # [40]

   70,    # [41]
   60,    # [42]
   48,    # [43]
   52,    # [44]
   56,    # [45]
   72,    # [46]
   70,    # [47]
   68,    # [48]
   78,    # [49]
   86,    # [50]

   82,    # [51]
   93,    # [52]
   98,    # [53]
   94,    # [54]
   110,   # [55]
   87,    # [56]
   95,    # [57]
   79,    # [58]
   88,    # [59]
   80,    # [60]

   69,    # [61]
   76,    # [62]
   91,    # [63]
   89,    # [64]
   93,    # [65]
   112,   # [66]
   108,   # [67]
   140,   # [68]
   129,   # [69]
   157,   # [70]

   138,   # [71]
   147,   # [72]
   129,   # [73]
   111,   # [74]
   101,   # [75]
   105,   # [76]
   98,    # [77]
   117,   # [78]
   106,   # [79]
   114,   # [80]

   131,   # [81]
   124,   # [82]
   132,   # [83]
   118,   # [84]
   128,   # [85]
   133,   # [86]
   128,   # [87]
   140,   # [88]
   129,   # [89]
   126,   # [90]

   140,   # [91]
   147,   # [92]
   168,   # [93]
   163,   # [94]
   174,   # [95]
   164,   # [96]
   170,   # [97]
   152,   # [98]
   150,   # [99]
   144,   # [100]
);

# Lidka cells after 100 ticks
my @lidka100 = (
   [ -29, 2 ],    # 1-10
   [ -28, 1 ],
   [ -28, 2 ],
   [ -28, 3 ],
   [ -27, 0 ],
   [ -27, 4 ],
   [ -26, 0 ],
   [ -26, 1 ],
   [ -26, 4 ],
   [ -26, 5 ],

   [ -25, 3 ],    # 10-20
   [ -25, 4 ],
   [ -24, 2 ],
   [ -24, 3 ],
   [ -23, 1 ],
   [ -23, 2 ],
   [ -22, 0 ],
   [ -21, 0 ],
   [ -21, 1 ],
   [ -17, -2 ],

   [ -17, -1 ],    # 20-30
   [ -16, -2 ],
   [ -16, -1 ],
   [ -12, 8 ],
   [ -11, 9 ],
   [ -11, 10 ],
   [ -11, 16 ],
   [ -11, 17 ],
   [ -10, 8 ],
   [ -10, 9 ],

   [ -10, 10 ],    # 30-40
   [ -10, 15 ],
   [ -10, 16 ],
   [ -10, 18 ],
   [ -10, 19 ],
   [ -9, 18 ],
   [ -9, 19 ],
   [ -8, 13 ],
   [ -8, 18 ],
   [ -8, 19 ],

   [ -7, 12 ],    # 40-50
   [ -6, -19 ],
   [ -6, -18 ],
   [ -6, 16 ],
   [ -6, 17 ],
   [ -5, -19 ],
   [ -5, -17 ],
   [ -5, 12 ],
   [ -5, 16 ],
   [ -5, 17 ],

   [ -4, -19 ],    # 50-60
   [ -4, 13 ],
   [ -4, 14 ],
   [ -4, 15 ],
   [ -4, 16 ],
   [ -4, 18 ],
   [ -3, 14 ],
   [ -3, 15 ],
   [ -3, 18 ],
   [ -2, 17 ],

   [ -2, 18 ],    # 60-70
   [ -1, -13 ],
   [ -1, 6 ],
   [ 0, -14 ],
   [ 0, -13 ],
   [ 0, -12 ],
   [ 0, 5 ],
   [ 0, 6 ],
   [ 0, 7 ],
   [ 1, -15 ],

   [ 1, -14 ],    # 70-80
   [ 1, -12 ],
   [ 1, 6 ],
   [ 1, 9 ],
   [ 1, 10 ],
   [ 2, -22 ],
   [ 2, -15 ],
   [ 2, 3 ],
   [ 2, 6 ],
   [ 3, -22 ],

   [ 3, -16 ],    # 80-90
   [ 3, 2 ],
   [ 3, 3 ],
   [ 3, 4 ],
   [ 3, 5 ],
   [ 3, 9 ],
   [ 4, -22 ],
   [ 4, -17 ],
   [ 4, -16 ],
   [ 4, -15 ],

   [ 4, -12 ],    # 90-100
   [ 4, -11 ],
   [ 4, 3 ],
   [ 4, 8 ],
   [ 5, -16 ],
   [ 5, -14 ],
   [ 5, -12 ],
   [ 5, 7 ],
   [ 6, -16 ],
   [ 6, -15 ],

   [ 6, -14 ],    # 100-110
   [ 6, 4 ],
   [ 6, 6 ],
   [ 7, -15 ],
   [ 7, -14 ],
   [ 7, 4 ],
   [ 7, 12 ],
   [ 7, 13 ],
   [ 8, -14 ],
   [ 8, 12 ],

   [ 8, 13 ],    # 110-120
   [ 9, -14 ],
   [ 9, -13 ],
   [ 9, -4 ],
   [ 9, -2 ],
   [ 9, 10 ],
   [ 9, 11 ],
   [ 10, -5 ],
   [ 10, -1 ],
   [ 10, 10 ],

   [ 10, 11 ],    # 120-130
   [ 11, -1 ],
   [ 11, 0 ],
   [ 12, -7 ],
   [ 12, -1 ],
   [ 12, 0 ],
   [ 13, -8 ],
   [ 13, -7 ],
   [ 13, -3 ],
   [ 13, -2 ],

   [ 13, -1 ],    # 130-140
   [ 14, -8 ],
   [ 14, -3 ],
   [ 14, -2 ],
   [ 14, 1 ],
   [ 15, -7 ],
   [ 15, -6 ],
   [ 15, -5 ],
   [ 15, -1 ],
   [ 15, 0 ],

   [ 15, 1 ],    # 140-144
   [ 16, -7 ],
   [ 16, -6 ],
   [ 16, -5 ],
);
my @slidka100 = sort { $a->[0] <=> $b->[0] || $a->[1] <=> $b->[1] } @l
+idka100;

# Lidka starting pattern
my @lidka0 = (
   [ -3, -7 ],
   [ -4, -6 ],
   [ -2, -6 ],
   [ -3, -5 ],
   [  4,  3 ],
   [  2,  4 ],
   [  4,  4 ],
   [  1,  5 ],
   [  2,  5 ],
   [  4,  5 ],
   [  0,  7 ],
   [  1,  7 ],
   [  2,  7 ],
);
my @slidka0 = sort { $a->[0] <=> $b->[0] || $a->[1] <=> $b->[1] } @lid
+ka0;

# Initial cell array
my $org = Organism->new();
$org->insert_cells(@lidka0);
{
   my @cells = $org->get_live_cells();
   is_deeply( \@cells, \@slidka0, "lidka initial cell array" );
   test_one( $org, "lidka 0", $lidka_ticks[0] );
}

# Test first 100 ticks
for my $i ( 1 .. $nticks ) {
   $org->tick();
   test_one( $org, "lidka $i", $lidka_ticks[$i] );
}

# Final cell array
{
   my @cells = $org->get_live_cells();
   cmp_ok( scalar(@cells), '==', $lidka_ticks[100], "lidka final array
+ count" );
   is_deeply( \@cells, \@slidka100, "lidka final cell array" );
}

        [download]

[reply]
[d/l]
[select]

In Section Meditations