File: rtest_dgeqrf.mac

package info (click to toggle)
maxima 5.44.0-2
  • links: PTS
  • area: main
  • in suites: bullseye, sid
  • size: 108,368 kB
  • sloc: lisp: 383,860; fortran: 14,665; perl: 14,369; tcl: 11,147; sh: 4,517; makefile: 2,580; ansic: 447; python: 262; xml: 59; awk: 37; sed: 17
file content (99 lines) | stat: -rw-r--r-- 2,159 bytes parent folder | download | duplicates (9)
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
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
/* rtest_dgeqrf.mac -- test script for dgeqrf
 * Copyright 2011 by Robert Dodier
 * I release this work under terms of the GNU General Public License.
 */

(if not ?fboundp ('dgeqrf) then load (lapack), 0);
0;

(random_matrix () := block ([nr : 5 + random (10), nc : 5 + random (10)],
  genmatrix (lambda ([i, j], random (100.0) - 50), nr, nc)), 0);
0;

set_random_state (make_random_state (12345));
done;

block ([a, q, r],
  a : random_matrix (),
  [q, r] : dgeqrf (a),
  mat_norm (q . r - a, 1),
  if %% < 1e-12 then true else [a, q, r, %%]);
true;

block ([a, q, r],
  a : random_matrix (),
  [q, r] : dgeqrf (a),
  mat_norm (q . r - a, 1),
  if %% < 1e-12 then true else [a, q, r, %%]);
true;

block ([a, q, r],
  a : random_matrix (),
  [q, r] : dgeqrf (a),
  mat_norm (q . r - a, 1),
  if %% < 1e-12 then true else [a, q, r, %%]);
true;

block ([a, q, r],
  a : random_matrix (),
  [q, r] : dgeqrf (a),
  mat_norm (q . r - a, 1),
  if %% < 1e-12 then true else [a, q, r, %%]);
true;

block ([a, q, r],
  a : random_matrix (),
  [q, r] : dgeqrf (a),
  mat_norm (q . r - a, 1),
  if %% < 1e-12 then true else [a, q, r, %%]);
true;

block ([a, q, r],
  a : random_matrix (),
  [q, r] : dgeqrf (a),
  mat_norm (q . r - a, 1),
  if %% < 1e-12 then true else [a, q, r, %%]);
true;

block ([a, q, r],
  a : random_matrix (),
  [q, r] : dgeqrf (a),
  mat_norm (q . r - a, 1),
  if %% < 1e-12 then true else [a, q, r, %%]);
true;

block ([a, q, r],
  a : random_matrix (),
  [q, r] : dgeqrf (a),
  mat_norm (q . r - a, 1),
  if %% < 1e-12 then true else [a, q, r, %%]);
true;

block ([a, q, r],
  a : random_matrix (),
  [q, r] : dgeqrf (a),
  mat_norm (q . r - a, 1),
  if %% < 1e-12 then true else [a, q, r, %%]);
true;

block ([a, q, r],
  a : random_matrix (),
  [q, r] : dgeqrf (a),
  mat_norm (q . r - a, 1),
  if %% < 1e-12 then true else [a, q, r, %%]);
true;

block ([a, q, r],
  a : random_matrix (),
  [q, r] : dgeqrf (a),
  mat_norm (q . r - a, 1),
  if %% < 1e-12 then true else [a, q, r, %%]);
true;

block ([a, q, r],
  a : random_matrix (),
  [q, r] : dgeqrf (a),
  mat_norm (q . r - a, 1),
  if %% < 1e-12 then true else [a, q, r, %%]);
true;