function Stereo = simulation(filename)
% The main function, it will call all necessary functions.

% Initialize simulation values
[General, Volume, Calibration] = initialize(filename);

% Generate the ground truth volumes and corresponding disparity maps.
Disparity = contruct_volumes(General, Volume);

% Calculate the parametric model based 
Model = parametric_model(Disparity);

% Generate random stereo image pairs.
Images = construct_images(General, Calibration);

% Reconstruct the three dimensional shape of a set of images based on a
% particular model.
Stereo = stereo_reconstruction(Images, Model);