[motoinzyniere]

Od rana z supportem Mathworks na telefonie. Problem dotyczy wersji 2019b (nowszych jeszcze nie sprawdzałem), 2017b działa OK. Mamy sobie dwie macierze.

a = [1, -8.392077756813356, 31.320231244770135, -68.227489240094530, 95.601443589819910, -89.357015966890490, 55.711712360873710, -22.341815453118137, 5.229278595243435, -0.544267372579095]
b = [1,1,1,1,1,1,1,1,1,1]

No i wykonujemy operacje:

c(1) = a*b(:)

c(2) = sum(a.*b)
W 2017b jest OK, w 2019b mamy dwa różne wyniki.
c(1) = 1.211586386773433e-09
c(2) = 1.211586164728828e-09


#matlab #programowanie

[Zendemion]

@motoinzyniere: nie żebym usprawiedliwiał ale co takiego liczysz że ósma cyfra znacząca jest dla Ciebie ważna?

[groman43]

@motoinzyniere: Prawdziwi mężczyźni nie używają floating point ( ͡° ͜ʖ ͡°)

[moonlisa]

@motoinzyniere: A co zwraca a*b' ?

[motoinzyniere]

@Zendemion: Automatyczne testy konwersji wyników z jednego softu do drugiego. Ten syntax a*b(:) jest uzywany przy funcji butter wbudowanej w matlaba. Wiec niech to naprawią.
@moonlisa: == c(1).

[P.....1]

@motoinzyniere: pamiętam jak na studiach używałem pirackiej wersji kiedyś to było.

[moonlisa]

@motoinzyniere: Sprawdziłem na 2021a. Jest jak u Ciebie. Śmieszna sprawa. xD

[Okcydent]

@motoinzyniere: W sumie ciekawe. Pewnie pod spodem rozchodzi się o kolejność dodawania - można by pobawić się z sortowaniem i podejrzeć wyniki.

[motoinzyniere]

@Positano1: Te czasy dawno za mną.
@Cryptograph: Mamy dedykowane hotline do problemów.

[moonlisa]

@motoinzyniere: Jest jeszcze śmieszniej, bo się przyjrzałem co ja dostałem xD

a = [1, -8.392077756813356, 31.320231244770135, -68.227489240094530, 95.601443589819910, -89.357015966890490, 55.711712360873710, -22.341815453118137, 5.229278595243435, -0.544267372579095]
b = [1,1,1,1,1,1,1,1,1,1]

a*b(:)
1.21158549859501e-09

sum(a.*b)
1.21158572063962e-09

Jeśli wszystko się dobrze przenosiło, bo kopiuję z telefonu na kompa przez pastebin i na odwrót, to mamy inne wyniki liczbowe xDDDD

[Okcydent]

@motoinzyniere: Na 2018b jest coś takiego:

sum(a)

ans =

1.211585720639619e-09

> sum(sort(a))

ans =

1.211603262163408e-09

> sum(sort(a,'descend'))

ans =

1.211574840453977e-09

> c(1)

ans =

1.211574840453977e-09

> c(2)

ans =

1.211585720639619e-09

[motoinzyniere]

@moonlisa: Grubo...

[Okcydent]

@motoinzyniere: Z tego co widzę to kolejność wykonywania działań. Eh, to nie matematyka, że dodawanie jest sobie łączne.

[moonlisa]

@motoinzyniere: Daj znać jak coś wymyślą, bo mam z tego bekę xD

Nie wiem oczywiście co się zmieniło w kodzie Matlaba, ale przetestowałem to sobie na pythonie. Jest różnica przy przechowywaniu danych jako double a long double. Podejrzewam, że tutaj jest grane coś podobnego, gdzie wkrada się błąd precyzji liczb zmiennoprzecinkowych do tablicy b z jedynkami, a pewnie też i tablicy a. Po drodze pewnie są różne konwersje danych. Mieli lepiej

[moonlisa]

@motoinzyniere: Jest już aktualizacja do 2021b. Mogę stestować jak chcesz. xD

[motoinzyniere]

@moonlisa: To prerelease - tez dostałem dzisiaj maila.