Users

users@lists.einsteintoolkit.org

3 participants
3984 discussions

ETK webpage tutorial not working
by Jose Fiestas Iquira 21 Feb '12

21 Feb '12

Dear colleagues, Is the ETK tutorial for new users not working currently? I can not open the link http://docs.einsteintoolkit.org/et-docs/Tutorial_for_New_Users I am trying to install ETK in a local cluster. Hope I am able to do it without much effort. Do you have some experiences in this task? Thanks, Jose

1 0

_No_ Einstein Toolkit Meeting Reminder
by Frank Loeffler 19 Feb '12

19 Feb '12

Hi There will be *no* Einstein Toolkit call tomorrow. Please consider joining the next Einstein Toolkit phone call on Feb 27th at 10 am US central time. 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#. Frank Loeffler

1 0

Meeting Minutes: 13. Feb 2012
by Tanja Bode 14 Feb '12

14 Feb '12

Present were: Peter, Erik, Ian, Elo, Josh, Yosef, Carl Rodriguez Spring Workshop in Atlanta: * Volunteers for Lectures (Josh, Bruno, Roland, Tanja, Frank, ...) please attach your name to sessions within the schedule [1], establishing what you *could* do, not just what you want to do. We will mete out final session assignments accordingly. [1] https://docs.einsteintoolkit.org/et-docs/Schedule_Proposal * Workshop Flyer suggestions (Am I missing some?) * Colorscheme should be ET colors * Should include ET logo * Font is too heavy, should be more inviting * Possible images to entice attendee interest * Text should be more enticing, not just bullet points * QR code for website URL link * Background grid interferes with readability * Some details (i.e. Hotel name/price, registration details) not necessary at the flyer level. * Laptop Requirement: Ian expressed reservations concerning the requirement to bring your own laptop. Possible options are hosting sessions in a computer lab or laptop loans from some source (GT? NSF?). Cheers, Tanja -- Dr. Tanja Bode, PhD Center for Relativistic Astrophysics tanja.bode(a)physics.gatech.edu Office: 2-31A Boggs Phone: (404) 385-8110

2 1

Einstein Toolkit Meeting Reminder
by Frank Loeffler 12 Feb '12

12 Feb '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#. Peter Diener will host this weeks call, and the main topic will be the organization of the Einstein Toolkit workshop in Atlanta [1] again, primarily who can present what material out of [2]. Also, a flyer was prepared by Dennis and, while not yet 100% finished, can be looked at [3] and should be discussed. Frank Loeffler [1] https://docs.einsteintoolkit.org/et-docs/ET_Workshop_Spring_2012 [2] https://docs.einsteintoolkit.org/et-docs/Schedule_Proposal [3] http://www.cct.lsu.edu/~knarf/flyer_test.png

1 0

Backporting ML_ADMConstraints fix
by Ian Hinder 08 Feb '12

08 Feb '12

Hi, I just had a question from someone concerning ML_ADMConstraints. There was an error in the equations, which was found and fixed in November, but this was after the release. I propose that this be back-ported to the release branch. The commits in question are: > commit 5fd1f73f9ca35f0852c0cb61f77abc8b053f5aff > Author: Ian Hinder <ian.hinder(a)aei.mpg.de> > Date: Thu Nov 17 17:47:43 2011 +0100 > > McLachlan_ADMConstraints.m: Fix index errors and > commit b3a25f307c9bb55370380b97e6074df1a30f645d > Author: Ian Hinder <ian.hinder(a)aei.mpg.de> > Date: Sun Nov 20 13:23:59 2011 +0100 > > Regenerate ML_ADMConstraints -- Ian Hinder http://numrel.aei.mpg.de/people/hinder

2 1

Re: [Users] problem with first example simulation with ETK
by Jose Fiestas Iquira 08 Feb '12

08 Feb '12

