1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
|
# modified from t/dgtsl.t in Photonic by W. Luis Mochán
use strict;
use warnings;
use PDL;
use PDL::NiceSlice;
use PDL::LinearAlgebra::Real;
use constant N=>10;
use Test::More tests => 2*N;
for my $D (3..N+2) { #first differences
#solve b_{n+1}-b_n=1 with homogeneous BCs
my $c=zeroes($D); $c->((-1)).=0;
my $d=-ones($D);
my $e=ones($D); $e->((-1)).=0;
my $b=ones($D); $b->((-1)).=1-$D;
my $info=pdl(short,0);
gtsv($c, $d, $e, $b, $info);
my $r=sequence($D);
ok($b->approx($r)->all, "1st diff in D=$D") or diag "info: ", $info, "\nGot: ", $b, "\nExpected: ", $r;
}
for my $D (3..N+2) { #second differences
#solve b_{n+1}-2{b_n}+b_{n-1}=1 with kinda homogeneous BCs
my $c=ones($D); $c->((-1)).=0;
my $d=-2*ones($D);
my $e=ones($D); $e->((-1)).=0;
my $b=ones($D);
my $info=pdl(short,0);
gtsv($c, $d, $e, $b, $info);
my $x=sequence($D);
my $r=-$D/2-($D-1)/2*$x+1/2*$x*$x;
ok($b->approx($r)->all, "2nd diff in D=$D") or diag "info: ", $info, "\nGot: ", $b, "\nExpected: ", $r;
}
|
