eschnett/cactusamrex - Change XFLAXC6NZLEO2YH34GCKBU5E5RSQPK3JI36SV42DOERSFVJ3CPLAC

BaikalX: remove support for Tmunu at least for now

Created by Roland Haas on December 17, 2019

XFLAXC6NZLEO2YH34GCKBU5E5RSQPK3JI36SV42DOERSFVJ3CPLAC

Dependencies

In channels

master

main

Change contents

Deletion in BaikalX/doc/BaikalX.ipynb at line 160 [4.16]

B:BD[4.8900] → [4.8900:8997]

    "# Enable stress-energy terms?\n",
    "add_stress_energy_source_terms = False\n",
    "\n",

Deletion in BaikalX/doc/BaikalX.ipynb at line 247 [4.16]

B:BD[4.14961] → [4.14961:15694]

    "if add_stress_energy_source_terms == True:\n",
    "    T4UU = ixp.register_gridfunctions_for_single_rank2(\"AUXEVOL\",\"T4UU\",\"sym01\",DIM=4)\n",
    "    import BSSN.BSSN_stress_energy_source_terms as Bsest\n",
    "    Bsest.BSSN_source_terms_for_BSSN_RHSs(T4UU)\n",
    "    rhs.trK_rhs += Bsest.sourceterm_trK_rhs\n",
    "    for i in range(DIM):\n",
    "        # Needed for Gamma-driving shift RHSs:\n",
    "        rhs.Lambdabar_rhsU[i] += Bsest.sourceterm_Lambdabar_rhsU[i]\n",
    "        # Needed for BSSN RHSs:\n",
    "        rhs.lambda_rhsU[i]    += Bsest.sourceterm_lambda_rhsU[i]\n",
    "        for j in range(DIM):\n",
    "            rhs.a_rhsDD[i][j] += Bsest.sourceterm_a_rhsDD[i][j]\n",
    "\n",

Deletion in BaikalX/doc/BaikalX.ipynb at line 377 [4.16]

B:BD[4.25850] → [4.25850:26300]

    "    if add_stress_energy_source_terms == True:\n",
    "#         T4UU = gri.register_gridfunctions_for_single_rank2(\"AUXEVOL\",\"T4UU\",\"sym01\",DIM=4)\n",
    "        import BSSN.BSSN_stress_energy_source_terms as Bsest\n",
    "        Bsest.BSSN_source_terms_for_BSSN_constraints(T4UU)\n",
    "        bssncon.H += Bsest.sourceterm_H\n",
    "        for i in range(DIM):\n",
    "            bssncon.MU[i] += Bsest.sourceterm_MU[i]\n",

Deletion in BaikalX/doc/BaikalX.ipynb at line 603 [4.16]
∅:D[2.48] → [4.39837:39948]
∅:D[5.64] → [4.39837:39948]
B:BD[4.39837] → [4.39837:39948]
```
    "    if add_stress_energy_source_terms == True:\n",
    "        file.write(\"inherits:   TmunuBase\")\n",
```

Deletion in BaikalX/doc/BaikalX.ipynb at line 680 [4.16]

B:BD[4.47185] → [4.47185:47293]

B:BD[4.47293] → [5.272:372]

∅:D[5.372] → [4.47390:47486]

B:BD[4.47390] → [4.47390:47486]

    "\"\"\")\n",
    "    if add_stress_energy_source_terms == True:\n",
    "        file.write(\"\"\"\n",
    "schedule driver_BSSN_T4UU in ODESolvers_RHS as BaikalETK_T4UU before BaikalETK_BSSN_to_ADM\n",
    "{\n",
    "  LANG: C\n",
    "} \"MoL: Compute T4UU, needed for BSSN RHSs.\"\n",
    "\n",

Replacement in BaikalX/doc/BaikalX.ipynb at line 684 [4.16]
B:BD[4.47693] → [4.47693:47741]
```
    "\"\"\")\n",
    "    file.write(\"\"\"\n",
```
[4.47693]
[5.476]
```
    "\n",
```

Insertion in BaikalX/doc/BaikalX.ipynb at line 706 [4.16]

[4.49139]

    "  WRITES: BaikalETK::evol_variables(everywhere)\n",
    "  READS: ADMBase::metric(everywhere)\n",
    "  READS: ADMBase::curv(everywhere)\n",
    "  READS: ADMBASE::shift(everywhere)\n",
    "  READS: ADMBASE::lapse(everywhere)\n",
    "  READS: ADMBASE::dtlapse(everywhere)\n",
    "  # these run everywhere so thy cannot set dtshift since vet and beta differ by an advection term\n",

Deletion in BaikalX/doc/BaikalX.ipynb at line 716 [4.16]

B:BD[4.49274] → [4.49274:49382]

    "\"\"\")\n",
    "    if add_stress_energy_source_terms == True:\n",
    "        file.write(\"\"\"\n",

Deletion in BaikalX/doc/BaikalX.ipynb at line 719 [4.16]

B:BD[5.853] → [5.853:863]

∅:D[5.863] → [4.49382:49608]

B:BD[4.49382] → [4.49382:49608]

    "\n",
    "schedule driver_BSSN_T4UU in MoL_PseudoEvolution before BaikalETK_BSSN_constraints\n",
    "{\n",
    "  LANG: C\n",
    "  OPTIONS: Local\n",
    "} \"MoL_PseudoEvolution: Compute T4UU, needed for BSSN constraints\"\n",

Deletion in BaikalX/doc/BaikalX.ipynb at line 1277 [4.16]

B:BD[4.80169] → [4.80169:80598]

B:BD[4.80598] → [6.676:704]

∅:D[6.704] → [4.80624:82550]

