Hello,
Hope all of you are doing well. I am trying to replicate some examples given in par directory of the Cactus for practice. I have run the parameter file kasner.par, have obtained the results but my graph is not correct according to the paper arXiv-1111.3344. I am attaching the parameter file which I have run, the resulting .asc file and the graph I have obtained. Kindly, anyone guide me what mistake I am making.
Thanks
Hello Nisa Amir,
(this is just a guess): You are including a file admbase-lapse.norm2.asc in the email which is the the 2-norm of the lapse (\sum_i \alpha(x_i)^2) while the paper in figure 20 shows the value of the metric components. So (assuming that attaching admbase-lapse instead of admbase-gxx was a mistake) I should try and plot something like admbase-gxx.maximum.asc which would be the maximum value on the grid which, since the spacetime is homogeneous, is identical to the value anywhere on the grid.
I am admittedly somewhat doubtful whether using the wrong norm is the only issue here though.
Yours, Roland
Hello,
Hope all of you are doing well. I am trying to replicate some examples given in par directory of the Cactus for practice. I have run the parameter file kasner.par, have obtained the results but my graph is not correct according to the paper arXiv-1111.3344. I am attaching the parameter file which I have run, the resulting .asc file and the graph I have obtained. Kindly, anyone guide me what mistake I am making.
Thanks
I have tried plotting the admbase-maximum.asc but again getting the same results. I am attaching my gnuplot commands, maybe I am making some mistake here. KIndly have a look. #set term png enhanced font "Vera,18" crop set term postscript enhanced color #dashed set size 0.6, 0.6 set xlabel 't' set style line 1 lt 1 lw 3 pt 1 set style line 2 lt 2 lw 3 pt 2 set style line 3 lt 3 lw 3 pt 3 set style line 4 lt 4 lw 3 pt 4 set style line 5 lt 5 lw 3 pt 5 set style line 6 lt rgb "blue" #lw 0.1 pt 6 set style line 7 lt rgb "red" #lw 0.1 pt 6 set key left Left reverse set yrange [:200] set log y set output "kasner.eps" plot 'admbase-metric.norm2.asc' ev 5 u 2:3 w p ls 7 title "g_{xx} (numerical)",\ 'admbase-metric.norm2.asc' ev 5 u 2:($2**(4.0/3.0)) w l ls 7 title "g_{xx} (exact)",\ 'admbase-metric.norm2.asc' ev 5 u 2:8 w p ls 6 title "g_{zz} (numerical)",\ 'admbase-metric.norm2.asc' ev 5 u 2:($2**(-2.0/3.0)) w l ls 6 title "g_{zz} (exact)"
On Fri, Apr 17, 2020 at 3:02 AM Roland Haas rhaas@illinois.edu wrote:
Hello Nisa Amir,
(this is just a guess): You are including a file admbase-lapse.norm2.asc in the email which is the the 2-norm of the lapse (\sum_i \alpha(x_i)^2) while the paper in figure 20 shows the value of the metric components. So (assuming that attaching admbase-lapse instead of admbase-gxx was a mistake) I should try and plot something like admbase-gxx.maximum.asc which would be the maximum value on the grid which, since the spacetime is homogeneous, is identical to the value anywhere on the grid.
I am admittedly somewhat doubtful whether using the wrong norm is the only issue here though.
Yours, Roland
Hello,
Hope all of you are doing well. I am trying to replicate some examples given in par directory of the Cactus for practice. I have run the
parameter
file kasner.par, have obtained the results but my graph is not correct according to the paper arXiv-1111.3344. I am attaching the parameter file which I have run, the resulting .asc file and the graph I have obtained. Kindly, anyone guide me what mistake I am making.
Thanks
-- My email is as private as my paper mail. I therefore support encrypting and signing email messages. Get my PGP key from http://pgp.mit.edu . _______________________________________________ Users mailing list Users@einsteintoolkit.org http://lists.einsteintoolkit.org/mailman/listinfo/users
Hello Nisa Amir,
from your gnuplot commands:
plot 'admbase-metric.norm2.asc' ev 5 u 2:3 w p ls 7 title "g_{xx}
you are still using the 2-norm since the filename contains "norm2". The file that contains the maximum values will be "admbase-metric.maximum.asc". If not such file is output, which, given your parameter file, seems likely, then you can either edit the parameter file by changing:
IOBasic::outInfo_reductions = "norm2"
to read
IOBasic::outInfo_reductions = "norm2 maximum"
or change:
IOASCII::out0D_every = 0 IOASCII::out0D_vars = " Carpet::timing "
to
IOASCII::out0D_every = 1 IOASCII::out0D_vars = " Carpet::timing ADMBase::metric "
and use the file admbase-metric..asc (two "." before "asc"),
or, alternatively, use the file
admbase-metric.x.asc
which is already output. That file contains the values of the metric along the x-axis and since the Kasner spacetime is homogeneous you can pick any of the values to plot them.
To pick the value for just one point you can try to use some inline filtering in gnuplot, something like:
plot "gawk '$10==0. && $11==0. && $12 == 0.' admbase-metric.x.asc" u 10:13 w p ls 7 title "g_{xx} (numerical)"
which select only those lines where columns 10, 11 and 12 (x,y,z coordinates) are 0.
Yours, Roland
I have tried plotting the admbase-maximum.asc but again getting the same results. I am attaching my gnuplot commands, maybe I am making some mistake here. KIndly have a look. #set term png enhanced font "Vera,18" crop set term postscript enhanced color #dashed set size 0.6, 0.6 set xlabel 't' set style line 1 lt 1 lw 3 pt 1 set style line 2 lt 2 lw 3 pt 2 set style line 3 lt 3 lw 3 pt 3 set style line 4 lt 4 lw 3 pt 4 set style line 5 lt 5 lw 3 pt 5 set style line 6 lt rgb "blue" #lw 0.1 pt 6 set style line 7 lt rgb "red" #lw 0.1 pt 6 set key left Left reverse set yrange [:200] set log y set output "kasner.eps" plot 'admbase-metric.norm2.asc' ev 5 u 2:3 w p ls 7 title "g_{xx} (numerical)",\ 'admbase-metric.norm2.asc' ev 5 u 2:($2**(4.0/3.0)) w l ls 7 title "g_{xx} (exact)",\ 'admbase-metric.norm2.asc' ev 5 u 2:8 w p ls 6 title "g_{zz} (numerical)",\ 'admbase-metric.norm2.asc' ev 5 u 2:($2**(-2.0/3.0)) w l ls 6 title "g_{zz} (exact)"
On Fri, Apr 17, 2020 at 3:02 AM Roland Haas rhaas@illinois.edu wrote:
Hello Nisa Amir,
(this is just a guess): You are including a file admbase-lapse.norm2.asc in the email which is the the 2-norm of the lapse (\sum_i \alpha(x_i)^2) while the paper in figure 20 shows the value of the metric components. So (assuming that attaching admbase-lapse instead of admbase-gxx was a mistake) I should try and plot something like admbase-gxx.maximum.asc which would be the maximum value on the grid which, since the spacetime is homogeneous, is identical to the value anywhere on the grid.
I am admittedly somewhat doubtful whether using the wrong norm is the only issue here though.
Yours, Roland
Hello,
Hope all of you are doing well. I am trying to replicate some examples given in par directory of the Cactus for practice. I have run the
parameter
file kasner.par, have obtained the results but my graph is not correct according to the paper arXiv-1111.3344. I am attaching the parameter file which I have run, the resulting .asc file and the graph I have obtained. Kindly, anyone guide me what mistake I am making.
Thanks
-- My email is as private as my paper mail. I therefore support encrypting and signing email messages. Get my PGP key from http://pgp.mit.edu . _______________________________________________ Users mailing list Users@einsteintoolkit.org http://lists.einsteintoolkit.org/mailman/listinfo/users
I have tried everything but again getting the same results. A day before Erick schentter send me an asc file they have used for plotting the figure. their value of gxx is growing slowly but in my file it grows rapidly. Attaching both files. is there something to do with the values of parameters p1 and p2?
On Fri, Apr 17, 2020 at 3:00 PM Roland Haas rhaas@illinois.edu wrote:
Hello Nisa Amir,
from your gnuplot commands:
plot 'admbase-metric.norm2.asc' ev 5 u 2:3 w p ls 7 title "g_{xx}
you are still using the 2-norm since the filename contains "norm2". The file that contains the maximum values will be "admbase-metric.maximum.asc". If not such file is output, which, given your parameter file, seems likely, then you can either edit the parameter file by changing:
IOBasic::outInfo_reductions = "norm2"
to read
IOBasic::outInfo_reductions = "norm2 maximum"
or change:
IOASCII::out0D_every = 0 IOASCII::out0D_vars = " Carpet::timing "
to
IOASCII::out0D_every = 1 IOASCII::out0D_vars = " Carpet::timing ADMBase::metric "
and use the file admbase-metric..asc (two "." before "asc"),
or, alternatively, use the file
admbase-metric.x.asc
which is already output. That file contains the values of the metric along the x-axis and since the Kasner spacetime is homogeneous you can pick any of the values to plot them.
To pick the value for just one point you can try to use some inline filtering in gnuplot, something like:
plot "gawk '$10==0. && $11==0. && $12 == 0.' admbase-metric.x.asc" u 10:13 w p ls 7 title "g_{xx} (numerical)"
which select only those lines where columns 10, 11 and 12 (x,y,z coordinates) are 0.
Yours, Roland
I have tried plotting the admbase-maximum.asc but again getting the same results. I am attaching my gnuplot commands, maybe I am making some
mistake
here. KIndly have a look. #set term png enhanced font "Vera,18" crop set term postscript enhanced color #dashed set size 0.6, 0.6 set xlabel 't' set style line 1 lt 1 lw 3 pt 1 set style line 2 lt 2 lw 3 pt 2 set style line 3 lt 3 lw 3 pt 3 set style line 4 lt 4 lw 3 pt 4 set style line 5 lt 5 lw 3 pt 5 set style line 6 lt rgb "blue" #lw 0.1 pt 6 set style line 7 lt rgb "red" #lw 0.1 pt 6 set key left Left reverse set yrange [:200] set log y set output "kasner.eps" plot 'admbase-metric.norm2.asc' ev 5 u 2:3 w p ls 7 title "g_{xx} (numerical)",\ 'admbase-metric.norm2.asc' ev 5 u 2:($2**(4.0/3.0)) w l ls 7 title "g_{xx} (exact)",\ 'admbase-metric.norm2.asc' ev 5 u 2:8 w p ls 6 title "g_{zz} (numerical)",\ 'admbase-metric.norm2.asc' ev 5 u 2:($2**(-2.0/3.0)) w l ls 6 title "g_{zz} (exact)"
On Fri, Apr 17, 2020 at 3:02 AM Roland Haas rhaas@illinois.edu wrote:
Hello Nisa Amir,
(this is just a guess): You are including a file admbase-lapse.norm2.asc in the email which is the the 2-norm of the lapse (\sum_i \alpha(x_i)^2) while the paper in figure 20 shows the value of the metric components. So (assuming that attaching admbase-lapse instead of admbase-gxx was a mistake) I should try and plot something like admbase-gxx.maximum.asc which would be the maximum value on the grid which, since the spacetime is homogeneous, is identical to the value anywhere on the grid.
I am admittedly somewhat doubtful whether using the wrong norm is the only issue here though.
Yours, Roland
Hello,
Hope all of you are doing well. I am trying to replicate some
examples
given in par directory of the Cactus for practice. I have run the
parameter
file kasner.par, have obtained the results but my graph is not
correct
according to the paper arXiv-1111.3344. I am attaching the parameter
file
which I have run, the resulting .asc file and the graph I have
obtained.
Kindly, anyone guide me what mistake I am making.
Thanks
-- My email is as private as my paper mail. I therefore support encrypting and signing email messages. Get my PGP key from http://pgp.mit.edu . _______________________________________________ Users mailing list Users@einsteintoolkit.org http://lists.einsteintoolkit.org/mailman/listinfo/users
-- My email is as private as my paper mail. I therefore support encrypting and signing email messages. Get my PGP key from http://pgp.mit.edu .
Hello Nisa Amir,
I have tried everything but again getting the same results. A day before Erick schentter send me an asc file they have used for plotting the figure. their value of gxx is growing slowly but in my file it grows rapidly. Attaching both files. is there something to do with the values of parameters p1 and p2?
I would not be able to tell. However I (finally) gave this a try see a similar very rapid growth in gxx (much more rapid than the growth in the paper). I will try and see if I can find out which version of the Einstein Toolkit first shows this behaviour.
Yours, Roland
Thanks alott
On Sun, 26 Apr 2020, 9:02 pm Roland Haas, rhaas@illinois.edu wrote:
Hello Nisa Amir,
I have tried everything but again getting the same results. A day
before
Erick schentter send me an asc file they have used for plotting the
figure.
their value of gxx is growing slowly but in my file it grows rapidly. Attaching both files. is there something to do with the values of parameters p1 and p2?
I would not be able to tell. However I (finally) gave this a try see a similar very rapid growth in gxx (much more rapid than the growth in the paper). I will try and see if I can find out which version of the Einstein Toolkit first shows this behaviour.
Yours, Roland
-- My email is as private as my paper mail. I therefore support encrypting and signing email messages. Get my PGP key from http://pgp.mit.edu .
Hello Nisa Amir,
Thanks alott
Alright. Here is how to make it work. Please add these lines to your parameters file (kasner.par or kasner-med.par or kasner-hi.par):
# ML_BSSN ignores gauge evolution method parameters, so manually force its RHS # to be all zero ML_BSSN::harmonicF = 0. ML_BSSN::LapseAdvectionCoeff = 0. ML_BSSN::ShiftGammaCoeff = 0. ML_BSSN::ShiftAdvectionCoeff = 0.
somewhere after the other ML_BSSN settings in that file.
The lets me run and produces the correct answer, at least for gxx.
I will add more detailed information to the ticket I created:
https://bitbucket.org/einsteintoolkit/tickets/issues/2373/kasner-example-fro...
Yours, Roland
users@lists.einsteintoolkit.org
-
Nisa Amir -
Roland Haas