-
Notifications
You must be signed in to change notification settings - Fork 19
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
DynVector dataset shapes sometimes change upon restart #41
Comments
What’s strange is that I think I was seeing the opposite issue when I made the postgkyl change. The Data* always had two dimensions (with 1 element in the second dimension), except after a restart there was a frame with only one dimension.
bpls rt-gk-sheath_esEnergy.bp
double TimeMesh0 {1}
double Data0 {1, 1}
double TimeMesh1 {110}
double Data1 {110, 1}
And after restart it would have given something like Data2 {100}. This is what the postgkyl work-around fixes.
… On Feb 2, 2021, at 10:31 AM, Manaure Francisquez ***@***.***> wrote:
The dataset inside a DynVector sometimes change shape after a restart. Noah tried to prepare postgkyl to handle this (postgkyl commit 89dbc7d5936bd31a4a26e0857bd5261f63cd996f), but it doesn't seem to always work.
This leads to errors like
pgkyl -f gk24-wham1x2v_esEnergy.bp pl
...
ValueError: all the input arrays must have same number of dimensions, but the array at index 0 has 2 dimension(s) and the array at index 1 has 1 dimension(s)
and inspecting the file with bpls we find
bpls gk24-wham1x2v_esEnergy.bp
...
double TimeMesh713 {2292}
double Data713 {2292}
double TimeMesh714 {2292}
double Data714 {2292, 1}
Also, I've also seen the insertion of empty datasets like:
bpls gk24-wham1x2v_elc_intM0.bp
...
double TimeMesh1019 {0}
double Data1019 {0, 1}
double TimeMesh1020 {1}
double Data1020 {1, 1}
double TimeMesh1021 {1}
double Data1021 {1, 1}
double TimeMesh1022 {0}
double Data1022 {0, 1}
but I'm not sure if this is related.
—
You are receiving this because you are subscribed to this thread.
Reply to this email directly, view it on GitHub, or unsubscribe.
|
It seems that the empty data sets can happen even at the beginning of a run, not just upon restart. For example
and more recently I've seen the appearance of scalars:
This, by way, makes the whole file useless from a |
Commit 12acfbd seems to fix the insertion of empty datasets in DynVectors. We were trying out output data more frequently than it was being calculated. So the solution is to use a max of calcIntQuantEvery and nFrame, so that we compute integrated quantities as frequently as we output data. |
I found a way to reproduce this problem of the appearance of the scalars, looping over a file with a bash script that performs many restarts. At around the 30th restart the scalar appears. I have also determined the following, although I don't see how any of these could happen:
|
… an Adios variable in DynVector:write. This shouldn't be necessary according to the ADIOS documentation, but I tried it and it seems to get rid of the scalar entries mentioned in issue #41 #41 (comment).
The dataset inside a DynVector sometimes change shape after a restart. Noah tried to prepare postgkyl to handle this (postgkyl commit 89dbc7d5936bd31a4a26e0857bd5261f63cd996f), but it doesn't seem to always work.
This leads to errors like
and inspecting the file with
bpls
we findAlso, I've also seen the insertion of empty datasets like:
but I'm not sure if this is related.
The text was updated successfully, but these errors were encountered: