2017 MHD Workshop - Revision history

192.16.204.70: /* Knights Landing on Stampede */

2017-03-29T07:10:43Z

‎Knights Landing on Stampede

Rhaas: strip down MHD workshopt page to remove non-essentials

2017-02-13T11:38:34Z

strip down MHD workshopt page to remove non-essentials

65.122.108.210 at 23:55, 5 February 2017

2017-02-05T23:55:26Z

Rhaas: update ENOG2 progress

2017-02-05T15:18:46Z

update ENOG2 progress

64.134.70.19: /* New Con2Prim Interface for GRHydro */

2017-02-04T12:55:50Z

‎New Con2Prim Interface for GRHydro

64.134.70.19: /* New Con2Prim Interface for GRHydro */

2017-02-04T12:54:18Z

‎New Con2Prim Interface for GRHydro

64.134.70.19: /* FunHPC parallelization progress */

2017-02-04T12:44:39Z

‎FunHPC parallelization progress

Eschnett: /* FunHPC parallelization progress */

2017-02-03T11:07:45Z

‎FunHPC parallelization progress

Eschnett: /* flux-CT implementation progress */

2017-02-03T11:04:39Z

‎flux-CT implementation progress

206.217.70.162 at 10:59, 3 February 2017

2017-02-03T10:59:13Z

@@ Line 39: / Line 39: @@
 * Stampede manual: https://portal.tacc.utexas.edu/user-guides/stampede
 * Login node: login-knl1.stampede.tacc.utexas.edu
-* option list (beta): https://dl.dropboxusercontent.com/u/17077801/stampede-knl.cfg
+* option list, submission scripts, and environment settings for both icc and gcc are available for batchtools: https://bitbucket.org/dradice/batchtools
 == flux-CT implementation progress ==

