Message: Re: Infinite loop with G4UniformElectricField Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Feedback Re: Infinite loop with G4UniformElectricField 

Forum: Fields: Magnetic and Otherwise
Re: Question Infinite loop with G4UniformElectricField (Matthew Dimmock)
Date: 24 Oct, 2013
From: Gumplinger Peter <Gumplinger Peter>

Hi Matthew,

You probably know that (faulty) programs can behave very differently, simply by adding an innocent line of code to the source. This causes the executable to be stored differently in memory and, bingo, if your program overwrites memory - completely unrelated to that extra code - it may now overwrite a part of the memory that suddenly crashes your progam or causes it to behave very differently. I believe, your issue is exactly that and has nothing to do with your field implementation.

You should ask, why:

> Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
>     0   -0.952   -0.618     3.74     0.174        0        0         0 ion_cham_layer_param_phys initStep
>     1   -0.756    -1.95        6     0.174 0.000382     2.64      2.64 ion_cham__sens_phys Transportation
>     2   -0.756    -1.95        6     0.174        0 -1.13e+25 -1.13e+25 ion_cham__sens_phys eIoni

is the step limiting process eIoni, yet the dE is 0, but the step length is completely unreasonable.

These are the worst kind of problems for a software developer. What happens if you change some of your em directives? What happens if you simplify your ElectricFieldSetup constructor? Only you can debug your progam.

Peter

 Add Message Add Message
to: "Re: Infinite loop with G4UniformElectricField"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]