Opened 4 years ago
Last modified 10 months ago
#426 assigned enhancement
add inequality constraints
Reported by: | pkienzle | Owned by: | awashington |
---|---|---|---|
Priority: | major | Milestone: | SasView Next Release +1 |
Component: | SasView | Keywords: | |
Cc: | Work Package: | SasView Fitting Redesign |
Description
Bumps supports constraints such as parameter a < parameter b by having a penalty cost function evaluated for each parameter set. In this case, the function would be something like:
((b-a)^{2 + penalty) * (a>b) }
That way, if the fit finds itself outside the valid range a⇐b , it will be pushed back into the range with a quadratic function; the penalty value forces a step at the boundary so that a lower chisq out of bounds is not preferred to a value at the bounds. The penalty value should be large enough to remove any minima outside the bounds, otherwise algorithms like DREAM may be able to stumble onto it. Using a faster growing function, such as e^{(b-a) will reduce the problem. If the boundary can be a little less strict, then penalty can be zero. }
The constraints panel needs an extended syntax to handle inequality constraints. From this syntax, a constraint function could be automatically generated for use with bumps.
Change History (3)
comment:1 Changed 4 years ago by butler
- Work Package changed from SasView Bug Fixing to SasView Fitting Redesign
comment:2 Changed 13 months ago by awashington
- Owner set to awashington
- Status changed from new to assigned
You can define a custom model for soft inequality constraints, returning "x - (C<0)*C". Use it with a data file containing the point (x,y,dy) = (1,1,1). For example, if you have M1.radius_core and M1.radius_shell, you could set
If the total radius is less than 100 then C < 0 is false and x=1 is returned. If it is greater than 100, then (C<0) is true and the function returns the theory value 1 - C. The fitter then computes (yth - ydata)^{2}/dy^{2} which is ((1-C)-1)^{2}/1^{2} = C^{2} and adds this as the residual to the fit. Note that we are using (x,y) = (1,1) so that it plots on a log-log graph.
You can extend this to allow multiple constraints by defining a file with N points (1,1,1), (1,1,1), … and fit it using the following plugin model:
The data file would be:
The parameters L1, …, Lm are available as underlying parameters that can be fitted rather than fitting the model parameters directly. They don't actually contribute to the model output, so they can have an arbitrary value. For example, you can fit the total radius and the core radius by setting
and fitting M2.L1.
You can add hard constraints Hk as well by returning (inf if Hk<0 else 0). This won't work for the Levenberg-Marquardt fitter since it is not differentiable, but DREAM, DE and Nelder-Mead should completely avoid points in the infeasible region.