#! /usr/bin/env perl
#
# Demo x25 for the PLplot PDL binding
#
# Copyright (C) 2006 Andrew Ross
#
# This file is part of PLplot.
#
# PLplot is free software; you can redistribute it and/or modify
# it under the terms of the GNU General Library Public License as published
# by the Free Software Foundation; either version 2 of the License, or
# (at your option) any later version.
#
# PLplot is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
# GNU Library General Public License for more details.
#
# You should have received a copy of the GNU Library General Public License
# along with PLplot; if not, write to the Free Software
# Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA

# SYNC: x25c.c 1.3

use PDL;
use PDL::Graphics::PLplot;
use Math::Trig qw [pi];

use Getopt::Long;

my ($i, $j);

# Parse and process command line arguments 

plParseOpts (\@ARGV, PL_PARSE_SKIP | PL_PARSE_NOPROGRAM);


# Initialize plplot

plssub(3,3);
plinit();

my $xextreme = pdl( [-120, -120, -120, -80, -220, -20, -20, -80, 20], [120, 120, 120, 80, -120, 20, 20, 80, 120]);
my $yextreme = pdl( [-120, 20, -20, -20, -120, -120, -20, -80, -120], [120, 120, 120, 120, 120, 120, 20, 80, 120]);
my ($x0, $y0);

for ( $j = 0; $j < 4; $j ++ ) {
    if ( $j == 0 ) {
# Polygon 1: a diamond 
	$x0 = pdl( 0, -100, 0, 100);
	$y0 = pdl( -100, 0, 100, 0);
    }
    if ( $j == 1 ) {
# Polygon 1: a diamond - reverse direction
	$x0 = pdl( 100, 0, -100, 0);
	$y0 = pdl( 0, 100, 0, -100);
    }
    if ( $j == 2 ) {
# Polygon 2: a square with punctures
	$x0 = pdl( -100, -100, 80, -100, -100, -80, 0, 80, 100, 100);
	$y0 = pdl( -100, -80, 0, 80, 100, 100, 80, 100, 100, -100);
    }
    if ( $j == 3 ) {
# Polygon 2: a square with punctures - reversed direction 
	$x0 = pdl( 100, 100, 80, 0, -80, -100, -100, 80, -100, -100);
	$y0 = pdl( -100, 100, 100, 80, 100, 100, 80, 0, -80, -100);
    }

    for ( $i = 0; $i < 9 ; $i ++ ) {
	pladv(0);
	plvsta();
	plwind($xextreme->at($i,0), $xextreme->at($i,1), $yextreme->at($i,0), $yextreme->at($i,1));

	plcol0(2) ;
	plbox(1.0, 0, 10.0, 0, "bc", "bcnv", );
	plcol0(1) ;
	plpsty(0) ;
	plfill($x0,$y0) ;
	plcol0(2) ;
	pllsty(1) ;
	plline($x0,$y0) ;
    }
}

# Don't forget to call plend() to finish off!

plend();