[FLASH-USERS] [EXT] Re: Segmentation fault after file modification

[Ryan Farber]

Yup, you currently calculate radius before i,j,k get assigned values(!)

! what you currently have:
! i,j,k aren't assigned any values until the do loops below!

radius = xCenter(i)**2+yCenter(j)**2+zCenter(k)**2


     ! Sum contributions from the indicated blkLimits of cells.

     do k = blkLimits(LOW,KAXIS), blkLimits(HIGH,KAXIS)

        do j = blkLimits(LOW,JAXIS), blkLimits(HIGH,JAXIS)

           do i = blkLimits(LOW,IAXIS), blkLimits(HIGH,IAXIS)

Best,
--------
Ryan


On Tue, Oct 31, 2023 at 11:43 AM Peter Manzella <peter.manzella at maine.edu>
wrote:

> I'm not sure what you mean by inside the i,j,k loops, is my radius on the
> wrong line?
>
> On Tue, Oct 31, 2023 at 10:32 AM Hansen, Eddie <ehansen at pas.rochester.edu>
> wrote:
>
>> The line calculating ‘radius’ needs to go inside the I,j,k do loops.
>> Also, keep in mind that his is radius^2 (maybe that’s what you want for
>> your modifications).
>>
>>
>>
>> --
>>
>> Eddie Hansen
>>
>> Applications Group Leader
>>
>> Flash Center for Computational Science
>>
>>
>>
>>
>>
>> *From: *Peter Manzella <peter.manzella at maine.edu>
>> *Date: *Tuesday, October 31, 2023 at 10:16 AM
>> *To: *Hansen, Eddie <ehansen at pas.rochester.edu>
>> *Cc: *flash-users at flash.rochester.edu <flash-users at flash.rochester.edu>
>> *Subject: *[EXT] Re: [FLASH-USERS] Segmentation fault after file
>> modification
>>
>> I can share that, and the debug message I get after I have added the
>> "-debug" flag as well. Here is the debug message,
>>
>>
>>
>> At line 119 of file IO_writeIntegralQuantities.F90
>> Fortran runtime error: Index '1997467855' of dimension 1 of array
>> 'xcenter' above upper bound of 24
>>
>> Error termination. Backtrace:
>> At line 119 of file IO_writeIntegralQuantities.F90
>> Fortran runtime error: Index '0' of dimension 1 of array 'xcenter' below
>> lower bound of 1
>>
>> Error termination. Backtrace:
>> Fortran runtime error: Index '0' of dimension 1 of array 'xcenter' below
>> lower bound of 1
>>
>> Error termination. Backtrace:
>> At line 119 of file IO_writeIntegralQuantities.F90
>> Fortran runtime error: Index '50696040' of dimension 1 of array 'xcenter'
>> above upper bound of 24
>>
>>
>>
>> And attached is the modified file.
>>
>>
>>
>> On Tue, Oct 31, 2023 at 10:05 AM Hansen, Eddie <ehansen at pas.rochester.edu>
>> wrote:
>>
>> Hello,
>>
>>
>>
>> It depends on how you modified it. If you share your modified routine,
>> someone might be able to find the bug.
>>
>>
>>
>> --
>>
>> Eddie Hansen
>>
>> Applications Group Leader
>>
>> Flash Center for Computational Science
>>
>>
>>
>>
>>
>> *From: *flash-users <flash-users-bounces at flash.rochester.edu> on behalf
>> of Peter Manzella <peter.manzella at maine.edu>
>> *Date: *Monday, October 30, 2023 at 3:01 PM
>> *To: *flash-users at flash.rochester.edu <flash-users at flash.rochester.edu>
>> *Subject: *[FLASH-USERS] Segmentation fault after file modification
>>
>> Hello,
>>
>>
>>
>> I recently modified IO_writeIntegralQuantities.F90, and in doing so, I've
>> gotten a segmentation fault issue. I am trying to track it down, but am
>> having some trouble, are there some ways I can go about doing this?
>> Any suggestions? Help is appreciated!
>>
>>
>>
>> --
>>
>> Peter Manzella
>>
>> PhD Candidate, Department of Physics and Astronomy
>>
>> 309 Bennett Hall
>>
>> University of Maine
>>
>> peter.manzella at maine.edu
>>
>>
>>
>> “If you run out of ideas follow the road; you'll get there.” – Edgar
>> Allan Poe
>>
>>
>>
>>
>> --
>>
>> Peter Manzella
>>
>> PhD Candidate, Department of Physics and Astronomy
>>
>> 309 Bennett Hall
>>
>> University of Maine
>>
>> peter.manzella at maine.edu
>>
>>
>>
>> “If you run out of ideas follow the road; you'll get there.” – Edgar
>> Allan Poe
>>
>
>
> --
> Peter Manzella
> PhD Candidate, Department of Physics and Astronomy
> 309 Bennett Hall
> University of Maine
> peter.manzella at maine.edu
>
> “If you run out of ideas follow the road; you'll get there.” – Edgar Allan
> Poe
> _______________________________________________
> flash-users mailing list
> flash-users at flash.rochester.edu
>
> For list info, including unsubscribe:
> https://flash.rochester.edu/mailman/listinfo/flash-users
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20231101/980d3f07/attachment.htm>