Skip to content

Commit

Permalink
Update READMEs with small fixes
Browse files Browse the repository at this point in the history
  • Loading branch information
sloede committed Aug 24, 2020
1 parent f341b60 commit bf5a2fc
Show file tree
Hide file tree
Showing 2 changed files with 15 additions and 15 deletions.
6 changes: 3 additions & 3 deletions section-4.2.1-jeans-instability/README.md
Original file line number Diff line number Diff line change
Expand Up @@ -39,11 +39,11 @@ This script accesses the precomputed exact energy profiles stored in `E_kin.txt`
## Sec. 4.2.1, Figures 5, Jeans test gravity sub-cycle distribution
**Fig. 5a:**
```julia
Trixi.run("parameters_jeans_instability.toml")
julia> Trixi.run("parameters_jeans_instability.toml")
```

**Fig. 5b:**
```julia
Trixi.run("parameters_jeans_instability.toml",
time_integration_scheme_gravity="timestep_gravity_erk52_3Sstar!", cfl_gravity=1.2)
julia> Trixi.run("parameters_jeans_instability.toml",
time_integration_scheme_gravity="timestep_gravity_erk52_3Sstar!", cfl_gravity=1.2)
```
24 changes: 12 additions & 12 deletions section-4.2.2-sedov-blast/README.md
Original file line number Diff line number Diff line change
Expand Up @@ -12,49 +12,48 @@ julia> Pkg.instantiate()

julia> using Trixi

julia> const TRIXI_DIR = pathof(Trixi) |> dirname |> dirname
julia> const TRIXI_DIR = pathof(Trixi) |> dirname |> dirname # Only required for postprocessing
```
After that, you can generate the Sedov self gravity data for the figures by running the following code in the Julia REPL.


## Sec. 4.2.2, Figure 6, T=0.5, AMR meshes for Sedov + gravity
**T = 0.0 and T = 0.5:**
```julia
Trixi.run("parameters_sedov_self_gravity.toml", t_end=0.5)
julia> Trixi.run("parameters_sedov_self_gravity.toml", t_end=0.5)
```

**T = 1.0:**
```julia
Trixi.run("parameters_sedov_self_gravity.toml")
julia> Trixi.run("parameters_sedov_self_gravity.toml")
```

## Sec. 4.2.2, Figure 7a, T=0.5, Sedov + gravity
**AMR mesh:**
```julia
Trixi.run("parameters_sedov_self_gravity.toml", t_end=0.5)
julia> Trixi.run("parameters_sedov_self_gravity.toml", t_end=0.5)
```

**Uniform mesh:**
```julia
Trixi.run("parameters_sedov_self_gravity.toml",
amr_interval=0, initial_refinement_level=8, t_end=0.5)
julia> Trixi.run("parameters_sedov_self_gravity.toml", amr_interval=0, initial_refinement_level=8, t_end=0.5)
```

## Sec. 4.2.2, Figure 7b, T=1.0, Sedov + gravity
**AMR mesh:**
```julia
Trixi.run("parameters_sedov_self_gravity.toml")
julia> Trixi.run("parameters_sedov_self_gravity.toml")
```

**Uniform mesh:**
```julia
Trixi.run("parameters_sedov_self_gravity.toml", amr_interval=0, initial_refinement_level=8)
julia> Trixi.run("parameters_sedov_self_gravity.toml", amr_interval=0, initial_refinement_level=8)
```

## Sec. 4.2.2, Table 6, Sedov + gravity, performance uniform vs. AMR
**AMR mesh:**
```julia
julia> Trixi.run("parameters_sedov_self_gravity.toml"))
julia> Trixi.run("parameters_sedov_self_gravity.toml")
```

**Uniform mesh:**
Expand All @@ -65,8 +64,9 @@ julia> Trixi.run(parameters_sedov_self_gravity.toml", amr_interval=0)
## Converting Trixi.jl's .h5 files to VTK/VTI
To postprocess the solution files use
```julia
shell> $TRIXI_DIR/postprocessing/trixi2vtk --nvisnodes 8 --format vti path/to/solution_euler_xxx.h5
shell> $TRIXI_DIR/postprocessing/trixi2vtk --nvisnodes 8 --format vti out/solution_euler_000000.h5
shell> $TRIXI_DIR/postprocessing/trixi2vtk --nvisnodes 8 --format vti path/to/solution_gravity_xxx.h5
shell> $TRIXI_DIR/postprocessing/trixi2vtk --nvisnodes 8 --format vti out/solution_gravity_000000.h5
```
Then one can generate the plots using a visualization tool like Paraview or VisIt.
Then one can generate the plots using a visualization tool like [ParaView](https://www.paraview.org)
or [VisIt](https://visit.llnl.gov).

0 comments on commit bf5a2fc

Please sign in to comment.