from simple import SIMPLE
from research import Research

model = SIMPLE((100, 200), (50, 50), 0.001, 100)

research = Research(model, f'{model.p.shape[0]}x{model.p.shape[1]}-{model.Re}')

is_complete = research.load()

if is_complete:
    research.inspect()
else:
    research.solve(preview=True, save_model=True, save_plot=True)