←Older revision		Revision as of 11:38, 13 February 2017
Line 1:		Line 1:
−	~~Workshop in New York City and at Columbia University.~~
−
−	~~== Chelsea Venue ==~~
−
−	* On the first afternoon, we'll meet off-campus to get the workshop started.
−	* We will start at 14:00.
−	* We have a conference room for 8 booked under Christian D. Ott, California Institute of Technology at <br>Select Office Suites Chelsea<br>116 W 23rd St, New York, NY 10011<br>https://goo.gl/maps/WDBNQpwadfE2 <br>From Columbia, the best way to get there is to take Subway Line 1 to 23rd Street Station.
−
−	~~== Columbia Venue==~~
−
−	* The workshop will take place in the Conference Room of the Center for Theoretical Physics at Columbia. The Conference Room (Room 907) is located on the 9th floor of the Physics building (Pupin Laboratories).
−	* Street address: 538 West 120th Street, New York, NY 10027
−	* Note that the main entrance is on the Campus Level (5th floor, south side) as indicated here:
−	~~[https://www.google.com/maps/place/40°48'35.9%22N+73°57'41.2%22W/@40.809983,-73.9619872,19z/data=!3m1!4b1!4m5!3m4!1s0x0:0x0!8m2!3d40.809982!4d-73.96144 map]~~
−	* You can enter Campus at 116th/Broadway. There is also the possibility of entering Pupin Hall from 120th/Broadway. In this case, enter through the North West Corner Building and take the stairs/escalators up to the Campus level, exit the building and enter Pupin through the main entrance.
−	* The conference room has a board and a big presentation screen -> There is the possibility of giving a black board talk or a slide presentation, or a combination of both.
−
	==Time and Agenda==		==Time and Agenda==

←Older revision		Revision as of 23:55, 5 February 2017
Line 99:		Line 99:

	As an example for an external Con2Prim Thorn that provides GRHydro_Con2PrimNG, Christian wrote Zelmani_Con2PrimNG_BBHGas, which implements a variant of the vanilla GRHydro Con2Prim routine and makes it available to GRHydro via the aliased function interface. Note that since the backend C2P routine is F90, an additional wrapper is required that massages the passed arguments in such a way so that the backend F90 routine can be called with CCTK_ARGUMENTS. In C or C++ such a wrapper would not be necessary.		As an example for an external Con2Prim Thorn that provides GRHydro_Con2PrimNG, Christian wrote Zelmani_Con2PrimNG_BBHGas, which implements a variant of the vanilla GRHydro Con2Prim routine and makes it available to GRHydro via the aliased function interface. Note that since the backend C2P routine is F90, an additional wrapper is required that massages the passed arguments in such a way so that the backend F90 routine can be called with CCTK_ARGUMENTS. In C or C++ such a wrapper would not be necessary.
		+
		+
		+	== Implementation of vector potential into GRHydro ==
		+
		+	Quick summary:
		+
		+	* Implemented staggered vector potential formalism according to Etienne et al. 2012 and the IllinoisGRMHD code
		+	* This is contained in the Zelmani branch vector_potential
		+	* Essentially isolated, rewrote or modified relevant code parts from IllinoisGRMHD and included them into GRHydro
		+	* Removed old (dead) vector potential parts in GRHydro and made several changes to GRHydro core routines to make it work with staggered Avec and Bvec fields
		+	* Currently compiles and only works with the cxx files/version of GRHydro
		+	* Added initial data routine to GRHydro_InitData to initialize magnetic field configurations for the staggered Avec and Bvec fields
		+
		+	Next steps:
		+	* First test runs for debugging
		+	* Modify GRHydro's reconstruction routines to also reconstruct staggered B and velocity fields for computation of RHS of staggered Avec evolution (currently using a modified routine from IllinoisGRMHD)
		+	* Modification of GRHydro's reconstruction routines to be coordinated with Philipp's work on new staggered CT scheme, as similar infrastructure required

@@ Line 74: / Line 74: @@
 * Continued to make FunHPC work, as alternative to OpenMP
 * Improve performance by taking cache alignment into account
-* * Many internal modifications to Cactus and Carpet
+** Many internal modifications to Cactus and Carpet
 * Benchmarking
-* * Need to measure performance details automatically, make available to users
+** Need to measure performance details automatically, make available to users
 == New Con2Prim Interface for GRHydro ==

@@ Line 89: / Line 89: @@
 * Introduced new C2P wrapper routine that is called from the schedule. This is in GRHydro/src/GRHydro_Con2Prim.cc. This routine provides full backward compatibility and calls previous C and F90 C2P routines if GRHydro_use_Con2PrimNG set to 'no' and no function GRHydro_Con2PrimNG is provided through the aliased function interface.
 * Cleaned up GRHydro/src/GRHydro_Con2Prim.F90 and removed code specific to the BBH+Gas project. The implementation in GRHydro_Con2Prim.F90 should not be back to the generic one in the ET version of GRHydro.
 As an example for an external Con2Prim Thorn that provides GRHydro_Con2PrimNG, Christian wrote Zelmani_Con2PrimNG_BBHGas, which implements a variant of the vanilla GRHydro Con2Prim routine and makes it available to GRHydro via the aliased function interface. Note that since the backend C2P routine is F90, an additional wrapper is required that massages the passed arguments in such a way so that the backend F90 routine can be called with CCTK_ARGUMENTS. In C or C++ such a wrapper would not be necessary.

←Older revision		Revision as of 12:54, 4 February 2017
Line 79:		Line 79:

	== New Con2Prim Interface for GRHydro ==		== New Con2Prim Interface for GRHydro ==
		+
		+	The goal of the Next Gen Con2Prim (Con2PrimNG) interface for GRHydro is to allow for customized Con2Prim routines provided by the user in a separate thorn. To enable this, GRHydro can now call an aliased function for Con2Prim provided by a dedicated user-provided Thorn.
		+
		+	The changes to GRHydro are as follows:
		+
		+	* In interface.ccl, introduced relevant ccl code so that GRHydro can use a function called GRHydro_Con2PrimNG for Con2Prim.
		+	* In param.ccl, introduced boolean 'GRHydro_use_Con2PrimNG'. It defaults to 'no', but if set to 'yes', GRHydro will attempt to call the aliased function GRHydro_Con2PrimNG for Con2Prim. Also removed two parameters that were used for the BBH+Gas project and that were recently introduced and have not been pushed upstream.
		+	* In schedule.ccl, commented out all calls to Con2Prim routines and replaced them with call to a new wrapper function defined in src/GRHydro_Con2Prim.cc.
		+	* Introduced new C2P wrapper routine that is called from the schedule. This is in GRHydro/src/GRHydro_Con2Prim.cc. This routine provides full backward compatibility and calls previous C and F90 C2P routines if GRHydro_use_Con2PrimNG set to 'no' and no function GRHydro_Con2PrimNG is provided through the aliased function interface.
		+	* Cleaned up GRHydro/src/GRHydro_Con2Prim.F90 and removed code specific to the BBH+Gas project. The implementation in GRHydro_Con2Prim.F90 should not be back to the generic one in the ET version of GRHydro.
		+
		+	As an example for an external Con2Prim Thorn that provides GRHydro_Con2PrimNG, Christian wrote Zelmani_Con2PrimNG_BBHGas, which implements a variant of the vanilla GRHydro Con2Prim routine and makes it available to GRHydro via the aliased function interface. Note that since the backend C2P routine is F90, an additional wrapper is required that massages the passed arguments in such a way so that the backend F90 routine can be called with CCTK_ARGUMENTS. In C or C++ such a wrapper would not be necessary.

←Older revision		Revision as of 12:44, 4 February 2017
Line 77:		Line 77:
	* Benchmarking		* Benchmarking
	* * Need to measure performance details automatically, make available to users		* * Need to measure performance details automatically, make available to users
		+
		+	== New Con2Prim Interface for GRHydro ==

←Older revision		Revision as of 11:07, 3 February 2017
Line 71:		Line 71:

	==FunHPC parallelization progress==		==FunHPC parallelization progress==
		+
		+	* Continued to make FunHPC work, as alternative to OpenMP
		+	* Improve performance by taking cache alignment into account
		+	* * Many internal modifications to Cactus and Carpet
		+	* Benchmarking
		+	* * Need to measure performance details automatically, make available to users

←Older revision		Revision as of 11:04, 3 February 2017
Line 69:		Line 69:
	* Started implementation for that and bookkeeping is going to be tedious but otherwise implementation should be straightforward		* Started implementation for that and bookkeeping is going to be tedious but otherwise implementation should be straightforward
	* Need to come up with good test cases to make sure this is implemented correctly		* Need to come up with good test cases to make sure this is implemented correctly
		+
		+	==FunHPC parallelization progress==

←Older revision		Revision as of 10:59, 3 February 2017
Line 61:		Line 61:

	David Radice's simfactory replacement: https://bitbucket.org/dradice/batchtools		David Radice's simfactory replacement: https://bitbucket.org/dradice/batchtools
		+
		+	== flux-CT implementation progress ==
		+
		+	* Identified which equations to solve and that there is indeed a Riemann problem to solve for the electric field
		+	* This is not a Riemann problem for the state vector (which is the magnetic field), but the electric field is used to update the magnetic field
		+	* The Riemann problem for the electric field requires reconstruction of the velocity and magnetic field to the cell edge, requiring a double reconstruction for the velocity
		+	* Started implementation for that and bookkeeping is going to be tedious but otherwise implementation should be straightforward
		+	* Need to come up with good test cases to make sure this is implemented correctly