[Gate-users] Can't run with 100 million particles
S. Saghamanesh
s.saghamanesh at yahoo.com
Wed Oct 15 15:09:15 CEST 2014
I read this in a Monte Carlo reference text: "It (the Jensen's inequality in Monte Carlo method) expresses that the end result of a measurement process that combines elementary measurement processes is different from the result of a measurement process of combined elementary measurement process results." That's why I am sort of doubtful about the ultimate uncertainty of this addition approach.
BestSoman
On Wednesday, October 15, 2014 2:10 PM, Grevillot Loic <loic.grevillot at medaustron.at> wrote:
#yiv7102759071 #yiv7102759071 -- _filtered #yiv7102759071 {font-family:Helvetica;panose-1:2 11 6 4 2 2 2 2 2 4;} _filtered #yiv7102759071 {font-family:Helvetica;panose-1:2 11 6 4 2 2 2 2 2 4;} _filtered #yiv7102759071 {font-family:Calibri;panose-1:2 15 5 2 2 2 4 3 2 4;} _filtered #yiv7102759071 {font-family:Tahoma;panose-1:2 11 6 4 3 5 4 4 2 4;} _filtered #yiv7102759071 {font-family:Verdana;panose-1:2 11 6 4 3 5 4 4 2 4;}#yiv7102759071 #yiv7102759071 p.yiv7102759071MsoNormal, #yiv7102759071 li.yiv7102759071MsoNormal, #yiv7102759071 div.yiv7102759071MsoNormal {margin:0cm;margin-bottom:.0001pt;font-size:12.0pt;}#yiv7102759071 a:link, #yiv7102759071 span.yiv7102759071MsoHyperlink {color:blue;text-decoration:underline;}#yiv7102759071 a:visited, #yiv7102759071 span.yiv7102759071MsoHyperlinkFollowed {color:purple;text-decoration:underline;}#yiv7102759071 p.yiv7102759071MsoAcetate, #yiv7102759071 li.yiv7102759071MsoAcetate, #yiv7102759071 div.yiv7102759071MsoAcetate {margin:0cm;margin-bottom:.0001pt;font-size:8.0pt;}#yiv7102759071 p.yiv7102759071msonormal, #yiv7102759071 li.yiv7102759071msonormal, #yiv7102759071 div.yiv7102759071msonormal {margin-right:0cm;margin-left:0cm;font-size:12.0pt;}#yiv7102759071 p.yiv7102759071msoacetate, #yiv7102759071 li.yiv7102759071msoacetate, #yiv7102759071 div.yiv7102759071msoacetate {margin-right:0cm;margin-left:0cm;font-size:12.0pt;}#yiv7102759071 p.yiv7102759071msochpdefault, #yiv7102759071 li.yiv7102759071msochpdefault, #yiv7102759071 div.yiv7102759071msochpdefault {margin-right:0cm;margin-left:0cm;font-size:12.0pt;}#yiv7102759071 p.yiv7102759071msonormal1, #yiv7102759071 li.yiv7102759071msonormal1, #yiv7102759071 div.yiv7102759071msonormal1 {margin-right:0cm;margin-left:0cm;font-size:12.0pt;}#yiv7102759071 p.yiv7102759071msochpdefault1, #yiv7102759071 li.yiv7102759071msochpdefault1, #yiv7102759071 div.yiv7102759071msochpdefault1 {margin-right:0cm;margin-left:0cm;font-size:12.0pt;}#yiv7102759071 p.yiv7102759071msonormal2, #yiv7102759071 li.yiv7102759071msonormal2, #yiv7102759071 div.yiv7102759071msonormal2 {margin-right:0cm;margin-left:0cm;font-size:12.0pt;}#yiv7102759071 p.yiv7102759071msoacetate1, #yiv7102759071 li.yiv7102759071msoacetate1, #yiv7102759071 div.yiv7102759071msoacetate1 {margin-right:0cm;margin-left:0cm;font-size:12.0pt;}#yiv7102759071 p.yiv7102759071msonormal3, #yiv7102759071 li.yiv7102759071msonormal3, #yiv7102759071 div.yiv7102759071msonormal3 {margin-right:0cm;margin-left:0cm;font-size:12.0pt;}#yiv7102759071 p.yiv7102759071msochpdefault2, #yiv7102759071 li.yiv7102759071msochpdefault2, #yiv7102759071 div.yiv7102759071msochpdefault2 {margin-right:0cm;margin-left:0cm;font-size:12.0pt;}#yiv7102759071 p.yiv7102759071msonormal11, #yiv7102759071 li.yiv7102759071msonormal11, #yiv7102759071 div.yiv7102759071msonormal11 {margin-right:0cm;margin-left:0cm;font-size:12.0pt;}#yiv7102759071 p.yiv7102759071msochpdefault11, #yiv7102759071 li.yiv7102759071msochpdefault11, #yiv7102759071 div.yiv7102759071msochpdefault11 {margin-right:0cm;margin-left:0cm;font-size:12.0pt;}#yiv7102759071 span.yiv7102759071msohyperlink {}#yiv7102759071 span.yiv7102759071msohyperlinkfollowed {}#yiv7102759071 span.yiv7102759071msohyperlink2 {}#yiv7102759071 span.yiv7102759071msohyperlinkfollowed2 {}#yiv7102759071 span.yiv7102759071msohyperlink11 {}#yiv7102759071 span.yiv7102759071msohyperlinkfollowed11 {}#yiv7102759071 span.yiv7102759071emailstyle1711 {}#yiv7102759071 span.yiv7102759071emailstyle271 {}#yiv7102759071 span.yiv7102759071emailstyle40 {}#yiv7102759071 p.yiv7102759071msonormal4, #yiv7102759071 li.yiv7102759071msonormal4, #yiv7102759071 div.yiv7102759071msonormal4 {margin:0cm;margin-bottom:.0001pt;font-size:12.0pt;}#yiv7102759071 span.yiv7102759071msohyperlink1 {color:blue;text-decoration:underline;}#yiv7102759071 span.yiv7102759071msohyperlinkfollowed1 {color:purple;text-decoration:underline;}#yiv7102759071 p.yiv7102759071msoacetate2, #yiv7102759071 li.yiv7102759071msoacetate2, #yiv7102759071 div.yiv7102759071msoacetate2 {margin:0cm;margin-bottom:.0001pt;font-size:12.0pt;}#yiv7102759071 p.yiv7102759071msochpdefault3, #yiv7102759071 li.yiv7102759071msochpdefault3, #yiv7102759071 div.yiv7102759071msochpdefault3 {margin-right:0cm;margin-left:0cm;font-size:12.0pt;}#yiv7102759071 p.yiv7102759071msonormal12, #yiv7102759071 li.yiv7102759071msonormal12, #yiv7102759071 div.yiv7102759071msonormal12 {margin-right:0cm;margin-left:0cm;font-size:12.0pt;}#yiv7102759071 p.yiv7102759071msochpdefault12, #yiv7102759071 li.yiv7102759071msochpdefault12, #yiv7102759071 div.yiv7102759071msochpdefault12 {margin-right:0cm;margin-left:0cm;font-size:12.0pt;}#yiv7102759071 p.yiv7102759071msonormal21, #yiv7102759071 li.yiv7102759071msonormal21, #yiv7102759071 div.yiv7102759071msonormal21 {margin:0cm;margin-bottom:.0001pt;font-size:12.0pt;}#yiv7102759071 span.yiv7102759071msohyperlink21 {color:blue;text-decoration:underline;}#yiv7102759071 span.yiv7102759071msohyperlinkfollowed21 {color:purple;text-decoration:underline;}#yiv7102759071 p.yiv7102759071msoacetate11, #yiv7102759071 li.yiv7102759071msoacetate11, #yiv7102759071 div.yiv7102759071msoacetate11 {margin:0cm;margin-bottom:.0001pt;font-size:8.0pt;}#yiv7102759071 p.yiv7102759071msonormal31, #yiv7102759071 li.yiv7102759071msonormal31, #yiv7102759071 div.yiv7102759071msonormal31 {margin-right:0cm;margin-left:0cm;font-size:12.0pt;}#yiv7102759071 p.yiv7102759071msochpdefault21, #yiv7102759071 li.yiv7102759071msochpdefault21, #yiv7102759071 div.yiv7102759071msochpdefault21 {margin-right:0cm;margin-left:0cm;font-size:12.0pt;}#yiv7102759071 p.yiv7102759071msonormal111, #yiv7102759071 li.yiv7102759071msonormal111, #yiv7102759071 div.yiv7102759071msonormal111 {margin:0cm;margin-bottom:.0001pt;font-size:12.0pt;}#yiv7102759071 span.yiv7102759071msohyperlink111 {color:blue;text-decoration:underline;}#yiv7102759071 span.yiv7102759071msohyperlinkfollowed111 {color:purple;text-decoration:underline;}#yiv7102759071 span.yiv7102759071emailstyle17111 {color:#1F497D;}#yiv7102759071 p.yiv7102759071msochpdefault111, #yiv7102759071 li.yiv7102759071msochpdefault111, #yiv7102759071 div.yiv7102759071msochpdefault111 {margin-right:0cm;margin-left:0cm;font-size:10.0pt;}#yiv7102759071 span.yiv7102759071emailstyle2711 {color:#1F497D;}#yiv7102759071 span.yiv7102759071emailstyle401 {color:#1F497D;}#yiv7102759071 span.yiv7102759071BalloonTextChar {}#yiv7102759071 span.yiv7102759071EmailStyle59 {color:#1F497D;}#yiv7102759071 .yiv7102759071MsoChpDefault {font-size:10.0pt;} _filtered #yiv7102759071 {margin:70.85pt 70.85pt 2.0cm 70.85pt;}#yiv7102759071 div.yiv7102759071WordSection1 {}#yiv7102759071 Of course not! If you buy 2 times 10 apples or 1 time 20 apples you get in both cases 20 apples. Reference to papers used for uncertainty calculations are available – please look at Gate reference papers! Please use the list, other people may benefit of discussion and my contributions are on free time which is very limited at the moment… Kind regards, Loïc Grevillot From: S. Saghamanesh [mailto:s.saghamanesh at yahoo.com]
Sent: Mittwoch, 15. Oktober 2014 12:34
To: Grevillot Loic
Subject: Re: [Gate-users] Can't run with 100 million particles Thanks for your note. However, due to the linear operation of this addition, I'll have a greater uncertainty! For example, I have two simulations with uncertainties of sigma1 = 0.7 (70%) and sigma2 = 0.6 (60%). After the addition, I'll have the total uncertainty of sqrt(sigma1^2 + sigma2^2) = 0.92 (92%). So I wanted to reduce my uncertainty by increasing my statistics (through several simulations), but ultimately I'd obtain a greater uncertainty. Am I right on this? Best Soman On Monday, October 13, 2014 8:41 AM, Grevillot Loic <loic.grevillot at medaustron.at> wrote: You can manually run 10 different simulations. Make sure to change output names for each simu and then make a simple addition. Sincerely, Loïc From: S. Saghamanesh [mailto:s.saghamanesh at yahoo.com]
Sent: Samstag, 11. Oktober 2014 13:53
To: Grevillot Loic; gate-users at lists.opengatecollaboration.org
Subject: Re: [Gate-users] Can't run with 100 million particles Dear Loïc, I've tried to run my main macro file through the following command: "/control/loop myMain.mac mycounter 1 10 10" which is supposed to run 10 simulations with 10 million particles with different seeds (by /gate/random/setEngineSeed auto). But It runs only the first 10 million particles and doesn't keep run the remains. Do you have any idea what's wrong with this command? I guess it's the only way I can proceed with 100 million primaries without getting error. Best regards Soman On Friday, October 3, 2014 9:10 AM, Grevillot Loic <loic.grevillot at medaustron.at> wrote: Dear Soman, You said you need to run 100 million primaries and artifacts appear after 10 million. My suggestion is therefore to run 10 simulations of 10 million particles (fully equivalent to a single simu of 100 million). You just have to make sure the seeds are different for every simulation, which can be easily achieved using the following command: /gate/random/setEngineSeed auto Best, Loïc From: S. Saghamanesh [mailto:s.saghamanesh at yahoo.com]
Sent: Freitag, 3. Oktober 2014 07:30
To: Grevillot Loic
Subject: Re: [Gate-users] Can't run with 100 million particles Hi, Thanks a lot for your note. I use of splitting 1000 for my emitted photons but I still have unacceptable uncertainties, so I have to run with 100 million particles. Do you tell me if I run 10 simulations with 10 million particles, it would lessen the uncertainty? I guess it wouldn't, because I have uncertain results with a single simulation of 10 million particles. Thanks in advance for your help - Soman On Wednesday, October 1, 2014 3:41 PM, Grevillot Loic <loic.grevillot at medaustron.at> wrote: Dear Soman, We noticed a similar simulation “artifact” when a too high number of particles was used. We also noticed that the effect is dependent on the particle type and simulation parameters (cut, step size). We will try correcting the bug for the next release. In the meantime I advise you to run 10 simulations of 10 million particles. Sincerely, Loïc From:gate-users-bounces at lists.opengatecollaboration.org [mailto:gate-users-bounces at lists.opengatecollaboration.org]On Behalf Of S. Saghamanesh
Sent: Mittwoch, 1. Oktober 2014 09:29
To: gate-users at lists.opengatecollaboration.org
Subject: [Gate-users] Can't run with 100 million particles Hi, Is there any limitation on the number of primary particles generated in Gate7.0 ? It seems running code with more than 10 million particles leads to false results. I need to run my code with 100 millions particles, but after checking (with a very simple example), I noticed that no more than 20 million particles would be generated. Thanks for your comments. - Soman Disclaimer:
Please notice our E-Mail Disclaimer http://www.medaustron.at/email-disclaimer/ Disclaimer:
Please notice our E-Mail Disclaimer http://www.medaustron.at/email-disclaimer/ Disclaimer:
Please notice our E-Mail Disclaimer http://www.medaustron.at/email-disclaimer/ Disclaimer:
Please notice our E-Mail Disclaimer http://www.medaustron.at/email-disclaimer/
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/mailman/private/gate-users/attachments/20141015/20e4b65c/attachment-0001.html>
More information about the Gate-users
mailing list