[Gate-users] Maximum Particle Number

[David Sarrut]

Hello,

yes currently limited to 32bits number, somewhere deep in the code ... You
probably need to split your simulation into several jobs.

David
PS: if someone wants to dig and investigate where is the limitation in the
code, he/she is welcome!

On Wed, Jan 29, 2020 at 3:07 PM Svetlana Shasharina <sveta at txcorp.com>
wrote:

> Dear Tom,
>
> I think you will need independent simulations with different random seeds.
>
> Best,
>
> Svetlana
>
> On 1/29/20 7:05 AM, Russ, Tom wrote:
> > Dear Dimitris,
> >
> > thank you very much for your help!
> > Unfortunately, I think this might not be possible in my case.
> > I want to simulate CT projections and use a general particle source
> (GPS) where I directly have to give it the number of primaries.
> > Besides I also think I cannot substitute the GPS with an activity source
> because I am using Gatelab, which (as I understood it) necessitates the use
> of the command 'setTotalNumberOfPrimaries'.
> >
> > Best regards,
> > Tom
> >
> > -----Ursprüngliche Nachricht-----
> > Von: Dimitrios Thanasas <dthanas at uoa.gr>
> > Gesendet: Mittwoch, 29. Januar 2020 14:21
> > An: Russ, Tom <Tom.Russ at medma.uni-heidelberg.de>
> > Cc: gate-users at lists.opengatecollaboration.org
> > Betreff: Re: [Gate-users] Maximum Particle Number
> >
> > Dear Tom,
> >
> > I had this problem too. I solved it by combining time and source
> activity.
> > If you want 3600000000 particles you can use 100 ΜΒq activity for 36
> seconds (or any other combination).
> > Comment the line (you don't need it):
> > #/gate/application/setTotalNumberOfPrimaries 3600000000 to run the
> simulation.
> >
> > Best Regards,
> > Dimitris
> >
> >> Dear Gate-Community,
> >>
> >> Gate seems to have a problem with primary particle numbers over
> >> exactly
> >> 2.147.483.648 (2^(32)/2) and simply throws an error ('Exiting with
> >> return value 6') and stops simulating for me.
> >> Can this be due to Gate having problems with numbers outside the
> >> 32bit-variable range?
> >> Did anybody else encounter this?
> >> This is how my code looks:
> >>
> >> '/gate/application/setTotalNumberOfPrimaries 3600000000
> >> /gate/application/setTimeSlice 1 s /gate/application/setTimeStart 0 s
> >> /gate/application/setTimeStop 90 s'
> >>
> >> I use this acquisition timing due to the rotation of a voxelized
> phantom.
> >>
> >> Any help will be much appreciated! :)
> >>
> >> Best regards,
> >> Tom
> >>
> >>
> >> _______________________________________________
> >> Gate-users mailing list
> >> Gate-users at lists.opengatecollaboration.org
> >> http://lists.opengatecollaboration.org/mailman/listinfo/gate-users
> >
> > _______________________________________________
> > Gate-users mailing list
> > Gate-users at lists.opengatecollaboration.org
> > http://lists.opengatecollaboration.org/mailman/listinfo/gate-users
> >
> _______________________________________________
> Gate-users mailing list
> Gate-users at lists.opengatecollaboration.org
> http://lists.opengatecollaboration.org/mailman/listinfo/gate-users



-- 
David Sarrut, Phd
Directeur de recherche CNRS
CREATIS, UMR CNRS 5220, Inserm U1206
Centre de lutte contre le cancer Léon Bérard
28 rue Laënnec, 69373 Lyon cedex 08
Tel : 04 78 78 51 51 / 06 74 72 05 42
http://www.creatis.insa-lyon.fr/~dsarrut
_________________________________
 "2 + 2 = 5,  for extremely large values of 2"
_________________________________
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20200129/ae97d8f3/attachment-0001.html>