200     unsigned int maxIteration;
1.0E-30 double 0.0 < epsilonStar;
1.0E4   double 0.0 < lambdaStar;
2.0     double 1.0 < omegaStar;
-1.0E5  double lowerBound;
1.0E5   double upperBound;
0.1     double 0.0 <= betaStar <  1.0;
0.3     double 0.0 <= betaBar  <  1.0, betaStar <= betaBar;
0.9     double 0.0 < gammaStar  <  1.0;
1.0E-30 double 0.0 < epsilonDash;
512     precision; only relevant for variable precision


The above are the default settings. For STABLE_BUT_SLOW use
this or similar

1000    unsigned int maxIteration;
1.0E-30 double 0.0 < epsilonStar;
1.0E2   double 0.0 < lambdaStar;
2.0     double 1.0 < omegaStar;
-1.0E5  double lowerBound;
1.0E5   double upperBound;
0.2     double 0.0 <= betaStar <  1.0;
0.4     double 0.0 <= betaBar  <  1.0, betaStar <= betaBar;
0.5     double 0.0 < gammaStar  <  1.0;
1.0E-30 double 0.0 < epsilonDash;
512     precision; only relevant for variable precision


For UNSTABLE_BUT_FAST use this or similar

100     unsigned int maxIteration;
1.0E-20	double 0.0 < epsilonStar;
1.0E2   double 0.0 < lambdaStar;
2.0   	double 1.0 < omegaStar;
-1.0E5  double lowerBound;
1.0E5   double upperBound;
0.01    double 0.0 <= betaStar <  1.0;
0.02    double 0.0 <= betaBar  <  1.0, betaStar <= betaBar;
0.98    double 0.0 < gammaStar  <  1.0;
1.0E-30	double 0.0 < epsilonDash;
512     precision; only relevant for variable precision