precisione di calcolo

[tony]

recentemente (problema "sistema solare in miniatura" della gara) ho cozzato contro i limiti dei numeri floating point da 80 bit ("extended precision" in c, pascal e alcuni basic).

per superarla ho un vecchio, lento, "UBASIC" un po' fuori dagli standard del linguaggio e un moderno, lento, "derive" non semplicissimo da programmare.

mi accontentavo, per quel problema, di una precisione di 24 cifre.

qualcuno dei lettori usa altri strumenti?

grazie. tony

[Nidhogg]

Un interessante articolo sul floating point

http://www.nvitalia.com/articoli/Floating_Point/1.htm

In questo articolo, come io pensavo, è considerato il FORTRAN il linguaggio più adatto per applicazioni matematiche.

http://www.vialattea.net/esperti/inform/precisione.htm

Ciao, Ermanno.

[GIOVANNI IL CHIMICO]

Ciao Tony, confermo che il fortran è ottimo per gli algoritmi di calcolo,ma anche Matlab fail suolavoro, con il vantaggio che è semplicissimo da programmare e contiene unalibreria di migliaia di subroutine di tutti i generi, oltre che una capacità grafica spettacolare...

[tony]

grazie, ragazzi,

(leonardo per i due link che mostrano le moderne tendenze dei calcolatori per scopi grafici, e Chimico che mi conferma che il fortran è ottimo per gli algoritmi di calcolo)

ma il problema che avevo posto è diverso:
<b>grosse</b> precisioni (avevo detto che mi servivano "solo" 24 cifre)

provate a fare in fortran qualche operazioncella aritmetica (non funzioni, per carità, non arrivo a pretendere tanto!) da 24 cifre,
per es.
X = 10^11 + 10^-12 e poi Y = X - 10^11
che vi aspettate? Y = 10^-12? campa, cavallo ...!


ecco, di questo parlavo.
(volete chiamarlo "unlimited precision arithmetics"?

c'e, anche qui su matematicamente, della teoria; ma a me interessano strumenti, e chiedevo quali strumenti (oltre quelli che ho citato: il ben noto "derive" e il forse ignoto "ubasic") usino i lettori per ottenere la precisione richiesta.
il Matlab che cita il Chimico che cosa riesce a fare?

tony

[tony]

nulla mi si dice sulla precisione ottenibile con questo matlab?

grazie
tony