Commits

commits@lists.einsteintoolkit.org

2 participants
6997 discussions

New changesets in Carpet Mercurial
by Carpet mercurial repositories 10 Sep '12

10 Sep '12

changeset: 3659:dc343eecda5f tag: tip user: Erik Schnetter <schnetter(a)gmail.com> date: Mon Sep 10 22:02:40 2012 -0400 summary: CarpetLib: Ensure that split/combined regions don't have a tree structure attached

1 0

New changesets in Carpet Mercurial
by Carpet mercurial repositories 08 Sep '12

08 Sep '12

changeset: 3656:2e33e718b5d5 user: Roland Haas <roland.haas(a)physics.gatech.edu> date: Fri Sep 07 23:17:05 2012 -0400 summary: Carpet: add aliased routines to obtain map, timelevel changeset: 3657:779e41b4def1 user: Roland Haas <roland.haas(a)physics.gatech.edu> date: Thu Sep 06 22:55:21 2012 -0400 summary: Carpet: add routines to record READS/WRITES at runtime changeset: 3658:dca1e61b47bd tag: tip user: Roland Haas <roland.haas(a)physics.gatech.edu> date: Fri Sep 07 23:18:34 2012 -0400 summary: Carpet: insert routines into Boundary group to capture boundary update

1 0

[svn:einsteintoolkit] incoming/Outflow/ (Rev. 73)
by rhaas＠tapir.caltech.edu 07 Sep '12

07 Sep '12

