Commits January 2012

commits@lists.einsteintoolkit.org

7 participants
25 discussions

[svn:einsteintoolkit] www/ (Rev. 774)
by knarf＠cct.lsu.edu 09 Jan '12

09 Jan '12

User: knarf Date: 2012/01/09 10:53 AM Added: /seminars/2012_01_16/ index.php Modified: /seminars/ index.php Log: add ians talk Directory Changes: Directory: /svn:executable/ =========================== + * File Changes: Directory: /seminars/2012_01_16/ ================================ File [added]: index.php Delta lines: +45 -0 =================================================================== --- seminars/2012_01_16/index.php (rev 0) +++ seminars/2012_01_16/index.php 2012-01-09 16:53:45 UTC (rev 774) @@ -0,0 +1,45 @@ +<?php $title='A conservation law formulation of nonlinear elasticity in general relativity'; +include_once($_SERVER['DOCUMENT_ROOT'].'/global/header.php');?> + +<h3>Monday, January 16, 2012, 10am US central time</h3> +<p>The <a href="http://evo.caltech.edu/evoNext/koala.jnlp?meeting=eIe9esvvv8anaiatII"><i>Einstein Toolkit weekly</i></a> EVO meeting will be used for audio and slides. The talk is based +on this <a href="http://lanl.arxiv.org/abs/1107.2551">paper</a>. +</p> + + +<h3>Abstract</h3> +<p> +We present a practical framework for ideal hyperelasticity in numerical +relativity. For this purpose, we recast the formalism of Carter and Quintana as +a set of Eulerian conservation laws in an arbitrary 3+1 split of spacetime. The +resulting equations are presented as an extension of the standard Valencia +formalism for a perfect fluid, with additional terms in the stress-energy +tensor, plus a set of kinematic conservation laws that evolve a configuration +gradient. We prove that the equations can be made symmetric hyperbolic by +suitable constraint additions, at least in a neighbourhood of the unsheared +state. We discuss the Newtonian limit of our formalism and its relation to a +second formalism also used in Newtonian elasticity. We validate our framework +by numerically solving a set of Riemann problems in Minkowski spacetime, as +well as Newtonian ones from the literature. +</p> + +<h3>Presenter: Ian Hawke</h3> +<img src="Zachariah_Etienne.jpg" width="20%" align="right" alt="Zachariah Etienne"> +<p> + Lecturer in the <a href="http://www.maths.soton.ac.uk/applied"> Applied Mathematics</a> + department of the <a href="http://www.maths.soton.ac.uk">School of Mathematics</a> + at the <a href="http://www.soton.ac.uk">University of Southampton</a>. +<p> + +<p>Research interests: + My primary interests are numerical simulations of Einstein's equations + coupled to hydrodynamic matter. I am one of the authors of the + <a href="http://www.whiskycode.org">Whisky</a> code which is used within the + <a href="http://www.eu-network.org">EU Network</a>. + I also contribute to the binary black hole simulations of the + <a href="http://numrel.aei.mpg.de">AEI Numerical Relativity group</a>. + I am also interested in mesh refinement codes such as + <a href="http://www.carpetcode.org">Carpet</a>. +</p> + +<?php include_once($_SERVER['DOCUMENT_ROOT'].'/global/footer.php');?> Property changes on: seminars/2012_01_16/index.php ___________________________________________________________________ Directory: /seminars/ ===================== File [modified]: index.php Delta lines: +8 -0 =================================================================== --- seminars/index.php 2012-01-06 03:29:25 UTC (rev 773) +++ seminars/index.php 2012-01-09 16:53:45 UTC (rev 774) @@ -17,6 +17,14 @@ <td><a href="2011_11_28/EM_Gauge_ET.pdf">pdf</a></td> <td><a href="2011_11_28/recording/2011_11_28.ogv">ogv</a> / <a href="2011_11_28/recording/2011_11_28_evo.tar.xz">evo.tar.xz</a></td> </tr> + <tr> + <td>2012/01/16</td> + <td><a href="2012_01_16/">A conservation law formulation of nonlinear elasticity in general relativity</a></td> + <td>Ian Hawke</td> + <td>University of Southampton</td> + <td> </td> + <td> </td> + </tr> </table> <?php include_once($_SERVER['DOCUMENT_ROOT'].'/global/footer.php');?>

1 0

