Users

users@lists.einsteintoolkit.org

3 participants
3984 discussions

link problems
by Mohror, Kathryn 26 Feb '11

26 Feb '11

Hi, I'm new at Einstein/Cactus and am having trouble at link time when following the tutorial instructions (https://docs.einsteintoolkit.org/et-docs/Tutorial_for_New_Users#Any_Problem…) I've gotten to the point where the output of: ./simfactory/sim build --optionlist=simfactory/optionlists/generic.cfg --thornlist=manifest/einsteintoolkit.th Says: Creating cactus_sim in /g/g23/mohror1/Einstein/Cactus/exe from ... (followed by lots of names of thorns) Then, I get a load of undefined reference errors, the first two of which are: /g/g23/mohror1/Einstein/Cactus/configs/sim/lib/libthorn_CactusBindings.a(EOSG_Polytrope_Functions.c.o): In function `CCTK_Wrapper_FtoC_EOS_GP_DPDIE': EOSG_Polytrope_Functions.c:(.text+0xe): undefined reference to `eos_gp_dpdie_' /g/g23/mohror1/Einstein/Cactus/configs/sim/lib/libthorn_CactusBindings.a(EOSG_Polytrope_Functions.c.o): In function `CCTK_Wrapper_FtoC_EOS_GP_IntEn': EOSG_Polytrope_Functions.c:(.text+0x2e): undefined reference to `eos_gp_inten_' In configs/sim/lib there is a libthorn_EOSG_Polytrope.a which I believe should contain these symbols. However, when I do an 'nm' of the library, the symbols are undefined: configs/sim/lib/libthorn_EOSG_Polytrope.a: EOS_GP_Setup.c.o: U CCTK_Warn U EOS_GP_Inverse_SetArray U EOS_GP_SetArray 0000000000000000 T EOS_GP_Startup U EOS_RegisterCall U Util_TableCreate U Util_TableSetInt U Util_TableSetString EOS_GP.c.o: 0000000000000000 T EOS_GP_Inverse_SetArray 00000000000000e0 T EOS_GP_SetArray U eos_gp_cs2_ U eos_gp_dpdie_ U eos_gp_dpdrho_ U eos_gp_inten_ U eos_gp_inv_cs2_ U eos_gp_inv_dpdie_ U eos_gp_inv_dpdrho_ U eos_gp_inv_inten_ U eos_gp_inv_rho_ U eos_gp_pressure_ Any ideas on what I'm doing wrong? (I'm using gcc/g++/gfortran on a Linux cluster) Thanks, -- Kathryn

2 3

segfault with kerr initial data using IDAnalyticBH
by Yosef Zlochower 24 Feb '11

24 Feb '11

Hi, If initial_data=kerr and metric_type=physical, IDAnalyticBH writes to psi[x-z] and psi[x-z][x-z] (it tries to set these GFs to zero) even if no storage is allocated. The problem is with the code IDAnalyticBH_zero_CCTK_REAL_array(npoints, psix); etc.., which is only executed when the metric is physical. I was able to use the Kerr data with the following patch, but just deleting the offending lines should be sufficient. --- Kerr.c 2011-02-10 09:08:58.000000000 -0500 +++ /tmp/Kerr.c 2011-02-23 18:15:38.000000000 -0500 @@ -103,6 +103,10 @@ conformal_storage); /*NOTREACHED*/ } } + else + { + *conformal_state = 0; + } /* printf("npoints: %i\n",npoints); */ @@ -311,15 +315,22 @@ } IDAnalyticBH_zero_CCTK_REAL_array(npoints, gxz); IDAnalyticBH_zero_CCTK_REAL_array(npoints, gyz); - IDAnalyticBH_zero_CCTK_REAL_array(npoints, psix); - IDAnalyticBH_zero_CCTK_REAL_array(npoints, psiy); - IDAnalyticBH_zero_CCTK_REAL_array(npoints, psiz); - IDAnalyticBH_zero_CCTK_REAL_array(npoints, psixx); - IDAnalyticBH_zero_CCTK_REAL_array(npoints, psixy); - IDAnalyticBH_zero_CCTK_REAL_array(npoints, psixz); - IDAnalyticBH_zero_CCTK_REAL_array(npoints, psiyy); - IDAnalyticBH_zero_CCTK_REAL_array(npoints, psiyz); - IDAnalyticBH_zero_CCTK_REAL_array(npoints, psizz); + + if (psix && *conformal_state < 2) + { + IDAnalyticBH_zero_CCTK_REAL_array(npoints, psix); + IDAnalyticBH_zero_CCTK_REAL_array(npoints, psiy); + IDAnalyticBH_zero_CCTK_REAL_array(npoints, psiz); + } + if (psixx && *conformal_state < 3) + { + IDAnalyticBH_zero_CCTK_REAL_array(npoints, psixx); + IDAnalyticBH_zero_CCTK_REAL_array(npoints, psixy); + IDAnalyticBH_zero_CCTK_REAL_array(npoints, psixz); + IDAnalyticBH_zero_CCTK_REAL_array(npoints, psiyy); + IDAnalyticBH_zero_CCTK_REAL_array(npoints, psiyz); + IDAnalyticBH_zero_CCTK_REAL_array(npoints, psizz); + } } return;

1 0

Output with carpet in more than 1 proc
by Francisco S. Guzman 24 Feb '11

24 Feb '11

Dear all, I just coded a simple wavetoy and have found that the output when using carpet (the number of time levels) depends on the number of processors used for both h5 or ascii. Am I doing something clearly worng and easily fixable? Bests, Francisco -- Francisco S. Guzman Murillo. Instituto de Fisica y Matematicas, Universidad Michoacana de San Nicolas de Hidalgo Phone: +52 (443) 322 3500 ext 4142 http://www.ifm.umich.mx/guzman/

3 2

Re: [Users] [Commits] [svn:einsteintoolkit] HydroBase/trunk/ (Rev. 38)
by Bruno Coutinho Mundim 21 Feb '11

21 Feb '11

Hi Christian: > File [modified]: interface.ccl > Delta lines: +2 -0 > =================================================================== > --- trunk/interface.ccl 2011-02-21 15:49:06 UTC (rev 37) > +++ trunk/interface.ccl 2011-02-21 17:14:39 UTC (rev 38) > @@ -35,6 +35,8 @@ > > CCTK_REAL vel[3] type = GF Timelevels = 3 tags='ProlongationParameter="HydroBase::prolongation_type" tensortypealias="U" interpolator="matter"' "velocity v^i" > > +CCTK_REAL w_lorentz type = GF Timelevels = 3 tags='ProlongationParameter="HydroBase::prolongation_type" tensortypealias="U" interpolator="matter"' "Lorentz Factor" > + > CCTK_REAL Y_e type = GF Timelevels = 3 tags='ProlongationParameter="HydroBase::prolongation_type" tensortypealias="Scalar" interpolator="matter"' "Electron Fraction" Shouldn't we have tensortypealias="Scalar" tensorweight=+1.0 instead of tensortypealias="U" for the Lorentz factor? Cheers, Bruno.

3 3

Cactus compilation issue
by Steven R. Brandt 21 Feb '11

21 Feb '11

I'm attempting to compile Einsteintoolkit on my laptop. I have gcc/g++/gfortran. The config looks like this: CPP = gcc FPP = gfortran CC = gcc CXX = g++ F77 = gfortran F90 = gfortran CPPFLAGS = -DMPICH_IGNORE_CXX_SEEK -DKRANC_DIFF_FUNCTIONS FPPFLAGS = -traditional # no -msse3 CFLAGS = -g3 -m128bit-long-double -march=native -std=gnu99 CXXFLAGS = -g3 -m128bit-long-double -march=native -std=gnu++0x F77FLAGS = -g3 -m128bit-long-double -march=native -fcray-pointer -ffixed-line-length-none F90FLAGS = -g3 -m128bit-long-double -march=native -fcray-pointer -ffixed-line-length-none LDFLAGS = My compiler gives a slew of errors, starting with this one. Checking status of thorn TwoPunctures ________________________________________________________________________ /home/sbrandt/Cactus/arrangements/LSUThorns/SummationByParts/src/Poisoning.F90:9.26: REAL*8 diss_fraction(3)&&REAL*8 epsdis&&REAL*8 h_scaling(3)&&REAL*8 poiso What's going wrong? How can I fix it? Cheers, Steve 1

2 2

Meetinq minutes
by Roland Haas 21 Feb '11

21 Feb '11

Present were: Frank, Tanja, Roland, Peter, Ian, Eloisa, Bruno, Erik, Josh, Christian HydroBase: * include w_lorentz as a grid function in HydroBase ET paper: * update references to use ET bibtex file (Peter, Frank, anyone else is welcome!) * add section on (Cactus) component structure/component interface (Erik) * Kranc section is finished (Ian) * ID data will be looked at during the next week (Josh), will likely have go over the 1/2 page limit * Evolution methods (Christian) * two level approach: ** produce a first "text only" paper that will be a CCT technical report and also posted to the arXive ** send to interested parties (LRR, Computational Science & Discovery) to ascertain if there is interest in a longer more detailed version (50-60 pages) showing basic equations, numerical implementation ** long version would be core reference for ET, would need more substance (more test cases) *** http://iopscience.iop.org/1749-4699/1/1/015005, http://arxiv.org/abs/0804.1394 *** http://iopscience.iop.org/1749-4699/1/1/015003 * add stellar core collapse to sample results spring ET workshop: * afternoon Mon April 25th - morning Wed April 27th * GT will find possible hotel venues 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

Meeting Reminder
by Frank Loeffler 21 Feb '11

21 Feb '11

Hi, This is a reminder about our usual Monday Einstein Toolkit phone call at 10 am US central time. As usual, you can find instructions how to join on the following web site: http://einsteintoolkit.org/info/support/ Some points of tomorrows discussion will be: - ET paper plans - ET workshops http://doodle.com/ayigys5kv7zzyxna -> looks like Apr Mon 25 - Wed 27? - "Poll" Visualization of simulation results https://docs.einsteintoolkit.org/et-docs/Visualization_of_simulation_results Connected to this, Ian Hinder will give an EVO-presentation on the mathematica scripts he uses. The EVO meeting is called "Einstein Toolkit weekly" and the link to it should be http://evo.caltech.edu/evoNext/koala.jnlp?meeting=eIe9esvvv8anaiatII We will use the phone line for audio, not evo. Evo will only be used for the slides/presentations. Frank Löffler

1 0

Compilation failure with xlc compiler
by Wolfgang Kastaun 19 Feb '11

19 Feb '11

Hi, I am trying to compile the Chandrasekhar release of the Einstein Toolkit on an AIX cluster (CINECA) using the native xlc compiler. It fails with the following message: "/sp6/usersissa/sists082/sim_build/pizza_dev/configs/pizza_dev/build/LoopControl/loopcontrol.c", line 46.8: 1506-1108 (I) The use of keyword '__inline__' is non-portable. "/sp6/usersissa/sists082/sim_build/pizza_dev/configs/pizza_dev/build/LoopControl/loopcontrol.c", line 946.1: 1506-343 (S) Redeclaration of lc_statmap_init differs from previous declaration on line 513 of "/sp6/usersissa/sists082/sim_build/pizza_dev/configs/pizza_dev/build/LoopControl/loopcontrol.c". "/sp6/usersissa/sists082/sim_build/pizza_dev/configs/pizza_dev/build/LoopControl/loopcontrol.c", line 946.1: 1506-376 (I) Redeclaration of lc_statmap_init has a different number of fixed parameters than the previous declaration. "/sp6/usersissa/sists082/sim_build/pizza_dev/configs/pizza_dev/build/LoopControl/loopcontrol.c", line 946.1: 1506-377 (I) The type "char*" of parameter 3 differs from the previous type "const char* restrict const". "/sp6/usersissa/sists082/sim_build/pizza_dev/configs/pizza_dev/build/LoopControl/loopcontrol.c", line 957.1: 1506-343 (S) Redeclaration of lc_control_init differs from previous declaration on line 546 of "/sp6/usersissa/sists082/sim_build/pizza_dev/configs/pizza_dev/build/LoopControl/loopcontrol.c". "/sp6/usersissa/sists082/sim_build/pizza_dev/configs/pizza_dev/build/LoopControl/loopcontrol.c", line 957.1: 1506-377 (I) The type "const int* restrict const" of parameter 3 differs from the previous type "const int". "/sp6/usersissa/sists082/sim_build/pizza_dev/configs/pizza_dev/build/LoopControl/loopcontrol.c", line 977.1: 1506-334 (S) Identifier lc_control_finish has already been defined on line 801 of "/sp6/usersissa/sists082/sim_build/pizza_dev/configs/pizza_dev/build/LoopControl/loopcontrol.c". gmake[3]: *** [loopcontrol.c.o] Error 1 gmake[2]: *** [make.checked] Error 2 gmake[1]: *** [/sp6/usersissa/sists082/sim_build/pizza_dev/configs/pizza_dev/lib/libthorn_LoopControl.a] Error 2 gmake: *** [pizza_dev] Error 2 -------------------------- my option file looks like this: CC = xlc_r CXX = xlC_r F90 = xlf90_r F77 = xlf_r CPP = cpp FPP = /opt/freeware/bin/cpp CPPFLAGS = CFLAGS = -qspill=100000 -qsuppress=1501-245 -qlanglvl=stdc99 CXXFLAGS = -qspill=100000 -qsuppress=1501-245 -qrtti=all F77FLAGS = -qspill=100000 -qsuppress=1501-245 F90FLAGS = -qspill=100000 -qsuppress=1501-245 DEBUG = no OPTIMISE = yes OPTIMISE = yes CPP_OPTIMISE_FLAGS = -DNDEBUG FPP_OPTIMISE_FLAGS = -DNDEBUG C_OPTIMISE_FLAGS = -O2 -qmaxmem=-1 -qarch=pwr6 -qtune=pwr6 CXX_OPTIMISE_FLAGS = -O2 -qmaxmem=-1 -qarch=pwr6 -qtune=pwr6 F77_OPTIMISE_FLAGS = -O2 -qmaxmem=-1 -qarch=pwr6 -qtune=pwr6 F90_OPTIMISE_FLAGS = -O2 -qmaxmem=-1 -qarch=pwr6 -qtune=pwr6 WARN = yes HDF5 = yes HDF5_DIR = /cineca/prod/libraries/hdf5/1.8.4_ser/xl--10.1/ LIBSZ_DIR = /cineca/prod/libraries/szlib/2.1/xl--10.1/ GSL = yes GSL_DIR = /cineca/prod/libraries/gsl/1.9/xl--10.1/ MPI = NATIVE PTHREADS = no -------------------------- The compiler versions are: IBM XL C/C++ for AIX, V10.1 Version: 10.01.0000.0008 IBM XL Fortran for AIX, V12.1 Version: 12.01.0000.0009 -------------------------- The same code compiled just fine on Ubuntu 10.10 using the gnu compiler. Any ideas ? Cheers, Wolfgang.

3 5

Configuring Cactus
by Steven R. Brandt 18 Feb '11

18 Feb '11

I am trying to use an MPI build on my laptop, but I'm getting an error during configuration. checking for crypt in library crypt... yes checking for finite... yes checking for isnan... yes checking for isinf... yes checking for mkstemp... yes checking for va_copy... yes Configuring with MPI MPI selected, but no known MPI method - what is "/usr/lib64/mpich2" ? In my configuration (through SimFactory) I have this: MPI = CUSTOM MPI_INC_DIRS = /usr/include/mpich2-x86_64 MPI_LIB_DIRS = /usr/lib64/mpich2/lib MPI_LIBS = mpich opa pthread rt Cheers, Steve

1 1

Lorentz factor in HydroBase?
by Christian D. Ott 17 Feb '11

17 Feb '11

Hi all, I find myself frequently in a situation where I have to re-calculate the Lorentz factor in various analysis thorns that are ment to be independent of a given GR hydro implementation. Would it make sense to store the Lorentz factor as part of HydroBase? - Christian

1 0

Jump to page:

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Users