ising_mf_t.py

ising_mf_t.py

from __future__ import print_function

import numpy as np

from ising_mf import find_solution, calc_entropy, calc_specific_heat


def main():
    h = 0
    tmin = 0
    tmax = 1
    n = 201
    for t in np.linspace(tmin, tmax, n):
        try:
            m = find_solution(1., t, h)
            s = calc_entropy(m)
            c = calc_specific_heat(m, 1., t)
        except ValueError:
            print(f"# {t:.4f}")
        else:
            print(f"{t:.4f} {m:.5e} {s:.5e} {c:.5e}")


if __name__ == '__main__':
    main()