User: rhaas Date: 2012/09/07 08:43 AM Modified: /Outflow/ interface.ccl Log: Outflow: collect all extra surface projections in a single group File Changes: Directory: /Outflow/ ==================== File [modified]: interface.ccl Delta lines: +10 -20 =================================================================== --- Outflow/interface.ccl 2012-09-07 13:36:06 UTC (rev 72) +++ Outflow/interface.ccl 2012-09-07 13:43:28 UTC (rev 73) @@ -20,23 +20,13 @@ REAL w_lorentz_projected[num_detectors] TYPE=ARRAY DIM=2 SIZE=SphericalSurface::maxntheta,SphericalSurface::maxnphi DISTRIB=CONSTANT tags='checkpoint="no"' "2D (theta,phi) grid arrays for Lorentz factor" # extra variables projected onto the spherical surfaces -REAL surface_projection_0[num_detectors] TYPE=ARRAY DIM=2 SIZE=SphericalSurface::maxntheta,SphericalSurface::maxnphi DISTRIB=CONSTANT tags='checkpoint="no"' "2D (theta,phi) grid arrays for points on the spherical surfaces" -REAL surface_projection_1[num_detectors] TYPE=ARRAY DIM=2 SIZE=SphericalSurface::maxntheta,SphericalSurface::maxnphi DISTRIB=CONSTANT tags='checkpoint="no"' "2D (theta,phi) grid arrays for points on the spherical surfaces" -REAL surface_projection_2[num_detectors] TYPE=ARRAY DIM=2 SIZE=SphericalSurface::maxntheta,SphericalSurface::maxnphi DISTRIB=CONSTANT tags='checkpoint="no"' "2D (theta,phi) grid arrays for points on the spherical surfaces" -REAL surface_projection_3[num_detectors] TYPE=ARRAY DIM=2 SIZE=SphericalSurface::maxntheta,SphericalSurface::maxnphi DISTRIB=CONSTANT tags='checkpoint="no"' "2D (theta,phi) grid arrays for points on the spherical surfaces" -REAL surface_projection_4[num_detectors] TYPE=ARRAY DIM=2 SIZE=SphericalSurface::maxntheta,SphericalSurface::maxnphi DISTRIB=CONSTANT tags='checkpoint="no"' "2D (theta,phi) grid arrays for points on the spherical surfaces" -REAL surface_projection_5[num_detectors] TYPE=ARRAY DIM=2 SIZE=SphericalSurface::maxntheta,SphericalSurface::maxnphi DISTRIB=CONSTANT tags='checkpoint="no"' "2D (theta,phi) grid arrays for points on the spherical surfaces" -REAL surface_projection_6[num_detectors] TYPE=ARRAY DIM=2 SIZE=SphericalSurface::maxntheta,SphericalSurface::maxnphi DISTRIB=CONSTANT tags='checkpoint="no"' "2D (theta,phi) grid arrays for points on the spherical surfaces" -REAL surface_projection_7[num_detectors] TYPE=ARRAY DIM=2 SIZE=SphericalSurface::maxntheta,SphericalSurface::maxnphi DISTRIB=CONSTANT tags='checkpoint="no"' "2D (theta,phi) grid arrays for points on the spherical surfaces" -REAL surface_projection_8[num_detectors] TYPE=ARRAY DIM=2 SIZE=SphericalSurface::maxntheta,SphericalSurface::maxnphi DISTRIB=CONSTANT tags='checkpoint="no"' "2D (theta,phi) grid arrays for points on the spherical surfaces" -REAL surface_projection_9[num_detectors] TYPE=ARRAY DIM=2 SIZE=SphericalSurface::maxntheta,SphericalSurface::maxnphi DISTRIB=CONSTANT tags='checkpoint="no"' "2D (theta,phi) grid arrays for points on the spherical surfaces" -REAL surface_projection_10[num_detectors] TYPE=ARRAY DIM=2 SIZE=SphericalSurface::maxntheta,SphericalSurface::maxnphi DISTRIB=CONSTANT tags='checkpoint="no"' "2D (theta,phi) grid arrays for points on the spherical surfaces" -REAL surface_projection_11[num_detectors] TYPE=ARRAY DIM=2 SIZE=SphericalSurface::maxntheta,SphericalSurface::maxnphi DISTRIB=CONSTANT tags='checkpoint="no"' "2D (theta,phi) grid arrays for points on the spherical surfaces" -REAL surface_projection_12[num_detectors] TYPE=ARRAY DIM=2 SIZE=SphericalSurface::maxntheta,SphericalSurface::maxnphi DISTRIB=CONSTANT tags='checkpoint="no"' "2D (theta,phi) grid arrays for points on the spherical surfaces" -REAL surface_projection_13[num_detectors] TYPE=ARRAY DIM=2 SIZE=SphericalSurface::maxntheta,SphericalSurface::maxnphi DISTRIB=CONSTANT tags='checkpoint="no"' "2D (theta,phi) grid arrays for points on the spherical surfaces" -REAL surface_projection_14[num_detectors] TYPE=ARRAY DIM=2 SIZE=SphericalSurface::maxntheta,SphericalSurface::maxnphi DISTRIB=CONSTANT tags='checkpoint="no"' "2D (theta,phi) grid arrays for points on the spherical surfaces" -REAL surface_projection_15[num_detectors] TYPE=ARRAY DIM=2 SIZE=SphericalSurface::maxntheta,SphericalSurface::maxnphi DISTRIB=CONSTANT tags='checkpoint="no"' "2D (theta,phi) grid arrays for points on the spherical surfaces" -REAL surface_projection_16[num_detectors] TYPE=ARRAY DIM=2 SIZE=SphericalSurface::maxntheta,SphericalSurface::maxnphi DISTRIB=CONSTANT tags='checkpoint="no"' "2D (theta,phi) grid arrays for points on the spherical surfaces" -REAL surface_projection_17[num_detectors] TYPE=ARRAY DIM=2 SIZE=SphericalSurface::maxntheta,SphericalSurface::maxnphi DISTRIB=CONSTANT tags='checkpoint="no"' "2D (theta,phi) grid arrays for points on the spherical surfaces" -REAL surface_projection_18[num_detectors] TYPE=ARRAY DIM=2 SIZE=SphericalSurface::maxntheta,SphericalSurface::maxnphi DISTRIB=CONSTANT tags='checkpoint="no"' "2D (theta,phi) grid arrays for points on the spherical surfaces" -REAL surface_projection_19[num_detectors] TYPE=ARRAY DIM=2 SIZE=SphericalSurface::maxntheta,SphericalSurface::maxnphi DISTRIB=CONSTANT tags='checkpoint="no"' "2D (theta,phi) grid arrays for points on the spherical surfaces" +REAL surface_projections[num_detectors] TYPE=ARRAY DIM=2 SIZE=SphericalSurface::maxntheta,SphericalSurface::maxnphi DISTRIB=CONSTANT tags='checkpoint="no"' +{ +surface_projection_0, surface_projection_1, surface_projection_2, +surface_projection_3, surface_projection_4, surface_projection_5, +surface_projection_6, surface_projection_7, surface_projection_8, +surface_projection_9, surface_projection_10 surface_projection_11 +surface_projection_12 surface_projection_13 surface_projection_14 +surface_projection_15 surface_projection_16, surface_projection_17, +surface_projection_18, surface_projection_19 +} "2D (theta,phi) grid arrays for points on the spherical surfaces"

1 0

[svn:einsteintoolkit] incoming/Outflow/ (Rev. 72)
by rhaas＠tapir.caltech.edu 07 Sep '12

07 Sep '12

User: rhaas Date: 2012/09/07 08:36 AM Modified: /Outflow/ interface.ccl /Outflow/src/ outflow.c Log: Outflow: inherit from HydroBase, remove StaticConformal this causes early aborts if HydroBase is not present, StaticConformal is deprecated. Code changes are only to avoid variable name conflicts. File Changes: Directory: /Outflow/ ==================== File [modified]: interface.ccl Delta lines: +1 -1 =================================================================== --- Outflow/interface.ccl 2012-05-27 07:41:08 UTC (rev 71) +++ Outflow/interface.ccl 2012-09-07 13:36:06 UTC (rev 72) @@ -2,7 +2,7 @@ # $Header:$ implements: outflow -inherits: ADMBase SphericalSurface StaticConformal +inherits: ADMBase HydroBase SphericalSurface CCTK_INT FUNCTION IO_TruncateOutputFiles \ (CCTK_POINTER_TO_CONST IN GH) Directory: /Outflow/src/ ======================== File [modified]: outflow.c Delta lines: +13 -13 =================================================================== --- Outflow/src/outflow.c 2012-05-27 07:41:08 UTC (rev 71) +++ Outflow/src/outflow.c 2012-09-07 13:36:06 UTC (rev 72) @@ -87,7 +87,7 @@ static char *sanitize_filename(const char *fn); /* write results to disk */ static int Outflow_write_output(CCTK_ARGUMENTS, CCTK_INT det, CCTK_REAL flux, - CCTK_REAL w_lorentz, const CCTK_REAL *threshold_fluxes); + CCTK_REAL avg_w_lorentz, const CCTK_REAL *threshold_fluxes); static int Outflow_write_2d_output(CCTK_ARGUMENTS, const char *varname, CCTK_INT det, const CCTK_REAL *data_det, const CCTK_REAL *w_det, const CCTK_REAL *surfaceelement_det, @@ -288,7 +288,7 @@ } static int Outflow_write_output(CCTK_ARGUMENTS, CCTK_INT det, CCTK_REAL flux, - CCTK_REAL w_lorentz, const CCTK_REAL *threshold_fluxes) + CCTK_REAL avg_w_lorentz, const CCTK_REAL *threshold_fluxes) { DECLARE_CCTK_PARAMETERS; DECLARE_CCTK_ARGUMENTS; @@ -356,7 +356,7 @@ sprintf (format_str_real, "%%d\t%%%s\t%%%s\t%%%s", out_format,out_format,out_format); - fprintf(file, format_str_real, cctk_iteration, cctk_time, flux, w_lorentz); + fprintf(file, format_str_real, cctk_iteration, cctk_time, flux, avg_w_lorentz); sprintf (format_str_real, "\t%%%s", out_format); for(thresh = 0 ; thresh < num_thresholds ; thresh++) { fprintf(file,format_str_real,threshold_fluxes[thresh]); @@ -428,7 +428,7 @@ CCTK_REAL th,ph,ct,st,cp,sp,rp; CCTK_INT ntheta,nphi,npoints; // auxilliary variables used in constructing j - static CCTK_REAL *rho = NULL, *velx = NULL, *vely = NULL, *velz = NULL; + static CCTK_REAL *rho0 = NULL, *velx = NULL, *vely = NULL, *velz = NULL; static CCTK_REAL *beta1 = NULL, *beta2 = NULL, *beta3 = NULL, *alpha = NULL; static CCTK_REAL *g11 = NULL, *g12 = NULL, *g13 = NULL, *g22 = NULL; static CCTK_REAL *g23 = NULL, *g33 = NULL; @@ -472,7 +472,7 @@ ALLOCATE_TEMP(g22); ALLOCATE_TEMP(g23); ALLOCATE_TEMP(g33); - ALLOCATE_TEMP(rho); + ALLOCATE_TEMP(rho0); ALLOCATE_TEMP(velx); ALLOCATE_TEMP(vely); ALLOCATE_TEMP(velz); @@ -569,7 +569,7 @@ (void *) velx, (void *) vely, (void *) velz, - (void *) rho, + (void *) rho0, (void *) beta1, (void *) beta2, @@ -646,7 +646,7 @@ // compute current from primitive values for(int i = 0 ; i < interp_npoints ; i++) { - CCTK_REAL detg, dens, v2, w_lorentz; + CCTK_REAL detg, dens, v2, my_w_lorentz; detg = 2*g12[i]*g13[i]*g23[i] + g33[i]*(g11[i]*g22[i] - pow2(g12[i])) - g22[i]*pow2(g13[i]) - g11[i]*pow2(g23[i]); @@ -671,8 +671,8 @@ 2*g12[i]*velx[i]*vely[i] + 2*g13[i]*velx[i]*velz[i] + 2*g23[i]*vely[i]*velz[i]; - w_lorentz = sqrt(1. / (1. - v2)); - if( w_lorentz < 1. || v2 > 1 ) + my_w_lorentz = sqrt(1. / (1. - v2)); + if( my_w_lorentz < 1. || v2 > 1 ) { static CCTK_INT last_warned = -1; @@ -682,20 +682,20 @@ "%s: Unphysical Lorentz factor %15.6g, v2 = %15.6g for data " "g = [%15.6g,%15.6g,%15.6g,%15.6g,%15.6g,%15.6g] " "vel = [%15.6g,%15.6g,%15.6g] occured in iteration %d at location [%15.6g,%15.6g,%15.6g]", - __func__, w_lorentz,v2, g11[i],g12[i],g13[i],g22[i],g23[i],g33[i], + __func__, my_w_lorentz,v2, g11[i],g12[i],g13[i],g22[i],g23[i],g33[i], velx[i],vely[i],velz[i], cctk_iteration, det_x[i],det_y[i],det_z[i]); last_warned = cctk_iteration; } - w_lorentz = 1.; + my_w_lorentz = 1.; } - dens = sqrt(detg)*rho[i]*w_lorentz; + dens = sqrt(detg)*rho0[i]*my_w_lorentz; jx[i] = dens * (alpha[i]*velx[i] - beta1[i]); jy[i] = dens * (alpha[i]*vely[i] - beta2[i]); jz[i] = dens * (alpha[i]*velz[i] - beta3[i]); - w[i] = w_lorentz; + w[i] = my_w_lorentz; } return interp_npoints;

