t/Matrix.t

   1 package Math::GSL::Matrix::Test;
   2 use base q{Test::Class};
   3 use Test::More tests => 209;
   4 use Math::GSL           qw/:all/;
   5 use Math::GSL::Test     qw/:all/;
   6 use Math::GSL::Matrix   qw/:all/;
   7 use Math::GSL::Vector   qw/:all/;
   8 use Math::GSL::Complex  qw/:all/;
   9 use Math::GSL::Errno    qw/:all/;
  10 use Data::Dumper;
  11 use strict;
  12
  13 BEGIN{ gsl_set_error_handler_off(); }
  14
  15 sub make_fixture : Test(setup) {
  16     my $self = shift;
  17     $self->{matrix} = gsl_matrix_alloc(5,5);
  18     $self->{obj}    = Math::GSL::Matrix->new(5,5);
  19     gsl_matrix_set_zero($self->{matrix});
  20 }
  21
  22 sub teardown : Test(teardown) {
  23     unlink 'matrix' if -f 'matrix';
  24 }
  25
  26 sub GSL_MATRIX_ALLOC : Tests {
  27     my $matrix = gsl_matrix_alloc(5,5);
  28     isa_ok($matrix, 'Math::GSL::Matrix');
  29 }
  30
  31 sub GSL_MATRIX_SET : Tests {
  32     my $self = shift;
  33     map { gsl_matrix_set($self->{matrix}, $_,$_, $_ ** 2) } (0..4);
  34
  35     my @got = map { gsl_matrix_get($self->{matrix}, $_, $_) } (0..4);
  36     ok_similar( [ @got ], [ map { $_**2 } (0..4) ] );
  37 }
  38
  39 sub GSL_MATRIX_CALLOC : Tests {
  40    my $matrix = gsl_matrix_calloc(5,5);
  41    isa_ok($matrix, 'Math::GSL::Matrix');
  42
  43    my @got = map { gsl_matrix_get($matrix, $_, $_) } (0..4);
  44    ok_similar( [ @got ], [ (0) x 5 ], 'gsl_matrix_calloc' );
  45 }
  46
  47 sub GSL_MATRIX_FREE : Tests {
  48    my $matrix = gsl_matrix_calloc(5,5);
  49    isa_ok($matrix, 'Math::GSL::Matrix');
  50
  51    is(gsl_matrix_get($matrix, 0, 0), 0);
  52    gsl_matrix_free($matrix);
  53 }
  54
  55 sub GSL_MATRIX_SUBMATRIX : Tests {
  56    my $matrix = gsl_matrix_alloc(5,5);
  57    map { gsl_matrix_set($matrix, $_,$_, $_) } (0..4);
  58    my $subMatrix = gsl_matrix_submatrix($matrix, 0, 0, 2, 2);
  59    my @got = map { gsl_matrix_get($matrix, $_, $_) } (0..2);
  60    ok_similar( [ @got ], [ 0..2 ] );
  61 }
  62
  63 sub GSL_MATRIX_ROW : Tests {
  64    my $self = shift;
  65    for my $line (0..4) {
  66         map { gsl_matrix_set($self->{matrix}, $_,$line, $_) } (0..4);
  67    }
  68
  69    my $vector_view = gsl_matrix_row($self->{matrix}, 2);
  70    my @got = map { gsl_vector_get($vector_view->{vector}, $_) } (0..4);
  71    ok_similar( [ @got ], [ (2)x 5], 'gsl_matrix_row' );
  72 }
  73
  74 sub GSL_MATRIX_COLUMN : Tests {
  75     my $self = shift;
  76     my $view = gsl_vector_alloc(5);
  77
  78     for my $line (0..4) {
  79         map { gsl_matrix_set($self->{matrix}, $line,$_, $line*$_) } (0..4);
  80     }
  81     $view = gsl_matrix_column($self->{matrix}, 2);
  82     my $vec = $view->swig_vector_get();
  83
  84     my @got = map { gsl_vector_get($vec, $_) } (0..4);
  85     ok_similar( [ @got ], [0,2,4,6,8 ], 'gsl_matrix_column' );
  86 }
  87
  88 sub GSL_MATRIX_DIAGONAL : Tests {
  89    my $matrix = gsl_matrix_alloc(4,4);
  90    map { gsl_matrix_set($matrix, $_,$_, $_) } (0..3);
  91    my $view = gsl_matrix_diagonal($matrix);
  92    # better interface is needed
  93    my $vec = $view->swig_vector_get();
  94
  95    my @got = map { gsl_vector_get($vec, $_) } (0..3);
  96    ok_similar( [ @got ], [ 0 .. 3 ], 'gsl_matrix_diagonal');
  97 }
  98
  99 sub GSL_MATRIX_SUBDIAGONAL : Tests {
 100    my $matrix = gsl_matrix_alloc(4,4);
 101
 102    map { gsl_matrix_set($matrix, $_,$_, $_)     } (0..3);
 103
 104    my $view = gsl_matrix_subdiagonal($matrix, 0);
 105    my $vec  = $view->swig_vector_get();
 106    my @got  = map { gsl_vector_get($vec, $_) } (0..3);
 107    ok_similar( [ @got ], [ 0 .. 3 ], 'gsl_matrix_subdiagonal');
 108 }
 109
 110 sub GSL_MATRIX_SWAP : Tests {
 111    my $self=shift;
 112    map { gsl_matrix_set($self->{matrix}, $_,$_, $_ ** 2) } (0..4);
 113    my $matrix = gsl_matrix_alloc(5,5);
 114    map { gsl_matrix_set($matrix, $_,$_, $_) } (0..4);
 115    ok_status(gsl_matrix_swap($self->{matrix}, $matrix));
 116
 117    my @got = map { gsl_matrix_get($self->{matrix}, $_, $_) } (0..4);
 118    map { is($got[$_], $_) } (0..4);
 119    @got = map { gsl_matrix_get($matrix, $_, $_) } (0..4);
 120    map { is($got[$_], $_** 2) } (0..4);
 121 }
 122
 123 sub GSL_MATRIX_MEMCPY : Tests {
 124    my $self = shift;
 125    my $matrix = gsl_matrix_alloc(5,5);
 126    map { gsl_matrix_set($self->{matrix}, $_,$_, $_ ** 2) } (0..4);
 127    ok_status(gsl_matrix_memcpy($matrix, $self->{matrix}));
 128    ok_similar( [ map { gsl_matrix_get($matrix, $_, $_) } (0..4) ],
 129                [ map {  $_** 2 } (0..4)                        ]
 130    );
 131 }
 132
 133 sub GSL_MATRIX_SWAP_ROWS : Tests {
 134    my $self = shift;
 135    map { gsl_matrix_set($self->{matrix}, 0,$_, $_) } (0..4);
 136    map { gsl_matrix_set($self->{matrix}, 1,$_, 3) } (0..4);
 137    ok_status(gsl_matrix_swap_rows($self->{matrix}, 0, 1));
 138    my @got = map { gsl_matrix_get($self->{matrix}, 1, $_) } (0..4);
 139    map { is($got[$_], $_) } (0..4);
 140    @got = map { gsl_matrix_get($self->{matrix}, 0, $_) } (0..4);
 141    map { is($got[$_], 3) } (0..4);
 142 }
 143
 144 sub GSL_MATRIX_SWAP_COLUMNS : Tests {
 145    my $self = shift;
 146    map { gsl_matrix_set($self->{matrix}, $_,0, $_) } (0..4);
 147    map { gsl_matrix_set($self->{matrix}, $_,1, 3) } (0..4);
 148    ok_status(gsl_matrix_swap_columns($self->{matrix}, 0, 1));
 149    my @got = map { gsl_matrix_get($self->{matrix}, $_, 1) } (0..4);
 150    map { is($got[$_], $_) } (0..4);
 151    @got = map { gsl_matrix_get($self->{matrix}, $_, 0) } (0..4);
 152    map { is($got[$_], 3) } (0..4);
 153 }
 154
 155 sub GSL_MATRIX_SWAP_ROWCOL : Tests {
 156    my $self = shift;
 157    map { gsl_matrix_set($self->{matrix}, 0,$_, $_) } (0..4);
 158    map { gsl_matrix_set($self->{matrix}, $_,2, 2) } (0..4);
 159    ok_status(gsl_matrix_swap_rowcol($self->{matrix}, 0, 2));
 160
 161    my @got = map { gsl_matrix_get($self->{matrix}, $_, 2) } (0..4);
 162    is_deeply( [ @got ], [ qw/2 1 0 3 4/  ] );
 163
 164    @got = map { gsl_matrix_get($self->{matrix}, 0, $_) } (0..4);
 165    is_deeply( [ @got ], [ qw/2 2 2 2 2/ ] );
 166 }
 167
 168 sub GSL_MATRIX_TRANSPOSE_MEMCPY : Tests {
 169    my $self = shift;
 170    map { gsl_matrix_set($self->{matrix}, 0,$_, $_) } (0..4);
 171    my $matrix = gsl_matrix_alloc(5,5);
 172    ok_status(gsl_matrix_transpose_memcpy($matrix, $self->{matrix}));
 173    my @got = map { gsl_matrix_get($matrix, $_, 0) } (0..4);
 174    map { is($got[$_], $_) } (0..4);
 175 }
 176
 177 sub GSL_MATRIX_TRANSPOSE : Tests {
 178    my $self = shift;
 179    map { gsl_matrix_set($self->{matrix}, 0,$_, $_) } (0..4);
 180    is(gsl_matrix_transpose($self->{matrix}), 0);
 181    my @got = map { gsl_matrix_get($self->{matrix}, $_, 0) } (0..4);
 182    map { is($got[$_], $_) } (0..4);
 183 }
 184
 185 sub GSL_MATRIX_ADD : Tests {
 186    my $self = shift;
 187    my $matrix = gsl_matrix_alloc(5, 5);
 188    map { gsl_matrix_set($self->{matrix}, $_, $_, $_) } (0..4);
 189    map { gsl_matrix_set($matrix, $_, $_, $_) } (0..4);
 190    is(gsl_matrix_add($self->{matrix}, $matrix), 0);
 191    my @got = map { gsl_matrix_get($self->{matrix}, $_, $_) } (0..4);
 192    map { is($got[$_], $_*2) } (0..4);
 193 }
 194
 195 sub GSL_MATRIX_SUB : Tests {
 196    my $self = shift;
 197    my $matrix = gsl_matrix_alloc(5, 5);
 198    map { gsl_matrix_set($self->{matrix}, $_, $_, $_) } (0..4);
 199    map { gsl_matrix_set($matrix, $_, $_, $_) } (0..4);
 200    is(gsl_matrix_sub($self->{matrix}, $matrix), 0);
 201    my @got = map { gsl_matrix_get($self->{matrix}, $_, $_) } (0..4);
 202    map { is($got[$_], 0) } (0..4);
 203 }
 204
 205 sub GSL_MATRIX_MUL_ELEMENTS : Tests {
 206    my $self = shift;
 207    my $matrix = gsl_matrix_alloc(5, 5);
 208    map { gsl_matrix_set($self->{matrix}, $_, $_, $_) } (0..4);
 209    map { gsl_matrix_set($matrix, $_, $_, $_) } (0..4);
 210    is(gsl_matrix_mul_elements($self->{matrix}, $matrix), 0);
 211    my @got = map { gsl_matrix_get($self->{matrix}, $_, $_) } (0..4);
 212    map { is($got[$_], $_**2) } (0..4);
 213 }
 214
 215 sub GSL_MATRIX_DIV_ELEMENTS : Tests {
 216    my $self = shift;
 217    my $matrix = gsl_matrix_alloc(5, 5);
 218    map { gsl_matrix_set($self->{matrix}, $_, $_, $_+1) } (0..4);
 219    map { gsl_matrix_set($matrix, $_, $_, $_+1) } (0..4);
 220    is(gsl_matrix_div_elements($self->{matrix}, $matrix), 0);
 221    my @got = map { gsl_matrix_get($self->{matrix}, $_, $_) } (0..4);
 222    map { is($got[$_], 1) } (0..4);
 223 }
 224
 225 sub GSL_MATRIX_SCALE : Tests {
 226    my $self = shift;
 227    map { gsl_matrix_set($self->{matrix}, $_, $_, $_) } (0..4);
 228    is(gsl_matrix_scale($self->{matrix}, 4), 0);
 229    my @got = map { gsl_matrix_get($self->{matrix}, $_, $_) } (0..4);
 230    map { is($got[$_], $_*4) } (0..4);
 231 }
 232
 233 sub GSL_MATRIX_ADD_CONSTANT : Tests {
 234    my $self = shift;
 235    map { gsl_matrix_set($self->{matrix}, $_, $_, $_) } (0..4);
 236    is(gsl_matrix_add_constant($self->{matrix}, 8), 0);
 237    my @got = map { gsl_matrix_get($self->{matrix}, $_, $_) } (0..4);
 238    map { is($got[$_], $_+8) } (0..4);
 239 }
 240
 241 sub GSL_MATRIX_MAX : Tests {
 242    my $self = shift;
 243    for my $row (0..4) {
 244         map { gsl_matrix_set($self->{matrix}, $row, $_, $_**2 ) } (0..4);
 245    }
 246    is(gsl_matrix_max($self->{matrix}), 16);
 247 }
 248
 249 sub GSL_MATRIX_MIN : Tests {
 250    my $self = shift;
 251    for my $row (0..4) {
 252         map { gsl_matrix_set($self->{matrix}, $row, $_, $_**2 ) } (0..4);
 253    }
 254    is(gsl_matrix_min($self->{matrix}), 0);
 255 }
 256
 257 sub GSL_MATRIX_MINMAX : Test {
 258    my $self = shift;
 259    map { gsl_matrix_set($self->{matrix}, $_, $_, $_**2) } (0..4);
 260    my ($min, $max) = gsl_matrix_minmax($self->{matrix});
 261    ok_similar( [ $min, $max ], [ 0, 16], 'gsl_matrix_minmax' );
 262 }
 263
 264 sub GSL_MATRIX_MAX_INDEX : Tests {
 265    my $self = shift;
 266    for my $row (0..3) {
 267         map { gsl_matrix_set($self->{matrix}, $row, $_, $_) } (0..4);
 268    }
 269    map { gsl_matrix_set($self->{matrix}, $_, $_, $_**2) } (0..4);
 270    my ($imax, $jmax) = gsl_matrix_max_index($self->{matrix});
 271    ok_similar( [ $imax, $jmax ], [ 4, 4 ], 'gsl_matrix_max_index' );
 272 }
 273
 274 sub GSL_MATRIX_MIN_INDEX : Tests {
 275    my $self = shift;
 276    map { gsl_matrix_set($self->{matrix}, $_, $_, $_**2) } (0..4);
 277    my ($imin, $jmin) = gsl_matrix_min_index($self->{matrix});
 278    ok_similar( [ $imin, $jmin ], [ 0, 0 ], 'gsl_matrix_min_index' );
 279
 280 }
 281
 282 sub GSL_MATRIX_ISNULL : Tests {
 283    my $self = shift;
 284    for my $line (0..4) {
 285         map { gsl_matrix_set($self->{matrix}, $line, $_, 0) } (0..4);
 286    }
 287    is(gsl_matrix_isnull($self->{matrix}), 1);
 288
 289    for my $line (0..4) {
 290         map { gsl_matrix_set($self->{matrix}, $line, $_, $_) } (0..4);
 291    }
 292
 293    is(gsl_matrix_isnull($self->{matrix}), 0);
 294 }
 295
 296 sub GSL_MATRIX_ISPOS : Tests {
 297    my $self = shift;
 298    my $line;
 299
 300    for($line=0; $line<5; $line++) {
 301    map { gsl_matrix_set($self->{matrix}, $line, $_, -1) } (0..4); }
 302    is(gsl_matrix_ispos($self->{matrix}), 0);
 303
 304    for($line=0; $line<5; $line++) {
 305    map { gsl_matrix_set($self->{matrix}, $line, $_, 1) } (0..4); }
 306    is(gsl_matrix_ispos($self->{matrix}), 1);
 307
 308    for($line=0; $line<5; $line++) {
 309    map { gsl_matrix_set($self->{matrix}, $line, $_, 0) } (0..4); }
 310    is(gsl_matrix_ispos($self->{matrix}), 0);
 311 }
 312
 313 sub GSL_MATRIX_ISNEG : Tests {
 314    my $self = shift;
 315    my $line;
 316
 317    for($line=0; $line<5; $line++) {
 318    map { gsl_matrix_set($self->{matrix}, $line, $_, -1) } (0..4); }
 319    is(gsl_matrix_isneg($self->{matrix}), 1);
 320
 321    for($line=0; $line<5; $line++) {
 322    map { gsl_matrix_set($self->{matrix}, $line, $_, 1) } (0..4); }
 323    is(gsl_matrix_isneg($self->{matrix}), 0);
 324
 325    for($line=0; $line<5; $line++) {
 326    map { gsl_matrix_set($self->{matrix}, $line, $_, 0) } (0..4); }
 327    is(gsl_matrix_isneg($self->{matrix}), 0);
 328 }
 329
 330 sub GSL_MATRIX_ISNONNEG : Tests {
 331     my $self = shift;
 332     for my $row (0..4) {
 333         map { gsl_matrix_set($self->{matrix}, $row, $_, -1) } (0..4);
 334     }
 335     is(gsl_matrix_isnonneg($self->{matrix}), 0);
 336
 337     for my $row (0..4) {
 338         map { gsl_matrix_set($self->{matrix}, $row, $_, 1) } (0..4);
 339     }
 340     is(gsl_matrix_isnonneg($self->{matrix}), 1);
 341
 342     for my $row (0..4) {
 343         map { gsl_matrix_set($self->{matrix}, $row, $_, 0) } (0..4);
 344     }
 345     is(gsl_matrix_isnonneg($self->{matrix}), 1);
 346
 347 }
 348
 349 sub GSL_MATRIX_GET_ROW : Tests {
 350    my $self = shift;
 351    my $vector = gsl_vector_alloc(5);
 352    map { gsl_matrix_set($self->{matrix}, 0, $_, $_) } (0..4);
 353    is(gsl_matrix_get_row($vector, $self->{matrix}, 0), 0);
 354    map { is(gsl_vector_get($vector, $_), $_) } (0..4);
 355 }
 356
 357 sub GSL_MATRIX_GET_COL : Tests {
 358    my $self = shift;
 359    my $vector = gsl_vector_alloc(5);
 360    map { gsl_matrix_set($self->{matrix}, $_, 0, $_) } (0..4);
 361    is(gsl_matrix_get_col($vector, $self->{matrix}, 0), 0);
 362    map { is(gsl_vector_get($vector, $_), $_) } (0..4);
 363 }
 364
 365 sub GSL_MATRIX_SET_ROW : Tests {
 366    my $self = shift;
 367    my $vector = gsl_vector_alloc(5);
 368    map { gsl_vector_set($vector, $_, $_**2) } (0..4);
 369    is(gsl_matrix_set_row($self->{matrix}, 1, $vector), 0);
 370
 371    ok_similar( [ map { gsl_matrix_get($self->{matrix}, 1, $_) } (0..4) ],
 372                [ map { $_ ** 2 } (0..4) ],
 373              );
 374 }
 375
 376 sub GSL_MATRIX_SET_COL : Tests {
 377    my $self = shift;
 378    my $vector = gsl_vector_alloc(5);
 379    map { gsl_vector_set($vector, $_, $_**2) } (0..4);
 380    is(gsl_matrix_set_col($self->{matrix}, 1, $vector), 0);
 381
 382    ok_similar( [ map { gsl_matrix_get($self->{matrix}, $_, 1) } (0..4) ],
 383                [ map { $_ ** 2 } (0..4) ],
 384              );
 385 }
 386
 387 sub GSL_MATRIX_FREAD_FWRITE : Tests {
 388    my $self = shift;
 389    for my $line (0..4) {
 390         map { gsl_matrix_set($self->{matrix}, $line, $_, $_**2) } (0..4);
 391    }
 392
 393    my $fh = gsl_fopen("matrix", 'w');
 394    is( gsl_matrix_fwrite($fh, $self->{matrix}), 0);
 395    gsl_fclose($fh);
 396
 397    for my $line (0..4) {
 398         map { gsl_matrix_set($self->{matrix}, $line, $_, $_**3) } (0..4);
 399    }
 400
 401    $fh = gsl_fopen("matrix", 'r');
 402
 403    is(gsl_matrix_fread($fh, $self->{matrix}), 0);
 404    for my $line (0..4) {
 405         ok_similar(
 406                     [ map { gsl_matrix_get($self->{matrix}, $line, $_) } (0..4) ],
 407                     [ map { $_**2 } (0..4) ],
 408                   );
 409    }
 410    gsl_fclose($fh);
 411 }
 412
 413 sub GSL_MATRIX_FPRINTF_FSCANF : Tests {
 414    my $self = shift;
 415
 416    for my $line (0..4) {
 417         map { gsl_matrix_set($self->{matrix}, $line, $_, $_**2) } (0..4);
 418    }
 419
 420    my $fh = gsl_fopen("matrix", 'w');
 421    is( gsl_matrix_fprintf($fh, $self->{matrix}, "%f"), 0);
 422    ok_status(gsl_fclose($fh));
 423
 424    for my $line (0..4) {
 425         map { gsl_matrix_set($self->{matrix}, $line, $_, $_**3) } (0..4);
 426    }
 427
 428    $fh = gsl_fopen("matrix", 'r');
 429    is(gsl_matrix_fscanf($fh, $self->{matrix}), 0);
 430    for my $line (0..4) {
 431         ok_similar(
 432                     [ map { gsl_matrix_get($self->{matrix}, $line, $_) } (0..4) ],
 433                     [ map { $_**2 } (0..4) ],
 434                   );
 435    }
 436    ok_status(gsl_fclose($fh));
 437 }
 438
 439 sub GSL_MATRIX_MINMAX_INDEX : Tests {
 440    my $self = shift;
 441    my $line;
 442    for ($line = 0; $line<4; $line ++) {
 443    map { gsl_matrix_set($self->{matrix}, $line, $_, $_) } (0..4); }
 444    map { gsl_matrix_set($self->{matrix}, 4, $_, $_**2) } (0..4);
 445    my ($imin, $jmin, $imax, $jmax) = gsl_matrix_minmax_index($self->{matrix});
 446    ok_similar( [ $imin, $jmin, $imax, $jmax ], [ 0, 0, 4, 4], 'gsl_matrix_minmax_index' );
 447 }
 448
 449 sub GSL_MATRIX_ADD_DIAGONAL : Tests {
 450    my $self = shift;
 451    map { gsl_matrix_set($self->{matrix}, $_, $_, $_) } (0..4);
 452    gsl_matrix_add_diagonal($self->{matrix}, 4);
 453    ok_similar( [ map { gsl_matrix_get($self->{matrix}, $_, $_)} (0..4) ],
 454                [ 4 .. 8 ],
 455              );
 456 }
 457
 458 sub GSL_MATRIX_NEW : Tests {
 459    my $self = shift;
 460    isa_ok( $self->{obj}, 'Math::GSL::Matrix' );
 461    isa_ok( $self->{obj}->raw, 'Math::GSL::Matrix::gsl_matrix' );
 462    ok( $self->{obj}->rows == 5, '->rows' );
 463    ok( $self->{obj}->cols == 5, '->cols' );
 464 }
 465
 466 sub AS_LIST_SQUARE : Tests {
 467     my $matrix = Math::GSL::Matrix->new(5,5);
 468     map { gsl_matrix_set($matrix->raw, $_, $_, 5 + $_**2) } (0..4);
 469     is_deeply( [
 470                 5, 0, 0, 0, 0,
 471                 0, 6, 0, 0, 0,
 472                 0, 0, 9, 0, 0,
 473                 0, 0, 0,14, 0,
 474                 0, 0, 0, 0, 21
 475                 ],
 476                [ $matrix->as_list],
 477                '$matrix->as_list',
 478     );
 479 }
 480
 481 sub AS_LIST_ROW : Tests {
 482     my $matrix = Math::GSL::Matrix->new(1,5);
 483     map { gsl_matrix_set($matrix->raw, 0 , $_, 5 + $_**2) } (0..4);
 484     is_deeply( [ 5, 6, 9, 14, 21, ],
 485                [ $matrix->as_list],
 486                '$matrix->as_list',
 487     );
 488 }
 489
 490 sub ROW : Tests {
 491     my $matrix = Math::GSL::Matrix->new(5,5);
 492     map { gsl_matrix_set($matrix->raw, $_, $_, 5 + $_**2) } (0..4);
 493     ok_similar( [qw/0 0 9 0 0/] , [$matrix->row(2)->as_list ] );
 494     map { gsl_matrix_set($matrix->raw, 0, $_, $_) } (0..4);
 495     ok_similar( [qw/0 1 2 3 4/] , [$matrix->row(0)->as_list ] );
 496 }
 497
 498 sub COL : Tests {
 499     my $matrix = Math::GSL::Matrix->new(3,3);
 500     map { gsl_matrix_set($matrix->raw, $_, $_, 7 + $_**2) } (0..2);
 501     ok_similar([7,0,0], [$matrix->col(0)->as_list]);
 502 }
 503
 504 sub NEW_SETS_VALUES_TO_ZERO : Tests {
 505     my $matrix = Math::GSL::Matrix->new(5,5);
 506     my $sum;
 507
 508     map { $sum += $_ } $matrix->as_list;
 509     ok( $sum == 0, 'new sets values to zero');
 510 }
 511
 512 sub HERMITIAN : Tests {
 513     my $matrix    = gsl_matrix_complex_alloc(2,2);
 514     my $transpose = gsl_matrix_complex_alloc(2,2);
 515     gsl_matrix_complex_set($matrix, 0, 0, gsl_complex_rect(3,0));
 516     gsl_matrix_complex_set($matrix, 0, 1, gsl_complex_rect(2,1));
 517     gsl_matrix_complex_set($matrix, 1, 0, gsl_complex_rect(2,-1));
 518     gsl_matrix_complex_set($matrix, 1, 1, gsl_complex_rect(1,0));
 519     gsl_matrix_complex_memcpy($transpose, $matrix);
 520     gsl_matrix_complex_transpose($transpose);
 521
 522     for my $row (0..1) {
 523         map { gsl_matrix_complex_set($transpose, $row, $_, gsl_complex_conjugate(gsl_matrix_complex_get($transpose, $row, $_))) } (0..1);
 524     }
 525
 526     my $upper_right = gsl_matrix_complex_get($matrix, 0, 1 );
 527     my $lower_left  = gsl_matrix_complex_get($matrix, 1, 0 );
 528
 529     ok( gsl_complex_eq( gsl_complex_conjugate($upper_right), $lower_left ), 'hermitian' );
 530 }
 531
 532 sub SET_ROW : Tests {
 533     my $m = Math::GSL::Matrix->new(3,3);
 534     $m->set_row(0, [1,2,3]);
 535     ok_similar([$m->row(0)->as_list], [1,2,3]);
 536 }
 537 sub SET_ROW_CHAINED : Tests {
 538     my $m = Math::GSL::Matrix->new(3,3);
 539     $m->set_row(1, [4,5,6])
 540       ->set_row(2, [9,8,7]);
 541     ok_similar([$m->row(1)->as_list], [4,5,6]);
 542     ok_similar([$m->row(2)->as_list], [9,8,7]);
 543 }
 544
 545 sub SET_COL_CHAINED : Tests {
 546     my $m = Math::GSL::Matrix->new(3,3);
 547     $m->set_col(1, [4,5,6])
 548       ->set_col(2, [9,8,7]);
 549     ok_similar([$m->col(1)->as_list], [4,5,6]);
 550     ok_similar([$m->col(2)->as_list], [9,8,7]);
 551 }
 552
 553 sub GSL_MATRIX_VIEW_ARRAY : Tests {
 554     my $array = [8,4,3,7];
 555     my $matrix_view = gsl_matrix_view_array ($array, 2,2);
 556     ok_similar([map { gsl_matrix_get($matrix_view->{matrix}, 0, $_) } 0..1], [8, 4]);
 557     ok_similar([map { gsl_matrix_get($matrix_view->{matrix}, 1, $_) } 0..1], [3, 7]);
 558 }
 559
 560 sub GSL_MATRIX_VIEW_ARRAY_WITH_TDA : Tests {
 561     my $array = [8,4,3,7,5];
 562     my $matrix_view = gsl_matrix_view_array_with_tda ($array, 2,2, 3);
 563     ok_similar([map { gsl_matrix_get($matrix_view->{matrix}, 0, $_) } 0..1], [8, 4]);
 564     ok_similar([map { gsl_matrix_get($matrix_view->{matrix}, 1, $_) } 0..1], [7, 5]);
 565 }
 566
 567 sub GSL_MATRIX_OO_ADDITION_CONSTANT : Tests {
 568     my $m = Math::GSL::Matrix->new(3,3);
 569     $m->set_col(1, [4,5,6])
 570       ->set_col(2, [9,8,7]);
 571     my $m2 = $m + 4;
 572     ok_similar([$m->col(2)->as_list], [9,8,7]);
 573     ok_similar([$m->col(1)->as_list], [4,5,6]);
 574     ok_similar([$m->col(0)->as_list], [0,0,0]);
 575
 576     ok_similar([$m2->col(1)->as_list], [8,9,10]);
 577     ok_similar([$m2->col(2)->as_list], [13,12,11]);
 578     ok_similar([$m2->col(0)->as_list], [4,4,4]);
 579
 580     my $m3 = 4 + $m;
 581     ok_similar([$m3->col(1)->as_list], [8,9,10]);
 582     ok_similar([$m3->col(2)->as_list], [13,12,11]);
 583     ok_similar([$m3->col(0)->as_list], [4,4,4]);
 584 }
 585
 586 sub GSL_MATRIX_OO_ADDITION_MATRICES : Tests {
 587     my $m = Math::GSL::Matrix->new(3,3);
 588     $m->set_col(1, [4,5,6])
 589       ->set_col(2, [9,8,7])
 590       ->set_col(0, [1,2,3]);
 591     my $m2 = $m + $m;
 592     ok_similar([$m->col(0)->as_list], [1,2,3]);
 593     ok_similar([$m->col(2)->as_list], [9,8,7]);
 594     ok_similar([$m->col(1)->as_list], [4,5,6]);
 595
 596     ok_similar([$m2->col(0)->as_list], [2,4,6]);
 597     ok_similar([$m2->col(1)->as_list], [8,10,12]);
 598     ok_similar([$m2->col(2)->as_list], [18,16,14]);
 599 }
 600
 601 sub GSL_MATRIX_OO_SUBSTRACTION_CONSTANT : Tests {
 602     my $m = Math::GSL::Matrix->new(3,3);
 603     $m->set_col(1, [4,5,6])
 604       ->set_col(2, [9,8,7]);
 605     my $m2 = $m - 4;
 606     ok_similar([$m->col(2)->as_list], [9,8,7]);
 607     ok_similar([$m->col(1)->as_list], [4,5,6]);
 608     ok_similar([$m->col(0)->as_list], [0,0,0]);
 609
 610     ok_similar([$m2->col(1)->as_list], [0,1,2]);
 611     ok_similar([$m2->col(2)->as_list], [5,4,3]);
 612     ok_similar([$m2->col(0)->as_list], [-4,-4,-4]);
 613 }
 614
 615 sub GSL_MATRIX_OO_SUBSTRACTION_MATRICES : Tests {
 616     my $m = Math::GSL::Matrix->new(3,3);
 617     my $m3 = Math::GSL::Matrix->new(3,3);
 618     $m->set_col(1, [4,5,6])
 619       ->set_col(2, [9,8,7])
 620       ->set_col(0, [1,2,3]);
 621     $m3->set_col(1, [1,2,3])
 622        ->set_col(2, [9,8,7])
 623        ->set_col(0, [1,2,3]);
 624     my $m2 = $m - $m3;
 625     ok_similar([$m->col(0)->as_list], [1,2,3]);
 626     ok_similar([$m->col(2)->as_list], [9,8,7]);
 627     ok_similar([$m->col(1)->as_list], [4,5,6]);
 628
 629     ok_similar([$m2->col(0)->as_list], [0,0,0]);
 630     ok_similar([$m2->col(1)->as_list], [3,3,3]);
 631     ok_similar([$m2->col(2)->as_list], [0,0,0]);
 632 }
 633
 634 sub GSL_MATRIX_OO_MULTIPLICATION_CONSTANT : Tests {
 635     my $m = Math::GSL::Matrix->new(3,3);
 636     $m->set_col(1, [4,5,6])
 637       ->set_col(2, [9,8,7]);
 638     my $m2 = $m * 4;
 639     ok_similar([$m->col(2)->as_list], [9,8,7]);
 640     ok_similar([$m->col(1)->as_list], [4,5,6]);
 641     ok_similar([$m->col(0)->as_list], [0,0,0]);
 642
 643     ok_similar([$m2->col(1)->as_list], [16,20,24]);
 644     ok_similar([$m2->col(2)->as_list], [36,32,28]);
 645     ok_similar([$m2->col(0)->as_list], [0,0,0]);
 646
 647     my $m3 = 4 * $m;
 648     ok_similar([$m3->col(1)->as_list], [16,20,24]);
 649     ok_similar([$m3->col(2)->as_list], [36,32,28]);
 650     ok_similar([$m3->col(0)->as_list], [0,0,0]);
 651 }
 652
 653 sub GSL_MATRIX_OO_MULTIPLICATION_MATRICES : Tests {
 654     my $m = Math::GSL::Matrix->new(3,3);
 655     my $m3 = Math::GSL::Matrix->new(3,3);
 656     $m->set_col(1, [4,5,6])
 657       ->set_col(2, [9,8,7])
 658       ->set_col(0, [1,2,3]);
 659     $m3->set_col(1, [1,2,3])
 660        ->set_col(2, [9,8,7])
 661        ->set_col(0, [1,2,3]);
 662     my $m2 = $m * $m3;
 663     ok_similar([$m->col(0)->as_list], [1,2,3]);
 664     ok_similar([$m->col(2)->as_list], [9,8,7]);
 665     ok_similar([$m->col(1)->as_list], [4,5,6]);
 666
 667     ok_similar([$m2->col(0)->as_list], [1,4,9]);
 668     ok_similar([$m2->col(1)->as_list], [4,10,18]);
 669     ok_similar([$m2->col(2)->as_list], [81,64,49]);
 670 }
 671 Test::Class->runtests;