Hello Prof. Giacomazzo,
Thank you very much for your kind reply.
//In the 10 point case it seems that you are using a lower number of points in Spritz than GRHydroX. In Spritz's plot we see indeed the boundary conditions (and the points start after x=-0.5). //
Thanks for pointing this out. There was an issue with the plots, so I am attaching the new plots with the correct x-axis. For the Spritz case, I considered the domain extending from -0.45 to 0.55 in all directions, and using 10 grid points in each direction. In the GRHydroX case, I was suggested at some point to use periodic boundary conditions because using the regular grid was giving failures (maybe this issue has already been resolved). When using periodic boundary conditions, I considered the domain -1 to 1 using 20 grid points.
Moreover, for the GRHydroX case, in the data output, the positions of the data values were shifted by half the grid spacing (unfortunately, I do not remember the reason for this, and will double-check with others in the next meeting). Therefore, for the Spritz case, I was suggested to shift the grid by half the grid spacing such that the data points coincide for the two cases.
The idea behind this comparison was to check differences between the two results in the first few iterations, which we discussed during the call. These plots, for the final time, were one of them, which Roland suggested to share.
//When using 100 points, GRHydroX seems to overestimate the density at the contact discontinuity, but otherwise it looks fine. //
It is likely that there is still a bug present somewhere in the GRHydroX code. In the next days, I will try to perform shock tube tests with different initial conditions to confirm if this problem persists in such cases too.
Thank you very much again for your comments. Cheers, Jay
On Thu, 24 Mar 2022 at 10:06, Bruno Giacomazzo bruno.giacomazzo@unimib.it wrote:
Hi Jay, thanks for the plots.
In the 10 point case it seems that you are using a lower number of points in Spritz than GRHydroX. In Spritz's plot we see indeed the boundary conditions (and the points start after x=-0.5).
When using 100 points, GRHydroX seems to overestimate the density at the contact discontinuity, but otherwise it looks fine.
Cheers, Bruno
Il giorno mer 23 mar 2022 alle ore 18:30 Jay Vijay Kalinani < jayvijay.kalinani@phd.unipd.it> ha scritto:
Hi all,
Roland had suggested sharing the plots of the shocktube test, which we earlier discussed during the Carpetx call today. I am attaching the plots for two cases: the first case considers a domain of 100 points, and the results are compared with exact solution at t=0.2, while the second case considers 10 points and results are compared with exact solution at t=0.4.
Thank you very much. Cheers, Jay
On Wed, 23 Mar 2022 at 18:00, < carpetx-developers-request@einsteintoolkit.org> wrote:
Send Carpetx-developers mailing list submissions to carpetx-developers@einsteintoolkit.org
To subscribe or unsubscribe via the World Wide Web, visit
http://lists.einsteintoolkit.org/mailman/listinfo/carpetx-developers or, via email, send a message with subject or body 'help' to carpetx-developers-request@einsteintoolkit.org
You can reach the person managing the list at carpetx-developers-owner@einsteintoolkit.org
When replying, please edit your Subject line so it is more specific than "Re: Contents of Carpetx-developers digest..."
Today's Topics:
- Weekly call at 10:00AM Eastern (rhaas@illinois.edu)
- Re: Weekly call at 10:00AM Eastern (Roland Haas)
- Re: Weekly call at 10:00AM Eastern (Roland Haas)
Message: 1 Date: Tue, 22 Mar 2022 16:08:01 -0500 From: rhaas@illinois.edu Subject: [Carpetx-developers] Weekly call at 10:00AM Eastern To: carpetx-developers@mail.einsteintoolkit.org Message-ID: E1nWljV-00F8ws-C1@ekohaes8 Content-Type: text/plain; charset="ANSI_X3.4-1968"
Title: CarpetX developer call When: Weekly from 10am to 11am Eastern Time on Wednesday
** DAYLIGHT SAVING TIME WARNING ** Please note that the US / EU has already / not yet transitioned to / from daylight saving time. The phone call will be at 15:00 Central EU time.
Join Zoom Meeting
https://illinois.zoom.us/j/86863054024?pwd=OXdBcTFSUVBoMkxWKzhXRld3d3RJUT09
Meeting ID: 868 6305 4024 Password: 876596
Message: 2 Date: Tue, 22 Mar 2022 16:16:36 -0500 From: Roland Haas rhaas@illinois.edu Subject: Re: [Carpetx-developers] Weekly call at 10:00AM Eastern To: carpetx-developers@mail.einsteintoolkit.org Message-ID: 20220322161636.37ae2a7e@fdea4908.spdns.org Content-Type: text/plain; charset="utf-8"
Hello all,
Erik will not be able to make it tomorrow. So we will use a different (see below) Zoom link.
Yours, Roland
Title: CarpetX developer call When: Weekly from 10am to 11am Eastern Time on Wednesday
** DAYLIGHT SAVING TIME WARNING ** Please note that the US / EU has already / not yet transitioned to / from daylight saving time. The phone call will be at 15:00 Central EU time.
Join Zoom Meeting
https://illinois.zoom.us/j/86863054024?pwd=OXdBcTFSUVBoMkxWKzhXRld3d3RJUT09
Meeting ID: 868 6305 4024 Password: 876596 _______________________________________________ Carpetx-developers mailing list Carpetx-developers@einsteintoolkit.org
https://urldefense.com/v3/__http://lists.einsteintoolkit.org/mailman/listinf...
-- My email is as private as my paper mail. I therefore support encrypting and signing email messages. Get my PGP key from http://keys.gnupg.net.