1 0

New changesets in Carpet Mercurial
by Carpet mercurial repositories 05 Sep '12

05 Sep '12

changeset: 3651:81389faf4401 parent: 3648:437389f2122e user: Roland Haas <roland.haas(a)physics.gatech.edu> date: Wed Sep 05 20:05:23 2012 -0400 summary: Carpet: add parameter to ignore variables during Requirements check changeset: 3652:c7027fc95d7f user: Roland Haas <roland.haas(a)physics.gatech.edu> date: Wed Sep 05 20:07:27 2012 -0400 summary: Carpet: match refinement level argument of Requirement::Recompose to changeset: 3653:8d95447fc475 user: Roland Haas <roland.haas(a)physics.gatech.edu> date: Wed Sep 05 20:09:21 2012 -0400 summary: Carpet: try to mimic enter_local_modes timelevel handling in changeset: 3654:c30137852913 user: Roland Haas <roland.haas(a)physics.gatech.edu> date: Wed Sep 05 20:12:12 2012 -0400 summary: Carpet: add output() routines to Requirements classes changeset: 3655:6f3e112bb184 tag: tip user: Roland Haas <roland.haas(a)physics.gatech.edu> date: Wed Sep 05 20:13:10 2012 -0400 summary: Carpet: output current state in Requirements::report_error

1 0

New changesets in Carpet Mercurial
by Carpet mercurial repositories 05 Sep '12

05 Sep '12

changeset: 3649:82e1c4b08aa7 branch: ET_2012_05 parent: 3556:b657ed971e64 user: Erik Schnetter <schnetter(a)gmail.com> date: Tue Aug 21 20:54:57 2012 -0400 summary: CarpetIOHDF5: Increase array size for last_output_iteration_slice changeset: 3650:32a217e2f8cc branch: ET_2012_05 tag: tip user: Erik Schnetter <schnetter(a)gmail.com> date: Tue Aug 21 20:57:25 2012 -0400 summary: CarpetIO*: Do not overwrite last_output_iteration after recovery

1 0

New changesets in Carpet Mercurial
by Carpet mercurial repositories 05 Sep '12

05 Sep '12

changeset: 3648:437389f2122e tag: tip user: Roland Haas <roland.haas(a)physics.gatech.edu> date: Wed Sep 05 10:13:52 2012 -0400 summary: CarpetIOHDF5: make index files sparse files with full datasets

1 0

New changesets in Carpet Mercurial
by Carpet mercurial repositories 05 Sep '12

05 Sep '12

changeset: 3647:e27bae01375c tag: tip user: Roland Haas <roland.haas(a)physics.gatech.edu> date: Wed Sep 05 00:44:46 2012 -0400 summary: Carpet: correct typos in Requirements handling

1 0

New changesets in Carpet Mercurial
by Carpet mercurial repositories 04 Sep '12

04 Sep '12

changeset: 3646:0a90dfca7b53 tag: tip user: Erik Schnetter <schnetter(a)gmail.com> date: Tue Sep 04 17:30:49 2012 -0400 summary: Carpet: Requirements: Allow specifying individual time levels

1 0

New changesets in Carpet Mercurial
by Carpet mercurial repositories 04 Sep '12

04 Sep '12

changeset: 3645:a88f0d03568c tag: tip user: Erik Schnetter <schnetter(a)gmail.com> date: Tue Sep 04 14:29:35 2012 -0400 summary: CarpetLib: Support prolongating faces

1 0

Jump to page:

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Commits