[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Octave-bug-tracker] [bug #60552] BIST for lscov.m tolerance failure wit
From: |
Mike Miller |
Subject: |
[Octave-bug-tracker] [bug #60552] BIST for lscov.m tolerance failure with OpenBLAS (v0.3.13–v0.3.15) |
Date: |
Fri, 7 May 2021 15:17:30 -0400 (EDT) |
User-agent: |
Mozilla/5.0 (X11; Linux x86_64) AppleWebKit/537.36 (KHTML, like Gecko) Chrome/90.0.4430.93 Safari/537.36 |
URL:
<https://savannah.gnu.org/bugs/?60552>
Summary: BIST for lscov.m tolerance failure with OpenBLAS
(v0.3.13–v0.3.15)
Project: GNU Octave
Submitted by: mtmiller
Submitted on: Fri 07 May 2021 12:17:28 PM PDT
Category: Test Suite
Severity: 3 - Normal
Priority: 5 - Normal
Item Group: Inaccurate Result
Status: None
Assigned to: None
Originator Name:
Originator Email:
Open/Closed: Open
Release: dev
Discussion Lock: Any
Operating System: GNU/Linux
_______________________________________________________
Details:
In Octave version 6.2.0 or version 7.0.0 development branch, with OpenBLAS
0.3.13 (on Debian native system):
>> test lscov.m
***** test
## Adapted from example in Matlab documentation
x1 = [.2 .5 .6 .8 1.0 1.1]';
x2 = [.1 .3 .4 .9 1.1 1.4]';
X = [ones(size(x1)) x1 x2];
y = [.17 .26 .28 .23 .27 .34]';
[b, se_b, mse, S] = lscov(X, y);
assert (b, [0.1203 0.3284 -0.1312]', 1E-4);
assert (se_b, [0.0643 0.2267 0.1488]', 1E-4);
assert (mse, 0.0015, 1E-4);
assert (S, [0.0041 -0.0130 0.0075; -0.0130 0.0514 -0.0328; 0.0075 -0.0328
0.0221], 1E-4);
w = [1 1 1 1 1 .1]';
[bw, sew_b, msew] = lscov (X, y, w);
assert (bw, [0.1046 0.4614 -0.2621]', 1E-4);
assert (sew_b, [0.0309 0.1152 0.0814]', 1E-4);
assert (msew, 3.4741e-004, -1E-4);
V = .2*ones (length (x1)) + .8*diag (ones (size (x1)));
[bg, sew_b, mseg] = lscov (X, y, V);
assert (bg, [0.1203 0.3284 -0.1312]', 1E-4);
assert (sew_b, [0.0672 0.2267 0.1488]', 1E-4);
assert (mseg, 0.0019, 1E-4);
y2 = [y 2*y];
[b2, se_b2, mse2, S2] = lscov (X, y2);
assert (b2, [b 2*b], 2*eps);
assert (se_b2, [se_b 2*se_b], 2*eps);
assert (mse2, [mse 4*mse], eps);
assert (S2(:, :, 1), S, eps);
assert (S2(:, :, 2), 4*S, eps);
!!!!! test failed
ASSERT errors for: assert (S2 (:, :, 2),4 * S,eps)
Location | Observed | Expected | Reason
(2,2) 0.20555 0.20555 Abs err 2.7756e-16 exceeds tol
2.2204e-16 by 6e-17
In Octave version 6.2.0 or version 7.0.0 development branch, with OpenBLAS
0.3.15, running in Flatpak container, same test, slightly different result:
!!!!! test failed
ASSERT errors for: assert (S2 (:, :, 2),4 * S,eps)
Location | Observed | Expected | Reason
(2,2) 0.20555 0.20555 Abs err 3.0531e-16 exceeds tol
2.2204e-16 by 8e-17
_______________________________________________________
Reply to this item at:
<https://savannah.gnu.org/bugs/?60552>
_______________________________________________
Message sent via Savannah
https://savannah.gnu.org/
- [Octave-bug-tracker] [bug #60552] BIST for lscov.m tolerance failure with OpenBLAS (v0.3.13–v0.3.15),
Mike Miller <=