Deploy solutions

Author: utkinis

slide-notebooks/deploy_notebooks.jl

@@ -1,6 +1,6 @@
 using Literate
 ## include Literate scripts starting with following 3 letters in the deploy
-incl = "l5_"
+incl = "l9_"
 ## Set `sol=true` to produce output with solutions contained and hints stripts. Otherwise the other way around.
 sol = true
 ##

slide-notebooks/notebooks/l7_1-xpu.ipynb

@@ -427,7 +427,7 @@
    "cell_type": "code",
    "source": [
     "@parallel function update_Pf!(Pf, qDx, qDy, _dx, _dy, _β_dτ)\n",
-    "    Pf = ...\n",
+    "    @all(Pf) = @all(Pf) - (@d_xa(qDx) * _dx + @d_ya(qDy) * _dy) * _β_dτ\n",
     "    return nothing\n",
     "end"
    ],

slide-notebooks/notebooks/l8_1-multi-xpu.ipynb

@@ -339,8 +339,8 @@
     "CL[2:end-1] .= CL[2:end-1] .+ dt*D*diff(diff(CL)/dx)/dx\n",
     "CR[2:end-1] .= CR[2:end-1] .+ dt*D*diff(diff(CR)/dx)/dx\n",
     "# Update boundaries (MPI)\n",
-    "CL[end] = ...\n",
-    "CR[1]   = ...\n",
+    "CL[end] = CR[2]\n",
+    "CR[1]   = CL[end-1]\n",
     "# Global picture\n",
     "C .= [CL[1:end-1]; CR[2:end]]\n",
     "```"
@@ -385,8 +385,8 @@
     "\n",
     "```julia\n",
     "# Update boundaries (MPI)\n",
-    "CL[end] = ...\n",
-    "CR[1]   = ...\n",
+    "CL[end] = CR[2]\n",
+    "CR[1]   = CL[end-1]\n",
     "```\n",
     "\n",
     "in order make the code work properly and run it again. Note what has changed in the visualisation."
@@ -422,7 +422,8 @@
     "    C[2:end-1,ip] .= C[2:end-1,ip] .+ dt*D*diff(diff(C[:,ip])/dxg)/dxg\n",
     "end\n",
     "for ip = 1:np-1 # update boundaries\n",
-    "   # ...\n",
+    "   C[end,ip  ] = C[    2,ip+1]\n",
+    "   C[  1,ip+1] = C[end-1,ip  ]\n",
     "end\n",
     "for ip = 1:np # global picture\n",
     "    i1 = 1 + (ip-1)*(nx-2)\n",
@@ -472,7 +473,7 @@
     "# Initial condition\n",
     "for ip = 1:np\n",
     "    for ix = 1:nx\n",
-    "        x[ix,ip] = ...\n",
+    "       x[ix,ip] = ( (ip-1)*(nx-2) + (ix-0.5) )*dxg - 0.5*lx\n",
     "        C[ix,ip] = exp(-x[ix,ip]^2)\n",
     "    end\n",
     "    i1 = 1 + (ip-1)*(nx-2)\n",
@@ -596,19 +597,19 @@
     "@views function update_halo!(A, neighbors_x, comm)\n",
     "    # Send to / receive from neighbour 1 (\"left neighbor\")\n",
     "    if neighbors_x[1] != MPI.PROC_NULL\n",
-    "        sendbuf = ??\n",
-    "        recvbuf = ??\n",
-    "        MPI.Send(??,  neighbors_x[?], 0, comm)\n",
-    "        MPI.Recv!(??, neighbors_x[?], 1, comm)\n",
-    "        A[1] = ??\n",
+    "        sendbuf = A[2]\n",
+    "        recvbuf = zeros(1)\n",
+    "        MPI.Send(sendbuf,  neighbors_x[1], 0, comm)\n",
+    "        MPI.Recv!(recvbuf, neighbors_x[1], 1, comm)\n",
+    "        A[1] = recvbuf[1]\n",
     "    end\n",
     "    # Send to / receive from neighbour 2 (\"right neighbor\")\n",
     "    if neighbors_x[2] != MPI.PROC_NULL\n",
-    "        sendbuf = ??\n",
-    "        recvbuf = ??\n",
-    "        MPI.Recv!(??, neighbors_x[?], 0, comm)\n",
-    "        MPI.Send(??,  neighbors_x[?], 1, comm)\n",
-    "        A[end] = ??\n",
+    "        sendbuf = A[end-1]\n",
+    "        recvbuf = zeros(1)\n",
+    "        MPI.Recv!(recvbuf, neighbors_x[2], 0, comm)\n",
+    "        MPI.Send(sendbuf,  neighbors_x[2], 1, comm)\n",
+    "        A[end] = recvbuf[1]\n",
     "    end\n",
     "    return\n",
     "end\n",
@@ -1077,11 +1078,11 @@
    "file_extension": ".jl",
    "mimetype": "application/julia",
    "name": "julia",