[svn:einsteintoolkit] WeylScal4/trunk/ (Rev. 101)
by schnetter＠cct.lsu.edu 06 Jan '12

06 Jan '12

User: eschnett Date: 2012/01/06 11:32 AM Modified: /trunk/ param.ccl, schedule.ccl /trunk/src/ psi4_calc_2nd.cc, psi4_calc_4th.cc, psi4_calc_Nth.cc, psis_calc_2nd.cc, psis_calc_4th.cc, psis_calc_Nth.cc Log: Regenerate code after changes to Kranc's vectorisation File Changes: Directory: /trunk/src/ ====================== File [modified]: psi4_calc_2nd.cc Delta lines: <Unavailable> File [modified]: psi4_calc_4th.cc Delta lines: <Unavailable> File [modified]: psi4_calc_Nth.cc Delta lines: <Unavailable> File [modified]: psis_calc_2nd.cc Delta lines: <Unavailable> File [modified]: psis_calc_4th.cc Delta lines: <Unavailable> File [modified]: psis_calc_Nth.cc Delta lines: <Unavailable> Directory: /trunk/ ================== File [modified]: param.ccl Delta lines: <Unavailable> File [modified]: schedule.ccl Delta lines: <Unavailable>

1 0

[svn:einsteintoolkit] www/about/members/ (Rev. 773)
by knarf＠cct.lsu.edu 05 Jan '12

05 Jan '12

