Instead of computing the unitarity bounds and oblique corrections analytically, in this work, we have
adopted another strategy and used the publicly available SARAH/SPheno framework [1]}, [2]}, [3]} which enables one to
compute them numerically in a given particle physics model for a particular parameter space point (for an earlier implementation of this procedure, see e.g. Ref. [4]}).
In our numerical analysis substantiated in detail below, we have implemented the 3HDM model under consideration into this framework where both the unitarity constraints and the electroweak precision bounds on \(S,T,U\)  parameters have been consistently incorporated in the parameter scans designed to search for the valid physical regions of the model. In particular, the oblique parameters are computed by SPheno-4.0.4 at one loop-order using tree-level masses and then are verified against
the experimental bounds for every parameter space point.
