[CIG-SEISMO] SPECFEM3D

[Dimitri Komatitsch]

Hi Dylan,

Yes, just use UNDO_ATTENUATION, which can undo any dissipative medium 
exactly, including PML. It is described in 
http://komatitsch.free.fr/preprints/GJI_undo_attenuation_2016.pdf .
We do that routinely here for FWI with PML, works beautifully.
Please cite that paper if you use it.

Thanks,
Best,
Dimitri.

On 10/30/2017 06:44 PM, James A. Smith wrote:
> Hi Dylan,
> 
> I have tested PML for the adjoint simulation. It seems to work, although 
> I am still playing with getting higher quality inversion results (as 
> compared to Stacey). Please see this issue: 
> https://github.com/geodynamics/specfem3d/issues/312
> 
> You will need to enable UNDO_ATTENUATION for PML to work in this case 
> (in setup/constanst.h). You will also have to comment out the lines in 
> the source code that stop the code if PML is turned on for an adjoint 
> simulation and recompile. I hope it works for you.
> 
> Best,
> James
> 
> ------------------------------------------------------------------------
> *From:* CIG-SEISMO [cig-seismo-bounces at geodynamics.org] on behalf of 
> Dylan Mikesell [dylanmikesell at boisestate.edu]
> *Sent:* Friday, October 27, 2017 2:17 PM
> *To:* cig-seismo at geodynamics.org
> *Subject:* [CIG-SEISMO] SPECFEM3D
> 
> Hi SPECFEM developers,
> 
> Is anyone currently working on the PML implementation for the adjoint. A 
> student and I are working on an adjoint tomography and planning to use 
> SPECFEM3D. We tried to save the forward simulations today and got the 
> following.
> 
> STOP PML_CONDITIONS is still under test for adjoint simulation
> 
> I checked the devel branch and this is still there. Is anyone actively 
> working on this topic? It would be nice to get the PML going with the 
> adjoint.
> 
> Best wishes,
> 
> Dylan
> 
> 
> _______________________________________________
> CIG-SEISMO mailing list
> CIG-SEISMO at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo
> 

-- 
Dimitri Komatitsch, CNRS Research Director (DR CNRS)
Laboratory of Mechanics and Acoustics, Marseille, France
http://komatitsch.free.fr