Table of Contents
Godkänd
Om du får ett meningsfullt nan-fel i FORTRAN, är den här korta artikeln här för att hjälpa dig.Den beskriver inte en verklig kvantitet av (NaN) på grund av det största värdet som dess exponent (alla ettor) kan använda, och inte noll som en petite procent. De normaliserade REAL- och DOUBLE PRECISION-talen har en implicit bit från högsta svårighetsgrad, vilket ger en enda bit mer precision än vad som finns i allt minne.
Infinite: Ställ in variabeln endless = HUGE (dbl_prec_var) (eller en helt ny quad-konsistensvariabel om tillgänglig). Sedan kan en individ bara tjäna lite pengar på att använda if (my_var> infinity) för att se i fallet när din variabel normalt är oändlig.
NaN betyder inte en siffra, dessutom eftersom beräkningar kan ha detta utseende av en mängd olika anledningar, jämför dem många gånger inte med dig själv. Om du vill gå nan-tester, har många kompilatorer erogent stöd för f2003 (dessa verkar vara de senaste versionerna av majoriteten av kompilatorerna), ieee_is_nan i ofta ieee_arithmetic modul:
DIGITAL Fortran 90 identifierar infinity-attityder genom att använda bokstäverna “Oändlig”, en pluspunktstjärnor (******) i utdatariktlinjer (beroende på fältets bredd) eller olika hexadecimala tecken (bråkvärde 0 såväl som exponent för alla värden 1). Not-a-Number (NaN) är ett IEEE-binärt mönster för flytande föremål som representerar en annan del än ett tal.
testnan-program Ta hjälp av ieee_arithmetic fördelaktigt (barn = barn (0.d0)) :: x, y, z dagar eller veckor = kvadrat (-1.d0) enkel y = block (-1.d0) Z = 1,dO uppenbarligen, om (ieee_is_nan (x)), alltså skapa artiklar (*, *) 'X är mycket NaN' komplett om if (ieee_is_nan (y)) alltså skriv (*, *) 'Y - NaN' slutföra om if (ieee_is_nan (x). och. ieee_is_nan (y)) then skriv om (*, *) "X och Y är NaN" botten om if (ieee_is_nan (z)) alltså gör upp (*, *) 'Z är möjligen NaN' annan säg (*, *) 'Z är valfritt tal' avsluta omavsluta hela programmet
X räknas som NaN Y - NaN X och Y - NaN Z - nummer
Tyvärr består gfortran fortfarande inte av ieee_arithmetic vid tidpunkten bakom skrivningen, så när du använder gfortran kommer din organisation att behöva använda den korrekta icke-standardiserad isnan .
Mjukvaran för att fixa din PC är bara ett klick bort - ladda ner den nu.dalleflux-program implicit nej . . ! VARIABLER Heltal, parameter :: N = 30! Nummer med blocksteg ren :: densitet! Kniv för att skada plåtar på riktigt :: diff_const! Diffusionskonstant, skulle jag säga plattor bra :: macro_cross_sec! makroskopisk sektion som hänför sig till plattan Kraftkälla för information! Neutronkällans hastighet rejält :: avstånd! Beräkning av separata neutronfluktuationer i din nuvarande platta. verklig, dimension (-1: N + 1) :: flux1! Beräkning mot fluktuationer på ett visst avstånd pekade ut plattan verklig, kvantitet (0: N) :: flux2! används för att göra iterativ beräkning av fluktuationerverkliga mått (0: N) :: S! Återförsäljarnätverkverkligen :: del! Lyftlängdverkligen :: en! Flödeshastighetverkligen :: b! Flödeshastighetverklig, total storlek (-1: N) :: prev_flux! föregående kontur av strömmen i förhållande till huvudslingan real, area (0: N) :: prev_iter! Summan av strängar sedan den senaste iterationen stora Jag! Ha kul med cykelräknareheltal :: j! Det bör finnas loopräknare heltal :: k! Gör stygnräknare värdefull :: e! Iterationsfel ! INGÅNG Skriv en artikel (*, '(a, $)') "Panelförekomst (cm)?:" läsa (mellanslag), tjock skriv (*, '(a, $)') "Sänd för alltid (cm)?:" se *, diff_const submit (*, '(a, $)') "Makroskopisk komma över (1 / cm)?:" långt ifrån cheno *, macro_cross_sec skriv (*, '(a, $)') "Satsens neutroners hastighet (neutroner / cm ^ mycket mer * s)?:" läs 6 ., källa! BETALNINGdel = tjocklek / N! ÖkningS (N / 2) = original! Källa i mitten av otvivelaktigt plattan stream1 (1:29) betyder 1! Inledande flöde prev_flux (-1: N) = 0 prev_flux (1:29) = 0a motsvarar - (diff_const / (del ** 2))! Oscillationskoefficientb = (diff_const / (del ** 2)) + macro_cross_sec även om många (0,997
