[Carpet] CarpetLib/src/vect.hh: Assertion `d>=0 && d<D' failed ??

Bruno Giacomazzo bgiacoma at sissa.it
Fri Jan 13 14:19:06 CET 2006


Hi,
 	I received the following error message from a run done with the 
development version of Carpet and Cactus:

cactus_whisky_carpet_devel: 
/home/people/bgiacoma/Cactus/arrangements/Carpet_devel/CarpetLib/src/vect.hh:168: 
Assertion `d>=0 && d<D' failed.

This run start from a checkpoint and the only difference with the previous 
run (the one that generated the checkpoint) is that now I activated 
CarpetIOBasic instead of CactusBase/IOBasic.

I have attached the par file, the Thornlist, the config-info, the stdout 
and stderr.

Thank you,
Bruno

-- 
Bruno Giacomazzo
PhD student - Astrophysics Sector
SISSA-ISAS (International School for Advanced Studies)
Via Beirut 4, 34014 Trieste, Italy

Tel.  : +39 0403787558
Fax   : +39 0403787528
email : bgiacoma at sissa.it

-------------------------------------------------
There are only 10 types of people in the world:
Those who understand binary, and those who don't
-------------------------------------------------
-------------- next part --------------
!REPOSITORY_LOCATION cvs.cactuscode.org
!REPOSITORY_NAME /cactusdevcvs
!REPOSITORY_TYPE pserver
!REPOSITORY_USER cvs_anon

CactusBase/MoL
CactusBase/Boundary
CactusBase/CartGrid3D
CactusBase/CoordBase
CactusBase/IOASCII
CactusBase/IOUtil
CactusBase/Time
CactusBase/LocalInterp
CactusBase/SymBase
CactusBase/LocalReduce
CactusEinstein/ADMBase
CactusEinstein/ADMMacros
CactusEinstein/ADMCoupling
CactusEinstein/ADMConstraints
CactusEinstein/CoordGauge
CactusEinstein/StaticConformal
CactusEinstein/SpaceMask
CactusEinstein/Extract
CactusEinstein/PsiKadelia
CactusElliptic/EllBase

Carpet_devel/Carpet
Carpet_devel/CarpetIOASCII
Carpet_devel/CarpetIOHDF5
Carpet_devel/CarpetIOScalar
Carpet_devel/CarpetIOBasic
Carpet_devel/CarpetInterp
Carpet_devel/CarpetLib
Carpet_devel/CarpetReduce
Carpet_devel/CarpetRegrid
Carpet_devel/CarpetSlab
Carpet_devel/CarpetTest


CactusUtils/NaNChecker
CactusExternal/jpeg6b
##CactusExternal/LAPACK
##CactusExternal/BLAS
CactusIO/IOJpeg
CactusConnect/HTTPD
CactusConnect/HTTPDExtra
CactusConnect/Socket


!REPOSITORY_NAME /arrangements
!REPOSITORY_USER cvs_anon
CactusEOS/EOS_Base
CactusEOS/EOS_Ideal_Fluid
TAT/Slab
TAT/TATelliptic
TAT/TATJacobi
##TAT/TGRtensor

!REPOSITORY_LOCATION cvs.aei.mpg.de
!REPOSITORY_NAME /numrelcvs
!REPOSITORY_TYPE pserver
!REPOSITORY_USER bgiacoma
Excision/LegoExcision

!REPOSITORY_LOCATION cvs.aei.mpg.de
!REPOSITORY_NAME /numrelcvs
!REPOSITORY_TYPE pserver
!REPOSITORY_USER bgiacoma
AEIThorns/AHFinderDirect
AEIThorns/BSSN_MoL
AEIThorns/AEILocalInterp
AEIThorns/SphericalSurface 
AEIThorns/Fortran
AEIDevelopment/RotatingSymmetry180
AEIDevelopment/WaveExtract
AEIDevelopment/ReflectionSymmetry
AEIDevelopment/ADMMass
##AEIDevelopment/IsolatedHorizon

!REPOSITORY_LOCATION cvs.eu-network.org
!REPOSITORY_NAME /EUNetwork
!REPOSITORY_TYPE pserver
!REPOSITORY_USER bgiacoma
EUHydro/EOS_Polytrope
EUHydro/Whisky_Init_Data
EUHydro/Whisky_RNSID
##EUHydro/Whisky_TOVSolver

!REPOSITORY_LOCATION cvs.aei.mpg.de
!REPOSITORY_NAME /numrelcvs
!REPOSITORY_MODULE HawkeCVS/EUHydro_Dev/Whisky
!REPOSITORY_TYPE pserver
!REPOSITORY_USER bgiacoma
!REPOSITORY_MODULE HawkeCVS/EUHydro_Dev/Whisky
EUHydro_Dev/Whisky
EUHydro_Dev/EOS_Table_Reader
EUHydro_Dev/Whisky_TATIVP
EUHydro_Dev/WhiskyModePower

