Create the linear constraint matrices A and b.
Standard notation
The parameter constraints for the motional amplitude parameters
are:
0 <= theta <= pi,
0 <= theta_x <= theta_y <= pi,
-0.125 <= S <= 1,
0 <= sigma_max <= pi,
The pivot point parameter, domain position parameters, and
eigenframe parameters are unconstrained.
Matrix notation
In the notation A.x >= b, where A is an matrix of coefficients, x
is an array of parameter values, and b is a vector of scalars, these
inequality constraints are:
| 1 0 0 0 0 | | 0 |
| | | |
|-1 0 0 0 0 | | -pi |
| | | |
| 0 1 0 0 0 | | 0 |
| | | |
| 0 -1 0 0 0 | | theta | | -pi |
| | | | | |
| 0 -1 1 0 0 | | theta_x | | 0 |
| | | | | |
| 0 0 1 0 0 | . | theta_y | >= | 0 |
| | | | | |
| 0 0 -1 0 0 | | S | | -pi |
| | | | | |
| 0 0 0 1 0 | | sigma_max | | -0.125 |
| | | |
| 0 0 0 -1 0 | | -1 |
| | | |
| 0 0 0 0 1 | | 0 |
| | | |
| 0 0 0 0 -1 | | -pi |
- Parameters:
scaling_matrix (numpy rank-2 square matrix) - The diagonal, square scaling matrix.
- Returns: numpy rank-2 NxM array, numpy rank-1 N array
- The matrices A and b.
|