[CIG-SEISMO] SPECFEM3D / PS

[Dimitri Komatitsch]

Hi all,

PS: we could thus consider renaming the UNDO_ATTENUATION flag to 
UNDO_ATTENUATION_AND_OR_PML. I know people do not like when we rename 
flags because Par_files then need to be updated, but this one has become 
a bit confusing.

Best,
Dimitri.

On 10/30/2017 08:59 PM, Dimitri Komatitsch wrote:
> 
> 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