1
+
2
+ <!doctype html>
3
+ < html lang ="en ">
4
+
5
+ < head >
6
+ < title > Code coverage report for lapack/base/dpttrf/lib/base.js</ title >
7
+ < meta charset ="utf-8 " />
8
+ < link rel ="stylesheet " href ="../../../../prettify.css " />
9
+ < link rel ="stylesheet " href ="../../../../base.css " />
10
+ < link rel ="shortcut icon " type ="image/x-icon " href ="../../../../favicon.png " />
11
+ < meta name ="viewport " content ="width=device-width, initial-scale=1 " />
12
+ < style type ='text/css '>
13
+ .coverage-summary .sorter {
14
+ background-image : url(../ ../ ../ ../ sor t-arrow-sprite.png);
15
+ }
16
+ </ style >
17
+ </ head >
18
+
19
+ < body >
20
+ < div class ='wrapper '>
21
+ < div class ='pad1 '>
22
+ < h1 > < a href ="../../../../index.html "> All files</ a > / < a href ="index.html "> lapack/base/dpttrf/lib</ a > base.js</ h1 >
23
+ < div class ='clearfix '>
24
+
25
+ < div class ='fl pad1y space-right2 '>
26
+ < span class ="strong "> 100% </ span >
27
+ < span class ="quiet "> Statements</ span >
28
+ < span class ='fraction '> 83/83</ span >
29
+ </ div >
30
+
31
+
32
+ < div class ='fl pad1y space-right2 '>
33
+ < span class ="strong "> 100% </ span >
34
+ < span class ="quiet "> Branches</ span >
35
+ < span class ='fraction '> 10/10</ span >
36
+ </ div >
37
+
38
+
39
+ < div class ='fl pad1y space-right2 '>
40
+ < span class ="strong "> 100% </ span >
41
+ < span class ="quiet "> Functions</ span >
42
+ < span class ='fraction '> 1/1</ span >
43
+ </ div >
44
+
45
+
46
+ < div class ='fl pad1y space-right2 '>
47
+ < span class ="strong "> 100% </ span >
48
+ < span class ="quiet "> Lines</ span >
49
+ < span class ='fraction '> 83/83</ span >
50
+ </ div >
51
+
52
+
53
+ </ div >
54
+ < p class ="quiet ">
55
+ Press < em > n</ em > or < em > j</ em > to go to the next uncovered block, < em > b</ em > , < em > p</ em > or < em > k</ em > for the previous block.
56
+ </ p >
57
+ < template id ="filterTemplate ">
58
+ < div class ="quiet ">
59
+ Filter:
60
+ < input type ="search " id ="fileSearch ">
61
+ </ div >
62
+ </ template >
63
+ </ div >
64
+ < div class ='status-line high '> </ div >
65
+ < pre > < table class ="coverage ">
66
+ < tr > < td class ="line-count quiet "> < a name ='L1 '> </ a > < a href ='#L1 '> 1</ a >
67
+ < a name ='L2 '> </ a > < a href ='#L2 '> 2</ a >
68
+ < a name ='L3 '> </ a > < a href ='#L3 '> 3</ a >
69
+ < a name ='L4 '> </ a > < a href ='#L4 '> 4</ a >
70
+ < a name ='L5 '> </ a > < a href ='#L5 '> 5</ a >
71
+ < a name ='L6 '> </ a > < a href ='#L6 '> 6</ a >
72
+ < a name ='L7 '> </ a > < a href ='#L7 '> 7</ a >
73
+ < a name ='L8 '> </ a > < a href ='#L8 '> 8</ a >
74
+ < a name ='L9 '> </ a > < a href ='#L9 '> 9</ a >
75
+ < a name ='L10 '> </ a > < a href ='#L10 '> 10</ a >
76
+ < a name ='L11 '> </ a > < a href ='#L11 '> 11</ a >
77
+ < a name ='L12 '> </ a > < a href ='#L12 '> 12</ a >
78
+ < a name ='L13 '> </ a > < a href ='#L13 '> 13</ a >
79
+ < a name ='L14 '> </ a > < a href ='#L14 '> 14</ a >
80
+ < a name ='L15 '> </ a > < a href ='#L15 '> 15</ a >
81
+ < a name ='L16 '> </ a > < a href ='#L16 '> 16</ a >
82
+ < a name ='L17 '> </ a > < a href ='#L17 '> 17</ a >
83
+ < a name ='L18 '> </ a > < a href ='#L18 '> 18</ a >
84
+ < a name ='L19 '> </ a > < a href ='#L19 '> 19</ a >
85
+ < a name ='L20 '> </ a > < a href ='#L20 '> 20</ a >
86
+ < a name ='L21 '> </ a > < a href ='#L21 '> 21</ a >
87
+ < a name ='L22 '> </ a > < a href ='#L22 '> 22</ a >
88
+ < a name ='L23 '> </ a > < a href ='#L23 '> 23</ a >
89
+ < a name ='L24 '> </ a > < a href ='#L24 '> 24</ a >
90
+ < a name ='L25 '> </ a > < a href ='#L25 '> 25</ a >
91
+ < a name ='L26 '> </ a > < a href ='#L26 '> 26</ a >
92
+ < a name ='L27 '> </ a > < a href ='#L27 '> 27</ a >
93
+ < a name ='L28 '> </ a > < a href ='#L28 '> 28</ a >
94
+ < a name ='L29 '> </ a > < a href ='#L29 '> 29</ a >
95
+ < a name ='L30 '> </ a > < a href ='#L30 '> 30</ a >
96
+ < a name ='L31 '> </ a > < a href ='#L31 '> 31</ a >
97
+ < a name ='L32 '> </ a > < a href ='#L32 '> 32</ a >
98
+ < a name ='L33 '> </ a > < a href ='#L33 '> 33</ a >
99
+ < a name ='L34 '> </ a > < a href ='#L34 '> 34</ a >
100
+ < a name ='L35 '> </ a > < a href ='#L35 '> 35</ a >
101
+ < a name ='L36 '> </ a > < a href ='#L36 '> 36</ a >
102
+ < a name ='L37 '> </ a > < a href ='#L37 '> 37</ a >
103
+ < a name ='L38 '> </ a > < a href ='#L38 '> 38</ a >
104
+ < a name ='L39 '> </ a > < a href ='#L39 '> 39</ a >
105
+ < a name ='L40 '> </ a > < a href ='#L40 '> 40</ a >
106
+ < a name ='L41 '> </ a > < a href ='#L41 '> 41</ a >
107
+ < a name ='L42 '> </ a > < a href ='#L42 '> 42</ a >
108
+ < a name ='L43 '> </ a > < a href ='#L43 '> 43</ a >
109
+ < a name ='L44 '> </ a > < a href ='#L44 '> 44</ a >
110
+ < a name ='L45 '> </ a > < a href ='#L45 '> 45</ a >
111
+ < a name ='L46 '> </ a > < a href ='#L46 '> 46</ a >
112
+ < a name ='L47 '> </ a > < a href ='#L47 '> 47</ a >
113
+ < a name ='L48 '> </ a > < a href ='#L48 '> 48</ a >
114
+ < a name ='L49 '> </ a > < a href ='#L49 '> 49</ a >
115
+ < a name ='L50 '> </ a > < a href ='#L50 '> 50</ a >
116
+ < a name ='L51 '> </ a > < a href ='#L51 '> 51</ a >
117
+ < a name ='L52 '> </ a > < a href ='#L52 '> 52</ a >
118
+ < a name ='L53 '> </ a > < a href ='#L53 '> 53</ a >
119
+ < a name ='L54 '> </ a > < a href ='#L54 '> 54</ a >
120
+ < a name ='L55 '> </ a > < a href ='#L55 '> 55</ a >
121
+ < a name ='L56 '> </ a > < a href ='#L56 '> 56</ a >
122
+ < a name ='L57 '> </ a > < a href ='#L57 '> 57</ a >
123
+ < a name ='L58 '> </ a > < a href ='#L58 '> 58</ a >
124
+ < a name ='L59 '> </ a > < a href ='#L59 '> 59</ a >
125
+ < a name ='L60 '> </ a > < a href ='#L60 '> 60</ a >
126
+ < a name ='L61 '> </ a > < a href ='#L61 '> 61</ a >
127
+ < a name ='L62 '> </ a > < a href ='#L62 '> 62</ a >
128
+ < a name ='L63 '> </ a > < a href ='#L63 '> 63</ a >
129
+ < a name ='L64 '> </ a > < a href ='#L64 '> 64</ a >
130
+ < a name ='L65 '> </ a > < a href ='#L65 '> 65</ a >
131
+ < a name ='L66 '> </ a > < a href ='#L66 '> 66</ a >
132
+ < a name ='L67 '> </ a > < a href ='#L67 '> 67</ a >
133
+ < a name ='L68 '> </ a > < a href ='#L68 '> 68</ a >
134
+ < a name ='L69 '> </ a > < a href ='#L69 '> 69</ a >
135
+ < a name ='L70 '> </ a > < a href ='#L70 '> 70</ a >
136
+ < a name ='L71 '> </ a > < a href ='#L71 '> 71</ a >
137
+ < a name ='L72 '> </ a > < a href ='#L72 '> 72</ a >
138
+ < a name ='L73 '> </ a > < a href ='#L73 '> 73</ a >
139
+ < a name ='L74 '> </ a > < a href ='#L74 '> 74</ a >
140
+ < a name ='L75 '> </ a > < a href ='#L75 '> 75</ a >
141
+ < a name ='L76 '> </ a > < a href ='#L76 '> 76</ a >
142
+ < a name ='L77 '> </ a > < a href ='#L77 '> 77</ a >
143
+ < a name ='L78 '> </ a > < a href ='#L78 '> 78</ a >
144
+ < a name ='L79 '> </ a > < a href ='#L79 '> 79</ a >
145
+ < a name ='L80 '> </ a > < a href ='#L80 '> 80</ a >
146
+ < a name ='L81 '> </ a > < a href ='#L81 '> 81</ a >
147
+ < a name ='L82 '> </ a > < a href ='#L82 '> 82</ a >
148
+ < a name ='L83 '> </ a > < a href ='#L83 '> 83</ a >
149
+ < a name ='L84 '> </ a > < a href ='#L84 '> 84</ a > </ td > < td class ="line-coverage quiet "> < span class ="cline-any cline-yes "> 3x</ span >
150
+ < span class ="cline-any cline-yes "> 3x</ span >
151
+ < span class ="cline-any cline-yes "> 3x</ span >
152
+ < span class ="cline-any cline-yes "> 3x</ span >
153
+ < span class ="cline-any cline-yes "> 3x</ span >
154
+ < span class ="cline-any cline-yes "> 3x</ span >
155
+ < span class ="cline-any cline-yes "> 3x</ span >
156
+ < span class ="cline-any cline-yes "> 3x</ span >
157
+ < span class ="cline-any cline-yes "> 3x</ span >
158
+ < span class ="cline-any cline-yes "> 3x</ span >
159
+ < span class ="cline-any cline-yes "> 3x</ span >
160
+ < span class ="cline-any cline-yes "> 3x</ span >
161
+ < span class ="cline-any cline-yes "> 3x</ span >
162
+ < span class ="cline-any cline-yes "> 3x</ span >
163
+ < span class ="cline-any cline-yes "> 3x</ span >
164
+ < span class ="cline-any cline-yes "> 3x</ span >
165
+ < span class ="cline-any cline-yes "> 3x</ span >
166
+ < span class ="cline-any cline-yes "> 3x</ span >
167
+ < span class ="cline-any cline-yes "> 3x</ span >
168
+ < span class ="cline-any cline-yes "> 3x</ span >
169
+ < span class ="cline-any cline-yes "> 3x</ span >
170
+ < span class ="cline-any cline-yes "> 3x</ span >
171
+ < span class ="cline-any cline-yes "> 3x</ span >
172
+ < span class ="cline-any cline-yes "> 3x</ span >
173
+ < span class ="cline-any cline-yes "> 3x</ span >
174
+ < span class ="cline-any cline-yes "> 3x</ span >
175
+ < span class ="cline-any cline-yes "> 3x</ span >
176
+ < span class ="cline-any cline-yes "> 3x</ span >
177
+ < span class ="cline-any cline-yes "> 3x</ span >
178
+ < span class ="cline-any cline-yes "> 3x</ span >
179
+ < span class ="cline-any cline-yes "> 3x</ span >
180
+ < span class ="cline-any cline-yes "> 3x</ span >
181
+ < span class ="cline-any cline-yes "> 3x</ span >
182
+ < span class ="cline-any cline-yes "> 3x</ span >
183
+ < span class ="cline-any cline-yes "> 3x</ span >
184
+ < span class ="cline-any cline-yes "> 3x</ span >
185
+ < span class ="cline-any cline-yes "> 3x</ span >
186
+ < span class ="cline-any cline-yes "> 3x</ span >
187
+ < span class ="cline-any cline-yes "> 3x</ span >
188
+ < span class ="cline-any cline-yes "> 3x</ span >
189
+ < span class ="cline-any cline-yes "> 3x</ span >
190
+ < span class ="cline-any cline-yes "> 3x</ span >
191
+ < span class ="cline-any cline-yes "> 3x</ span >
192
+ < span class ="cline-any cline-yes "> 3x</ span >
193
+ < span class ="cline-any cline-yes "> 3x</ span >
194
+ < span class ="cline-any cline-yes "> 16x</ span >
195
+ < span class ="cline-any cline-yes "> 16x</ span >
196
+ < span class ="cline-any cline-yes "> 16x</ span >
197
+ < span class ="cline-any cline-yes "> 16x</ span >
198
+ < span class ="cline-any cline-yes "> 16x</ span >
199
+ < span class ="cline-any cline-yes "> 16x</ span >
200
+ < span class ="cline-any cline-yes "> 16x</ span >
201
+ < span class ="cline-any cline-yes "> 2x</ span >
202
+ < span class ="cline-any cline-yes "> 2x</ span >
203
+ < span class ="cline-any cline-yes "> 14x</ span >
204
+ < span class ="cline-any cline-yes "> 14x</ span >
205
+ < span class ="cline-any cline-yes "> 14x</ span >
206
+ < span class ="cline-any cline-yes "> 14x</ span >
207
+ < span class ="cline-any cline-yes "> 16x</ span >
208
+ < span class ="cline-any cline-yes "> 34x</ span >
209
+ < span class ="cline-any cline-yes "> 34x</ span >
210
+ < span class ="cline-any cline-yes "> 4x</ span >
211
+ < span class ="cline-any cline-yes "> 4x</ span >
212
+ < span class ="cline-any cline-yes "> 30x</ span >
213
+ < span class ="cline-any cline-yes "> 30x</ span >
214
+ < span class ="cline-any cline-yes "> 30x</ span >
215
+ < span class ="cline-any cline-yes "> 30x</ span >
216
+ < span class ="cline-any cline-yes "> 30x</ span >
217
+ < span class ="cline-any cline-yes "> 30x</ span >
218
+ < span class ="cline-any cline-yes "> 30x</ span >
219
+ < span class ="cline-any cline-yes "> 30x</ span >
220
+ < span class ="cline-any cline-yes "> 30x</ span >
221
+ < span class ="cline-any cline-yes "> 10x</ span >
222
+ < span class ="cline-any cline-yes "> 16x</ span >
223
+ < span class ="cline-any cline-yes "> 2x</ span >
224
+ < span class ="cline-any cline-yes "> 2x</ span >
225
+ < span class ="cline-any cline-yes "> 8x</ span >
226
+ < span class ="cline-any cline-yes "> 16x</ span >
227
+ < span class ="cline-any cline-yes "> 3x</ span >
228
+ < span class ="cline-any cline-yes "> 3x</ span >
229
+ < span class ="cline-any cline-yes "> 3x</ span >
230
+ < span class ="cline-any cline-yes "> 3x</ span >
231
+ < span class ="cline-any cline-yes "> 3x</ span >
232
+ < span class ="cline-any cline-neutral "> </ span > </ td > < td class ="text "> < pre class ="prettyprint lang-js "> /**
233
+ * @license Apache-2.0
234
+ *
235
+ * Copyright (c) 2024 The Stdlib Authors.
236
+ *
237
+ * Licensed under the Apache License, Version 2.0 (the "License");
238
+ * you may not use this file except in compliance with the License.
239
+ * You may obtain a copy of the License at
240
+ *
241
+ * http://www.apache.org/licenses/LICENSE-2.0
242
+ *
243
+ * Unless required by applicable law or agreed to in writing, software
244
+ * distributed under the License is distributed on an "AS IS" BASIS,
245
+ * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
246
+ * See the License for the specific language governing permissions and
247
+ * limitations under the License.
248
+ */
249
+
250
+ 'use strict';
251
+
252
+ // MAIN //
253
+
254
+ /**
255
+ * Computes the `L * D * L^T` factorization of a real symmetric positive definite tridiagonal matrix `A`.
256
+ *
257
+ * @private
258
+ * @param {NonNegativeInteger} N - order of matrix `A`
259
+ * @param {Float64Array} D - the `N` diagonal elements of `A`
260
+ * @param {integer} strideD - stride length for `D`
261
+ * @param {NonNegativeInteger} offsetD - starting index of `D`
262
+ * @param {Float64Array} E - the `N-1` subdiagonal elements of `A`
263
+ * @param {integer} strideE - stride length for `E`
264
+ * @param {NonNegativeInteger} offsetE - starting index of `E`
265
+ * @returns {integer} status code
266
+ *
267
+ * @example
268
+ * var Float64Array = require( '@stdlib/array/float64' );
269
+ *
270
+ * var D = new Float64Array( [ 4.0, 5.0, 6.0 ] );
271
+ * var E = new Float64Array( [ 1.0, 2.0 ] );
272
+ *
273
+ * dpttrf( 3, D, 1, 0, E, 1, 0 );
274
+ * // D => <Float64Array>[ 4, 4.75, ~5.15789 ]
275
+ * // E => <Float64Array>[ 0.25, ~0.4210 ]
276
+ */
277
+ function dpttrf( N, D, strideD, offsetD, E, strideE, offsetE ) {
278
+ var id;
279
+ var ie;
280
+ var v;
281
+ var i;
282
+
283
+ if ( N === 0 ) {
284
+ return 0;
285
+ }
286
+ ie = offsetE;
287
+ id = offsetD;
288
+
289
+ // Compute the `L * D * L^T` factorization of `A`...
290
+ for ( i = 0; i < N-1; i++ ) {
291
+ // If `D[k] <= 0`, then the matrix is not positive definite...
292
+ if ( D[ id ] <= 0.0 ) {
293
+ return i+1;
294
+ }
295
+ // Solve for E[k] and D[k+1]...
296
+ v = E[ ie ];
297
+ E[ ie ] = v / D[ id ];
298
+
299
+ id += strideD;
300
+ D[ id ] -= E[ ie ] * v;
301
+
302
+ ie += strideE;
303
+ }
304
+ // Check `D[k]` for positive definiteness...
305
+ if ( D[ id ] <= 0.0 ) {
306
+ return N;
307
+ }
308
+ return 0;
309
+ }
310
+
311
+
312
+ // EXPORTS //
313
+
314
+ module.exports = dpttrf;
315
+ </ pre > </ td > </ tr > </ table > </ pre >
316
+
317
+ < div class ='push '> </ div > <!-- for sticky footer -->
318
+ </ div > <!-- /wrapper -->
319
+ < div class ='footer quiet pad2 space-top1 center small '>
320
+ Code coverage generated by
321
+ < a href ="https://istanbul.js.org/ " target ="_blank " rel ="noopener noreferrer "> istanbul</ a >
322
+ at 2024-07-28T11:31:54.082Z
323
+ </ div >
324
+ < script src ="../../../../prettify.js "> </ script >
325
+ < script >
326
+ window . onload = function ( ) {
327
+ prettyPrint ( ) ;
328
+ } ;
329
+ </ script >
330
+ < script src ="../../../../sorter.js "> </ script >
331
+ < script src ="../../../../block-navigation.js "> </ script >
332
+ </ body >
333
+ </ html >
334
+
0 commit comments