Test (latest update 05, 2019)
1. 3D FWI with...2. Internal mesh + external tomo in 2D/3D 3. Lasif: Python3 + LBFGS4. Sp and Ps Kernels 5. Adjoint tomography: 2D and 3D 6. 2D/3D SEM-FK Model and Inversion%%%6. Discontinuous Galerkin Spectral-Element Method (SEM) vs. continuous Galerkin SEM (x)%%%7. Spectral-Element Method (SEM, by Specfem3D)and Finite-Difference Model (FDM, by SW4): PI-Lab ,Test (latest update 05, 2019)
to do:
- Lasif: adjoint source's formula and set the source location
- %%% second derivatives of the misfit function (see IFOS3D), compared with Tromp et al. 2005. How to compute the Hessian in IFO3D, i.e., how to compute the complex conjugate (du/dm)* , if we know du/dm. see (eq. 5.2 of Butzer's PhD thesis, 2015). see paper: Hessian kernels of seismic data functions based upon adjoint techniques, 2011.
interface model (common used): 1) use interface model to get converted waves that is used to compute the adjoin source (it includes time delay). 2) use the adjoint source as source time function to compute the adjoint wavefield.interface model (for Sp only): 1) use interface model to get converted waves that is used to compute the adjoin source (it includes time delay). 2) use the background model (without interface) to compute forward wavefield. 3) use the same background model as input, and use the adjoint determined at step 1) to compute the adjoint wavefield. 4) Kernels are computed based on the forward and adjoint wavefield.
5. separation:
- https://scikit-image.org/docs/dev/auto_examples/segmentation/plot_random_walker_segmentation.html#sphx-glr-auto-examples-segmentation-plot-random-walker-segmentation-py
- https://www.bogotobogo.com/python/OpenCV_Python/python_opencv3_Image_Watershed_Algorithm_Marker_Based_Segmentation.php
- https://www.pyimagesearch.com/2015/11/02/watershed-opencv/ (only for circle shape detection, e.g. coin)
- (Topological watershed: programs):https://perso.esiee.fr/~dpt-it/tw/TW_programs.html
- http://mathmistakes.info/facts/CalculusFacts/learn/doi/doi.html
- https://www.khanacademy.org/math/multivariable-calculus/multivariable-derivatives
A update of the specfem2D and specfem3D package:
- 1. Specfem2D: reading several external tomography files (fixed).
- 2.Specfem2D: reading several external tomography files and coupling acoustic and elastic simulation (fixed).
- 3. Specfem2D: full Hessian kernel calculation (developed).
- 4.Specfem2D: SEM-FK for external tomography files with coupling acoustic and elastic simulation (developed).
- 5.Specfem3D: reading several external tomography files (fixed).
- 6.Specfem3D: reading several external tomography files and coupling acoustic and elastic simulation (fixed).
- 7.Specfem3D: full Hessian kernel calculation (easy to do)
- 8.Specfem3D:SEM-FK for external tomography files with coupling acoustic and elastic simulation (can be done)
usr/bin/bashcd ./OUTPUT_FILES/DATABASES_MPIscp proc000000_vp.bin proc000000_vp_tomo.binscp proc000001_vp.bin proc000001_vp_tomo.binscp proc000002_vp.bin proc000002_vp_tomo.binscp proc000003_vp.bin proc000003_vp_tomo.binscp proc000004_vp.bin proc000004_vp_tomo.binscp proc000005_vp.bin proc000005_vp_tomo.binscp proc000006_vp.bin proc000006_vp_tomo.binscp proc000007_vp.bin proc000007_vp_tomo.binscp proc000008_vp.bin proc000008_vp_tomo.binscp proc000009_vp.bin proc000009_vp_tomo.binfor (( i=10; i<=99; i++ ))do scp proc0000${i}_vp.bin proc0000${i}_vp_tomo.bindonefor (( i=100; i<=119; i++))do scp proc000${i}_vp.bin proc000${i}_vp_tomo.bindone# copy to 1D to make vp_tomo - vp_1Dscp ./*vp_tomo.bin ../../../many_interfaces_PiLab_1D/OUTPUT_FILES/DATABASES_MPI/-------- submit job----qsub.