[Carpetx-developers] Carpetx-developers Digest, Vol 7, Issue 4

Thu Mar 24 12:55:56 CDT 2022

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 at 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 at 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 at einsteintoolkit.org> wrote:
>>
>>> Send Carpetx-developers mailing list submissions to
>>>         carpetx-developers at 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 at einsteintoolkit.org
>>>
>>> You can reach the person managing the list at
>>>         carpetx-developers-owner at einsteintoolkit.org
>>>
>>> When replying, please edit your Subject line so it is more specific
>>> than "Re: Contents of Carpetx-developers digest..."
>>>
>>>
>>> Today's Topics:
>>>
>>>    1. Weekly call at 10:00AM Eastern (rhaas at illinois.edu)
>>>    2. Re: Weekly call at 10:00AM Eastern (Roland Haas)
>>>    3. Re: Weekly call at 10:00AM Eastern (Roland Haas)
>>>
>>>
>>> ----------------------------------------------------------------------
>>>
>>> Message: 1
>>> Date: Tue, 22 Mar 2022 16:08:01 -0500
>>> From: rhaas at illinois.edu
>>> Subject: [Carpetx-developers] Weekly call at 10:00AM Eastern
>>> To: carpetx-developers at mail.einsteintoolkit.org
>>> Message-ID: <E1nWljV-00F8ws-C1 at 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 at illinois.edu>
>>> Subject: Re: [Carpetx-developers] Weekly call at 10:00AM Eastern
>>> To: carpetx-developers at mail.einsteintoolkit.org
>>> Message-ID: <20220322161636.37ae2a7e at 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 at einsteintoolkit.org
>>> >
>>> https://urldefense.com/v3/__http://lists.einsteintoolkit.org/mailman/listinfo/carpetx-developers__;!!DZ3fjg!oxcHskS1z9r4YyKPUBR3thucel9FU1sGJNd1lM8XoN3Uc0Jei_c6dwnxU5gZX3Ha$
>>>
>>>
>>>
>>> --
>>> 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.
>>> -------------- next part --------------
>>> A non-text attachment was scrubbed...
>>> Name: not available
>>> Type: application/pgp-signature
>>> Size: 833 bytes
>>> Desc: OpenPGP digital signature
>>> Url :
>>> http://lists.einsteintoolkit.org/pipermail/carpetx-developers/attachments/20220322/50a66f29/attachment-0001.bin
>>>
>>> ------------------------------
>>>
>>> Message: 3
>>> Date: Wed, 23 Mar 2022 10:18:24 -0500
>>> From: Roland Haas <rhaas at illinois.edu>
>>> Subject: Re: [Carpetx-developers] Weekly call at 10:00AM Eastern
>>> To: carpetx-developers at mail.einsteintoolkit.org
>>> Message-ID: <20220323101824.40438d93 at fdea4908.spdns.org>
>>> Content-Type: text/plain; charset="utf-8"
>>>
>>> present: Roland, Lorenzo, Jay, Yosef, Federico
>>>
>>> hydro toy code
>>> * Jay ran a trial for shock tube using low resolution (100 points
>>>   total) comparing HydroToy and Spritz. Both do similarly well, if
>>>   anything Spritz does worse since Spritz has the wrong value on the
>>>   left hand side of the shock, while HydroToy's was more correct.
>>> * Jay will run some more of the Balsara tests
>>> * Jay to send around plots with results to mailing list
>>>
>>> ET summer school talk
>>> * Jay, and Federico and Lorenzo (in order of involvemet), will contact
>>>   Leo who is the executive organizer
>>> * target audience would be ET developers who want to use GPUs in the
>>>   ET, which is facilitated by CarpetX
>>>
>>> late TODO items:
>>> * Roland has not yet worked on the unit test for HydroToy
>>> * no decision on new (if any) name for GRHydroToyGPU yet, "Primavera"
>>>   or "GRHydroX" were suggested.
>>>
>>> > 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 at einsteintoolkit.org
>>> > >
>>> https://urldefense.com/v3/__http://lists.einsteintoolkit.org/mailman/listinfo/carpetx-developers__;!!DZ3fjg!oxcHskS1z9r4YyKPUBR3thucel9FU1sGJNd1lM8XoN3Uc0Jei_c6dwnxU5gZX3Ha$
>>>
>>> >
>>> >
>>> >
>>>
>>>
>>>
>>> --
>>> 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.
>>> -------------- next part --------------
>>> A non-text attachment was scrubbed...
>>> Name: not available
>>> Type: application/pgp-signature
>>> Size: 833 bytes
>>> Desc: OpenPGP digital signature
>>> Url :
>>> http://lists.einsteintoolkit.org/pipermail/carpetx-developers/attachments/20220323/01122914/attachment-0001.bin
>>>
>>> ------------------------------
>>>
>>> _______________________________________________
>>> Carpetx-developers mailing list
>>> Carpetx-developers at einsteintoolkit.org
>>> http://lists.einsteintoolkit.org/mailman/listinfo/carpetx-developers
>>>
>>>
>>> End of Carpetx-developers Digest, Vol 7, Issue 4
>>> ************************************************
>>>
>> _______________________________________________
>> Carpetx-developers mailing list
>> Carpetx-developers at einsteintoolkit.org
>> http://lists.einsteintoolkit.org/mailman/listinfo/carpetx-developers
>>
>
>
> --
>
> Prof. Bruno Giacomazzo
> Department of Physics
> University of Milano-Bicocca
> Piazza della Scienza 3
> 20126 Milano
> Italy
>
> email: bruno.giacomazzo at unimib.it
> phone: (+39) 02 6448 2321
> web: http://www.brunogiacomazzo.org
>
> ----------------------------------------------------------------------
> There are only 10 types of people in the world:
> Those who understand binary, and those who don't
> ----------------------------------------------------------------------
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.einsteintoolkit.org/pipermail/carpetx-developers/attachments/20220324/a8c2056a/attachment-0001.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: shocktube_10pts.pdf
Type: application/pdf
Size: 74203 bytes
Desc: not available
Url : http://lists.einsteintoolkit.org/pipermail/carpetx-developers/attachments/20220324/a8c2056a/attachment-0001.pdf 