from simple import SIMPLE
from research import Research

model = SIMPLE((50, 100), (40, 40), 0.001, 700)

research = Research(model, f'{model.p.shape[0] * model.step}x{model.p.shape[1] * model.step}_{model.bfs_shape[0] * model.step}x{model.bfs_shape[1] * model.step}_Re{model.Re}/{model.step}')

is_complete = research.load()

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