Re: snapshot file problem

From: Weiguang Cui <cuiweiguang_at_gmail.com>
Date: Mon, 11 Oct 2021 09:35:53 +0100

Hi Volker,



> How did you specify your output times? Have you done this with a list of
> output times? In this case, have you used a second column with a flag
> variable for each output time (0/1)? If you specify the flag, it needs to
> be 1 to get a full dump at the corresponding time, but if you specify a 0,
> you only get the snapshot with the most bound particles.
>

Apologize for this late reply, cosma was down for some time. I didn't have
time to take more tests.
I am using the output file, which only contains a list of the scale factor.
There is no second column in this output list file.
First I thought that can be caused by some compiling options, this is
because I can dump the snapshots by simply commenting out these two lines
in the src/main/run.cc:
```
diff --git a/src/main/run.cc b/src/main/run.cc
index ce4c231..268d5bf 100644
--- a/src/main/run.cc
+++ b/src/main/run.cc
_at__at_ -580,14 +580,15 @@ void sim::create_snapshot_if_desired(void)
 #endif
 #endif

- if(All.DumpFlag_nextoutput)
- {
+ //if(All.DumpFlag_nextoutput)
+ // {
             snap_io Snap(&Sp, Communicator, All.SnapFormat);
/* get an I/O object */
             Snap.write_snapshot(All.SnapshotFileCount, NORMAL_SNAPSHOT);
/* write snapshot file */
- }
+ // }
```
However, I have tried to exclude different compiling options (currently
with -std=c++11 -ggdb -O3 -march=native) which did not work.
I think this could be caused by some compiling parameter which forces
All.DumpFlag_nextoutput (always<1) to 0 in this if condition. But the weird
part is the code compiled before worked fine until snapshot_027, only after
that (recompiled the code for certain), it fails to dump the snapshots. I
also checked if there are other modifications of the code, but not besides
the Makefile.
Please let me know your thoughts.

Best,
Weiguang

>
>
Received on 2021-10-11 10:36:39

This archive was generated by hypermail 2.3.0 : 2023-01-10 10:01:33 CET