macierze - kompilator szaleje?, błędne wyniki dla int, poprawne dla double

macierze - kompilator szaleje?, błędne wyniki dla int, poprawne dla double
PP
  • Rejestracja:ponad 13 lat
  • Ostatnio:ponad 12 lat
  • Postów:9
0

czy te kompilatory do reszty oszalały? :)

mam szablon macierzowy, algorytm wyznaczania wyznacznika macierzy kwadratowej. inicjuje zmienną m1(3,3) jakimiś poprawnymi wartościami: o coś takiego...

Kopiuj
    matrix2D<double> m10(3,3);
    m10(0,0)=1;
    m10(0,1)=2;
    m10(0,2)=5;
    m10(1,0)=73;
    m10(1,1)=4;
    m10(1,2)=754;
    m10(2,0)=4;
    m10(2,1)=54;
    m10(2,2)=4;
    cout<<"wyznacznik macierzy jest rowny: "<<m10.det()<<endl;
    //ZWRACA WYNIK POPRAWNY

teraz gdy zamiast double zrobie <int> m10(3,3); to zwraca błędny wynik. To się kupy d**y nie trzyma. Nie mogę wymyślić co kompilator robi z tym kodem, że zwraca on błędne wartości... na google specjalnie nie wiem nawet co wpisać. To taki błąd z którym nigdy się nie spotkałem i nie mam pojęcia skąd on się bierze. Zazwyczaj jak jest coś nie halo to conajmniej kojarzę gdzie może być problem. Tutaj nic, według logiki to powinno chodzić. Może jakiś spec się wypowie.

Endrju
  • Rejestracja:około 22 lata
  • Ostatnio:ponad rok
0

Kto Ci powiedział, że wyznacznik macierzy z intów będzie intem? Wychodzi źle bo wszystko jest pewnie zaokrąglane albo ucinane, ciężko powiedzieć, bo nie ma kodu.


"(...) otherwise, the behavior is undefined".
PP
  • Rejestracja:ponad 13 lat
  • Ostatnio:ponad 12 lat
  • Postów:9
0

no tak, przecież tam jest dzielenie... do zamknięcia :)

hauleth
Moderator
  • Rejestracja:ponad 17 lat
  • Ostatnio:10 dni
0

Proszę bardzo, sesja z Octave, która powinna Ci mniej więcej powiedzieć czego możesz oczekiwać:

Kopiuj
octave:13> x
x =

     1     2     5
    73     4   754
     4    54     4

octave:14> ceil(det(x))
ans = -15621
octave:15> floor(det(x))
ans = -15622
octave:16> det(x)
ans = -1.5622e+04

MarekR22
ewidentnie Octave dała ciała (no chyba, że ceil lub floor w Octave dla argumentu całkowitego nie daje tej samej liczby).
PP
  • Rejestracja:ponad 13 lat
  • Ostatnio:ponad 12 lat
  • Postów:9
0

okazuje się, że to nie kompilator szaleje tylko ja :D
dzięki winefresh, znałem poprawny wynik bo sobie go policzyłem ale na to, że int/int!=int to już nie wpadłem... :D

MarekR22
Moderator C/C++
  • Rejestracja:około 17 lat
  • Ostatnio:7 minut
0

definicja wyznacznika nie zawiera dzielenia jedynie sumę po permutacjach z iloczynów odpowiednich elementów i parzystości permutacji.
\det M=\sum_{\sigma \in S_n} (-1)^{\rm{Inv}(\sigma)}~a_{1\sigma(1)}\cdot a_{2\sigma(2)}\cdot ... \cdot a_{n\sigma(n)}
Ty stosujesz jakaś metodę, do liczenia wyznacznika. która używa dzielenia.
Gdybyś zastosował rozwinięcie Laplace'a nie miałbyś tego problemu.


offtop: LOL jakiś totalny dyslektyk edytował wiki na temat wyznacznika: "wyznacznik ogulny" albo "pżykładowo"

Jeśli chcesz pomocy, NIE pisz na priva, ale zadaj dobre pytanie na forum.
edytowany 4x, ostatnio: MarekR22

Zarejestruj się i dołącz do największej społeczności programistów w Polsce.

Otrzymaj wsparcie, dziel się wiedzą i rozwijaj swoje umiejętności z najlepszymi.