Bring back VectorField

PTsolvers · May 23, 2024 · dda8709 · dda8709
1 parent cdb9cf1
commit dda8709
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/examples/stokes_3d_inc_ve_T_mpi_perf.jl b/examples/stokes_3d_inc_ve_T_mpi_perf.jl
@@ -135,7 +135,7 @@ end
     τ_old = TensorField(backend, grid)
     T     = Field(backend, grid, Center())
     T_old = Field(backend, grid, Center())
-    # qT    = VectorField(backend, grid)
+    qT    = VectorField(backend, grid)
     # Pr_v  = (me==0) ? KernelAbstractions.zeros(CPU(), Float64, size(interior(Pr)) .* dims(topo)) : nothing
     # initial conditions
     init_incl(x, y, z, x0, y0, z0, r, in, out) = ifelse((x - x0)^2 + (y - y0)^2 + (z - z0)^2 < r^2, in, out)