Users May 2012

users@lists.einsteintoolkit.org

21 participants
44 discussions

Mclachan, CCZ4: Wrong "phi" term in Ricci tensor?
by Erik Schnetter 21 May '12

21 May '12

The current implementation of the CCZ4 formulation in McLachlan contains the expression (2/phi) in the definition of R[li,lj]. However, the meaning of phi depends on the conformalMethod, and this is not explicitly taken into account here. Is this wrong? -erik -- Erik Schnetter <schnetter(a)cct.lsu.edu> http://www.perimeterinstitute.ca/personal/eschnetter/

3 3

Re: [Users] [Einstein Toolkit] #176: Test parameter files without running them
by Steven R. Brandt 21 May '12

21 May '12

On 05/16/2012 09:25 PM, Einstein Toolkit wrote: > #176: Test parameter files without running them > --------------------------+------------------------------------------------- > Reporter: hinder | Owner: sbrandt > Type: enhancement | Status: review > Priority: minor | Milestone: > Component: Cactus | Version: > Resolution: | Keywords: > --------------------------+------------------------------------------------- > > Comment (by knarf): > > Is the CCTK_Exit(NULL,0); within CactusInitialiseGH intentional? Doesn't > this always abort the run? > Yes. The purpose of the ticket was to cause Cactus to exit right after checking the parameters (if the flag was used). Cheers, Steve

2 1

Einstein Toolkit Meeting
by Frank Loeffler 20 May '12

20 May '12

Hi Please consider joining the weekly Einstein Toolkit phone call at 10 am US central time on Mondays. As usual, you can find instructions how to join on the following web site: http://einsteintoolkit.org/community/support/ I short: the number is (+1) 225-578-4942 or (+1) 866-573-0359 and the conference id is 118682#. The main topic will be the organization of the next ET release. In particular, this means we will discuss: - Release-relevant tickets (18) [1] - Schedule: - soft freeze since May 14th - no additions (thorns ect.), major changes, new features - hard freeze date of May 21th (evening) - only 'important' bug fixes or documentation improvements ^--> to be decided by two maintainers - release date of May 28th. - Now most important: running testsuites, write release notes Not necessarily connected to the release: - What happend during last week: 27 issue were opened or changed and are still open, and 18 tickets could be closed. Some of them were: - Some column headers in Carpet ASCII output were wrong - ADMMass is not part of the ET - GRHydro changes from Zelmani are now (mostly) in the ET - HDF5 output works in AHFinderDirect - Various testsuites have been fixed, and/or made work with Mclachlan instead of BSSN_MoL - EinsteinExact has been added to the ET As usual, feel free to add to this agenda. Frank Loeffler [1] https://trac.einsteintoolkit.org/query?status=accepted&status=assigned&stat…

1 0

weyl scalar symmetries
by Yosef Zlochower 18 May '12

18 May '12

Hi, I noticed that WeylScal4 sets the symmetries of the Weyl scalars by defining the real and imaginary parts of each of the psi's in a different group. RotatingSymmetry and the ReflectionSymmetry thorns define a tensor alias "weylscalars_real", but the symmetries given there do not correspond to those of the psi's. The Weyl symmetries as provided in RotatingSymmetry180 static int const weylparities[10][3] = {{+1,+1,+1}, {-1,-1,-1}, {+1,+1,+1}, {-1,-1,-1}, {+1,+1,+1}, {-1,-1,-1}, {+1,+1,+1}, {-1,-1,-1}, {+1,+1,+1}, {-1,-1,-1}}; The symmetries for rpsi0, ipsi0, rpsi1, ipsi1, ..., rpsi4, ipsi4 static int const weylparities[10][3] = {{+1,+1,+1}, /* rpsi0 */ {-1,-1,-1}, /* ipsi0 */ {+1,+1,-1}, /* rpsi1 */ {-1,-1,+1}, /* ipsi1 */ {+1,+1,+1}, /* rpsi2 */ {-1,-1,-1}, /* ipsi2 */ {+1,+1,-1}, /* rpsi3 */ {-1,-1,+1}, /* ipsi3 */ {+1,+1,+1}, /* rpsi4 */ {-1,-1,-1}}; /* ipsi4 */ -- Dr. Yosef Zlochower Center for Computational Relativity and Gravitation Assistant Professor School of Mathematical Sciences Rochester Institute of Technology 85 Lomb Memorial Drive Rochester, NY 14623 Office:74-2067 Phone: +1 585-475-6103 yosef(a)astro.rit.edu CONFIDENTIALITY NOTE: The information transmitted, including attachments, is intended only for the person(s) or entity to which it is addressed and may contain confidential and/or privileged material. Any review, retransmission, dissemination or other use of, or taking of any action in reliance upon this information by persons or entities other than the intended recipient is prohibited. If you received this in error, please contact the sender and destroy any copies of this information.

2 2