Hi, The first problem was solved as I reported, by disabling the QuasiLocalMeasures. Some MPI problem. Btw, which information am I missing without QuasiLocalMeasures? (the output file was the wrong one, of a successful simulation without WeyScal4) About the .par file. I found out (please correct me if necessary) that I need to change the name of the run by the simfactory submission (here static_tov): ./simfactory/bin/sim submit *static_tov* --parfile=par/static_tov.par --procs=32 --walltime=8:0:0 Otherwise the run works with the original .par (first used with this run) , and not to the newly specified .par file. So, it worked by doing ./simfactory/bin/sim submit *static_tov01* --parfile=par/*static_tov_new.par * --procs=32 --walltime=8:0:0 Best, Jose On Wed, Feb 8, 2012 at 10:05 AM, Ian Hinder <ian.hinder(a)aei.mpg.de> wrote: > > On 23 Jan 2012, at 22:36, Jose Fiestas Iquira wrote: > > Hi, > > In order to use the wave generator, I am probably not submitting the > job properly. > > I had modified the .par file located at: > /home/ettest24/Cactus/par/qc0-mclachlan.par > since I run the example like this: > ./simfactory/bin/sim submit qc0-mclachlan > --parfile=par/qc0-mclachlan.par --procs=32 --walltime=8:0:0 > > Here the run stopped with errors, as previously reported. > > > Hi Jose, > > Could you update me on your status? Did you solve this problem? In the > standard output file that you sent to the ET mailing list on 16-Jan-2012, I > do not see any problem. The run reached 100 M of evolution time. The only > problem is a non-fatal warning message which is caused by some problem with > the way that Carpet keeps track of points which should or should not be > included in norm operations. You can ignore this warning. > > Could be possible that the proper .par is located somewhere else? I > found it also here: > > /home/ettest24/Cactus/simfactory/etc/par/qc0-mclachlan.par > /home/ettest24/Cactus/arrangements/McLachlan/par/qc0-mclachlan.par > > Btw, in the last one WeylScal4 is activated, in a similar way as you told > me. > > > SimFactory takes the parameter file that you specify on the command line - > it does not look anywhere else. Can you be more specific about the problem > you are having? > > PS: if you ask these questions on the mailing list, other people can help > and benefit from the discussion. > > > > Best, > Jose > > > > On Mon, Jan 23, 2012 at 1:11 PM, Ian Hinder <ian.hinder(a)aei.mpg.de> wrote: > > > On 23 Jan 2012, at 14:57, Jose Fiestas Iquira wrote: > > > Hello, > > > Regarding the Binary Black Hole coalescence example from the ETK. > > > In the BH diagnostics files are columns with quadrupole terms, which > > are not constant. Do you know of some tool in ETK to calculate the GW > > emission from these terms? > > > > No, not that I am aware of. These are the quadrupole moments of the > > individual apparent horizons, rather than the quadrupole moment of the > > system as a whole. The gravitational wave emission can be computed by > using > > the WeylScal4 and Multipole thorns. These give you the spherical harmonic > > multipoles of Psi4 (e.g. l = 2, m = 2 is the dominant one), which can then > > be converted to strain by "two integrations in time" [this is actually > quite > > subtle and hard to get right - see Reisswig and Pollney, arXiv:1006.1632, > > for an overview and solution). > > > > > Thanks, > > Jose > > > On Tue, Jan 17, 2012 at 5:18 AM, Ian Hinder <ian.hinder(a)aei.mpg.de> wrote: > > > > On 17 Jan 2012, at 13:37, Jose Fiestas Iquira wrote: > > > > The Einstein Analysis guide is very useful. I found here > > > > http://einsteintoolkit.org/documentation/ThornGuide/ThornGuidech142.html#x1… > > > the units used in 'Hydrobase'. > > > > If I assume these are the units for all thorns, I obtain from data in > > > BH_diagnostics.ah1.gp a unit of distance of about 1 Km. To estimate > > > the time I am using a mass of ~0.5 M_sun for both BHs (m_irreducible > > > ?) , and I read a final time (cctk_time) ~ 1microsec. Are these > > > quantities in the usual range? > > > > > The length unit is G M / c^2 and the time unit is G M / c^3. For M = > M_sun, > > > I get the length unit as 1.5 km and the time unit as 4.9 us. The final > time > > > in the simulation is 100 M (i.e. 100 G M / c^3), which would be 0.5 ms. > > > > > > > > > Jose > > > > > On Tue, Jan 17, 2012 at 4:04 AM, Jose Fiestas Iquira <jafiestas(a)lbl.gov> > > > wrote: > > > > I sent the wrong file, here is the right one. > > > > > $ hg tip > > > > gives me: > > > > ----------------------------------------------------- > > > > changeset: 3421:dff4df2b3687 > > > > tag: tip > > > > user: Erik Schnetter <schnetter(a)cct.lsu.edu> > > > > date: Thu Jan 12 14:11:27 2012 -0500 > > > > summary: CarpetLib: Allocate space for dstbox before > > > > restricting/prolongating > > > > > [40514 refs] > > > > ----------------------------------------------------- > > > > > I copied this version from: > > > > > ./GetComponents -a > > > > http://svn.einsteintoolkit.org/manifest/trunk/einsteintoolkit.th > > > > > Thanks, > > > > Jose > > > > > > > > On Tue, Jan 17, 2012 at 1:13 AM, Ian Hinder <ian.hinder(a)aei.mpg.de> wrote: > > > > > On 16 Jan 2012, at 23:21, Jose Fiestas Iquira wrote: > > > > > Thanks, I am attaching the output file. > > > > > > The output file looks fine - I can't see any problem with the simulation. > > > > According to the output file, it reached t = 100 M, and in the parameter > > > > file that is listed as cctk_final_time, so all looks well! > > > > > I did copy the development version. As I can see, it was already > > > > fixed, so I should not get this warning, since I copied the > > > > development version. > > > > > > Yes, if you are using the current development version, then I don't know > why > > > > you get the warning. Can you check which version of Carpet you have by > > > > doing > > > > > cd repos/carpet > > > > hg tip > > > > > to be sure you are using the development version? > > > > > Is the puncture tracker not a better way to plot BH trajectories. > > > > > > It's cheaper to use the puncture tracker, because computing the horizon > > > > centroid is more expensive, but it has to be done every iteration, whereas > > > > the horizon can be found less frequently. For comparison with > > > > point-particle results, I don't think it makes much difference; the results > > > > are gauge dependent anyway so differences on the order of the horizon size > > > > cannot be interpreted physically. Since the parameter file already has > > > > horizon-finding built in, it is easier for you to use that. > > > > > If you want to figure out how to add the puncture tracker to the > simulation, > > > > you can start with these parameters: > > > > > ActiveThorns = "PunctureTracker" > > > > > PunctureTracker::track [0] = yes > > > > PunctureTracker::initial_x [0] = $xp > > > > PunctureTracker::which_surface_to_store_info[0] = 0 > > > > PunctureTracker::track [1] = yes > > > > PunctureTracker::initial_x [1] = $xm > > > > PunctureTracker::which_surface_to_store_info[1] = 1 > > > > > IOASCII::out0D_vars = "PunctureTracker::pt_loc" > > > > > Replace $xp and $xm with the initial x coordinates of the +x and -x > > > > punctures (this will be \pm TwoPunctures::par_b). You then need to choose > > > > two spherical surfaces > > > > (https://docs.einsteintoolkit.org/et-docs/Einstein_Toolkit_standards) > which > > > > are not being used currently to store the puncture locations in. If there > > > > is already 0D output, just add the puncturetracker variable separated with > a > > > > space. > > > > > I > > > > did some plots, what are really the combinations xx, xy, xz, yy, zz ? > > > > Are they related to the trajectories of one BH (projections)? What > > > > about the second BH? > > > > > > You should have BH_diagnostics.ah1.gp and BH_diagnostics.ah2.gp. These > > > > correspond to the first and second BHs. If you look at the top of these > > > > files, you will see an explanation of the columns. More information about > > > > each column is available in the AHFinderDirect documentation > > > > ( > http://einsteintoolkit.org/documentation/ThornGuide/ThornGuidech128.html#x1… > ). > > > > Specifically, information about the BH_diagnostics output files is > > > > at > http://einsteintoolkit.org/documentation/ThornGuide/ThornGuidech128.html#x1… > > > > > To plot the trajectory of one of the BHs, you want to plot column 3 and > > > > column 4 in a parametric plot. > > > > > > Cheers, > > > > Jose > > > > > On Mon, Jan 16, 2012 at 4:39 AM, Ian Hinder <ian.hinder(a)aei.mpg.de> wrote: > > > > > > On 14 Jan 2012, at 00:41, Jose Fiestas Iquira wrote: > > > > > > Hi Ian, > > > > > I did the changes but it started and finished early. I am attaching > > > > > the qc0-mclachlan.err file. At the end there is a warning in one of > > > > > the Carpet files (mask_test.c) > > > > > > > The warning you see, > > > > > > WARNING level 1 in thorn CarpetReduce processor 0 host qb517 > > > > > (line 120 of > > > > > /home/ettest24/Cactus/arrangements/Carpet/CarpetReduce/src/mask_test.c): > > > > > -> Simulation domain volume and reduction weight sum differ > > > > > > is a known issue (https://trac.einsteintoolkit.org/ticket/434) which was > > > > > fixed after the last release and the fix is only available in the > > > > > development version. I believe that the warning is not fatal, and it > should > > > > > not have terminated the simulation. Can you also post > the qc0-mclachlan.out > > > > > file? > > > > > > Where can I find the puncture tracker? > > > > > > > The PunctureTracker thorn is part of the Einstein Toolkit thornlist, so you > > > > > should have it already. However, I think the best thing for you to do > > > > > initially would be to look at the BH_diagnostics files for the apparent > > > > > horizon centroid locations for tracking the locations of the BHs. > > > > > > > > > Thanks, > > > > > Jose > > > > > > On Fri, Jan 13, 2012 at 4:38 AM, Ian Hinder <ian.hinder(a)aei.mpg.de> wrote: > > > > > > > On 13 Jan 2012, at 12:32, Jose Fiestas Iquira wrote: > > > > > > > Yes, trajectories and wave forms would be enough for now. Which are > > > > > > the names of this files? Many of the .asc files have data like kxx, > > > > > > kxy, kxz, etc. Is there any description of them? > > > > > > > Not that I am aware of. Files with extension .x.asc are 1D ascii output of > > > > > grid functions from CarpetIOASCII. kxy is the variable name for the x,y > > > > > component of the extrinsic curvature tensor K. The waveform file should be > > > > > in mp_psi4_l2_m2_r*.asc where * corresponds to the radius of the extraction > > > > > sphere. However, looking at the qc0-mclachlan.par file, I see that there > is > > > > > no wave extraction switched on. You should be able to enable wave > > > > > extraction by adding something like the following to the parameter file: > > > > > > > ActiveThorns = "WeylScal4 CarpetInterp AEILocalInterp Multipole" > > > > > > > Multipole::nradii = 3 > > > > > > Multipole::radius[0] = 30 > > > > > > Multipole::radius[1] = 40 > > > > > > Multipole::radius[2] = 50 > > > > > > Multipole::variables = "WeylScal4::Psi4r{sw=-2 cmplx=’WeylScal4::Psi4i’}" > > > > > > Multipole::l_max = 4 > > > > > > WeylScal4::fdOrder = 4 > > > > > > > (I haven't tested this, so let me know if it doesn't work) > > > > > > > I've just found that the usual PunctureTracker thorn is not enabled in this > > > > > example parameter file. So to look at trajectories you will have to use > the > > > > > apparent horizon centroids. These are in the files BH_diagnostics.ah*.gp. > > > > > The meanings of the columns are at the top of the files. You want the > > > > > centroid x and y columns to plot in GnuPlot. > > > > > > > > > > > Thanks, > > > > > > Jose > > > > > > > On Fri, Jan 13, 2012 at 12:41 AM, Ian Hinder <ian.hinder(a)aei.mpg.de> > wrote: > > > > > > > On 12 Jan 2012, at 23:46, Jose Fiestas Iquira wrote: > > > > > > > Dear all, > > > > > > > I run the tutorial and found other interesting examples like the > > > > > > McLachlan tool for Binary Black Holes. > > > > > > > Regarding the McLachlan tool, > > > > > > Could somebody help me finding a description of the output? Is all > > > > > > contained in the .asc files? I would like to visualize it. > > > > > > > > Hi Jose, > > > > > > > What sort of output do you want to visualise? Cactus thorns typically > > > > > > output "lightweight" information such as BH trajectories and waveforms as > > > > > > ASCII, and more heavyweight data -- like gridfunctions -- as HDF5. For the > > > > > > trajectories (from thorn Puncture Tracker) and waveforms (from Multipole) > > > > > > you can use gnuplot or anything that can read standard ASCII output to > > > > > > visualise them. For gridfunctions in HDF5 format, I think the easiest way > > > > > > to get started is to use VisIt, for which there is a Carpet plugin (though > > > > > > installing VisIt is sometimes not so easy!). > > > > > > > Does anyone know if the Carpet plugin is distributed with VisIt yet? I > > > > > > believe that was the intention. > > > > > > > > > Cheers, > > > > > > Jose > > > > > > > > On Wed, Jan 11, 2012 at 7:06 PM, Jose Fiestas Iquira <jafiestas(a)lbl.gov> > > > > > > wrote: > > > > > > > Thanks to Erik for advise. ETK example simulation run successfully. > > > > > > > > > On Wed, Jan 11, 2012 at 11:49 AM, Jose Fiestas Iquira <jafiestas(a)lbl.gov> > > > > > > wrote: > > > > > > > Right, it looks good now: > > > > > > > > ./simfactory/bin/sim list-simulations > > > > > > > static_tov [ACTIVE (RUNNING), restart 0000, job id 595490] > > > > > > > > Let see how it looks when it finishes. The example was set to 8 hours. > > > > > > > > Thanks, > > > > > > > Jose > > > > > > > > > On Wed, Jan 11, 2012 at 11:45 AM, Erik Schnetter <schnetter(a)cct.lsu.edu> > > > > > > wrote: > > > > > > > On Wed, Jan 11, 2012 at 2:42 PM, Jose Fiestas Iquira <jafiestas(a)lbl.gov> > > > > > > wrote: > > > > > > > > The only thing is that now > > > > > > > $ ./simfactory/bin/sim list-simulations > > > > > > > returns > > > > > > > > Error: unknown command list-simulation > > > > > > > Aborting Simfactory. > > > > > > > > Did you really type "list-simulations", and the error message says > > > > > > > "list-simulation" (without the final "s")? > > > > > > > > -erik > > > > > > > > -- > > > > > > > Erik Schnetter <schnetter(a)cct.lsu.edu> http://www.cct.lsu.edu/~eschnett/ > > > > > > > _______________________________________________ > > > > > > > Users mailing list > > > > > > Users(a)einsteintoolkit.org > > > > > > http://lists.einsteintoolkit.org/mailman/listinfo/users > > > > > > > > -- > > > > > > Ian Hinder > > > > > > http://numrel.aei.mpg.de/people/hinder > > > > > > > > -- > > > > > > Ian Hinder > > > > > > http://numrel.aei.mpg.de/people/hinder > > > > > > > <qc0-mclachlan.err> > > > > > > > -- > > > > > Ian Hinder > > > > > http://numrel.aei.mpg.de/people/hinder > > > > > > <qc0-mclachlan.out> > > > > > > -- > > > > Ian Hinder > > > > http://numrel.aei.mpg.de/people/hinder > > > > > > -- > > > Ian Hinder > > > http://numrel.aei.mpg.de/people/hinder > > > > > -- > > Ian Hinder > > http://numrel.aei.mpg.de/people/hinder > > > > -- > Ian Hinder > http://numrel.aei.mpg.de/people/hinder > >

3 2

How to output ADM Mass in ET
by Petr Tsatsin 08 Feb '12

08 Feb '12

Hello, I was wondering how to output an ADM Mass in ET? -- Petr Tsatsin Graduate Student - Department of Physics Charles E. Schmidt College of Science Florida Atlantic University (561)-297-3386 cell: 561-692-2040 ptsatin(a)fau.edu, ptsatsin(a)gmail.com

3 2

"signbit" not defined during compilation of ETK on XSEDE:Kraken?
by Kelly, Bernard J. (GSFC-660.0)[UNIVERSITY OF MARYLAND BALTIMORE COUNTY] 07 Feb '12

07 Feb '12

Hi all. I recently checked out ETK Maxwell on the machine XSEDE machine Kraken, and am having trouble compiling and running a simple test case. The problem is with LSUThorns/Vectors/src/test.cc, where the "signbit" function is (apparently) undefined. No doubt the problem recurs with later source files, but I can't get past this one without abandoning Vectors and hence all of Carpet. Now I can compile a "hello world" program in C invoking signbit.cc without a problem, as long as I include <math.h>, so presumably the ETK process is looking in the wrong location for this library or its headers. Since Kraken is a "known" machine in ETK's simfactory, and lots of Numerical Relativity groups must have used it during the recently expired "NR-AR" allocation, I assume the problem's at my end. Any ideas? Thanks, Bernard Some details: (1) I initially used the simfactory "build" mechanism, but I have this problem with a straightforward Cactus "make" as well. (2) Here's a list of my loaded modules on Kraken. I'm using the PGI compilers (item 18 below): 1) modules/3.1.6.5 2) torque/2.4.14 3) moab/5.4.3.s16991 4) /opt/cray/xt-asyncpe/default/modulefiles/xtpe-istanbul 5) tgusage/3.0-r2 6) altd/1.0 7) DefApps 8) xtpe-target-cnl 9) xt-service/2.2.74 10) xt-os/2.2.74 11) xt-boot/2.2.74 12) xt-lustre-ss/2.2.74_1.6.5 13) cray/job/1.5.5-0.1_2.0202.21413.56.7 14) cray/csa/3.0.0-1_2.0202.21426.77.7 15) cray/account/1.0.0-2.0202.19482.49.18 16) cray/projdb/1.0.0-1.0202.19483.52.1 17) Base-opts/2.2.74 18) pgi/11.4.0 19) xt-libsci/10.5.02 20) pmi/2.1.4-1.0000.8596.15.1.ss 21) xt-mpt/5.2.3 22) xt-pe/2.2.74 23) xt-asyncpe/4.9 24) PrgEnv-pgi/2.2.74 25) cray/MySQL/5.0.64-1.0202.2899.21.1 26) git/1.7.4.2 27) python/2.6.4 28) mercurial/1.6.4 29) subversion/1.6.9 30) hdf5-parallel/1.6.10 --------------------------------------------------------------------------- -- Bernard Kelly -- CRESST Research Associate, NASA/GSFC Phone: +1 (301) 286-7243 E-Mail: bernard.j.kelly(a)nasa.gov Web: http://science.gsfc.nasa.gov/sed/index.cfm?fuseAction=people.jumpBio&iphone bookid=13052 --------------------------------------------------------------------------- --

4 9

output from WeylScal
by Jose Fiestas Iquira 06 Feb '12

06 Feb '12

Dear all, I am using WeylScal for BH binary coalescence, and I am obtaining 'weylscal4::psi4r...' and 'weylscal4::psi4i ...' files, with which I could plot the wave against distance/time. Could somebody advise me how to get the energy emission due to grav. waves from the WeylScal output? I would like to construct a plot: GW Energy vs. distance/time. Thanks, Jose

4 4

Meeting minutes
by Roland Haas 06 Feb '12

06 Feb '12

Hello all, present were: Frank, Roland, Bruno, Tanja, Yosef, Josh, Peter, Christian, Christine Corbett Moran Spring workshop in Atlanta: * schedule at https://docs.einsteintoolkit.org/et-docs/Schedule_Proposal * hours: have to end at 5:30 sharp * shorten lunch hour by 30min, since we are in the same building as the food * distance from hotel: 5min walking * topics for invited talks are not yet fixed * announce soon so that people can take it into account when planning travel. APS and topical group mailing lists, hyperspace, prepare and send around 1 page flyer * there is a suggestion to turn any leftover invited talks into hands-on-sessions Elliptic solver: * Christine Corbett Moran is interested in an elliptic solver for solving a Poisson type equation * willing to spearhead effort * time estimate to make work (for general solver): >= 1/2 year * interest in general community: for initial data in GR, interest in CFD community * see discussion at last workshop: https://docs.einsteintoolkit.org/et-docs/Generic_elliptic_solver 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://keys.gnupg.net.

1 0

Jump to page:

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Users