@@ -54,18 +54,34 @@ function test_cat()
54
54
B = simple_matrix (T, 100 , 10 )
55
55
C = simple_matrix (T, 100 , 90 )
56
56
D = [A B C]
57
- Ao = LinearOperator (T, 100 , 100 , false , false ,
58
- (res, v) -> mul! (res, A, v),
59
- (res, u) -> mul! (res, transpose (A), u),
60
- (res, w) -> mul! (res, adjoint (A), w))
61
- Bo = LinearOperator (T, size (B)... , false , false ,
62
- (res, v) -> mul! (res, B, v),
63
- (res, u) -> mul! (res, transpose (B), u),
64
- (res, w) -> mul! (res, adjoint (B), w))
65
- Co = LinearOperator (T, size (C)... , false , false ,
66
- (res, v) -> mul! (res, C, v),
67
- (res, u) -> mul! (res, transpose (C), u),
68
- (res, w) -> mul! (res, adjoint (C), w))
57
+ Ao = LinearOperator (
58
+ T,
59
+ 100 ,
60
+ 100 ,
61
+ false ,
62
+ false ,
63
+ (res, v) -> mul! (res, A, v),
64
+ (res, u) -> mul! (res, transpose (A), u),
65
+ (res, w) -> mul! (res, adjoint (A), w),
66
+ )
67
+ Bo = LinearOperator (
68
+ T,
69
+ size (B)... ,
70
+ false ,
71
+ false ,
72
+ (res, v) -> mul! (res, B, v),
73
+ (res, u) -> mul! (res, transpose (B), u),
74
+ (res, w) -> mul! (res, adjoint (B), w),
75
+ )
76
+ Co = LinearOperator (
77
+ T,
78
+ size (C)... ,
79
+ false ,
80
+ false ,
81
+ (res, v) -> mul! (res, C, v),
82
+ (res, u) -> mul! (res, transpose (C), u),
83
+ (res, w) -> mul! (res, adjoint (C), w),
84
+ )
69
85
Do = [Ao Bo Co]
70
86
Do2 = LinearOperator (D)
71
87
rhs = simple_vector (T, 200 )
@@ -106,18 +122,34 @@ function test_cat()
106
122
B = simple_matrix (T, 10 , 100 )
107
123
C = simple_matrix (T, 90 , 100 )
108
124
D = [A; B; C]
109
- Ao = LinearOperator (T, 100 , 100 , false , false ,
110
- (res, v) -> mul! (res, A, v),
111
- (res, u) -> mul! (res, transpose (A), u),
112
- (res, w) -> mul! (res, adjoint (A), w))
113
- Bo = LinearOperator (T, size (B)... , false , false ,
114
- (res, v) -> mul! (res, B, v),
115
- (res, u) -> mul! (res, transpose (B), u),
116
- (res, w) -> mul! (res, adjoint (B), w))
117
- Co = LinearOperator (T, size (C)... , false , false ,
118
- (res, v) -> mul! (res, C, v),
119
- (res, u) -> mul! (res, transpose (C), u),
120
- (res, w) -> mul! (res, adjoint (C), w))
125
+ Ao = LinearOperator (
126
+ T,
127
+ 100 ,
128
+ 100 ,
129
+ false ,
130
+ false ,
131
+ (res, v) -> mul! (res, A, v),
132
+ (res, u) -> mul! (res, transpose (A), u),
133
+ (res, w) -> mul! (res, adjoint (A), w),
134
+ )
135
+ Bo = LinearOperator (
136
+ T,
137
+ size (B)... ,
138
+ false ,
139
+ false ,
140
+ (res, v) -> mul! (res, B, v),
141
+ (res, u) -> mul! (res, transpose (B), u),
142
+ (res, w) -> mul! (res, adjoint (B), w),
143
+ )
144
+ Co = LinearOperator (
145
+ T,
146
+ size (C)... ,
147
+ false ,
148
+ false ,
149
+ (res, v) -> mul! (res, C, v),
150
+ (res, u) -> mul! (res, transpose (C), u),
151
+ (res, w) -> mul! (res, adjoint (C), w),
152
+ )
121
153
Do = [Ao; Bo; Co]
122
154
Do2 = LinearOperator (D)
123
155
rhs = simple_vector (T, 100 )
0 commit comments