McLachlan/CCZ4: Z terms missing in Rt[li,lj]?
by Erik Schnetter 18 May '12

18 May '12

The definition of the conformally rescaled Ricci tensor Rt[li,lj] has the comment "(* Adding Z term by changing Xtn to Xt *)". However, all the Xtn terms in the definition are still there. Does this mean that there is an omission in this expression? -erik -- Erik Schnetter <schnetter(a)cct.lsu.edu> http://www.perimeterinstitute.ca/personal/eschnetter/

3 2

OpenCL support for release?
by Erik Schnetter 17 May '12

17 May '12

I forgot to bring up yesterday the thorns ExternalLibraries/OpenCL, CactusUtils/Accelerator, CactusUtils/OpenCLRunTime, and CactusExamples/WaveToyOpenCL. These thorns work, OpenCL has documentation, and WaveToyOpenCL has examples and test cases. Is it still okay to add these thorns to the Einstein Toolkit thorn list? I will add documentation to thorns Accelerator and OpenCLRunTime as well. -erik -- Erik Schnetter <schnetter(a)cct.lsu.edu> http://www.perimeterinstitute.ca/personal/eschnetter/

2 2

Re: [Users] [Einstein Toolkit] #905: NaNs in ADMBase::gxx using Development Version of ET
by Yosef Zlochower 16 May '12

16 May '12

This was something I stumbled on when trying to test CCE with the development version of ET. There were nans in the output files for particular observer radii, and only at early times. On 05/16/2012 05:09 PM, Einstein Toolkit wrote: > #905: NaNs in ADMBase::gxx using Development Version of ET > ---------------------------------+------------------------------------------ > Reporter: yosef@… | Type: defect > Status: new | Priority: major > Milestone: | Component: Other > Version: | Keywords: > ---------------------------------+------------------------------------------ > While testing the CCE in the ET, I stumbled across an issue > interpolating the admbase::metric at CCTK_POSTSTEP in GLOBAL mode > at particular spatial points. > > > The "bug" (or parfile/scheduling error) seems to be sensitive to > the exact grid setup. Basically, when I use the > qc0-mclachlan-CCE_Cauchy.par parfile, I get nans > for gxx. I made a small test thorn that only interpolates > gxx at a particular point. Nans show up at iteration 1. > > attached is the parfile I used and the test thorn. > > > Here is the output from the above test thorn: > gxx(3.966899,0.000000,-17.630662) = 1.1160523626974641 > at iter 0 on proc 1 > gxx(3.966899,0.000000,-17.630662) = 1.1160523626974641 > at iter 0 on proc 0 > gxx(3.966899,0.000000,-17.630662) = -nan > at iter 8 on proc 1 > gxx(3.966899,0.000000,-17.630662) = -nan > at iter 8 on proc 0 > gxx(3.966899,0.000000,-17.630662) = -nan > at iter 16 on proc 1 > gxx(3.966899,0.000000,-17.630662) = -nan > at iter 16 on proc 0 > > etc.... > -- Dr. Yosef Zlochower Center for Computational Relativity and Gravitation Assistant Professor School of Mathematical Sciences Rochester Institute of Technology 85 Lomb Memorial Drive Rochester, NY 14623 Office:74-2067 Phone: +1 585-475-6103 yosef(a)astro.rit.edu CONFIDENTIALITY NOTE: The information transmitted, including attachments, is intended only for the person(s) or entity to which it is addressed and may contain confidential and/or privileged material. Any review, retransmission, dissemination or other use of, or taking of any action in reliance upon this information by persons or entities other than the intended recipient is prohibited. If you received this in error, please contact the sender and destroy any copies of this information.

1 0

Einstein Toolkit on Mac OS
by Babiuc-Hamilton, Maria 16 May '12

16 May '12

Hello, I was trying to get the new version, and I had this error: babiuc$ ./GetComponents -p -a https://svn.einsteintoolkit.org/manifest/trunk/einsteintoolkit.th Can't locate object method "new" via package "Thread::Semaphore" at ./GetComponents line 600. Apropos of Mac OS, is there a still maintained wiki tutorial on how to get the Toolkit working? I could not find the path to it. Maria

6 15

Re: [Users] How I got Visit 2.4.2 & CarpetHDF5 to work on Mac OS X Lion
by Scott Hawley 16 May '12

16 May '12

