eschnett/cactusamrex - Change H6GIS2UO4G25J7MFBEEC5TV4LVM66TRPJFHXXVXAMJTZIDVAWA7QC

BaikalETK: make compile as CarpetX using thorn

Created by Roland Haas on December 16, 2019

H6GIS2UO4G25J7MFBEEC5TV4LVM66TRPJFHXXVXAMJTZIDVAWA7QC

Dependencies

In channels

main

master

Change contents

Insertion in BaikalX/doc/BaikalX.ipynb at line 89 [3.16]

[3.7256]

   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Step 0 define my own looping constructs\n",
    "def loop_int(kernel):\n",
    "    return \"\"\"Loop::loop_int<0,0,0>(cctkGH, [&](const Loop::PointDesc &p){\n",
    "    const int i0 = p.i;\n",
    "    const int i1 = p.j;\n",
    "    const int i2 = p.k;\n",
    "\"\"\"+kernel+\"\"\"\n",
    "});\\n\"\"\"\n",
    "def loop_all(kernel):\n",
    "    return \"\"\"Loop::loop_all<0,0,0>(cctkGH, [&](const Loop::PointDesc &p){\n",
    "    const int i0 = p.i;\n",
    "    const int i1 = p.j;\n",
    "    const int i2 = p.k;\n",
    "\"\"\"+kernel+\"\"\"\n",
    "});\\n\"\"\""
   ]
  },
  {

Deletion in BaikalX/doc/BaikalX.ipynb at line 178 [3.16]

B:BD[3.10546] → [3.10546:10778]

    "# Copy SIMD/SIMD_intrinsics.h to $outdir/SIMD/SIMD_intrinsics.h\n",
    "cmd.mkdir(os.path.join(outdir,\"SIMD\"))\n",
    "shutil.copy(os.path.join(\"SIMD/\")+\"SIMD_intrinsics.h\",os.path.join(outdir,\"SIMD/\"))\n",
    "\n",

Replacement in BaikalX/doc/BaikalX.ipynb at line 180 [3.16]

B:BD[3.10940] → [3.10940:11007]

    "par.set_parval_from_str(\"grid::GridFuncMemAccess\",\"ETK\")"

[3.10940]

[3.11007]

    "par.set_parval_from_str(\"grid::GridFuncMemAccess\",\"CarpetX\")"

Replacement in BaikalX/doc/BaikalX.ipynb at line 330 [3.16]

B:BD[3.21315] → [3.21315:21442]

    "    BSSN_RHSs_string = fin.FD_outputC(\"returnstring\",BSSN_evol_rhss, params=\"outCverbose=False,SIMD_enable=True\",\n",

[3.21315]

[3.21442]

    "    BSSN_RHSs_string = fin.FD_outputC(\"returnstring\",BSSN_evol_rhss, params=\"outCverbose=False,SIMD_enable=False\",\n",

Replacement in BaikalX/doc/BaikalX.ipynb at line 334 [3.16]

B:BD[3.21600] → [3.21600:22361]

    "        file.write(lp.loop([\"i2\",\"i1\",\"i0\"],[\"cctk_nghostzones[2]\",\"cctk_nghostzones[1]\",\"cctk_nghostzones[0]\"],\n",
    "       [\"cctk_lsh[2]-cctk_nghostzones[2]\",\"cctk_lsh[1]-cctk_nghostzones[1]\",\"cctk_lsh[0]-cctk_nghostzones[0]\"],\n",
    "                           [\"1\",\"1\",\"SIMD_width\"],\n",
    "                            [\"#pragma omp parallel for\",\n",
    "                                         \"#include \\\"rfm_files/rfm_struct__SIMD_outer_read2.h\\\"\",\n",
    "                                         \"#include \\\"rfm_files/rfm_struct__SIMD_outer_read1.h\\\"\"],\"\",\n",
    "                                         \"#include \\\"rfm_files/rfm_struct__SIMD_inner_read0.h\\\"\\n\"+BSSN_RHSs_string))\n",

[3.21600]

[3.22361]

    "        file.write(loop_int(BSSN_RHSs_string))\n",

Replacement in BaikalX/doc/BaikalX.ipynb at line 348 [3.16]

B:BD[3.23454] → [3.23454:23545]

    "                                   params=\"outCverbose=False,SIMD_enable=True\")\n",

[3.23454]

[3.23545]

    "                                   params=\"outCverbose=False,SIMD_enable=False\")\n",

Replacement in BaikalX/doc/BaikalX.ipynb at line 350 [3.16]

B:BD[3.23623] → [3.23623:24380]

    "        file.write(lp.loop([\"i2\",\"i1\",\"i0\"],[\"cctk_nghostzones[2]\",\"cctk_nghostzones[1]\",\"cctk_nghostzones[0]\"],\n",
    "       [\"cctk_lsh[2]-cctk_nghostzones[2]\",\"cctk_lsh[1]-cctk_nghostzones[1]\",\"cctk_lsh[0]-cctk_nghostzones[0]\"],\n",
    "                           [\"1\",\"1\",\"SIMD_width\"],\n",
    "                            [\"#pragma omp parallel for\",\n",
    "                                         \"#include \\\"rfm_files/rfm_struct__SIMD_outer_read2.h\\\"\",\n",
    "                                         \"#include \\\"rfm_files/rfm_struct__SIMD_outer_read1.h\\\"\"],\"\",\n",
    "                                         \"#include \\\"rfm_files/rfm_struct__SIMD_inner_read0.h\\\"\\n\"+Ricci_string))\n",

[3.23623]

[3.24380]

    "        file.write(loop_int(Ricci_string))\n",

Replacement in BaikalX/doc/BaikalX.ipynb at line 400 [3.16]

B:BD[3.27150] → [3.27150:27532]

    "        file.write(lp.loop([\"i2\",\"i1\",\"i0\"],[\"cctk_nghostzones[2]\",\"cctk_nghostzones[1]\",\"cctk_nghostzones[0]\"],\n",
    "       [\"cctk_lsh[2]-cctk_nghostzones[2]\",\"cctk_lsh[1]-cctk_nghostzones[1]\",\"cctk_lsh[0]-cctk_nghostzones[0]\"],\n",
    "                           [\"1\",\"1\",\"1\"],[\"#pragma omp parallel for\",\"\",\"\"], \"\", Ham_mom_string))\n",

[3.27150]

[3.27532]

    "        file.write(loop_int(Ham_mom_string))\n",

Replacement in BaikalX/doc/BaikalX.ipynb at line 432 [3.16]

B:BD[3.29405] → [3.29405:30039]

    "        file.write(lp.loop([\"i2\",\"i1\",\"i0\"],[\"0\", \"0\", \"0\"],\n",
    "                           [\"cctk_lsh[2]\",\"cctk_lsh[1]\",\"cctk_lsh[0]\"],\n",
    "                           [\"1\",\"1\",\"1\"],\n",
    "                            [\"#pragma omp parallel for\",\n",
    "                                         \"#include \\\"rfm_files/rfm_struct__read2.h\\\"\",\n",
    "                                         \"#include \\\"rfm_files/rfm_struct__read1.h\\\"\"],\"\",\n",
    "                                         \"#include \\\"rfm_files/rfm_struct__read0.h\\\"\\n\"+enforce_gammadet_string))\n",

[3.29405]

[3.30039]

    "        file.write(loop_all(enforce_gammadet_string))\n",

Replacement in BaikalX/doc/BaikalX.ipynb at line 796 [3.16]

B:BD[3.53114] → [3.53114:53193]

    "with open(append_to_make_code_defn_list(\"Banner.c\"),\"w\") as file:\n",

[3.53114]

[3.53193]

    "with open(append_to_make_code_defn_list(\"Banner.cc\"),\"w\") as file:\n",

Replacement in BaikalX/doc/BaikalX.ipynb at line 800 [3.16]

B:BD[3.53262] → [3.53262:53296]

    "void BaikalETK_Banner() \n",

[3.53262]

[3.53296]

    "#include \"cctk.h\"\n",
    "\n",
    "int BaikalETK_Banner(void) \n",

Insertion in BaikalX/doc/BaikalX.ipynb at line 809 [3.16]
[3.53619]
[3.56879]
```
    "    file.write(\"return 0;\\n\")\n",
```

Replacement in BaikalX/doc/BaikalX.ipynb at line 820 [3.16]

B:BD[3.57062] → [3.57062:57144]

    "with open(append_to_make_code_defn_list(\"zero_rhss.c\"),\"w\") as file:\n",

[3.57062]

[3.57144]

    "with open(append_to_make_code_defn_list(\"zero_rhss.cc\"),\"w\") as file:\n",

Insertion in BaikalX/doc/BaikalX.ipynb at line 822 [3.16]
[3.57175]
[3.57175]
```
    "#include <loop.hxx>    \n",
    "\n",
```

Insertion in BaikalX/doc/BaikalX.ipynb at line 834 [3.16]

[3.57475]

    "    rhss_loopers = \"\"\n",
    "    for gf in rhs_list:\n",
    "        rhss_loopers += \"Loop::GF3D<CCTK_REAL, 0, 0, 0> \"+gf+\"_(cctkGH, \"+gf+\");\\n\"\n",
    "    file.write(rhss_loopers)\n",
    "    \n",

Replacement in BaikalX/doc/BaikalX.ipynb at line 841 [3.16]

B:BD[3.57545] → [3.57545:57633]

    "        set_rhss_to_zero += gf+\"[CCTK_GFINDEX3D(cctkGH,i0,i1,i2)] = 0.0;\\n\"\n",

[3.57545]

[3.57633]

    "        set_rhss_to_zero += gf+\"_(i0, i1, i2) = 0.0;\\n\"\n",

Replacement in BaikalX/doc/BaikalX.ipynb at line 843 [3.16]

B:BD[3.57647] → [3.57647:57960]

    "    file.write(lp.loop([\"i2\",\"i1\",\"i0\"],[\"0\", \"0\", \"0\"],\n",
    "                       [\"cctk_lsh[2]\",\"cctk_lsh[1]\",\"cctk_lsh[0]\"],\n",
    "                       [\"1\",\"1\",\"1\"],\n",
    "                        [\"#pragma omp parallel for\",\"\",\"\",],\"\",set_rhss_to_zero))\n",

[3.57647]

[3.63684]

    "    file.write(loop_all(set_rhss_to_zero))\n",

Replacement in BaikalX/doc/BaikalX.ipynb at line 883 [3.16]

B:BD[3.65560] → [3.65560:65645]

    "with open(append_to_make_code_defn_list(\"ADM_to_BSSN.c\"), \"w\") as file:\n",

[3.65560]

[3.65645]

    "with open(append_to_make_code_defn_list(\"ADM_to_BSSN.cc\"), \"w\") as file:\n",

Insertion in BaikalX/doc/BaikalX.ipynb at line 887 [3.16]
[3.65713]
[3.65713]
```
    "#include <loop.hxx>\n",
    "\n",
```

Insertion in BaikalX/doc/BaikalX.ipynb at line 930 [3.16]

[3.66901]

    "\n",
    "    adm_to_bssn_gfs = []\n",
    "    adm_to_bssn_gfs.append(\"alphaSphorCart\"+\"GF\")\n",
    "    adm_to_bssn_gfs.append(\"trK\"+\"GF\")\n",
    "    adm_to_bssn_gfs.append(\"alpha\"+\"GF\")\n",
    "    adm_to_bssn_gfs.append(\"cf\"+\"GF\")\n",
    "    for i in range(DIM):\n",
    "        adm_to_bssn_gfs.append(\"betaSphorCartU\"+str(i)+\"GF\")\n",
    "        adm_to_bssn_gfs.append(\"BSphorCartU\"+str(i)+\"GF\")\n",
    "        adm_to_bssn_gfs.append(\"vetU\"+str(i)+\"GF\")\n",
    "        adm_to_bssn_gfs.append(\"betU\"+str(i)+\"GF\")\n",
    "        adm_to_bssn_gfs.append(\"lambdaU\"+str(i)+\"GF\")\n",
    "        for j in range(i,DIM):\n",
    "            adm_to_bssn_gfs.append(\"gammaSphorCartDD\"+str(i)+str(j)+\"GF\")\n",
    "            adm_to_bssn_gfs.append(\"KSphorCartDD\"+str(i)+str(j)+\"GF\")\n",
    "            adm_to_bssn_gfs.append(\"hDD\"+str(i)+str(j)+\"GF\")\n",
    "            adm_to_bssn_gfs.append(\"aDD\"+str(i)+str(j)+\"GF\")\n",
    "    for gf in adm_to_bssn_gfs:\n",
    "        file.write(\"Loop::GF3D<CCTK_REAL,0,0,0> \"+gf+\"_(cctkGH,\"+gf+\");\\n\")\n",
    "    file.write(\"\\n\")\n",
    "\n",

Replacement in BaikalX/doc/BaikalX.ipynb at line 982 [3.16]

B:BD[3.69069] → [3.69069:69316]

    "    file.write(lp.loop([\"i2\",\"i1\",\"i0\"],[\"0\",\"0\",\"0\"],[\"cctk_lsh[2]\",\"cctk_lsh[1]\",\"cctk_lsh[0]\"],\n",
    "                       [\"1\",\"1\",\"1\"],[\"#pragma omp parallel for\",\"\",\"\"],\"\",all_but_lambdaU_outC))\n",

[3.69069]

[3.69316]

    "    file.write(loop_all(all_but_lambdaU_outC))\n",

Replacement in BaikalX/doc/BaikalX.ipynb at line 990 [3.16]

B:BD[3.69845] → [3.69845:70212]

    "    file.write(lp.loop([\"i2\",\"i1\",\"i0\"],[\"cctk_nghostzones[2]\",\"cctk_nghostzones[1]\",\"cctk_nghostzones[0]\"],\n",
    "       [\"cctk_lsh[2]-cctk_nghostzones[2]\",\"cctk_lsh[1]-cctk_nghostzones[1]\",\"cctk_lsh[0]-cctk_nghostzones[0]\"],\n",
    "       [\"1\",\"1\",\"1\"],[\"#pragma omp parallel for\",\"\",\"\"],\"\",lambdaU_expressions_FDout))\n",

[3.69845]

[3.70212]

    "    file.write(loop_all(lambdaU_expressions_FDout))\n",

Replacement in BaikalX/doc/BaikalX.ipynb at line 993 [3.16]

B:BD[3.70257] → [3.70257:70410]

    "    ExtrapolateGammas(cctkGH,lambdaU0GF);\n",
    "    ExtrapolateGammas(cctkGH,lambdaU1GF);\n",
    "    ExtrapolateGammas(cctkGH,lambdaU2GF);\n",

[3.70257]

[3.70410]

    "    //ExtrapolateGammas(cctkGH,lambdaU0GF);\n",
    "    //ExtrapolateGammas(cctkGH,lambdaU1GF);\n",
    "    //ExtrapolateGammas(cctkGH,lambdaU2GF);\n",

Replacement in BaikalX/doc/BaikalX.ipynb at line 1022 [3.16]

B:BD[3.71255] → [3.71255:71340]

    "with open(append_to_make_code_defn_list(\"BSSN_to_ADM.c\"), \"w\") as file:\n",

[3.71255]

[3.71340]

    "with open(append_to_make_code_defn_list(\"BSSN_to_ADM.cc\"), \"w\") as file:\n",

Insertion in BaikalX/doc/BaikalX.ipynb at line 1024 [3.16]
[3.71371]
[3.71371]
```
    "#include <loop.hxx>\n",
    "\n",
```

Insertion in BaikalX/doc/BaikalX.ipynb at line 1037 [3.16]

[3.71686]

    "    bssn_to_adm_gfs = []\n",
    "    bssn_to_adm_gfs.append(\"trK\"+\"GF\")\n",
    "    bssn_to_adm_gfs.append(\"alpha\"+\"GF\")\n",
    "    bssn_to_adm_gfs.append(\"alp\")\n",
    "    bssn_to_adm_gfs.append(\"cf\"+\"GF\")\n",
    "    for i in range(DIM):\n",
    "        bssn_to_adm_gfs.append(\"beta\"+chr(ord('x')+i))\n",
    "        bssn_to_adm_gfs.append(\"dtbeta\"+chr(ord('x')+i))\n",
    "        bssn_to_adm_gfs.append(\"vetU\"+str(i)+\"GF\")\n",
    "        bssn_to_adm_gfs.append(\"betU\"+str(i)+\"GF\")\n",
    "        for j in range(i,DIM):\n",
    "            bssn_to_adm_gfs.append(\"g\"+chr(ord('x')+i)+chr(ord('x')+j))\n",
    "            bssn_to_adm_gfs.append(\"k\"+chr(ord('x')+i)+chr(ord('x')+j))\n",
    "            bssn_to_adm_gfs.append(\"hDD\"+str(i)+str(j)+\"GF\")\n",
    "            bssn_to_adm_gfs.append(\"aDD\"+str(i)+str(j)+\"GF\")\n",
    "    for gf in bssn_to_adm_gfs:\n",
    "        file.write(\"Loop::GF3D<CCTK_REAL,0,0,0> \"+gf+\"_(cctkGH,\"+gf+\");\\n\")\n",
    "    file.write(\"\\n\")    \n",
    "\n",

Replacement in BaikalX/doc/BaikalX.ipynb at line 1059 [3.16]

B:BD[3.71788] → [3.71788:71912]

    "            btoa_lhrh.append(lhrh(lhs=\"g\"+chr(ord('x')+i)+chr(ord('x')+j)+\"[CCTK_GFINDEX3D(cctkGH,i0,i1,i2)]\",\n",

[3.71788]

[3.71912]

    "            btoa_lhrh.append(lhrh(lhs=\"g\"+chr(ord('x')+i)+chr(ord('x')+j)+\"_(p.I)\",\n",

Replacement in BaikalX/doc/BaikalX.ipynb at line 1063 [3.16]

B:BD[3.72054] → [3.72054:72178]

    "            btoa_lhrh.append(lhrh(lhs=\"k\"+chr(ord('x')+i)+chr(ord('x')+j)+\"[CCTK_GFINDEX3D(cctkGH,i0,i1,i2)]\",\n",

[3.72054]

[3.72178]

    "            btoa_lhrh.append(lhrh(lhs=\"k\"+chr(ord('x')+i)+chr(ord('x')+j)+\"_(p.I)\",\n",

Replacement in BaikalX/doc/BaikalX.ipynb at line 1065 [3.16]

B:BD[3.72242] → [3.72242:72337]

    "    btoa_lhrh.append(lhrh(lhs=\"alp[CCTK_GFINDEX3D(cctkGH,i0,i1,i2)]\",rhs=Bq.alpha))\n",

[3.72242]

[3.72337]

    "    btoa_lhrh.append(lhrh(lhs=\"alp_(p.I)\",rhs=Bq.alpha))\n",

Replacement in BaikalX/doc/BaikalX.ipynb at line 1068 [3.16]

B:BD[3.72393] → [3.72393:72500]

    "        btoa_lhrh.append(lhrh(lhs=\"beta\"+chr(ord('x')+i)+\"[CCTK_GFINDEX3D(cctkGH,i0,i1,i2)]\",\n",

[3.72393]

[3.72500]

    "        btoa_lhrh.append(lhrh(lhs=\"beta\"+chr(ord('x')+i)+\"_(p.I)\",\n",

Replacement in BaikalX/doc/BaikalX.ipynb at line 1072 [3.16]

B:BD[3.72605] → [3.72605:72714]

    "        btoa_lhrh.append(lhrh(lhs=\"dtbeta\"+chr(ord('x')+i)+\"[CCTK_GFINDEX3D(cctkGH,i0,i1,i2)]\",\n",

[3.72605]

[3.72714]

    "        btoa_lhrh.append(lhrh(lhs=\"dtbeta\"+chr(ord('x')+i)+\"_(p.I)\",\n",

Replacement in BaikalX/doc/BaikalX.ipynb at line 1077 [3.16]

B:BD[3.72971] → [3.72971:73215]

    "    file.write(lp.loop([\"i2\",\"i1\",\"i0\"],[\"0\",\"0\",\"0\"],[\"cctk_lsh[2]\",\"cctk_lsh[1]\",\"cctk_lsh[0]\"],\n",
    "                       [\"1\",\"1\",\"1\"],[\"#pragma omp parallel for\",\"\",\"\"],\"\",bssn_to_adm_Ccode))\n",

[3.72971]

[3.73215]

    "    file.write(loop_all(bssn_to_adm_Ccode))\n",

Replacement in BaikalX/doc/BaikalX.ipynb at line 1105 [3.16]

B:BD[3.74079] → [3.74079:74174]

    "with open(append_to_make_code_defn_list(\"driver_pt1_BSSN_Ricci.c\"), \"w\") as file:\n",

[3.74079]

[3.74174]

    "with open(append_to_make_code_defn_list(\"driver_pt1_BSSN_Ricci.cc\"), \"w\") as file:\n",

Replacement in BaikalX/doc/BaikalX.ipynb at line 1109 [3.16]
B:BD[3.74242] → [3.74242:74287]
```
    "#include \"SIMD/SIMD_intrinsics.h\"\n",
```
[3.74242]
[3.74287]
```
    "#include <loop.hxx>\n",
```

Replacement in BaikalX/doc/BaikalX.ipynb at line 1118 [3.16]

B:BD[3.74520] → [3.74520:74969]

    "    const CCTK_REAL NOSIMDinvdx0 = 1.0/CCTK_DELTA_SPACE(0);\n",
    "    const REAL_SIMD_ARRAY invdx0 = ConstSIMD(NOSIMDinvdx0);\n",
    "    const CCTK_REAL NOSIMDinvdx1 = 1.0/CCTK_DELTA_SPACE(1);\n",
    "    const REAL_SIMD_ARRAY invdx1 = ConstSIMD(NOSIMDinvdx1);\n",
    "    const CCTK_REAL NOSIMDinvdx2 = 1.0/CCTK_DELTA_SPACE(2);\n",
    "    const REAL_SIMD_ARRAY invdx2 = ConstSIMD(NOSIMDinvdx2);\n",
    "#include \"BSSN_Ricci.h\"\n",

[3.74520]

[3.74969]

    "    const CCTK_REAL invdx0 = 1.0/CCTK_DELTA_SPACE(0);\n",
    "    const CCTK_REAL invdx1 = 1.0/CCTK_DELTA_SPACE(1);\n",
    "    const CCTK_REAL invdx2 = 1.0/CCTK_DELTA_SPACE(2);\n",
    "\"\"\")\n",
    "    pt1_bssn_ricci_gfs = []\n",
    "    for i in range(DIM):\n",
    "        pt1_bssn_ricci_gfs.append(\"lambdaU\"+str(i)+\"GF\")\n",
    "        for j in range(i,DIM):\n",
    "            pt1_bssn_ricci_gfs.append(\"hDD\"+str(i)+str(j)+\"GF\")\n",
    "            pt1_bssn_ricci_gfs.append(\"RbarDD\"+str(i)+str(j)+\"GF\")\n",
    "    for gf in pt1_bssn_ricci_gfs:\n",
    "        file.write(\"Loop::GF3D<CCTK_REAL,0,0,0> \"+gf+\"_(cctkGH,\"+gf+\");\\n\")\n",
    "    file.write(\"\\n\")    \n",
    "\n",
    "    file.write(\"\"\"#include \"BSSN_Ricci.h\"\n",

Deletion in BaikalX/doc/BaikalX.ipynb at line 1146 [3.16]

B:BD[3.75421] → [3.75421:75459]

B:BD[3.75459] → [2.560:588]

∅:D[2.588] → [3.75485:76543]

B:BD[3.75485] → [3.75485:76543]

   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "def SIMD_declare_C_params():\n",
    "    SIMD_declare_C_params_str = \"\"\n",
    "    for i in range(len(par.glb_Cparams_list)):\n",
    "        # keep_param is a boolean indicating whether we should accept or reject\n",
    "        #    the parameter. singleparstring will contain the string indicating\n",
    "        #    the variable type.\n",
    "        keep_param, singleparstring = keep_param__return_type(par.glb_Cparams_list[i])\n",
    "\n",
    "        if (keep_param) and (\"CCTK_REAL\" in singleparstring):\n",
    "            parname = par.glb_Cparams_list[i].parname\n",
    "            SIMD_declare_C_params_str += \"    const \"+singleparstring + \"*NOSIMD\"+parname+\\\n",
    "                                         \" = CCTK_ParameterGet(\\\"\"+parname+\"\\\",\\\"BaikalETK\\\",NULL);\\n\"\n",
    "            SIMD_declare_C_params_str += \"    const REAL_SIMD_ARRAY \"+parname+\" = ConstSIMD(*NOSIMD\"+parname+\");\\n\"\n",
    "    return SIMD_declare_C_params_str"

Replacement in BaikalX/doc/BaikalX.ipynb at line 1154 [3.16]

B:BD[3.76659] → [3.76659:76753]

    "with open(append_to_make_code_defn_list(\"driver_pt2_BSSN_RHSs.c\"), \"w\") as file:\n",

[3.76659]

[3.76753]

    "with open(append_to_make_code_defn_list(\"driver_pt2_BSSN_RHSs.cc\"), \"w\") as file:\n",

Replacement in BaikalX/doc/BaikalX.ipynb at line 1158 [3.16]
B:BD[3.76821] → [3.76821:76866]
```
    "#include \"SIMD/SIMD_intrinsics.h\"\n",
```
[3.76821]
[3.76866]
```
    "#include <loop.hxx>\n",
```

Insertion in BaikalX/doc/BaikalX.ipynb at line 1166 [3.16]

[3.77032]

    "#define UPWIND_ALG(UpwindVecU) UpwindVecU > 0.0 ? 1.0 : 0.0\n",
    "\n",

Insertion in BaikalX/doc/BaikalX.ipynb at line 1170 [3.16]
[3.77122]
[3.77122]
```
    "    DECLARE_CCTK_PARAMETERS;\n",
```

Replacement in BaikalX/doc/BaikalX.ipynb at line 1172 [3.16]

B:BD[3.77132] → [3.77132:77593]

    "    const CCTK_REAL NOSIMDinvdx0 = 1.0/CCTK_DELTA_SPACE(0);\n",
    "    const REAL_SIMD_ARRAY invdx0 = ConstSIMD(NOSIMDinvdx0);\n",
    "    const CCTK_REAL NOSIMDinvdx1 = 1.0/CCTK_DELTA_SPACE(1);\n",
    "    const REAL_SIMD_ARRAY invdx1 = ConstSIMD(NOSIMDinvdx1);\n",
    "    const CCTK_REAL NOSIMDinvdx2 = 1.0/CCTK_DELTA_SPACE(2);\n",
    "    const REAL_SIMD_ARRAY invdx2 = ConstSIMD(NOSIMDinvdx2);\n",
    "\"\"\"+SIMD_declare_C_params()+\"\"\"\n",

[3.77132]

[3.77593]

    "    const CCTK_REAL invdx0 = 1.0/CCTK_DELTA_SPACE(0);\n",
    "    const CCTK_REAL invdx1 = 1.0/CCTK_DELTA_SPACE(1);\n",
    "    const CCTK_REAL invdx2 = 1.0/CCTK_DELTA_SPACE(2);\n",
    "\"\"\")\n",
    "    pt2_bssn_rhs_gfs = []\n",
    "    pt2_bssn_rhs_gfs.append(\"trK\"+\"GF\")\n",
    "    pt2_bssn_rhs_gfs.append(\"trK\"+\"_rhsGF\")\n",
    "    pt2_bssn_rhs_gfs.append(\"alpha\"+\"GF\")\n",
    "    pt2_bssn_rhs_gfs.append(\"alpha\"+\"_rhsGF\")\n",
    "    pt2_bssn_rhs_gfs.append(\"cf\"+\"GF\")\n",
    "    pt2_bssn_rhs_gfs.append(\"cf\"+\"_rhsGF\")\n",
    "    for i in range(DIM):\n",
    "        pt2_bssn_rhs_gfs.append(\"lambdaU\"+str(i)+\"GF\")\n",
    "        pt2_bssn_rhs_gfs.append(\"lambdaU\"+str(i)+\"_rhsGF\")\n",
    "        pt2_bssn_rhs_gfs.append(\"vetU\"+str(i)+\"GF\")\n",
    "        pt2_bssn_rhs_gfs.append(\"vetU\"+str(i)+\"_rhsGF\")\n",
    "        pt2_bssn_rhs_gfs.append(\"betU\"+str(i)+\"GF\")\n",
    "        pt2_bssn_rhs_gfs.append(\"betU\"+str(i)+\"_rhsGF\")\n",
    "        for j in range(i,DIM):\n",
    "            pt2_bssn_rhs_gfs.append(\"hDD\"+str(i)+str(j)+\"GF\")\n",
    "            pt2_bssn_rhs_gfs.append(\"aDD\"+str(i)+str(j)+\"GF\")\n",
    "            pt2_bssn_rhs_gfs.append(\"hDD\"+str(i)+str(j)+\"_rhsGF\")\n",
    "            pt2_bssn_rhs_gfs.append(\"aDD\"+str(i)+str(j)+\"_rhsGF\")\n",
    "            pt2_bssn_rhs_gfs.append(\"RbarDD\"+str(i)+str(j)+\"GF\")\n",
    "    for gf in pt2_bssn_rhs_gfs:\n",
    "        file.write(\"Loop::GF3D<CCTK_REAL,0,0,0> \"+gf+\"_(cctkGH,\"+gf+\");\\n\")\n",
    "    file.write(\"\\n\")    \n",
    "\n",
    "    file.write(\"\"\"\n",

Replacement in BaikalX/doc/BaikalX.ipynb at line 1223 [3.16]

B:BD[3.78382] → [3.78382:78486]

    "with open(append_to_make_code_defn_list(\"enforce_detgammabar_constraint.c\"), \"w\") as file:\n",

[3.78382]

[3.78486]

    "with open(append_to_make_code_defn_list(\"enforce_detgammabar_constraint.cc\"), \"w\") as file:\n",

Insertion in BaikalX/doc/BaikalX.ipynb at line 1226 [3.16]
[3.78544]
[3.78544]
```
    "\n",
    "#include <loop.hxx>\n",
```

Insertion in BaikalX/doc/BaikalX.ipynb at line 1236 [3.16]

[3.78810]

    "\"\"\")\n",
    "    enforce_detgammabar_constraint_gfs = []\n",
    "    for i in range(DIM):\n",
    "        for j in range(i,DIM):\n",
    "            enforce_detgammabar_constraint_gfs.append(\"hDD\"+str(i)+str(j)+\"GF\")\n",
    "    for gf in enforce_detgammabar_constraint_gfs:\n",
    "        file.write(\"Loop::GF3D<CCTK_REAL,0,0,0> \"+gf+\"_(cctkGH,\"+gf+\");\\n\")\n",
    "    file.write(\"\\n\")    \n",

Insertion in BaikalX/doc/BaikalX.ipynb at line 1245 [3.16]
[3.78820]
[3.78820]
```
    "    file.write(\"\"\"\n",
```

Replacement in BaikalX/doc/BaikalX.ipynb at line 1268 [3.16]

B:BD[3.79509] → [3.79509:79599]

    "with open(append_to_make_code_defn_list(\"BSSN_constraints.c\"), \"w\") as file:\n",

[3.79509]

[3.79599]

    "with open(append_to_make_code_defn_list(\"BSSN_constraints.cc\"), \"w\") as file:\n",

Insertion in BaikalX/doc/BaikalX.ipynb at line 1272 [3.16]
[3.79667]
[3.79667]
```
    "#include <loop.hxx>\n",
    "\n",
```

Insertion in BaikalX/doc/BaikalX.ipynb at line 1285 [3.16]

[3.80118]

    "\"\"\")\n",
    "    BSSN_constraints_gfs = []\n",
    "    BSSN_constraints_gfs.append(\"trK\"+\"GF\")\n",
    "    BSSN_constraints_gfs.append(\"cf\"+\"GF\")\n",
    "    BSSN_constraints_gfs.append(\"H\"+\"GF\")\n",
    "    for i in range(DIM):\n",
    "        BSSN_constraints_gfs.append(\"lambdaU\"+str(i)+\"GF\")\n",
    "        BSSN_constraints_gfs.append(\"MU\"+str(i)+\"GF\")\n",
    "        for j in range(i,DIM):\n",
    "            BSSN_constraints_gfs.append(\"hDD\"+str(i)+str(j)+\"GF\")\n",
    "            BSSN_constraints_gfs.append(\"aDD\"+str(i)+str(j)+\"GF\")\n",
    "    for gf in BSSN_constraints_gfs:\n",
    "        file.write(\"Loop::GF3D<CCTK_REAL,0,0,0> \"+gf+\"_(cctkGH,\"+gf+\");\\n\")\n",
    "    file.write(\"\\n\")    \n",

Insertion in BaikalX/doc/BaikalX.ipynb at line 1300 [3.16]
[3.80128]
[3.80128]
```
    "    file.write(\"\"\"\n",
```

Replacement in BaikalX/doc/BaikalX.ipynb at line 1353 [3.16]

B:BD[3.82550] → [3.82550:82647]

    "    outCparams = \"outCverbose=False,includebraces=False,preindent=2,SIMD_enable=True\"\n",

[3.82550]

[3.82647]

    "    outCparams = \"outCverbose=False,includebraces=False,preindent=2,SIMD_enable=False\"\n",

Replacement in BaikalX/doc/BaikalX.ipynb at line 1355 [3.16]

B:BD[3.82732] → [3.82732:82982]

    "    T4UUstr_loop = lp.loop([\"i2\",\"i1\",\"i0\"],[\"0\",\"0\",\"0\"],[\"cctk_lsh[2]\",\"cctk_lsh[1]\",\"cctk_lsh[0]\"],\n",
    "                           [\"1\",\"1\",\"SIMD_width\"],[\"#pragma omp parallel for\",\"\",\"\"],\"\",T4UUstr)\n",

[3.82732]

[3.82982]

    "    T4UUstr_loop = loop_all(T4UUstr)\n",

Replacement in BaikalX/doc/BaikalX.ipynb at line 1357 [3.16]

B:BD[3.82992] → [3.82992:83086]

    "    with open(append_to_make_code_defn_list(\"driver_BSSN_T4UU.c\"), \"w\") as file:\n",

[3.82992]

[3.83086]

    "    with open(append_to_make_code_defn_list(\"driver_BSSN_T4UU.cc\"), \"w\") as file:\n",

Deletion in BaikalX/doc/BaikalX.ipynb at line 1364 [3.16]
B:BD[3.83266] → [3.83266:83321]
```
    "\n",
    "#include \"SIMD/SIMD_intrinsics.h\"\n",
```