For more details, refer to the [theory overview](docs/source/documentation/mf_notes.md) and [algorithm description](docs/source/documentation/algorithm.md).
## How to use `pymf`?
The calculation of a mean-field Hamiltonian is a simple 3-step process:
1.**Define**
To specify the interacting problem, use a `Model` object which collects:
- Non-interacting Hamiltonian as a tight-binding dictionary.
- Interaction Hamiltonian as a tight-binding dictionary.
- Particle filling number in the unit cell.
2.**Guess**
Construct a starting guess for the mean-field correction.
3.**Solve**
Solve for the mean-field correction using the `solver` function and add it to the non-interacting part to obtain the total mean-field Hamiltonian.