Skip to content

Commit

Permalink
update
Browse files Browse the repository at this point in the history
  • Loading branch information
PharmCat committed Feb 28, 2024
1 parent f7eda56 commit 1068a32
Show file tree
Hide file tree
Showing 3 changed files with 3 additions and 3 deletions.
2 changes: 1 addition & 1 deletion examples/animate.jl
Original file line number Diff line number Diff line change
Expand Up @@ -48,6 +48,6 @@ sphprob = SPHProblem(system, h, H, sphkernel, ρ, v, ptype, ρ₀, m₀, Δt,
sphprob.dpc_l₀ = 0.01
sphprob.dpc_λ = 0.01
sphprob.dpc_pmax = 36000
timesolve!(sphprob; batch = 100, timeframe = 1.1, writetime = 0.01, path = "D:/vtk/", pvc = true, anim = true)
timesolve!(sphprob; batch = 100, timeframe = 1.1, writetime = 0.001, path = "D:/vtk/", pvc = true, anim = true)

#makedf(sphprob)
2 changes: 1 addition & 1 deletion src/celllist.jl
Original file line number Diff line number Diff line change
Expand Up @@ -157,7 +157,7 @@ end
Update only distance
"""
@noinline function partialupdate!(c::GPUCellList)
@noinline function partialupdate!(c::GPUCellList, fillzero::Bool = true)
fill!(c.cnt, zero(Int32))
if fillzero fill!(c.pairs, (zero(Int32), zero(Int32))) end
neib_internal_2d!(c.pairs, c.cnt, c.cellpnum, c.points, c.celllist, c.dist)
Expand Down
2 changes: 1 addition & 1 deletion src/gpukernels2d.jl
Original file line number Diff line number Diff line change
Expand Up @@ -751,7 +751,7 @@ Compute ∂Π∂t - artificial viscosity. Add to `∑∂Π∂t`
```math
\\Pi_{ij} = \\begin{cases} \\frac{- \\alpha \\overline{c}_{ij} \\mu_{ij} + \\beta \\mu_{ij}^2 }{\\overline{\rho}_{ij}} & \\textbf{v}_{ij}\\cdot \\textbf{r}_{ij} < 0 \\\\ 0 & otherwise \\end{cases}
\\Pi_{ij} = \\begin{cases} \\frac{- \\alpha \\overline{c}_{ij} \\mu_{ij} + \\beta \\mu_{ij}^2 }{\\overline{\\rho}_{ij}} & \\textbf{v}_{ij}\\cdot \\textbf{r}_{ij} < 0 \\\\ 0 & otherwise \\end{cases}
\\\\
\\\\
Expand Down

0 comments on commit 1068a32

Please sign in to comment.