for consistency with other interface classes

    t = re.sub(r'(\W)(Simulation)(\W)', r'\1I\2\3', t)
    t = re.sub(r'(/)I(Simulation)', r'\1\2', t)

As recommended by
https://google.github.io/styleguide/cppguide.html#The__define_Guard

Previous implementation provided wrong number of elements in
the intensity map for off-spec simulation. In the case of GISAS
it has skipped masked areas

Simulation::numberOfSimulationElements was previously used as a synonym
to the number of values in the intensity map. However it is different
in the case of SpecularSimulation with beam divergence and could lead
to malfunctioning in the future. To avoid confusion,
numberOfSimulationElements method was placed in the
protected scope of Simulation class. Simulation::intensityMapSize returns
the actual size of intensity map.

Redmine: #2259

Redmine: #1917