-   "version": "1.10.5"
+   "version": "1.11.1"
   },
   "kernelspec": {
-   "name": "julia-1.10",
-   "display_name": "Julia 1.10.5",
+   "name": "julia-1.11",
+   "display_name": "Julia 1.11.1",
    "language": "julia"
   }
  },

slide-notebooks/notebooks/l9_1-projects.ipynb

@@ -469,11 +469,11 @@
    "file_extension": ".jl",
    "mimetype": "application/julia",
    "name": "julia",
-   "version": "1.10.5"
+   "version": "1.11.1"
   },
   "kernelspec": {
-   "name": "julia-1.10",
-   "display_name": "Julia 1.10.5",
+   "name": "julia-1.11",
+   "display_name": "Julia 1.11.1",
    "language": "julia"
   }
  },

slide-notebooks/notebooks/l9_2-doc.ipynb

@@ -324,11 +324,11 @@
    "file_extension": ".jl",
    "mimetype": "application/julia",
    "name": "julia",
-   "version": "1.10.5"
+   "version": "1.11.1"
   },
   "kernelspec": {
-   "name": "julia-1.10",
-   "display_name": "Julia 1.10.5",
+   "name": "julia-1.11",
+   "display_name": "Julia 1.11.1",
    "language": "julia"
   }
  },

website/_literate/l7_1-xpu_web.jl

@@ -231,7 +231,7 @@ Note that currently the shorthand `-=` notation is not supported and we need to
 By analogy, update `update_Pf!`.
 """
 @parallel function update_Pf!(Pf, qDx, qDy, _dx, _dy, _β_dτ)
-    Pf = ...
+    @all(Pf) = @all(Pf) - (@d_xa(qDx) * _dx + @d_ya(qDy) * _dy) * _β_dτ
     return nothing
 end
 

website/_literate/l8_1-multi-xpu_web.jl

@@ -182,8 +182,8 @@ The idea of this fake parallelisation approach is the following:
 CL[2:end-1] .= CL[2:end-1] .+ dt*D*diff(diff(CL)/dx)/dx
 CR[2:end-1] .= CR[2:end-1] .+ dt*D*diff(diff(CR)/dx)/dx
 # Update boundaries (MPI)
-CL[end] = ...
-CR[1]   = ...
+CL[end] = CR[2]
+CR[1]   = CL[end-1]
 # Global picture
 C .= [CL[1:end-1]; CR[2:end]]
 ```
@@ -209,8 +209,8 @@ Then, add the required boundary update:
 
 ```julia
 # Update boundaries (MPI)
-CL[end] = ...
-CR[1]   = ...
+CL[end] = CR[2]
+CR[1]   = CL[end-1]
 ```
 
 in order make the code work properly and run it again. Note what has changed in the visualisation.
@@ -234,7 +234,8 @@ for ip = 1:np # compute physics locally
     C[2:end-1,ip] .= C[2:end-1,ip] .+ dt*D*diff(diff(C[:,ip])/dxg)/dxg
 end
 for ip = 1:np-1 # update boundaries
-   # ...
+   C[end,ip  ] = C[    2,ip+1]
+   C[  1,ip+1] = C[end-1,ip  ]
 end
 for ip = 1:np # global picture
     i1 = 1 + (ip-1)*(nx-2)
@@ -266,7 +267,7 @@ md"""
 # Initial condition
 for ip = 1:np
     for ix = 1:nx
-        x[ix,ip] = ...
+       x[ix,ip] = ( (ip-1)*(nx-2) + (ix-0.5) )*dxg - 0.5*lx
         C[ix,ip] = exp(-x[ix,ip]^2)
     end
     i1 = 1 + (ip-1)*(nx-2)
@@ -351,19 +352,19 @@ Then, we need to (2.) implement a boundary update routine, which can have the fo
 @views function update_halo!(A, neighbors_x, comm)
     # Send to / receive from neighbour 1 ("left neighbor")
     if neighbors_x[1] != MPI.PROC_NULL
-        sendbuf = ??
-        recvbuf = ??
-        MPI.Send(??,  neighbors_x[?], 0, comm)
-        MPI.Recv!(??, neighbors_x[?], 1, comm)
-        A[1] = ??
+        sendbuf = A[2]
+        recvbuf = zeros(1)
+        MPI.Send(sendbuf,  neighbors_x[1], 0, comm)
+        MPI.Recv!(recvbuf, neighbors_x[1], 1, comm)
+        A[1] = recvbuf[1]
     end
     # Send to / receive from neighbour 2 ("right neighbor")
     if neighbors_x[2] != MPI.PROC_NULL
-        sendbuf = ??
-        recvbuf = ??
-        MPI.Recv!(??, neighbors_x[?], 0, comm)
-        MPI.Send(??,  neighbors_x[?], 1, comm)
-        A[end] = ??
+        sendbuf = A[end-1]
+        recvbuf = zeros(1)
+        MPI.Recv!(recvbuf, neighbors_x[2], 0, comm)
+        MPI.Send(sendbuf,  neighbors_x[2], 1, comm)
+        A[end] = recvbuf[1]
     end
     return
 end