Re: Index '100001' of dimension 1 of array 'tx' above upper bound of 100000

From: Sam Tygier <sam.tygier_at_hep.manchester.ac.uk>
Date: Fri, 08 Feb 2013 11:20:01 +0000

2013 12:20:09 +0100 (CET)
X-Virus-Scanned: clamav-milter 0.97.6 at smtp1.mi.infn.it
X-Virus-Status: Clean
Sender: owner-fluka-discuss_at_smtp2.mi.infn.it

Hi Stefan,

Do you mean I need to find a simpler expression of the SUR region into
my inputfile, or that I should be using the WHAT(5) option in GEOBEGIN?
or is it the GEOBEGIN optimisation that has expanded my expression from
31 to 100k terms?

Thanks,
Sam

On 06/02/13 15:46, Stefan Roesler wrote:
> Dear Sam
>
> I am afraid that you have to optimize your parenthesis. The code
> arrived at 100k terms in the expansion that would slow down a lot the
> calculation if the array dimension was increased.
>
> There is a protection missing in the code that may need to be fixed.
> But for your case it would only give you a 'stop' instead of an
> overflow...
>
> Cheers
> Stefan
>
>
> On Tue, 5 Feb 2013, Sam Tygier wrote:
>
>> Hi,
>>
>> I am hitting the error:
>> At line 84 of file rpnorm.f
>> Fortran runtime error: Index '100001' of dimension 1 of array 'tx' above
>> upper bound of 100000
>> With a fairly simple geometry in FLUKA2011 Version 2.17 Dec-12. I have
>> striped out most of my geometry to just the problematic elements.
>>
>> I have 7 rings/coils, defined using 14 planes, and 8 cylinders (ZCC). A
>> vacuum region surrounds them. This surrounding region needs to have the
>> 7 rings cut out of it, so it has a couple of lines in its definition.
>> This still seems to me like a simple region, and the 100001 vs 100000
>> would make me suspicious of an out-by-one error.
>>
>> The 'Handling of errors and crashes' tutorial suggests that infinite
>> bodies are are problem, but the fluka documentation recommends infinite
>> bodies over finite ones.
>>
>> What should I do? Thanks,
>>
>> Sam
>>
>> ===============
>>
>> TITLE
>> piontarget
>> FREE
>> GLOBAL, 10000, , , 3 , 1 ,
>> DEFAULTS , , , , , , ,NEW-DEFA
>> BEAM, -8, 0, 0, 0, 0, 0, PROTON
>> BEAMPOS, -15.643446504, 0, 98.7688340595, 0.15643446504, 0, , NEGATIVE
>> FREE
>> GEOBEGIN , , , , , , , COMBNAME
>> 0, 0, pt
>> XYP r0pzs -95.0
>> XYP r0pze -77.5
>> XYP r1pzs -72.5
>> XYP r1pze 62.5
>> XYP r2pzs 67.5
>> XYP r2pze 247.5
>> XYP r3pzs 254.5
>> XYP r3pze 292.5
>> XYP r4pzs 284.5
>> XYP r4pze 305.5
>> XYP r5pzs 318.0
>> XYP r5pze 357.0
>> XYP r6pzs 367.0
>> XYP r6pze 398.0
>> ZCC zcc74_2 0 0 74.2
>> ZCC zcc66_2 0 0 66.2
>> ZCC zcc29_8 0 0 29.8
>> ZCC zcc26_6 0 0 26.6
>> ZCC zcc25_0 0 0 25.0
>> ZCC zcc80_6 0 0 80.6
>> ZCC zcc28_2 0 0 28.2
>> ZCC zcc79_0 0 0 79.0
>> ZCC sur 0 0 100
>> XYP surps -105
>> XYP surpe 408.0
>> SPH bbs 0.0 0.0 0.0 1000
>> END
>> R0 5 -r0pzs +r0pze -zcc66_2 +zcc80_6
>> R1 5 -r1pzs +r1pze -zcc66_2 +zcc80_6
>> R2 5 -r2pzs +r2pze -zcc66_2 +zcc74_2
>> R3 5 -r3pzs +r3pze -zcc66_2 +zcc79_0
>> R4 5 -r4pzs +r4pze -zcc25_0 +zcc26_6
>> R5 5 -r5pzs +r5pze -zcc25_0 +zcc28_2
>> R6 5 -r6pzs +r6pze -zcc25_0 +zcc29_8
>> SUR 5 - surps + surpe + sur -( -r0pzs +r0pze -zcc66_2 +zcc80_6) -(
>> -r1pzs +r1pze -zcc66_2 +zcc80_6) -( -r2pzs +r2pze -zcc66_2 +zcc74_2) -(
>> -r3pzs +r3pze -zcc66_2 +zcc79_0) -( -r4pzs +r4pze -zcc25_0 +zcc26_6) -(
>> -r5pzs +r5pze -zcc25_0 +zcc28_2) -( -r6pzs +r6pze -zcc25_0 +zcc29_8)
>> BB 5 +bbs - ( - surps + surpe + sur)
>> END
>> GEOEND
>> ASSIGNMAT, COPPER, R0, , ,1
>> ASSIGNMAT, COPPER, R1, , ,1
>> ASSIGNMAT, COPPER, R2, , ,1
>> ASSIGNMAT, COPPER, R3, , ,1
>> ASSIGNMAT, COPPER, R4, , ,1
>> ASSIGNMAT, COPPER, R5, , ,1
>> ASSIGNMAT, COPPER, R6, , ,1
>> ASSIGNMAT, VACUUM, SUR, , ,1
>> ASSIGNMAT, BLCKHOLE, BB, , ,0
>> SCORE, 208
>> RANDOMIZE, 1, 0
>> START , 10.0
>> STOP
>>
>>
>>
Received on Mon Feb 11 2013 - 10:02:49 CET