Re: Gadget crashes at IC readin

From: Volker Springel <volker_at_MPA-Garching.MPG.DE>
Date: Thu, 26 Nov 2015 11:17:19 +0100

Hi Sylvia,

Presumably you have an IC file with 64-bit IDs but use gadget with 32-bit IDs. This could be resolved by activating LONGIDS in the code. Alternatively you recreate the ICs with 32-bit IDs.

Alternatively the IC file has an incorrect structure, but itís hard to diagnose this without information where this file comes from.

Volker



On Nov 26, 2015, at 10:42 AM, Sylvia Ploeckinger <ploeckinger_at_strw.leidenuniv.nl> wrote:

> Hi Michael,
>
> thanks a lot for your input!
> The compiler on my system should use 64bit as default,
> but to be sure I included the -m64 flag for both the IC generator
> as well as the code itself. Unfortunately it still doesn't work.
>
> Do you remember which flags you used to get it to work again after this error occurred?
> Was it also -m64 or did you use something else?
>
> Greets,
>
> Sylvia
>
>
> On 25/11/15 19:13, Michael Hansen wrote:
>> HI Sylvia.
>>
>> I had the same problem some while ago. Check that you are compiling both code and IC generator with or without 64 bit compatibility.
>>
>> Cheers,
>> Michael
>>
>> Sent from my iPhone
>>
>> On 25 Nov 2015, at 17:32, Sylvia Ploeckinger <ploeckinger_at_strw.leidenuniv.nl> wrote:
>>
>>> Hi all,
>>>
>>> as the subject line already gives away, Gadget crashes when I try to read in an input file.
>>>
>>> The error is:
>>> I/O error (fread) on task=0 has occured: end of file
>>> task 0: endrun called with an error level of 778
>>>
>>> I found a posting in this mailing list from 2006(!) with the same error:
>>> http://wwwmpa.mpa-garching.mpg.de/gadget/gadget-list/0098.html
>>>
>>> but the solution seems to include some magic:
>>> "So, the good news is that Gadget2 runs with the LCDM gas simulation, but I don't have a clue why it didn't before!"
>>> (from the above posting)
>>>
>>> Did anyone else in the last 9 years have the same problem and knows how to fix it?
>>> I included the following extra lines in the routine my_fread to get some extra information:
>>> printf("nread = %zd, nmemb = %zd\n", nread, nmemb);
>>> printf("fread = %zd\n", fread(ptr, size, nmemb, stream));
>>> printf("size = %zd\n", size);
>>>
>>> and in the file read_ic.c, I included this line (as suggested by Volker in the 2006 discussion)
>>>
>>> printf("task=%d blocknr=%d bytes_per_blockelement=%d npart=%d\n",
>>> ThisTask, blocknr, bytes_per_blockelement, npart);
>>> fflush(stdout);
>>>
>>>
>>> Please find below the corresponding output. It comes from a one core test run, but it crashes with the same error message when I run it on 4 cores.
>>> The IC file seems to be okay, but is there an easy way to find out if the file is damaged to exclude this possibility?
>>>
>>> Thanks in advance!
>>>
>>> Regards,
>>>
>>> Sylvia
>>>
>>>
>>> Allocated 129 MByte for FFT kernel(s).
>>>
>>>
>>> Allocated 84.4627 MByte for particle storage.
>>>
>>> Allocated 1.47507 MByte for storage of SPH data.
>>>
>>>
>>> reading file `/net/para35/data2/ploeckinger/IC_Gadget/output/m1e11f30.ics' on task=0 (contains 288018 particles.)
>>> distributing this file to tasks 0-0
>>> Type 0 (gas): 6445 (tot= 0000006445) masstab=5.28564e-06
>>> Type 1 (halo): 262773 (tot= 0000262773) masstab=2.6309e-05
>>> Type 2 (disk): 15040 (tot= 0000015040) masstab=5.28505e-06
>>> Type 3 (bulge): 3760 (tot= 0000003760) masstab=5.28505e-06
>>> Type 4 (stars): 0 (tot= 0000000000) masstab=0
>>> Type 5 (bndry): 0 (tot= 0000000000) masstab=0
>>>
>>> task=0 blocknr=0 bytes_per_blockelement=12 npart=288018
>>> task=0 blocknr=2 bytes_per_blockelement=12 npart=288018
>>> task=0 blocknr=3 bytes_per_blockelement=4 npart=288018
>>> task=0 blocknr=4 bytes_per_blockelement=4 npart=0
>>> task=0 blocknr=6 bytes_per_blockelement=4 npart=6445
>>> task=0 blocknr=7 bytes_per_blockelement=4 npart=6445
>>> I/O error (fread) on task=0 has occured: end of file
>>> nread = 0, nmemb = 1
>>> fread = 0
>>> size = 4
>>> task 0: endrun called with an error level of 778
>>>
>>>
>>> application called MPI_Abort(MPI_COMM_WORLD, 778) - process 0
>>>
>>> ===================================================================================
>>> = BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES
>>> = PID 16587 RUNNING AT para37.strw.leidenuniv.nl
>>> = EXIT CODE: 10
>>> = CLEANING UP REMAINING PROCESSES
>>> = YOU CAN IGNORE THE BELOW CLEANUP MESSAGES
>>> ===================================================================================
>>>
>>> --
>>> Sylvia Ploeckinger
>>> Sterrewacht Leiden
>>> The Netherlands
>>>
>>> Room Nr. 438
>>> email:
>>> ploeckinger_at_strw.leidenuniv.nl
>>>
>>> -----------------------------------------------------------
>>>
>>> If you wish to unsubscribe from this mailing, send mail to
>>> minimalist_at_MPA-Garching.MPG.de with a subject of: unsubscribe gadget-list
>>> A web-archive of this mailing list is available here:
>>> http://www.mpa-garching.mpg.de/gadget/gadget-list
>
> --
> Sylvia Ploeckinger
> Sterrewacht Leiden
> The Netherlands
>
> Room Nr. 438
> email:
> ploeckinger_at_strw.leidenuniv.nl
>
> -----------------------------------------------------------
>
> If you wish to unsubscribe from this mailing, send mail to
> minimalist_at_MPA-Garching.MPG.de with a subject of: unsubscribe gadget-list
> A web-archive of this mailing list is available here:
> http://www.mpa-garching.mpg.de/gadget/gadget-list
Received on 2015-11-26 11:17:21

This archive was generated by hypermail 2.3.0 : 2022-09-01 14:03:42 CEST