This is a code for finding the Lefschetz thimbles. The Lefschetz thimbles are optimal integration contours for evaluating oscillatory integrals based on the Picard-Lefschetz theory, which has been applied to fields including radio astronomy, condensed matter physics, and quantum cosmology.
- This particular code features the constant phase contour method which can locate the Lefschetz thimbles efficiently for 1D integrals and can be extended for integrating 2D integrals.
- See Shi (2024) for a full description of the concept of the code.
- The main file is "find_thimble.py"
- Prerequirements are Numpy, Scipy, Matplotlib, and Networkx
- Contact Xun Shi for further questions relating to usage: [email protected].
If your work makes use of this code, please cite Shi (2024) and provide a url link to this github page.