Parma_Dev/Whisky_Hydro_Analysis

-------------- next part --------------
# CONFIGURATION  : whisky_carpet_devel
# CONFIG-DATE    : Thu Jan 12 20:37:26 2006 (GMT)
# CONFIG-HOST    : albert.pr.infn.it
# CONFIG-STATUS  : 0
# CONFIG-OPTIONS :
CC=icc
CXX=icc
F77=ifc
F90=ifc
FPP=${PERL} ${CCTK_HOME}/lib/sbin/cpp.pl
HDF5=yes
HDF5_DIR=/opt/intel/compiler70/hdf5
MPI=CUSTOM
MPI_INC_DIRS=/opt/intel/compiler70/lam-7.1/include
MPI_LIBS=pthread pbs nsl lammpio lamf77mpi mpi lam util dl
MPI_LIB_DIRS=/opt/intel/compiler70/lam-7.1/lib
THORNLIST=Whisky_Carpet_devel.th
THORNLIST_DIR=thornlists
-------------- next part --------------
ActiveThorns = "Carpet LocalReduce CarpetLib CarpetRegrid CarpetInterp CarpetReduce CarpetSlab cartgrid3d boundary admbase staticconformal spacemask admcoupling coordgauge admmacros whisky whisky_init_data eos_base eos_polytrope eos_ideal_fluid time mol nanchecker whisky_rnsid CarpetIOBasic CarpetIOScalar CarpetIOASCII ioutil admconstraints CarpetIOHDF5 bssn_mol coordbase aeilocalinterp legoexcision localinterp sphericalsurface symbase slab ahfinderdirect whiskymodepower Whisky_Hydro_Analysis ADMMass"

ActiveThorns = "waveextract"

ActiveThorns = "reflectionsymmetry"

##ActiveThorns = "TATJacobi TATElliptic Whisky_TATIVP"


IOScalar::outScalar_vars = " admbase::metric admbase::curv admbase::lapse admbase::shift adm_bssn::adm_bssn_phi adm_bssn::adm_bssn_metric adm_bssn::adm_bssn_curv adm_bssn::adm_bssn_k adm_bssn::adm_bssn_gamma admconstraints::momentum admconstraints::ham whisky::whisky_prim_vector_var whisky::whisky_prim_scalar_var whisky::whisky_cons_vector_var whisky::whisky_cons_scalar_var whisky::whisky_quadrupole whisky::whisky_masses whisky::whisky_angular_momentum whisky_hydro_analysis::hydro_quadrupole whisky_hydro_analysis::hydro_angular_momentum whisky_hydro_analysis::hydro_center_mass whisky_hydro_analysis::hydro_modes whisky::whisky_total_rest_mass whisky::whisky_masses_fractions"

IOBasic::outInfo_vars = " whisky::rho admbase::alp "


IOASCII::out1D_vars = " admbase::metric admbase::curv admbase::lapse admbase::shift adm_bssn::adm_bssn_phi adm_bssn::adm_bssn_metric{out_every=1} adm_bssn::adm_bssn_curv adm_bssn::adm_bssn_k adm_bssn::adm_bssn_gamma admconstraints::momentum admconstraints::ham whisky::whisky_prim_vector_var whisky::whisky_prim_scalar_var whisky::whisky_cons_vector_var whisky::whisky_cons_scalar_var whisky::whisky_quadrupole_temps"

IO::strict_io_parameter_check="yes"

IO::parfile_write = "generate"

#======================================
# IVP Solver (TATIVP)
#======================================
#whisky_tativp::solve_whisky_ivp		=	"yes"
#whisky_tativp::whisky_ivp_tolerance	=	1.e-10


#======================================
# Whisky with PPM parameters
#======================================
whisky::whisky_rho_central		=	9.02479e-01
whisky::riemann_solver			=	"Marquina"
whisky::whisky_eos_type			=	"General"
whisky::whisky_eos_table		=	"Ideal_Fluid"
whisky::recon_method	            	=       "ppm"
whisky::ppm_detect			=	"yes"
whisky::whisky_stencil          	=       3
whisky::bound		                =       "none"


#======================================
# MoL parameters with Whisky
#======================================
mol::ode_method				=	"icn"
adm_bssn::bssn_maxnumsandrvars          = 	0
whisky::whisky_maxnumsandrvars 		= 	0
adm_bssn::timelevels                    =	3
whisky::timelevels                      =	3


#======================================
# Whisky atmosphere parameters
#======================================
whisky::rho_rel_min                     =       3.e-7
whisky::initial_atmosphere_factor       =       0.9
whisky::initial_rho_abs_min             =       5.e-13

