Fix some factor function corner cases
Ensure initial value for Gibbs solver continous variables are within domain bounds if only one bound is specified