[WetProg] 2e zit

[burnie]

even voor alle duidelijkheid; morgen (maandag dus) om 8.30
u kunne we anders best is ff de bekomen resultaten vergelijken!

iedereen die zig geroepen voelt, 8.30u @ Campus Middelheim!

P.S. als ge't nog ni gemaakt hebt, ge moet ni vr de opl'n komen!

[Lesly]

ok... ik kom... 8u30 (ergens in het G-gebouw)

[Lesly]

voor vraag 1 heb ik nu:

Code: Select all

    Het conditie getal van A is 524.06    (berekend m.b.v. Octave)

    

    Oplossing voor i = 1:   1.59   b(1) - Ax~:   0

                           -1.14                 0.937

                            5.94                 1.811

    

    Oplossing voor i = 2:  -0.9   b(2) - Ax~:    0

                            9.6                  1.766666667

                           -3                   -7.666666667

    

    Oplossing voor i = 3:   5.1    b(3) - Ax~:   4.440892099e-16

                          -18                   -0.2666666667

                           21                   16.96666667

    

    Oplossing voor i = 4:  -3     b(4) - Ax~:    8.881784197e-16

                           24                    2.6

                          -18                  -21.2

    

    Oplossing voor i = 5:   1      b(5) - Ax~:   0

                            2                    1.15

                            3                   -1.05

[Shinta]

Heeft er iemand toevallig een cursus van wetenschappelijk programmeren online staan ?

[Lesly]

staat op de website van cant

[Dave]

Mijn resultaten op oef 1 zijn ongeveer indentiek

[Lesly]

mmm... raar... ofwel zijn ze identiek ofwel ni

ik ben trouwens bijna klaar... nu nog beetje uitleg in verslag schrijve en klaar

dus als er iemand nog wil vergelijken: post output hier (GEEN code... dat gaan ze verdacht en/of niet leuk vinden) en dan zulle we wel is kijke

[Dave]

Verschillende precisies/architecturen kunnen andere resultaten opleveren he
Alsik double precision naboots op mijn intel cpu, krijgik ook licht andere resultaten dan wanneer ik native precision gebruikt. Meer bepaald resultaten als 4.440892099e-16 worden 0, doordat extended precision wel toerijkend is om hier een exact resultaat te geven blijkbaar.

[Lesly]

blijkbaar nog een foutje bij mijn controle berekening:


hier is mijn nieuwe output

Code: Select all

/* Output:

    

    Het conditie getal van A is 524.06    (berekend m.b.v. Octave)

    

    Oplossing voor i = 1:   1.59   b(1) - Ax~:   0

                           -1.14                 0

                            5.94                 0

    

    Oplossing voor i = 2:  -0.9    b(2) - Ax~:   0

                            9.6                  2.220446049e-16

                           -3                    0

    

    Oplossing voor i = 3:   5.1    b(3) - Ax~:   4.440892099e-16

                          -18                    2.220446049e-16

                           21                   -4.440892099e-16

    

    Oplossing voor i = 4:  -3      b(4) - Ax~:   8.881784197e-16

                           24                    0

                          -18                    4.440892099e-16

    

    Oplossing voor i = 5:   1      b(5) - Ax~:   0

                            2                    0

                            3                    2.220446049e-16

    

    

*/

architectuur = intel core duo... double precision

[Dave]

Euh, wat hebt ge nu gedaan? want ik had dezelfde resultaten als voor dat ge der dat foutje had uitgehaald.
intel core duo is standaard extended precision he, ik mimic double precision met GSL_IEEE_MODE="double-precision,round-to-nearest".

[Lesly]

in de linker kolom staan mijn oplossingen voor het stelsel... die zijn onveranderd gebleven... rechts is de controle die ik uitvoer en daar was eerst een fout.

Je zou links zeker hetzelfde moeten hebben...



ik maakte volgende fout:

1. ik definieerde A

2. ik splitste A op in LU (maar dat noemt ook A)

3. ik rekende met A (eigenlijk LU nu) de controle uit


oplossing: zorg ervoor dat je een back-up hebt van A

[Dave]

Ah, daar zegt ge zoiets, die fout heb ik dan eerder ook al gemaakt. Thx!

[Lesly]

Dave, je zit nu toch al verder dan de 1ste opgave he?

anyway... als je problemen hebt, output wilt vergelijken ofzo, en er reageert ni direct iemand op dit forum stuur dan maar ne email naar mij... ik heb alles af (morge afgeve... )