Commit b9e017b
committed
Add Phase 2 content: High-Order Methods and Elastic Wave Equations
New chapter:
- Chapter 8: High-Order Methods
- Dispersion Analysis (numerical dispersion, phase/group velocity)
- Fornberg Algorithm (computing FD weights for arbitrary order)
- Dispersion-Relation-Preserving (DRP) Schemes (Tam-Webb optimization)
- Custom weights in Devito: u.dx2(weights=weights)
New section in Chapter 7:
- 7.3 Elastic Wave Equations
- Velocity-stress formulation
- VectorTimeFunction and TensorTimeFunction
- Vector operators: div, grad, diag
- Staggered grid discretization
- Varying Lame parameters for heterogeneous media
Key Devito patterns introduced:
- Custom FD weights for optimized dispersion
- VectorTimeFunction for velocity fields
- TensorTimeFunction for stress tensors (symmetric=True)
- Vector calculus operators (div, grad, diag)
- Staggered grid handling for first-order systems
New solvers: dispersion.py, drp_devito.py, elastic_devito.py
New tests: 73 tests (39 highorder, 34 elastic)
Added references for Fornberg, Tam-Webb DRP schemes1 parent ab9b38c commit b9e017b
13 files changed
Lines changed: 4488 additions & 11 deletions
File tree
- chapters
- highorder
- systems
- src
- highorder
- systems
- tests
| Original file line number | Diff line number | Diff line change | |
|---|---|---|---|
| |||
8 | 8 | | |
9 | 9 | | |
10 | 10 | | |
| 11 | + | |
| 12 | + | |
| Original file line number | Diff line number | Diff line change | |
|---|---|---|---|
| |||
21 | 21 | | |
22 | 22 | | |
23 | 23 | | |
| 24 | + | |
24 | 25 | | |
25 | 26 | | |
26 | 27 | | |
| |||
173 | 174 | | |
174 | 175 | | |
175 | 176 | | |
| 177 | + | |
176 | 178 | | |
177 | 179 | | |
178 | 180 | | |
| |||
0 commit comments