Metadata-Version: 2.4
Name: doublef
Version: 0.1.0
Summary: DoubleF: A fast and flexible phase association and earthquake location method
Author-email: Longfei Duan <duanlongfei20@mails.ucas.ac.cn>
Requires-Python: >=3.9
Description-Content-Type: text/markdown
Requires-Dist: matplotlib==3.9.4
Requires-Dist: numpy==1.26.3
Requires-Dist: obspy==1.4.2
Requires-Dist: pandas==2.3.0
Requires-Dist: pyrocko==2025.1.21
Requires-Dist: scipy==1.13.1
Requires-Dist: sobol-seq==0.2.0
Requires-Dist: torch==2.3.1
Requires-Dist: torchaudio==2.3.1
Requires-Dist: torchvision==0.18.1
Requires-Dist: tqdm==4.67.1
Requires-Dist: natsort==8.4.0
Requires-Dist: pillow==10.4.0


# DoubleF

**DoubleF** is a fast and flexible algorithm for seismic phase association and earthquake location. It uses an adaptive Sobol sampling strategy to efficiently search the source parameter space. The algorithm initializes a quasi-uniform sample set using a Sobol sequence. At each iteration, the objective function is evaluated, and the search region is refined based on the top-performing samples. This iterative narrowing continues until convergence. The highest-scoring sample is then selected, and associated phases within a residual threshold are included as the optimal event.

![Iterative search process](src/doublef/img/plot.png)

## Environment Setup and Installation

A dedicated Conda environment is recommended.

### Create a Conda environment

```bash
conda create -n doublef python=3.9
conda activate doublef
pip install doublef
```

If you prefer to install from source inside the environment:

```bash
conda create -n doublef python=3.9
conda activate doublef
git clone https://github.com/Lonngfei/DoubleF.git
cd DoubleF
pip install .
```


### GPU support

If you plan to run DoubleF on a GPU, make sure that your **PyTorch** version matches your **CUDA** version.

Please refer to the official PyTorch installation guide: [https://pytorch.org/get-started/previous-versions](https://pytorch.org/get-started/previous-versions).

### Confirm CUDA availability

Run the following in Python:

```python
import torch
print(torch.cuda.is_available())
```

* `True`: CUDA is available and PyTorch can use your GPU.
* `False`: Check the NVIDIA driver and CUDA-PyTorch compatibility.

### Show program information

```bash
doublef
```

This prints the program name,  version, and basic usage information.

## Input Files

Before running DoubleF, make sure the required input files are correctly prepared.

Typical inputs include:

```text
Picks/YYYYMMDD.csv       # Daily pick files
TravelTime/mymodel.nd    # Velocity model used for travel-time calculation
example.config           # Configuration file
```

The exact directory structure can be adjusted in the configuration file.

## Configuration

DoubleF is controlled through a configuration file such as `example.config`.

Most parameters do not need frequent modification. In most cases, only a few settings require special attention.

### 1. Velocity model and travel-time table

Set `cal_tt = True` when using a new velocity model.

This is usually required only once to generate the travel-time tables.

After the tables have been generated, set: `cal_tt = False`, so that the program loads the existing tables and skips recalculation.

### 2. Sampling parameters

In most applications, the default sampling settings are sufficient.

If the nearest-station distance is larger than **0.6°**, increasing the number of samples may improve the results.

### 3. Score calculation

DoubleF provides several alternative objective functions.

In most cases, the choice among these objective functions does not significantly affect the final results.

Users who are familiar with the method may further customize the scoring strategy if needed.

A custom objective function can be implemented by modifying: ``weight.py``, ``batch_weight.py`` in the source code.

### 4. Output settings

Set the output directory and related options according to your needs.

DoubleF automatically writes:

* logs
* configuration records
* phase association results

to the specified output path.

### 5. Memory and speed

#### `max_batch_size`

This parameter only affects computational efficiency and does **not** affect the final results. 

In general, a larger value may improve speed, but this is not always the case.  

Once the computation reaches saturation, further increasing `max_batch_size` may provide little or no additional speedup, while leading to higher memory usage.

### 6. Visualization

Visualization is usually recommended to be turned off during normal runs. 

It should only be enabled when intermediate inspection, debugging, or result checking is needed.


## Running the Program

Once the input files and configuration file are ready, run:

```bash
doublef example.config
```

If the installation and configuration are correct, DoubleF will start processing and write logs and results to the output directory.


## Output Phase File Format

A typical output phase file has the following format:

```text
# Year Month Day Hour Minute Second Latitude Longitude Depth Magnitude ErrHorizontal ErrVertical ErrTime RMS NumP NumS NumBoth NumSum ID
NET Station Distance PhaseTime Probability PhaseType Residual ML Mag Amplitude
```

### Notes

* `ErrHorizontal`, `ErrVertical`, and `ErrTime` do not necessarily represent the true location uncertainty.
*  These values describe the spatial dispersion of candidate solutions that are associated with the same nearby location. Specifically, they quantify the statistical deviation (e.g., mean or standard deviation) of these candidate locations relative to the final solution.
*  If only a single candidate solution exists, the dispersion cannot be computed and the value is reported as **NaN**.
*  When these values are unusually large or reported as NaN, the corresponding results should be interpreted with caution.


## Recommended Citation

If you use DoubleF in your research, please cite the related paper(s).

## Contact

For questions, suggestions, or bug reports, please open an issue on GitHub or contact:

**Longfei Duan**
Email: [duanlongfei20@mails.ucas.ac.cn](mailto:duanlongfei20@mails.ucas.ac.cn)
