Fix coeffs in trivial dimensions to work with BLAS.

This commit is contained in:
Alex Shinn 2023-06-04 06:05:29 +09:00
parent 6ea80c5ea4
commit 4b5e885f31

View file

@ -328,6 +328,9 @@
res)) res))
((= (+ 1 (interval-lower-bound domain (- i 1))) ((= (+ 1 (interval-lower-bound domain (- i 1)))
(interval-upper-bound domain (- i 1))) (interval-upper-bound domain (- i 1)))
(vector-set! res i (if (< i (interval-dimension domain))
(interval-width domain i)
1))
(lp (+ i 1) (cdr ls) offset count)) (lp (+ i 1) (cdr ls) offset count))
(else (else
(let ((dir (if (and (> count 0) (let ((dir (if (and (> count 0)
@ -396,7 +399,9 @@
res) res)
((= (+ 1 (interval-lower-bound domain i)) ((= (+ 1 (interval-lower-bound domain i))
(interval-upper-bound domain i)) (interval-upper-bound domain i))
(vector-set! res (+ i 1) 0) (vector-set! res (+ i 1) (if (< (+ i 1) dim)
(interval-width domain (+ i 1))
1))
(lp (- i 1) scale)) (lp (- i 1) scale))
(else (else
(let ((coeff (* scale (- (interval-upper-bound domain i) (let ((coeff (* scale (- (interval-upper-bound domain i)