my $c = cos($rotateby);
my $transform = Math::GSL::Matrix->new(2,2);
$transform->set_row(0, [ $c, -$s ]);
- $transform->set_row(1, [ -$s, $c ]);
+ $transform->set_row(1, [ $s, $c ]);
print STDERR "rotate was=",prvec($dir)," theta=$theta",
" rotateby=$rotateby s=$s c=$c\n";
my @lims;
my $v = $topend ? -$GSL_DBL_MAX : $GSL_DBL_MAX;
$lims[$topend] = Math::GSL::Vector->new([$v,$v]);
}
- # xxx rotation is not happening
for_each_pos(sub {
my ($pr, $rr, $why) = @_;
my $y = Math::GSL::Vector->new(2);
$$pr->raw,
0, $y->raw)
and confess;
- gsl_blas_dcopy($$pr->raw, $y->raw)
+ #print STDERR "gsl_blas_dgemv ",prvec($$pr)," => ",prvec($y),"\n";
+ gsl_blas_dcopy($y->raw, $$pr->raw)
and confess;
foreach my $topend (qw(0 1)) {
foreach my $xy (qw(0 1)) {