whisky::EoS_Change			=	"no"
whisky::EoS_Change_type                 =       "K"
whisky::initial_k                       =       1.0
whisky::c2p_reset_pressure		=	"yes"

#================================================
# Whisky: For computing masses within given radii
#================================================
whisky::compute_angular_momentum        =       "yes"
whisky::compute_masses                  =       "yes"

whisky::masses_nr                       =       4
whisky::reference_radius[1]             =       1.5
whisky::reference_radius[2]             =       0.5
whisky::reference_radius[3]             =       0.2
whisky::reference_radius[4]             =       0.05

#================================================
# Whisky_Hydro_Analysis
#================================================
Whisky_Hydro_Analysis::compute_hydro_quadrupole = "yes"
Whisky_Hydro_Analysis::mass_radius_1    =       0.05
Whisky_Hydro_Analysis::mass_radius_2    =       1.0
Whisky_Hydro_Analysis::mass_radius_3    =       1.5
Whisky_Hydro_Analysis::modes_radius_min =       0.0
Whisky_Hydro_Analysis::modes_radius_max =       1.5

#======================================
# RNSID initial data with the EOS changed
#======================================
admbase::initial_data			=	"rnsid"
admbase::initial_lapse			=	"rnsid"
admbase::initial_shift			=	"rnsid"
admbase::metric_type			=	"physical"
rnsid::rotation_type			=	"diff"
rnsid::A_diff				=	1.4
rnsid::axes_ratio			=	0.52
rnsid::accuracy				=	1e-6
rnsid::cf                               =       0.4
rnsid::consistent_boundaries		=	"yes"
rnsid::MDIV                             =       301
rnsid::SDIV                             =       601

rnsid::RNS_Gamma			=	3.0
rnsid::RNS_K				=	1.0
eos_polytrope::eos_gamma		=	3.0
eos_polytrope::eos_k			=	0.98
eos_ideal_fluid::eos_ideal_fluid_gamma	=	3.0

#======================================
# Whisky output
#======================================
IO::out_dir             = "data"
IOScalar::outScalar_dir = "data/Scalar"
IOASCII::out1D_dir      = "data/1d"
IO::out_unchunked    = "no"
IO::out_mode         = "proc"

# 1D output
IOScalar::outScalar_every 		=	8
IOASCII::out1D_every	     		= 	16
IOBasic::outInfo_every   		= 	2

#====================================================================
# Whisky output
#====================================================================
# 3D output
IOHDF5::out3D_every	     		= 	32
IOHDF5::out_dir      			= 	"data/hdf5"
# save disk space by using single precision and downsampling
IO::out_single_precision	      	= 	"yes"
IOHDF5::out3D_vars	      		=	"whisky::whisky_prim_vector_var whisky::rho whisky::eps"


#======================================
# BSSN as run with Whisky
#======================================
admmacros::spatial_order		=	2
admbase::evolution_method		=	"adm_bssn"
admbase::lapse_evolution_method		=	"1+log"
admbase::shift_evolution_method		=	"gamma2"
adm_bssn::lapsesource			=	"modified"
adm_bssn::bound				=	"radiative"
boundary::radpower			=	2
adm_bssn::harmonic_f			=	2.0
adm_bssn::ShiftGammaCoeff		=	0.75
adm_bssn::BetaDriver			=	3


#==============================================================================
# GRID used by Whisky
#==============================================================================
Carpet::domain_from_coordbase		= "yes"
Carpet::verbose				= "yes"
#Carpet::enable_all_storage        	= "no"

time::dtfac                             = 0.375

grid::type				= "coordbase"

grid::avoid_origin			= "no"

coordbase::boundary_size_z_lower 	= 3 #ghost size
coordbase::boundary_shiftout_z_lower 	= 1 #due to avoid_origin="no"

carpetregrid::smart_outer_boundaries 	= "yes"

reflectionsymmetry::reflection_z	= "yes"
reflectionsymmetry::avoid_origin_z	= "no"

coordbase::xmin =   -1.60
coordbase::xmax =    1.60
coordbase::ymin =   -1.60
coordbase::ymax =    1.60
coordbase::zmin =    0.0 
coordbase::zmax =    1.92

coordbase::dx = 0.04
coordbase::dy = 0.04
coordbase::dz = 0.012

carpetregrid::verbose			= "yes"
carpetregrid::refined_regions		= "manual-coordinate-list"
carpetregrid::coordinates		= "[
[ ([-0.8,-0.8,0.0 ]:[0.8,0.8,1.5 ]:[0.02,0.02,0.006])]]"

#level1 80*80*250
#level0 80*80*160

carpet::max_refinement_levels		= 2

carpetregrid::refinement_levels		= 2 

carpetlib::save_memory_during_regridding = "yes"
#carpetlib::max_allowed_memory_MB	 = 750
CarpetLib::print_memstats_every          = 16

