| 623 | } |
| 624 | |
| 625 | static inline double pml_x(int i, double dx, double bloc, double a) { |
| 626 | double here = i * 0.5 / a; |
| 627 | return (0.5 / a * ((int)(dx * (2 * a) + 0.5) - (int)(fabs(bloc - here) * (2 * a) + 0.5))); |
| 628 | } |
| 629 | |
| 630 | void structure_chunk::use_pml(direction d, double dx, double bloc, double Rasymptotic, |
| 631 | double mean_stretch, pml_profile_func pml_profile, |