From ebc10e2ef46f0fd21fbebd036179bc7adf3b337a Mon Sep 17 00:00:00 2001 From: Erich Eckner Date: Wed, 14 Oct 2015 16:17:57 +0200 Subject: dumpErhaltungsgroeszen vertrÃĪgt jetzt auch, wenn keine Teilchen da sind MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- Physikunit.pas | 20 +++++++++++--------- 1 file changed, 11 insertions(+), 9 deletions(-) diff --git a/Physikunit.pas b/Physikunit.pas index 254c8d4..d546621 100644 --- a/Physikunit.pas +++ b/Physikunit.pas @@ -939,16 +939,18 @@ var i,j: longint; dens: double; begin - dens:=0; - for i:=0 to length(massen)-1 do - dens:=dens+massen[i]/teilchen[i].eigenschaften[tsgMasse]; - gitter.prot.schreibe('Gesamtdefizit: '+floattostr(gesamtDefizit)+' (Anteil '+floattostr(gesamtDefizit/dens*dX*dP)+')',true); - for i:=0 to length(massen)-1 do begin + if length(massen)>0 then begin dens:=0; - for j:=0 to aX-1 do - dens:=dens+impulsIntegral(j,i,msN); - dens:=dens*teilchen[i].eigenschaften[tsgMasse]*gitter.dX; - gitter.prot.schreibe('n['+inttostr(i+1)+'] = '+floattostr(dens)+' (relative Abweichung: '+floattostr(dens/massen[i]-1)+')',true); + for i:=0 to length(massen)-1 do + dens:=dens+massen[i]/teilchen[i].eigenschaften[tsgMasse]; + gitter.prot.schreibe('Gesamtdefizit: '+floattostr(gesamtDefizit)+' (Anteil '+floattostr(gesamtDefizit/dens*dX*dP)+')',true); + for i:=0 to length(massen)-1 do begin + dens:=0; + for j:=0 to aX-1 do + dens:=dens+impulsIntegral(j,i,msN); + dens:=dens*teilchen[i].eigenschaften[tsgMasse]*gitter.dX; + gitter.prot.schreibe('n['+inttostr(i+1)+'] = '+floattostr(dens)+' (relative Abweichung: '+floattostr(dens/massen[i]-1)+')',true); + end; end; end; -- cgit v1.2.3-54-g00ecf