[JModelica] Singularities in Buildings.BoundaryConditions.SolarGeometry.Examples.ProjectedShadowLength and Buildings.HeatTransfer.Windows.Examples.BeamDepthInRoom

[agnesramle]

Since tan(acos(x) = sqrt(1-x^2)/x and 0.7433783508418161 * cos(proShaLenSou.decAng.decAng) * cos(proShaLenSou.solHouAng.solHouAng) + 0.6688711591178841 * sin(proShaLenSou.decAng.decAng) in the assignment below becomes 0 at two time points, you can get an arbitrary big spike in the trajectory for proSheLenSou.shaLen.u:

proShaLenSou.lessThreshold.u := acos(0.7433783508418161 * cos(proShaLenSou.decAng.decAng) * cos(proShaLenSou.solHouAng.solHouAng) + 0.6688711591178841 * sin(proShaLenSou.decAng.decAng))
--- Solved equation ---
proShaLenSou.shaLen.u := tan(proShaLenSou.lessThreshold.u)
(Equations from Buildings.BoundaryConditions.SolarGeometry.Examples.ProjectedShadowLength but similar equations in Buildings.HeatTransfer.Windows.Examples.BeamDepthInRoom)

This division by zero leads to verification problems between JModelica and Dymola since very small differences close to these points leads to big differences in the trajectory.