User: knarf Date: 2012/01/05 09:29 PM Modified: /about/members/ send-registration.php Log: insert space File Changes: Directory: /about/members/ ========================== File [modified]: send-registration.php Delta lines: +1 -1 =================================================================== --- about/members/send-registration.php 2011-12-19 22:02:13 UTC (rev 772) +++ about/members/send-registration.php 2012-01-06 03:29:25 UTC (rev 773) @@ -24,7 +24,7 @@ $addtolist = "This person's email address is ".$email."; however, this person does not wish to be added to the ET users mailing list."; } -$message = "Einstein Toolkit maintainers: \n\n".$name.$frominstitution." has submitted a request to register with the Einstein Toolkit. ".$addtolist."\n\n Thanks,\n Einstein Toolkit Registration Bot\n"; +$message = "Einstein Toolkit maintainers: \n\n".$name." ".$frominstitution." has submitted a request to register with the Einstein Toolkit. ".$addtolist."\n\n Thanks,\n Einstein Toolkit Registration Bot\n"; /* PHP form validation: the script checks that the Email field contains a valid email address and the Subject field isn't empty. preg_match performs a regular expression match. It's a very powerful PHP function to validate form fields and other strings - see PHP manual for details. */ if (empty($name)) {

1 0

[svn:einsteintoolkit] GRHydro/trunk/ (Rev. 310)
by bcmsma＠astro.rit.edu 03 Jan '12

03 Jan '12

User: bmundim Date: 2012/01/03 04:13 PM Added: /trunk/src/ GRHydro_CalcBcom.F90 Modified: /trunk/ schedule.ccl /trunk/src/ make.code.defn Log: Add a routine to calculate mag pressure/comoving field * When Tmunu is not calculated, ie when no storage is selected, the comoving magnetic field/magnetic pressure wasn't calculated at all. This commit allows that now. File Changes: Directory: /trunk/src/ ====================== File [added]: GRHydro_CalcBcom.F90 Delta lines: +66 -0 =================================================================== --- trunk/src/GRHydro_CalcBcom.F90 (rev 0) +++ trunk/src/GRHydro_CalcBcom.F90 2012-01-03 22:13:04 UTC (rev 310) @@ -0,0 +1,66 @@ + /*@@ + @file GRHydro_CalcBcom.F90 + @date Jan 01, 2012 + @author Bruno Mundim + @histpry + Based on GRHydro_TmunuM.F90 file. + @desc + The calculation of the magnetic pressure and the comoving magnetic + field. + @enddesc + @@*/ +#include "cctk.h" +#include "cctk_Arguments.h" + +#define velx(i,j,k) vel(i,j,k,1) +#define vely(i,j,k) vel(i,j,k,2) +#define velz(i,j,k) vel(i,j,k,3) +#define Bvecx(i,j,k) Bvec(i,j,k,1) +#define Bvecy(i,j,k) Bvec(i,j,k,2) +#define Bvecz(i,j,k) Bvec(i,j,k,3) +#define bcomx(i,j,k) bcom(i,j,k,1) +#define bcomy(i,j,k) bcom(i,j,k,2) +#define bcomz(i,j,k) bcom(i,j,k,3) + + subroutine GRHydro_CalcBcom(CCTK_ARGUMENTS) + + implicit none + + DECLARE_CCTK_ARGUMENTS + + CCTK_REAL :: velxlow, velylow, velzlow + CCTK_REAL :: Bvecxlow,Bvecylow,Bveczlow + CCTK_REAL :: bdotv,b2,bxlow,bylow,bzlow,btlow,dum + CCTK_INT :: i,j,k + + if(stress_energy_state.ne.0) then + return + end if + + !$OMP PARALLEL DO PRIVATE(i,j,k,velxlow, velylow, velzlow,& + !$OMP Bvecxlow,Bvecylow,Bveczlow, bdotv,dum,b2,bxlow,bylow,bzlow) + + do k = 1, cctk_lsh(3) + do j = 1, cctk_lsh(2) + do i = 1, cctk_lsh(1) + + call calc_vlow_blow(gxx(i,j,k),gxy(i,j,k),gxz(i,j,k),& + gyy(i,j,k),gyz(i,j,k),gzz(i,j,k), & + velx(i,j,k),vely(i,j,k),velz(i,j,k),Bvecx(i,j,k),Bvecy(i,j,k),Bvecz(i,j,k), & + velxlow,velylow,velzlow,Bvecxlow,Bvecylow,Bveczlow, & + bdotv,b2,dum,dum,bxlow,bylow,bzlow) + + bcom_sq(i,j,k) = b2 + bcom0(i,j,k) = w_lorentz(i,j,k)*bdotv/alp(i,j,k) + bcomx(i,j,k) = Bvecx(i,j,k)/w_lorentz(i,j,k) + bcom0(i,j,k)*(alp(i,j,k)*velx(i,j,k)-betax(i,j,k)) + bcomy(i,j,k) = Bvecy(i,j,k)/w_lorentz(i,j,k) + bcom0(i,j,k)*(alp(i,j,k)*vely(i,j,k)-betay(i,j,k)) + bcomz(i,j,k) = Bvecz(i,j,k)/w_lorentz(i,j,k) + bcom0(i,j,k)*(alp(i,j,k)*velz(i,j,k)-betaz(i,j,k)) + + end do + end do + end do + !$OMP END PARALLEL DO + + return + + end subroutine GRHydro_CalcBcom File [modified]: make.code.defn Delta lines: +1 -0 =================================================================== --- trunk/src/make.code.defn 2011-12-07 22:33:23 UTC (rev 309) +++ trunk/src/make.code.defn 2012-01-03 22:13:04 UTC (rev 310) @@ -5,6 +5,7 @@ SRCS = Utils.F90 \ GRHydro_Boundaries.F90 \ + GRHydro_CalcBcom.F90 \ GRHydro_CalcUpdate.F90 \ GRHydro_Con2Prim.F90 \ GRHydro_DivergenceClean.F90 \ Directory: /trunk/ ================== File [modified]: schedule.ccl Delta lines: +10 -0 =================================================================== --- trunk/schedule.ccl 2011-12-07 22:33:23 UTC (rev 309) +++ trunk/schedule.ccl 2012-01-03 22:13:04 UTC (rev 310) @@ -1139,6 +1139,16 @@ { LANG: Fortran } "Compute the energy-momentum tensor - MHD version" +# bcom is already calculated in GRHydro_TmunuM. +# When Tmunu is not calculated, then the following +# routine is used to evaluate bcom, if requested. + if (calculate_bcom){ + schedule GRHydro_CalcBcom IN HydroBase_PostStep + { + LANG: Fortran + } "Compute comoving magnetic field, pressure, etc..." + } + } else { schedule GRHydro_Tmunu IN AddToTmunu {

1 0

New changesets in Carpet Mercurial
by Carpet mercurial repositories 02 Jan '12

02 Jan '12

changeset: 3411:8dee1dd409db tag: tip user: Erik Schnetter <schnetter(a)cct.lsu.edu> date: Mon Jan 02 17:10:03 2012 -0500 summary: CarpetInterp2: Don't define dim, use Carpet's dim instead

1 0

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Commits January 2012