double mh_ref = 117.;
// Create SM and set top mass to 175 GeV                                                                                            
SM sm;  sm.set_qmass_pole(6, 175.);
THDM model;
model.set_SM(sm);
model.set_yukawas_type(type);
double m12_2 = 0.;
  bool pset = model.set_param_phys
    (mh,mH,mA,mC,sba,
     /*lambda_6*/ 0.,
     /*lambda_7*/ 0., m12_2,tb);
Constraints constr(model);
double S,T,U,V,W,X;
constr.oblique_param(mh_ref,S,T,U,V,W,X);
status = 0;
if ( !constr.check_stability()) { status += 1; }
if ( !constr.check_unitarity()) { status += 2; }
if ( !constr.check_perturbativity()) { status += 4; }
if ( !constr.check_masses()) { status +=8; }