carpetregrid::keep_same_grid_structure	= "yes"

carpet::prolongation_order_space	= 3
carpet::prolongation_order_time		= 2
carpet::buffer_width			= 6

carpet::init_3_timelevels		= "no" #"yes"

driver::ghost_size                      = 3
cactus::cctk_itlast                     = 2000000


#==============================================================================
# AHFinderDirect parameters for a single star collapse
#==============================================================================
ahfinderdirect::N_horizons				= 3
ahfinderdirect::find_every				= 32

ahfinderdirect::max_allowable_horizon_radius[1]		= 0.7
ahfinderdirect::max_allowable_horizon_radius[2]		= 0.7
ahfinderdirect::max_allowable_horizon_radius[3]		= 0.7

ahfinderdirect::initial_guess__coord_sphere__radius[1]	= 0.01
ahfinderdirect::initial_guess__coord_sphere__radius[2]	= 0.1
ahfinderdirect::initial_guess__coord_sphere__radius[3]	= 0.3


#==============================================================================
# Excise hydro and spacetime
#==============================================================================
spacemask::use_mask			= 	"true"
adm_bssn::excise			=	"yes"
adm_bssn::excisionType			=	"lego"
ahfinderdirect::set_mask_for_all_horizons	= 	"true"
ahfinderdirect::set_old_style_mask	= 	"true"
ahfinderdirect::set_new_style_mask	= 	"true"
ahfinderdirect::mask_radius_multiplier	=	0.6
ahfinderdirect::mask_radius_offset	= 	0.0
ahfinderdirect::mask_buffer_thickness	= 	0.0
ahfinderdirect::new_style_mask_bitfield_name	= 	"Hydro_Excision"
ahfinderdirect::new_style_mask_inside_value	= 	"excised"
ahfinderdirect::new_style_mask_buffer_value	= 	"normal"
ahfinderdirect::new_style_mask_outside_value	= 	"normal"


# Checkpoint parameters
iohdf5::checkpoint   = "yes"

io::recover          = "auto"
io::recover_dir      = "/home/people/bgiacoma/C-N05-A14.4492.CHKPT"

io::checkpoint_every = 1024
IO::checkpoint_dir = "data/CHECKPOINT"
io::checkpoint_keep=1

# wave extraction
waveextract::active			= "yes"
waveextract::out_every			= 32

waveextract::interpolation_operator 	= "Lagrange polynomial interpolation"
waveextract::interpolation_order 	= 2
waveextract::maximum_detector_number 	= 5
waveextract::switch_output_format 	= 20
waveextract::rsch2_computation       	= "average Schwarzschild metric"
waveextract::verbose       		= 1
waveextract::l_mode        		= 5
waveextract::m_mode        		= 5

waveextract::detector_radius[0]     	= 0.5
waveextract::detector_radius[1]     	= 0.8
waveextract::detector_radius[2]     	= 1.1
waveextract::detector_radius[3]     	= 1.4
waveextract::detector_radius[4]     	= 1.5

#============================================================
# Mode power
#============================================================
whiskymodepower::n_rings                = 9
whiskymodepower::ring_radius[0]         = 0.01
whiskymodepower::ring_radius[1]         = 0.02
whiskymodepower::ring_radius[2]         = 0.1
whiskymodepower::ring_radius[3]         = 0.2
whiskymodepower::ring_radius[4]         = 0.4
whiskymodepower::ring_radius[5]         = 0.6
whiskymodepower::ring_radius[6]         = 0.8
whiskymodepower::ring_radius[7]         = 1.0
whiskymodepower::ring_radius[8]         = 1.5
#whiskymodepower::ring_radius[9]         = 2.0

whiskymodepower::n_m_modes              = 6
whiskymodepower::ring_rin               = 1
whiskymodepower::ring_rout              = 7
whiskymodepower::n_angles               = 80

ADMMass::ADMMass_use_surface_distance_as_volume_radius = "no"
ADMMass::ADMMass_use_all_volume_as_volume_radius = "yes"
ADMMass::ADMMass_distance_from_grid_boundary[0] = 0.024
CarpetIOASCII::out0D_every = 8
CarpetIOASCII::out0D_vars = "admmass::ADMMass_Masses"
-------------- next part --------------
A non-text attachment was scrubbed...
Name: C-N05-A14.e4497.bz2
Type: application/x-bzip2
Size: 729 bytes
Desc: 
Url : /archives/developers/attachments/20060113/2c6dc533/attachment.bin 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: C-N05-A14.o4497.bz2
Type: application/x-bzip2
Size: 12189 bytes
Desc: 
Url : /archives/developers/attachments/20060113/2c6dc533/attachment-0001.bin 


More information about the developers mailing list