groupeLIAMG / ttcr

Codes to do raytracing for geophysical applications
GNU General Public License v3.0
85 stars 34 forks source link

Overload the method projectOnFace when the gradient fall on an edge #41

Closed maher-nasr closed 3 years ago

maher-nasr commented 3 years ago

-the gradient projection is no longer done by taking the minimum angle between gradient vector and its projection. The local gradient at the cell scale is used to determine which face is to consider. Some tests are removed to speed up the method. Use the head wave traveltimes when the gradient projection points outside the tested edges.