@@ -211,11 +211,11 @@ import MultivariateOrthogonalPolynomials: tri_forwardrecurrence, grid, TriangleR
211
211
if N > 1
212
212
@test TriangleRecurrenceC (N, X, Y) ≈ B⁺ (N)* [Cˣ; Cʸ]
213
213
end
214
-
214
+
215
215
x,y = 0.1 ,0.2
216
216
v = randn (N)
217
217
w = randn (N+ 1 )
218
- @test xy_muladd! ((x,y), TriangleRecurrenceA (N,X,Y), v, 2.0 , copy (w)) ≈
218
+ @test xy_muladd! ((x,y), TriangleRecurrenceA (N,X,Y), v, 2.0 , copy (w)) ≈
219
219
B⁺ (N) * [x* Eye (N); y* Eye (N)]* v + 2 w
220
220
@test xy_muladd! ((x,y), TriangleRecurrenceA (N,X,Y)' , w, 2.0 , copy (v)) ≈
221
221
(B⁺ (N) * [x* Eye (N); y* Eye (N)])' * w + 2 v
@@ -228,9 +228,9 @@ import MultivariateOrthogonalPolynomials: tri_forwardrecurrence, grid, TriangleR
228
228
229
229
if N > 1
230
230
u = randn (N- 1 )
231
- @test muladd! (3.0 , TriangleRecurrenceC (N, X, Y), u, 2.0 , copy (w)) ≈
231
+ @test muladd! (3.0 , TriangleRecurrenceC (N, X, Y), u, 2.0 , copy (w)) ≈
232
232
3 B⁺ (N)* [Cˣ; Cʸ]* u + 2 w
233
- @test muladd! (3.0 , TriangleRecurrenceC (N, X, Y)' , w, 2.0 , copy (u)) ≈
233
+ @test muladd! (3.0 , TriangleRecurrenceC (N, X, Y)' , w, 2.0 , copy (u)) ≈
234
234
3 (B⁺ (N)* [Cˣ; Cʸ])' * w + 2 u
235
235
236
236
# need in-place to minimise buffers in Clenshaw
0 commit comments