Skip to content

Commit

Permalink
examples: enforce so>=8 for self adjoint operator for stability
Browse files Browse the repository at this point in the history
  • Loading branch information
mloubout committed Sep 3, 2021
1 parent 18c842a commit 465c67d
Showing 1 changed file with 3 additions and 0 deletions.
3 changes: 3 additions & 0 deletions examples/seismic/self_adjoint/example_iso.py
Original file line number Diff line number Diff line change
Expand Up @@ -11,6 +11,9 @@
def acoustic_sa_setup(shape=(50, 50, 50), spacing=(10.0, 10.0, 10.0),
tn=500., space_order=8, nbl=10, **kwargs):
# SA parameters
if space_order < 8:
info("Low space order not supported, running space_order=8")
space_order = 8
qmin = 0.1
qmax = 1000.0
fpeak = 0.010
Expand Down

0 comments on commit 465c67d

Please sign in to comment.