main.py

from create_initial_matrix import create_initial_matrix
from copying import turing_pattern
from create_data import create_data
from optimising import optimise
from MCMC import MCMC_main, MCMC_save_plot
from turing_class import NumericalSolver


a = NumericalSolver()
a.plot_fig()

size = 100 # size of matrix dimensions
time_lst = [0.0, 1.0, 2.0, 3.0, 4.0, 5.0]  # times that we might want to solve for

create_initial_matrix(size)  # create initial matrix
turing_pattern(a=2.8e-4, b=5e-3, tau=2, k=0)  # making clean data
create_data(time_lst)  # making noisy data
a_est, b_est = optimise(time=1.0, a_initial=0.0, a_final=1e-3, b_initial=0.0, b_final=1e-2, iters=10)  # rough parameter recovery of a and b
res, final_guess = MCMC_main(a_init=a_est, b_init=b_est, a_max=1e-3, b_max=1e-2, iters=1000)  # this is quite slow
MCMC_save_plot(res)