B:BD[4.80624] → [4.80624:82550]

B:BD[4.82550] → [3.9498:9596]

∅:D[3.9596] → [4.82647:82732]

B:BD[4.82647] → [4.82647:82732]

B:BD[4.82732] → [3.9597:9643]

∅:D[3.9643] → [4.82982:82992]

B:BD[4.82982] → [4.82982:82992]

B:BD[4.82992] → [3.9644:9739]

∅:D[3.9739] → [4.83086:83266]

B:BD[4.83086] → [4.83086:83266]

B:BD[4.83321] → [4.83321:84061]

    "}\\n\"\"\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<a id='t4uu'></a>\n",
    "\n",
    "## Step 4.f: `driver_BSSN_T4UU()`: Compute $T^{\\mu\\nu}$ from `TmunuBase`'s $T_{\\mu\\nu}$ \\[Back to [top](#toc)\\]\n",
    "$$\\label{t4uu}$$\n",
    "\n",
    "Here we implement $T^{\\mu\\nu} = g^{\\mu \\delta} g^{\\nu \\gamma} T_{\\delta \\gamma}.$"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "if add_stress_energy_source_terms == True:\n",
    "    # Declare T4DD as a set of gridfunctions. These won't \n",
    "    #    actually appear in interface.ccl, as interface.ccl \n",
    "    #    was set above. Thus before calling the code output\n",
    "    #    by FD_outputC(), we'll have to set pointers\n",
    "    #    to the actual gridfunctions they reference.\n",
    "    #    (In this case the eTab's.)\n",
    "    T4DD = ixp.register_gridfunctions_for_single_rank2(\"AUXEVOL\",\"T4DD\",\"sym01\",DIM=4)\n",
    "    import BSSN.ADMBSSN_tofrom_4metric as AB4m\n",
    "    AB4m.g4UU_ito_BSSN_or_ADM(\"BSSN\")\n",
    "\n",
    "    T4UUraised = ixp.zerorank2(DIM=4)\n",
    "    for mu in range(4):\n",
    "        for nu in range(4):\n",
    "            for delta in range(4):\n",
    "                for gamma in range(4):\n",
    "                    T4UUraised[mu][nu] += AB4m.g4UU[mu][delta]*AB4m.g4UU[nu][gamma]*T4DD[delta][gamma]\n",
    "    \n",
    "    T4UU_expressions = [\n",
    "        lhrh(lhs=gri.gfaccess(\"in_gfs\",\"T4UU00\"),rhs=T4UUraised[0][0]),\n",
    "        lhrh(lhs=gri.gfaccess(\"in_gfs\",\"T4UU01\"),rhs=T4UUraised[0][1]),\n",
    "        lhrh(lhs=gri.gfaccess(\"in_gfs\",\"T4UU02\"),rhs=T4UUraised[0][2]),\n",
    "        lhrh(lhs=gri.gfaccess(\"in_gfs\",\"T4UU03\"),rhs=T4UUraised[0][3]),\n",
    "        lhrh(lhs=gri.gfaccess(\"in_gfs\",\"T4UU11\"),rhs=T4UUraised[1][1]),\n",
    "        lhrh(lhs=gri.gfaccess(\"in_gfs\",\"T4UU12\"),rhs=T4UUraised[1][2]),\n",
    "        lhrh(lhs=gri.gfaccess(\"in_gfs\",\"T4UU13\"),rhs=T4UUraised[1][3]),\n",
    "        lhrh(lhs=gri.gfaccess(\"in_gfs\",\"T4UU22\"),rhs=T4UUraised[2][2]),\n",
    "        lhrh(lhs=gri.gfaccess(\"in_gfs\",\"T4UU23\"),rhs=T4UUraised[2][3]),\n",
    "        lhrh(lhs=gri.gfaccess(\"in_gfs\",\"T4UU33\"),rhs=T4UUraised[3][3])]\n",
    "\n",
    "    outCparams = \"outCverbose=False,includebraces=False,preindent=2,SIMD_enable=False\"\n",
    "    T4UUstr = fin.FD_outputC(\"returnstring\",T4UU_expressions, outCparams)\n",
    "    T4UUstr_loop = loop_all(T4UUstr)\n",
    "\n",
    "    with open(append_to_make_code_defn_list(\"driver_BSSN_T4UU.cc\"), \"w\") as file:\n",
    "        file.write(\"\"\"\n",
    "#include <math.h>\n",
    "\n",
    "#include \"cctk.h\"\n",
    "#include \"cctk_Arguments.h\"\n",
    "#include \"cctk_Parameters.h\"\n",
    "\n",
    "void driver_BSSN_T4UU(CCTK_ARGUMENTS) {\n",
    "    DECLARE_CCTK_ARGUMENTS;\n",
    "    DECLARE_CCTK_PARAMETERS;\n",
    "\n",
    "    const CCTK_REAL *restrict T4DD00GF = eTtt;\n",
    "    const CCTK_REAL *restrict T4DD01GF = eTtx;\n",
    "    const CCTK_REAL *restrict T4DD02GF = eTty;\n",
    "    const CCTK_REAL *restrict T4DD03GF = eTtz;\n",
    "    const CCTK_REAL *restrict T4DD11GF = eTxx;\n",
    "    const CCTK_REAL *restrict T4DD12GF = eTxy;\n",
    "    const CCTK_REAL *restrict T4DD13GF = eTxz;\n",
    "    const CCTK_REAL *restrict T4DD22GF = eTyy;\n",
    "    const CCTK_REAL *restrict T4DD23GF = eTyz;\n",
    "    const CCTK_REAL *restrict T4DD33GF = eTzz;\n",
    "\"\"\"+T4UUstr_loop+\"\"\"\n",