@@ -39,16 +39,24 @@ var Float64Array = require( '@stdlib/array/float64' );
39
39
var Int32Array = require ( ' @stdlib/array/int32' );
40
40
var dgttrf = require ( ' @stdlib/lapack/base/dgttrf' );
41
41
42
- var DL = new Float64Array ( [ 1 .0 , 1 .0 ] );
43
- var D = new Float64Array ( [ 2 .0 , 3 .0 , 1 .0 ] );
44
- var DU = new Float64Array ( [ 1 .0 , 1 .0 ] );
42
+ var DL = new Float64Array ( [ 6 .0 , 6 .0 ] );
43
+ var D = new Float64Array ( [ 20 .0 , 30 .0 , 10 .0 ] );
44
+ var DU = new Float64Array ( [ 8 .0 , 8 .0 ] );
45
45
var DU2 = new Float64Array ( [ 0.0 ] );
46
46
var IPIV = new Int32Array ( [ 0 , 0 , 0 ] );
47
47
48
+ /*
49
+ A = [
50
+ [ 20.0, 8.0, 0.0 ],
51
+ [ 6.0, 30.0, 8.0 ],
52
+ [ 0.0, 6.0, 10.0 ]
53
+ ]
54
+ */
55
+
48
56
dgttrf ( 3 , DL , D , DU , DU2 , IPIV );
49
- // DL => <Float64Array>[ 0.5, 0.4 ]
50
- // D => <Float64Array>[ 2 .0, 2.5, 0.6 ]
51
- // DU => <Float64Array>[ 1 .0, 1 .0 ]
57
+ // DL => <Float64Array>[ 0.3, ~0.22 ]
58
+ // D => <Float64Array>[ 20 .0, 27.6, ~8.26 ]
59
+ // DU => <Float64Array>[ 8 .0, 8 .0 ]
52
60
// DU2 => <Float64Array>[ 0.0 ]
53
61
// IPIV => <Int32Array>[ 0, 1, 2 ]
54
62
```
@@ -72,12 +80,20 @@ var Int32Array = require( '@stdlib/array/int32' );
72
80
var dgttrf = require ( ' @stdlib/lapack/base/dgttrf' );
73
81
74
82
// Initial arrays...
75
- var DL0 = new Float64Array ( [ 0.0 , 1 .0 , 1 .0 ] );
76
- var D0 = new Float64Array ( [ 0.0 , 2 .0 , 3 .0 , 1 .0 ] );
77
- var DU0 = new Float64Array ( [ 0.0 , 1 .0 , 1 .0 ] );
83
+ var DL0 = new Float64Array ( [ 0.0 , 6 .0 , 6 .0 ] );
84
+ var D0 = new Float64Array ( [ 0.0 , 20 .0 , 30 .0 , 10 .0 ] );
85
+ var DU0 = new Float64Array ( [ 0.0 , 8 .0 , 8 .0 ] );
78
86
var DU20 = new Float64Array ( [ 0.0 , 0.0 ] );
79
87
var IPIV0 = new Int32Array ( [ 0 , 0 , 0 , 0 ] );
80
88
89
+ /*
90
+ A = [
91
+ [ 20.0, 8.0, 0.0 ],
92
+ [ 6.0, 30.0, 8.0 ],
93
+ [ 0.0, 6.0, 10.0 ]
94
+ ]
95
+ */
96
+
81
97
// Create offset views...
82
98
var DL = new Float64Array ( DL0 .buffer , DL0 .BYTES_PER_ELEMENT * 1 ); // start at 2nd element
83
99
var D = new Float64Array ( D0 .buffer , D0 .BYTES_PER_ELEMENT * 1 ); // start at 2nd element
@@ -86,9 +102,9 @@ var DU2 = new Float64Array( DU20.buffer, DU20.BYTES_PER_ELEMENT*1 ); // start at
86
102
var IPIV = new Int32Array ( IPIV0 .buffer , IPIV0 .BYTES_PER_ELEMENT * 1 ); // start at 2nd element
87
103
88
104
dgttrf ( 3 , DL , D , DU , DU2 , IPIV );
89
- // DL0 => <Float64Array>[ 0.0, 0.5, 0.4 ]
90
- // D0 => <Float64Array>[ 0.0, 2 .0, 2.5, 0.6 ]
91
- // DU0 => <Float64Array>[ 0.0, 1 .0, 1 .0 ]
105
+ // DL0 => <Float64Array>[ 0.0, 0.3, ~0.22 ]
106
+ // D0 => <Float64Array>[ 0.0, 20 .0, 27.6, ~8.26 ]
107
+ // DU0 => <Float64Array>[ 0.0, 8 .0, 8 .0 ]
92
108
// DU20 => <Float64Array>[ 0.0, 0.0 ]
93
109
// IPIV0 => <Int32Array>[ 0, 0, 1, 2 ]
94
110
```
@@ -104,16 +120,24 @@ var Float64Array = require( '@stdlib/array/float64' );
104
120
var Int32Array = require ( ' @stdlib/array/int32' );
105
121
var dgttrf = require ( ' @stdlib/lapack/base/dgttrf' );
106
122
107
- var DL = new Float64Array ( [ 1 .0 , 1 .0 ] );
108
- var D = new Float64Array ( [ 2 .0 , 3 .0 , 1 .0 ] );
109
- var DU = new Float64Array ( [ 1 .0 , 1 .0 ] );
123
+ var DL = new Float64Array ( [ 6 .0 , 6 .0 ] );
124
+ var D = new Float64Array ( [ 20 .0 , 30 .0 , 10 .0 ] );
125
+ var DU = new Float64Array ( [ 8 .0 , 8 .0 ] );
110
126
var DU2 = new Float64Array ( [ 0.0 ] );
111
127
var IPIV = new Int32Array ( [ 0 , 0 , 0 ] );
112
128
129
+ /*
130
+ A = [
131
+ [ 20.0, 8.0, 0.0 ],
132
+ [ 6.0, 30.0, 8.0 ],
133
+ [ 0.0, 6.0, 10.0 ]
134
+ ]
135
+ */
136
+
113
137
dgttrf .ndarray ( 3 , DL , 1 , 0 , D , 1 , 0 , DU , 1 , 0 , DU2 , 1 , 0 , IPIV , 1 , 0 );
114
- // DL => <Float64Array>[ 0.5, 0.4 ]
115
- // D => <Float64Array>[ 2 .0, 2.5, 0.6 ]
116
- // DU => <Float64Array>[ 1 .0, 1 .0 ]
138
+ // DL => <Float64Array>[ 0.3, ~0.22 ]
139
+ // D => <Float64Array>[ 20 .0, 27.6, ~8.26 ]
140
+ // DU => <Float64Array>[ 8 .0, 8 .0 ]
117
141
// DU2 => <Float64Array>[ 0.0 ]
118
142
// IPIV => <Int32Array>[ 0, 1, 2 ]
119
143
```
@@ -140,16 +164,24 @@ var Float64Array = require( '@stdlib/array/float64' );
140
164
var Int32Array = require ( ' @stdlib/array/int32' );
141
165
var dgttrf = require ( ' @stdlib/lapack/base/dgttrf' );
142
166
143
- var DL = new Float64Array ( [ 0.0 , 1 .0 , 1 .0 ] );
144
- var D = new Float64Array ( [ 0.0 , 2 .0 , 3 .0 , 1 .0 ] );
145
- var DU = new Float64Array ( [ 0.0 , 1 .0 , 1 .0 ] );
167
+ var DL = new Float64Array ( [ 0.0 , 6 .0 , 6 .0 ] );
168
+ var D = new Float64Array ( [ 0.0 , 20 .0 , 30 .0 , 10 .0 ] );
169
+ var DU = new Float64Array ( [ 0.0 , 8 .0 , 8 .0 ] );
146
170
var DU2 = new Float64Array ( [ 0.0 , 0.0 ] );
147
171
var IPIV = new Int32Array ( [ 0 , 0 , 0 , 0 ] );
148
172
173
+ /*
174
+ A = [
175
+ [ 20.0, 8.0, 0.0 ],
176
+ [ 6.0, 30.0, 8.0 ],
177
+ [ 0.0, 6.0, 10.0 ]
178
+ ]
179
+ */
180
+
149
181
dgttrf .ndarray ( 3 , DL , 1 , 1 , D , 1 , 1 , DU , 1 , 1 , DU2 , 1 , 1 , IPIV , 1 , 1 );
150
- // DL => <Float64Array>[ 0.0, 0.5, 0.4 ]
151
- // D => <Float64Array>[ 0.0, 2 .0, 2.5, 0.6 ]
152
- // DU => <Float64Array>[ 0.0, 1 .0, 1 .0 ]
182
+ // DL => <Float64Array>[ 0.0, 0.3, ~0.22 ]
183
+ // D => <Float64Array>[ 0.0, 20 .0, 27.6, ~8.26 ]
184
+ // DU => <Float64Array>[ 0.0, 8 .0, 8 .0 ]
153
185
// DU2 => <Float64Array>[ 0.0, 0.0 ]
154
186
// IPIV => <Int32Array>[ 0, 0, 1, 2 ]
155
187
```
0 commit comments