Skip to content

Commit 9904810

Browse files
committed
non-Euclid use $range2
1 parent 4136610 commit 9904810

1 file changed

Lines changed: 8 additions & 12 deletions

File tree

lib/PDL/Graphics/TriD/Graph.pm

Lines changed: 8 additions & 12 deletions
Original file line numberDiff line numberDiff line change
@@ -297,17 +297,15 @@ sub transform {
297297
my ($this,$point,$data,$inds) = @_;
298298
barf "no \$inds given" if !defined $inds;
299299
barf "Wrong number of arguments to transform $this\n" if @$inds != 3;
300-
my ($longrange, $latrange) = $this->{Scale}->slice('0:1')->t->diff2->dog;
300+
my $range2 = $this->{Scale}->t->diff2->t->slice('0:1');
301301
my $pressure_max = $this->{Scale}->slice('2,1');
302302
$data = $data->dice_axis(0, $inds);
303-
my $data01_ctr = $data->slice("0:1")-$this->{Center};
303+
my $data01_ctr = ($data->slice("0:1")-$this->{Center}) / $range2;
304304
$point->slice("(0)") +=
305-
0.5+$data01_ctr->slice("(0)") /
306-
$longrange
305+
0.5+$data01_ctr->slice("(0)")
307306
*cos($data->slice("(1)")*DEG2RAD);
308307
$point->slice("(1)") +=
309-
0.5+$data01_ctr->slice("(1)") /
310-
$latrange;
308+
0.5+$data01_ctr->slice("(1)");
311309
$point->slice("(2)") .=
312310
log($data->slice("(2)")/1012.5)/log($pressure_max/1012.5);
313311
$point;
@@ -388,16 +386,14 @@ sub transform {
388386
my $i = 0;
389387
barf "Wrong number of arguments to transform $this\n" if @$inds != 3;
390388
$data = $data->dice_axis(0, $inds);
391-
my ($longrange, $latrange) = $this->{Scale}->slice('0:1')->t->diff2->dog;
389+
my $range2 = $this->{Scale}->t->diff2->t->slice('0:1');
392390
my $pressure_max = $this->{Scale}->slice('2,1');
393-
my $data01_ctr = $data->slice("0:1")-$this->{Center};
391+
my $data01_ctr = ($data->slice("0:1")-$this->{Center}) / $range2;
394392
$point->slice("(0)") +=
395-
0.5+$data01_ctr->slice("(0)") /
396-
$longrange
393+
0.5+$data01_ctr->slice("(0)")
397394
*cos($data->slice("(1)")*DEG2RAD);
398395
$point->slice("(1)") +=
399-
0.5+$data01_ctr->slice("(1)") /
400-
$latrange
396+
0.5+$data01_ctr->slice("(1)")
401397
*cos($data->slice("(1)")*DEG2RAD);
402398
# Longitude transformation
403399
# $point->slice("(0)") =

0 commit comments

Comments
 (0)