Permalink
Switch branches/tags
Commits on May 3, 2018
  1. Merge pull request #13 from anpar/fix-HorizontalReflector2D

    rhewett committed May 3, 2018
    Fix segy.core import from obspy
Commits on May 2, 2018
Commits on Oct 28, 2016
Commits on Oct 26, 2016
Commits on Oct 25, 2016
  1. Merge pull request #9 from rhewett/petsc4py_fix

    rhewett committed Oct 25, 2016
    Fixed petsc wrapper so that petsc is not required for install.
  2. Fixed petsc wrapper so that petsc is not required for install.

    rhewett committed Oct 25, 2016
    This should kill any request to use the petsc wrapper in any solver.
Commits on Oct 17, 2016
  1. Merge pull request #6 from lawillem/variable_density_higher_orders_ac…

    rhewett committed Oct 17, 2016
    …curacy
    
    Variable density higher orders accuracy
Commits on Mar 18, 2016
  1. Merge pull request #5 from lawillem/fix_slow_perm_mat_building

    rhewett committed Mar 18, 2016
    Improving 'build_permutation_matrix' and fixing a bug in 'image_processing'
Commits on Feb 2, 2016
  1. Removed the old build_heterogeneous_laplacian routine, no longer used…

    lawillem committed Feb 2, 2016
    …. It was hardcoded to second order accuracy. The other 2nd order derivative functions Ian wrote is still there because it is used in the adjoint routines. But when those are updated by using the Dx and Dz from ' build_derivative_matrix_VDA' those functions can also be removed
  2. cleaned up the code a bit. Added the new Laplacian to the VDA frequen…

    lawillem committed Feb 2, 2016
    …cy domain solver as well. Added warnings to the adjoint calculations of the VDA routines. They currently use (2nd order) operators from Ian still even though the solver has been modified and also works with arbitrary spatial order of accuracy. The adjoint field calculation and the linear forward model will have to be updated by someone in the future. I added an example demonstrating how the VDA solver agrees with CDA on a uniform model and 2nd order spatial accuracy. At higher orders of accuracy there is a small difference due to the different structure of the Laplacian. I also show the solver on a variable density medium.
  3. early results show that the VDA solver is self adjoint. It also gives…

    lawillem committed Feb 2, 2016
    … exactly the same wavefield as the CDA solver on a constant density model. Running the unchanged adjoint tests in temporal modeling shows that the adjoint test is passed. My formulation is in agreement with what Ian wrote. Some optimization is needed, I modify a sparse matrix which is not efficient of course. I have also not yet verified that the derivatives of the density around the (PML) Dirichlet boundaries are correctly implemented, as the derivative entries may no longer add up to 0. TODO
Commits on Jan 29, 2016
  1. fixed some bugs in temporal modeling which prevented the CDA adjoint …

    lawillem committed Jan 29, 2016
    …test from running. They are the result of the introduction of the argument 'imaging_period'. Found out that the VDA solver I implemented is not self adjoint...'
Commits on Jan 28, 2016
Commits on Jan 26, 2016
  1. Drastically improved efficiency, speed and memory usage of the functi…

    lawillem committed Jan 26, 2016
    …on build_permutation_matrix. This function was slowing down variable density simulations significantly. I also fixed a bug in the gallery where an array of float values was passed to np.pad, causing a crash on the new numpy versions
Commits on Jan 22, 2016
  1. Merge pull request #3 from mtcli/fix_layered_medium

    rhewett committed Jan 22, 2016
    Fixed small bugs in layered_medium.py
  2. Merge pull request #4 from rhewett/lebrat_pr_merged

    rhewett committed Jan 22, 2016
    Lebrat pr merged
  3. Add Compact Helmoltz Operator for PML

    lebrat authored and rhewett committed Sep 29, 2015
  4. Add PETSc solver (i.e. Mumps SUPERLU_Dist ...) and the multi right ha…

    lebrat authored and rhewett committed Sep 28, 2015
    …nd side solve over shots
Commits on Jan 16, 2016
  1. Open MP for multithreading

    lebrat authored and rhewett committed Sep 21, 2015
  2. Adding the save option to reload shots

    lebrat authored and rhewett committed Sep 20, 2015
Commits on Jan 13, 2016
Commits on Oct 8, 2015
  1. Merge pull request #1 from imcbrear/master

    rhewett committed Oct 8, 2015
    VariableDensity
Commits on Aug 19, 2015
  1. fixed sign of gradient in temporal _gradient_helper

    Ian McBrearty Ian McBrearty
    Ian McBrearty authored and Ian McBrearty committed Aug 19, 2015
  2. Merge branch 'Imaging_Period_Temporal'

    Ian McBrearty Ian McBrearty
    Ian McBrearty authored and Ian McBrearty committed Aug 19, 2015
  3. bug with Imaging Period for variable density fixed.

    Ian McBrearty Ian McBrearty
    Ian McBrearty authored and Ian McBrearty committed Aug 19, 2015
  4. Merge branch 'VariableDensity'

    Ian McBrearty Ian McBrearty
    Ian McBrearty authored and Ian McBrearty committed Aug 19, 2015
  5. Added Solvers for 2D time/frequency Variable Density

    Ian McBrearty Ian McBrearty
    Ian McBrearty authored and Ian McBrearty committed Aug 19, 2015
Commits on Jun 9, 2015
  1. Moved hgignore to gitignore

    rhewett committed Jun 9, 2015
Commits on Jun 4, 2015
  1. gradient accuracy is preserved with this version. Can run new test in…

    lawillem committed Jun 4, 2015
    … tests folder. Not tested in inversion yet
  2. first implementation of saving the forward wavefield at a reduced rat…

    lawillem committed Jun 4, 2015
    …e. Not implemented in optimization yet
Commits on Apr 27, 2015
Commits on Feb 18, 2015
  1. merge in model parameter fix

    lawillem committed Feb 18, 2015