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
|
;; -*-theme-d-*-
;; Copyright (C) 2014, 2024 Tommi Höynälänmaa
;; Distributed under GNU General Public License version 3,
;; see file doc/GPL-3.
;; Expected results: translation and running OK
(define-proper-program (tests test709)
(import (standard-library core)
(standard-library list-utilities)
(standard-library matrix)
(tests numerical-test-env))
(define-param-proc my-matrix (%number)
(((elements (:uniform-list (:uniform-list %number))))
(:matrix %number)
force-pure)
(let ((rows (length elements)))
(assert (> rows 0))
(let ((columns (length (gen-car elements))))
(assert (for-all?
(lambda (((row (:uniform-list %number))) <boolean> (pure))
(= (length row) columns))
(gen-cdr elements)))
(let ((result (make-matrix rows columns (zero %number))))
(do ((i-row <integer> 0 (+ i-row 1))
(cur-elements (:uniform-list (:uniform-list %number))
elements (gen-cdr cur-elements)))
((>= i-row rows))
(do ((i-column <integer> 0 (+ i-column 1))
(cur-elements2 (:uniform-list %number)
(gen-car cur-elements) (gen-cdr cur-elements2)))
((>= i-column columns))
(matrix-set! result i-row i-column (gen-car cur-elements2))))
result))))
(define-param-proc my-diagonal-matrix (%number)
(((elements (:uniform-list %number)))
(:diagonal-matrix %number)
force-pure)
(let* ((i-len (length elements))
(mx (make-diagonal-matrix i-len (zero %number))))
(do ((i <integer> 0 (+ i 1))
(l-cur (:uniform-list %number) elements (gen-cdr l-cur)))
((>= i i-len))
(diagonal-matrix-set! mx i (gen-car l-cur)))
mx))
(define-main-proc (() <none> nonpure)
(let ((mx1 (my-matrix
(static-cast (:uniform-list (:uniform-list <real>))
(list
(list 1.0 0.0 2.0 0.0)
(list -1.5 2.5 0.0 5.0)
(list -3.0 -2.0 1.0 1.0)))))
(mx2 (my-matrix
(static-cast (:uniform-list (:uniform-list <real>))
(list
(list 1.0 0.0 2.0 0.0)
(list -1.5 2.5 0.0 5.0)
(list -3.0 -2.0 1.0 1.0)))))
(mx3 (my-matrix
(static-cast (:uniform-list (:uniform-list <real>))
(list
(list 1.0 0.0 2.0 0.0)
(list -1.5 2.5 0.0 5.0)
(list -3.0 -2.0 -1.0 1.0)))))
(mx4 (my-matrix
(static-cast (:uniform-list (:uniform-list <real>))
(list
(list 5.0 0.0 0.0)
(list 0.0 2.0 0.0)
(list 0.0 0.0 -1.0)))))
(mx5 (my-diagonal-matrix
(static-cast (:uniform-list <real>)
(list 5.0 2.0 -1.0))))
(mx6 (my-diagonal-matrix
(static-cast (:uniform-list <real>)
(list 5.0 2.0 -1.0))))
(mx7 (my-diagonal-matrix
(static-cast (:uniform-list <real>)
(list -15.0 3.0 -2.1))))
(mx8 (my-diagonal-matrix
(static-cast (:uniform-list <real>)
(list -15.0 3.0 -2.1 1.2)))))
(report-boolean-test (= mx1 mx1) #t)
(report-boolean-test (= mx1 mx2) #t)
(report-boolean-test (= mx1 mx3) #f)
(report-boolean-test (= mx1 mx4) #f)
(report-boolean-test (= mx5 mx5) #t)
(report-boolean-test (= mx5 mx6) #t)
(report-boolean-test (= mx5 mx7) #f)
(report-boolean-test (= mx7 mx8) #f))))
|