One more thing: I set the environment variable HDF5_DISABLE_VERSION_CHECK=2 On 5/15/12 11:01 AM, "Kelly, Bernard J. (GSFC-660.0)[UNIVERSITY OF MARYLAND BALTIMORE COUNTY]" <bernard.j.kelly(a)nasa.gov> wrote: >Hi Scott. > >Thanks for the detailed walk-through, which of course, I'm *not* following >to the letter. Unfortunately the results of my fiddlings are the same as >before. Details below. > >I really don't think it can be my HDF version per se. Perhaps it's my MPI? >Any reason why MPICH shouldn't work? > > >Bernard > >0th step: I'm using Snow Leopard, not Lion. > > >On 5/14/12 5:23 PM, "Scott Hawley" <scott.hawley(a)belmont.edu> wrote: > >> >>Here's what I did: >> >>1. Downloaded the 64-bit binary of Visit 2.4.2 from the Visit web site > >check > >>2. Downloaded the CarpetHDF directory via svn. > >check > >>3. Downloaded hdf5-1.8.9 source via ftp. I just put it in the CarpetHDF5 >>directory. Did a make and make installŠ. >> "HDF5 Directory" is then .../hdf5-1.8.9/hdf5/ (see below) > >I've stuck with my MacPorts-installed 1.8.8 > >>4. Made sure "Remote Login" is enabled under "Sharing" in "System >>Preferences". Visit uses this. >> Went ahead and did a "ssh localhost" just to make sure. I was >>actually >>prompted to confirm an address because I'm on LSU's network right now. >> (Presumably Visit needs a clean login with no pesky questions or >>password prompts.) > >check > >>5. Created symbolic links in /usr/lib to whichever MPI libraries I >>planned >>to use. >> Older versions of OS X had native MPI support, and this is hard coded >>in the binary you get from the visit site. >> I use OpenMPI and MacPorts, so I went to /usr/lib and did "sudo ln -s >>/opt/local/lib/libmpi* ." > > >check. But I'm using MPICH (MacPorts again), so my MPI libraries are >libmpich, libmpichcxx, libmpichf90. There is an old libmpi in /usr/local, >which seems to have been installed with the OS. I don't want to delete >this without a good reason. > >>6. Deleted my ~/.visit directory entirely. This may not be necessary. > >check > >>7. Went in the CarpetHDF5 and ran ./install and Šbasically just did what >>was in the readme, that is: >> in the xmledit window, clicked on "Cmake" and then updated the flags, >>as >>in >> >> CXXFLAGS = >>-I/Applications/VisIt.app/Contents/Resources/2.4.2/darwin-x86_64/include >>-I/Users/shawley/apps/VizTools/CarpetHDF5/hdf5-1.8.9/hdf5/include >> LDFLAGS = >>-L/Applications/VisIt.app/Contents/Resources/2.4.2/darwin-x86_64/lib >>-L/Users/shawley/apps/VizTools/CarpetHDF5/hdf5-1.8.9/hdf5/lib >> Saved and closed the window, let it do its thing... > >check (though my HDF installation is in /opt/local/include, >/opt/local/lib) > >>8. Ran visit from the command line, using the -noconfig option. The >>latest version of Visit has a problem with "localhost" on Mac, for some >>reason. >>so >>visit -noconfig >> >>..and then visit came up, I opened a .h5 file and was able to do all the >>usual things with it! > > >Sadly, UNcheck. > >

1 0

How I got Visit 2.4.2 & CarpetHDF5 to work on Mac OS X Lion
by Scott Hawley 14 May '12

14 May '12

Here's what I did: 1. Downloaded the 64-bit binary of Visit 2.4.2 from the Visit web site 2. Downloaded the CarpetHDF directory via svn. 3. Downloaded hdf5-1.8.9 source via ftp. I just put it in the CarpetHDF5 directory. Did a make and make installŠ. "HDF5 Directory" is then .../hdf5-1.8.9/hdf5/ (see below) 4. Made sure "Remote Login" is enabled under "Sharing" in "System Preferences". Visit uses this. Went ahead and did a "ssh localhost" just to make sure. I was actually prompted to confirm an address because I'm on LSU's network right now. (Presumably Visit needs a clean login with no pesky questions or password prompts.) 5. Created symbolic links in /usr/lib to whichever MPI libraries I planned to use. Older versions of OS X had native MPI support, and this is hard coded in the binary you get from the visit site. I use OpenMPI and MacPorts, so I went to /usr/lib and did "sudo ln -s /opt/local/lib/libmpi* ." 6. Deleted my ~/.visit directory entirely. This may not be necessary. 7. Went in the CarpetHDF5 and ran ./install and Šbasically just did what was in the readme, that is: in the xmledit window, clicked on "Cmake" and then updated the flags, as in CXXFLAGS = -I/Applications/VisIt.app/Contents/Resources/2.4.2/darwin-x86_64/include -I/Users/shawley/apps/VizTools/CarpetHDF5/hdf5-1.8.9/hdf5/include LDFLAGS = -L/Applications/VisIt.app/Contents/Resources/2.4.2/darwin-x86_64/lib -L/Users/shawley/apps/VizTools/CarpetHDF5/hdf5-1.8.9/hdf5/lib Saved and closed the window, let it do its thing... 8. Ran visit from the command line, using the -noconfig option. The latest version of Visit has a problem with "localhost" on Mac, for some reason. so visit -noconfig ..and then visit came up, I opened a .h5 file and was able to do all the usual things with it! Good luck! -Scott

1 0

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Users May 2012