In this section we describe each operator, function, set, and command that you are likely to encounter in CVX. In some cases, limitations of the underlying solver place certain restrictions or caveats on their use:

- Functions marked with a dagger (†) are not supported natively by manu
solvers that CVX uses. Models incorporating these functions will be solved
most effectively by Mosek, the only bundled solver with support for the
exponential cone upon which these functions are constructed. For all
other solvers, these functions are handled using a successive
approximation method which makes multiple calls to the underlying
solver, achieving the same final precision. If your problem requires the use
of the successive approximation method, a warning will be issued.
This technique is discussed further in
*The successive approximation method*. As this section discusses, this is an experimental approach that works well in many cases, but cannot be guaranteed. - Functions involving powers (
*e.g.*,`x^p`) and \(p\)-norms (*e.g.*,`norm(x,p)`) are marked with a double dagger (‡). CVX represents these functions exactly when \(p\) is a rational number. For irrational values of`p`, a nearby rational is selected instead. See*Power functions and p-norms*for details on how both cases are handled.

Matlab’s standard arithmetic operations for addition `+`, subtraction `-`,
multiplication, `*` `.*`, division `/` `./` `\` `.\`, and
exponentiation `^` `.^` have been overloaded to work in
CVX whenever appropriate—that is, whenever their use is consistent
with both standard mathematical and Matlab conventions *and* the DCP
ruleset. For example:

- Two CVX expressions can be added together if they are of the same
dimension (or one is scalar) and have the same curvature (
*i.e.*, both are convex, concave, or affine). - A CVX expression can be multiplied or divided by a scalar constant. If the constant is positive, the curvature is preserved; if it is negative, curvature is reversed.
- An affine column vector CVX expression can be multiplied by a constant matrix of appropriate dimensions; or it can be left-divided by a non-singular constant matrix of appropriate dimension.

Numerous other combinations are possible, of course. The use of the exponentiation
operators `^` `.^` are somewhat limited;
see the definitions of `power` in *Nonlinear* below.

Matlab’s basic matrix manipulation and arithmetic operations have been extended to work with CVX expressions as well, including:

- Concatenation:
`[ A, B ; C, D ]` - Indexing:
`x(n+1:end)`,`X([3,4],:)`,*etc.* - Indexed assignment, including deletion:
`y(2:4) = 1`,`Z(1:4,:) = []`,*etc.* - Transpose and conjugate transpose:
`Z.'`,`y'`

A number of Matlab’s basic linear and bilinear functions either work automatically
with `cvx` expressions or have been extended to do so, including:
`conj`, `conv`, `cumsum`, `diag`, `dot`,
`find`, `fliplr`, `flipud`, `flipdim`,
`horzcat`, `hankel`, `ipermute`, `kron`, `mean`,
`permute`, `repmat`, `reshape`, `rot90`,
`sparse`, `sum`, `trace`, `tril`, `triu`,
`toeplitz`, `vertcat`.

Most should behave identically with CVX expressions as they do with
numeric expressions. Those that perform some sort of summation, such as
`cumsum`, `sum`, or multiplication, such as `conv`, `dot` or
`kron`, can only be used in accordance with the disciplined convex
programming rules. For example, `kron(X,Y)` is valid only if either
`X` or `Y` is constant; and `trace(Z)` is valid only if the
elements along the diagonal have the same curvature.

`abs`- absolute value for real and complex arrays. Convex.
- †
`exp` - exponential. Convex and nondecreasing.
- †
`log` - logarithm. Concave and nondecreasing.
`max`- maximum. Convex and nondecreasing.
`min`- minimum. Concave and nondecreasing.
`norm`norms for real and complex vectors and matrices. Convex. Thus function follows the Matlab conventions closely. Thus the one-argument version

`norm(x)`computes the 2-norm for vectors, and the 2-norm (maximum singular value) for matrices. The two-argument version`norm(x,p)`is supported as follows:- ‡ For vectors, all values \(p\geq 1\) are accepted.
- For matrices,
`p`must be`1`,`2`,`Inf`, or`'Fro'`.

`polyval`polynomial evaluation.

`polyval(p,x)`, where`p`is a vector of length`n`, computesp(1) * x.^(n-1) + p(2) * x.^(n-2) + ... + p(n-1) * x + p(n)

This function can be used in CVX in two ways:

- If
`p`is a variable and`x`is a constant, then`polyval(x,p)`computes a linear combination of the elements of`p`. The combination must satisfy the DCP rules for addition and scaling. - If
`p`is a constant and`x`is a variable, then`polyval(x,p)`constructs a polynomial function of the variable`x`. The polynomial must be affine, convex, or concave, and`x`must be real and affine.

- If
- ‡
`power(x,p)` `x^p`and`x.^p`, where`x`is a real variable and and`p`is a real constant. For`x^p`, both`x`and`p`must be scalars. Only those values of`p`which can reasonably and unambiguously interpreted as convex or concave are accepted:- \(p=0\). Constant.
`x.^p`is treated as identically 1. - \(0 < p < 1\). Concave. The argument \(x\) must be concave (or affine), and is implicitly constrained to be nonnegative.
- \(p = 1\). Affine.
`x.^p`is simply`x`. - \(p \in \{2,4,6,8,...\}\). Convex. Argument \(x\) must be affine.
- \(p > 1\), \(p\not\in\{2,3,4,5,...\}\). Convex. Argument \(x\) must be affine, and is implicitly constrained to be nonnegative.

Negative and odd integral values of \(p\) are not permitted, but see the functions

`pow_p`,`pow_pos`, and`pow_abs`in the next section for useful alternatives.- \(p=0\). Constant.
- †
`power(p,x)` `p.^x`and`p^x`, where`p`is a real constant and`x`is a real variable. For`p^x`, both`p`and`x`must be scalars. Valid values of`p`include:- \(p \in \{0,1\}\). Constant.
- \(0 < p < 1\). Convex and nonincreasing;
`x`must be concave. - \(p > 1\). Convex and nondecreasing;
`x`must be convex.

Negative values of

`p`are not permitted.`std`- Standard deviation. Convex.
`sqrt`- Square root. Implicitly constrains its argument to be nonnegative. Concave and nondecreasing.
`var`- Variance. Convex.

Even though these functions were developed specifically for CVX, they work outside of a CVX specification as well, when supplied with numeric arguments.

`avg_abs_dev`- The average absolute deviation about the mean \(\mu(x)\) of \(x\). Convex.
- \[f_{\text{aad}}(x) = \frac{1}{n} \sum_{i=1}^n |x_i-\mu(x)| = \frac{1}{n} \sum_{i=1}^n \left| x_i - {\textstyle\frac{1}{n}\sum_i x_i}\right| = \frac{1}{n}\left\| (I-\tfrac{1}{n}\textbf{1}\textbf{1}^T)x \right\|_1.\]

`avg_abs_dev_med`- The average absolute deviation about the median \(\mathop{\textrm{m}}(x)\) of \(x\). Convex.
- \[f_{\text{aadm}}(x) = \frac{1}{n} \sum_{i=1}^n |x_i-\mathop{\textrm{m}}(x)| = \inf_y \frac{1}{n} \sum_{i=1}^n |x_i-y|\]

`berhu(x,M)`- The reversed Huber function (hence, Berhu), defined as
- \[\begin{split}f_{\text{berhu}}(x,M) \triangleq \begin{cases} |x| & |x| \leq M \\ (|x|^2+M^2)/2M & |x| \geq M \end{cases}\end{split}\]

Convex. If \(M\) is omitted, \(M=1\) is assumed; but if supplied, it must be a positive constant. Also callable with three arguments as

`berhu(x,M,t)`, which computes`t+t*berhu(x/t,M)`, useful for concomitant scale estimation (see [Owen06]).`det_inv`- determinant of inverse of a symmetric (or Hermitian) positive
definite matrix, \(\det X^{-1}\), which is the same as the
product of the inverses of the eigenvalues. When used inside a
CVX specification,
`det_inv`constrains the matrix to be symmetric (if real) or Hermitian (if complex) and positive semidefinite. When used with numerical arguments,`det_inv`returns`+Inf`if these constraints are not met. Convex. `det_rootn`- \(n\)-th root of the determinant of a semidefinite matrix,
\((\det X)^{1/n}\). When used inside a CVX specification,
`det_rootn`constrains the matrix to be symmetric (if real) or Hermitian (if complex) and positive semidefinite. When used with numerical arguments,`det_rootn`returns`-Inf`if these constraints are not met. Concave. `det_root2n`- the \(2n\)-th root of the determinant of a semidefinite matrix;
*i.e.*,`det_root2n(X)=sqrt(det_rootn(X))`. Concave. Maintained solely for back-compatibility purposes. - †
`entr` - the elementwise entropy function:
`entr(x)=-x.*log(x)`. Concave. Returns`-Inf`when called with a constant argument that has a negative entry. `geo_mean`- the geometric mean of a vector,
\(\left( \prod_{k=1}^n x_k \right)^{1/n}\). When used inside a
CVX specification,
`geo_mean`constrains the elements of the vector to be nonnegative. When used with numerical arguments,`geo_mean`returns`-Inf`if any element is negative. Concave and increasing. `huber(x,M)`- The Huber function, defined as
- \[\begin{split}f_{\text{huber}}(x,M) \triangleq \begin{cases} |x|^2 & |x| \leq M \\ 2M|x|-M^2 & |x| \geq M \end{cases}\end{split}\]

Convex. If $x$ is a vector or array, the function is applied on an elementwise basis. If $M$ is omitted, then $M=1$ is assumed; but if it supplied, it must be a positive constant. Also callable as

`huber(x,M,t)`, which computes`t+t*huber(x/t,M)`, useful for concomitant scale estimation (see [Owen06]).`huber_circ(x,M)`- The circularly symmetric Huber function, defined as
- \[\begin{split}f_{\text{huber\_circ}}(x,M) \triangleq \begin{cases} \|x\|_2^2 & \|x\|_2 \leq M \\ 2M\|x\|_2-M^2 & \|x\|_2 \geq M \end{cases}\end{split}\]

Convex. Same (and implemented) as

`huber_pos(norm(x),M)`.`huber_pos(x,M)`- The same as the Huber function for nonnegative
`x`; zero for negative`x`. Convex and nondecreasing. `inv_pos`- The inverse of the positive portion, \(1/\max\{x,0\}\). Inside CVX specification, imposes constraint that its argument is positive. Outside CVX specification, returns \(+\infty\) if \(x\leq 0\). Convex and decreasing.
- †
`kl_div` Kullback-Leibler distance:

\[\begin{split}f_{\text{kl}}(x,y) \triangleq \begin{cases} x\log(x/y)-x+y & x,y>0 \\ 0 & x=y=0 \\ +\infty & \text{otherwise} \end{cases}\end{split}\]Convex. Outside CVX specification, returns \(+\infty\) if arguments aren’t in the domain.

`lambda_max`- maximum eigenvalue of a real symmetric or complex Hermitian matrix. Inside CVX, imposes constraint that its argument is symmetric (if real) or Hermitian (if complex). Convex.
`lambda_min`- minimum eigenvalue of a real symmetric or complex Hermitian matrix. Inside CVX, imposes constraint that its argument is symmetric (if real) or Hermitian (if complex). Concave.
`lambda_sum_largest(X,k)`- sum of the largest \(k\) values of a real symmetric or complex Hermitian matrix. Inside CVX, imposes constraint that its argument is symmetric (if real) or Hermitian (if complex). Convex.
`lambda_sum_smallest(X,k)`- sum of the smallest \(k\) values of a real symmetric or complex Hermitian matrix. Inside CVX, imposes constraint that its argument is symmetric (if real) or Hermitian (if complex). Concave.
`log_det`- log of determinant of a positive definite matrix,
\(\log \det(X)\). When used inside a CVX specification,
`log_det`constrains its argument to be symmetric (if real) or Hermitian (if complex) and positive definite. With numerical argument,`log_det`returns`-Inf`if these constraints are not met. Concave. - ‡
`log_normcdf(x)` - logarithm of cumulative distribution function of standard normal random variable. Concave and increasing. The current implementation is a fairly crude SDP-representable approximation, with modest accuracy over the interval \([-4,4]\); we intend to replace it with a much better approximation at some point.
- †
`log_prod(x)` - \(\log\prod_i x_i\) if when \(x\) is positive; \(-\infty\) otherwise.
Concave and nonincreasing. Equivalent to
`sum_log(x)`. - †
`log_sum_exp(x)` - the logarithm of the sum of the elementwise exponentials of
`x`. Convex and nondecreasing. `logsumexp_sdp`- a polynomial approximation to the log-sum-exp function with global absolute accuracy. This can be used to estimate the log-sum-exp function without using the successive approximation method.
`matrix_frac(x,Y)`- matrix fractional function, \(x^TY^{-1}x\). In CVX, imposes constraint that \(Y\) is symmetric (or Hermitian) and positive definite; outside CVX, returns \(+\infty\) unless \(Y=Y^T\succ 0\). Convex.
`norm_largest(x,k)`- For real and complex vectors, returns the sum of the largest
`k`*magnitudes*in the vector`x`. Convex. `norm_nuc(X)`- The sum of the singular values of a real or complex matrix
`X`. (This is the dual of the usual spectral matrix norm,*i.e.*, the largest singular value.) Convex. - ‡
`norms(x,p,dim)`,`norms_largest(x,k,dim)` - Computes
*vector*norms along a specified dimension of a matrix or N-d array. Useful for sum-of-norms and max-of-norms problems. Convex. `poly_env(p,x)`Computes the value of the

*convex or concave envelope*of the polynomial described by`p`(in the`polyval`sense).`p`must be a real constant vector whose length`n`is 0, 1, 2, 3, or some other*odd*length; and`x`must be real and affine. The sign of the first nonzero element of`p`determines whether a convex (positive) or concave (negative) envelope is constructed. For example, consider the function \(p(x)\triangleq (x^2-1)^2=x^4-2x^2+1\), depicted along with its convex envelope in the figure below.The two coincide when \(|x|\geq 1\), but deviate when \(|x|<1\). Attempting to call

`polyval([1,0,2,0,1],x)`in a CVX model would yield an error, but a call to`poly_env([1,0,2,0,1],x)`yields a valid representation of the envelope. For convex or concave polynomials, this function produces the same result as`polyval`.`pos(x)`- \(\max\{x,0\}\), for real \(x\). Convex and increasing.
- ‡
`pow_abs(x,p)` - \(|x|^p\) for \(x\in\mathbf{R}\) or \(x\in\mathbf{C}\) and \(p\geq 1\). Convex.
- ‡
`pow_pos(x,p)` - \(\max\{x,0\}^p\) for \(x\in\mathbf{R}\) and \(p\geq 1\). Convex and nondecreasing.
- ‡
`pow_p(x,p)` for \(x\in\mathbf{R}\) and real constant \(p\), computes nonnegative convex and concave branches of the power function:

\[\begin{split}\begin{array}{ccl} p\leq 0 & f_p(x) \triangleq \begin{cases} x^p & x > 0 \\ +\infty & x \leq 0 \end{cases} & \text{convex, nonincreasing} \\ 0 < p \leq 1 & f_p(x) \triangleq \begin{cases} x^p & x \geq 0 \\ -\infty & x < 0 \end{cases} & \text{concave, nondecreasing} \\ p \geq 1 & f_p(x) \triangleq \begin{cases} x^p & x \geq 0 \\ +\infty & x < 0 \end{cases} & \text{convex, nonmonotonic} \end{array}\end{split}\]`prod_inv(x)`- \(\prod_i x_i^{-1}\) when \(x\) is positive; \(+\infty\) otherwise. Convex and nonincreasing.
`quad_form(x,P)`- \(x^TPx\) for real \(x\) and symmetric \(P\), and \(x^HPx\) for complex \(x\) and Hermitian \(P\). Convex in \(x\) for \(P\) constant and positive semidefinite; concave in \(x\) for \(P\) constant and negative semidefinite.

Note

Quadratic functions such as `quad_form`, `sum_square` can often be replaced
by the `norm` function without sacrificing equivalence. For numerical reasons,
this alternate formulation is *preferred*. Please see *Eliminating quadratic forms* for
more information.

`quad_over_lin(x,y)`- \(x^Tx/y\) for \(x \in \mathbf{R}^n\), \(y >0\); for \(x \in \mathbf{C}^n\), \(y>0\), \(x^*x/y\). In CVX specification, adds constraint that \(y>0\). Outside CVX specification, returns \(+\infty\) if \(y\leq 0\). Convex, and decreasing in \(y\).
`quad_pos_over_lin(x,y)``sum_square_pos( x )/y`for \(x\in\mathbf{R}^n\), \(y>0\). Convex, increasing in \(x\), and decreasing in \(y\).- †
`rel_entr(x)` - Scalar relative entropy;
`rel_entr(x,y)=x.*log(x/y)`. Convex. `sigma_max`- maximum singular value of real or complex matrix. Same as
`norm`. Convex. `square`- \(x^2\) for \(x \in \mathbf{R}\). Convex.
`square_abs`- \(|x|^2\) for \(x\in\mathbf{R}\) or \(x\in\mathbf{C}\).
`square_pos`- \(\max\{x,0\}^2\) for \(x\in\mathbf{R}\). Convex and increasing.
`sum_largest(x,k)`- sum of the largest \(k\) values, for real vector \(x\). Convex and nondecreasing.
- †
`sum_log(x)` - \(\sum_i\log(x_i)\) when \(x\) is positive; \(-\infty\) otherwise. Concave and nondecreasing.
`sum_smallest(x,k)`- sum of the smallest \(k\) values,
*i.e.*, equivalent to`-sum_largest(-x,k)`. Concave and nondecreasing. `sum_square`- Equivalent to
`sum(square(x))`, but more efficient. Convex. Works only for real values. `sum_square_abs`- Equivalent to
`sum(square_abs(x))`, but more efficient. Convex. `sum_square_pos`- Equivalent to
`sum(square_pos(x))`, but more efficient. Works only for real values. Convex and increasing. `trace_inv(X)`- trace of the inverse of an SPD matrix
`X`, which is the same as the sum of the inverses of the eigenvalues. Convex. Outside of CVX, returns`+Inf`if argument is not positive definite. `trace_sqrtm(X)`- trace of the matrix square root of a positive semidefinite matrix
`X`. which is the same as the sum of the squareroots of the eigenvalues. Concave. Outside of CVX, returns`+Inf`if argument is not positive semidefinite.

CVX currently supports the following sets; in each case, `n` is a
positive integer constant.

`nonnegative(n)`- \[R^n_+ \triangleq \left\{\,x\in\mathbf{R}^n\,~|~\,x_i\geq 0,~i=1,2,\dots,n\,\right\}\]
`simplex(n)`- \[R^n_{1+} \triangleq \left\{\,x\in\mathbf{R}^n\,~|~\,x_i\geq 0,~i=1,2,\dots,n,~\textstyle\sum_ix_i=1\,\right\}\]
`lorentz(n)`- \[\mathbf{Q}^n \triangleq \left\{\,(x,y)\in\mathbf{R}^n\times\mathbf{R}\,~|~\,\|x\|_2\leq y\,\right\}\]
`rotated_lorentz(n)`- \[\mathbf{Q}^n_r \triangleq \left\{\,(x,y,z)\in\mathbf{R}^n\times\mathbf{R}\times\mathbf{R}\,~|~\,\|x\|_2\leq \sqrt{yz},~y,z\geq 0\,\right\}\]
`complex_lorentz(n)`- \[\mathbf{Q}^n_c \triangleq \left\{\,(x,y)\in\mathbf{C}^n\times\mathbf{R}\,~|~\,\|x\|_2\leq y\,\right\}\]
`rotated_complex_lorentz(n)`- \[\mathbf{Q}^n_{rc} \triangleq \left\{\,(x,y,z)\in\mathbf{C}^n\times\mathbf{R}\times\mathbf{R}\,~|~\,\|x\|_2\leq \sqrt{yz},~y,z\geq 0\,\right\}\]
`semidefinite(n)`- \[\mathbf{S}^n_+ \triangleq \left\{\,X\in\mathbf{R}^{n\times n}\,~|~\,X=X^T,~X\succeq 0\,\right\}\]
`hermitian_semidefinite(n)`- \[\mathbf{H}^n_+ \triangleq \left\{\,Z\in\mathbf{C}^{n\times n}\,~|~\,Z=Z^H,~Z\succeq 0\,\right\}\]
`nonneg_poly_coeffs(n)`The cone of all coefficients of nonnegative polynomials of degree \(n\); \(n\) must be even:

\[\mathbf{P}_{+,n} \triangleq \left\{\,p\in\mathbf{R}^n[n+1]\,~|~\,\sum_{i=0}^n p_{i+1} x^{n-i} \geq 0 ~ \forall x\in\mathbf{R}\,\right\}\]`convex_poly_coeffs(n)`The cone of all coefficients of convex polynomials of degree \(n\); \(n\) must be even:

\[\mathbf{P}_{+,n} \triangleq \left\{\,p\in\mathbf{R}^n[n+1]\,~|~\,\sum_{i=0}^{n-2} (n-i)(n-i-1) p_{i+1} x^{n-i-2} \geq 0 ~ \forall x\in\mathbf{R}\,\right\}\]`exp_cone`- \[\begin{split}\mathbf{E} \triangleq \text{cl}\left\{\,(x,y,z)\in\mathbf{R}\times\mathbf{R}\times\mathbf{R}\,~|~\,y>0,~ye^{x/y}\leq z\,\right\}\end{split}\]
`geo_mean_cone(n)`- \[\mathbf{G}_n \triangleq \text{cl}\left\{\,(x,y)\in\mathbf{R}^n\times\mathbf{R}^n\times\mathbf{R}^n\,~|~\,x\geq 0,~(\prod_{i=1}^n x_i)^{1/n} \geq y\,\right\}\]

`cvx_begin`- Begins a new CVX model. If a model is already in progress, it will issue a warning
and clear it. See
*cvx_begin and cvx_end*for a full description, including the modifying keywords that control solver output, SDP mode, GDP mode, etc. `cvx_clear`- Clears any model being constructed. Useful when an error has been made and it is
necessary to start from the beginning. Whereas
`cvx_begin`issues a warning if called with a model in progress,`cvx_clear`is silent. `cvx_end`- Signals the end of a CVX model. In typical use, this instructs CVX to begin the solution process.
See
*cvx_begin and cvx_end*. `cvx_expert`- Controls the issuance of warnings when models requiring the use of successive
approximation are employed; see
*The successive approximation method*more details. `cvx_power_warning`- Controls if and when CVX issues warnings during the construction of models involving
rational power functions (i.e.,
`x^p`, where`x`is a variable and`p`is a constant); see*Power functions and p-norms*. `cvx_precision`- Controls solver precision; see
*Controlling precision*. `cvx_quiet`- Enables or disables screen output during the solution process; see
*Controlling screen output*. Also see*cvx_begin and cvx_end*for the newer, preferred syntax`cvx_begin quiet`. `cvx_save_prefs`- Saves the current states for
`cvx_expert`,`cvx_power_warning`,`cvx_precision`, and`cvx_solver`to disk, so that their values are retained when quitting and re-starting MATLAB. The file is saved in MATLAB’s preference directory, which can be located by typing the`prefdir`command. `cvx_setup`- The setup script used to install and configure CVX; see
*Installation*. `cvx_solver`- Selects the solver to be employed when solving CVX models; see
*Selecting a solver*. `cvx_solver_settings`- Allows the user to deliver advanced, solver-specific settings to the solver that CVX
does not otherwise support; see
*Advanced solver settings*. `cvx_version`- Prints information about the current versions of CVX, Matlab, and the operating system. When submitting bug reports, please include the output of this command.
`cvx_where`- Returns the directory where CVX is installed.
`dual variable`,`dual variables`- Creates one or more dual variables to be connected to constraints in the current model;
see
*Dual variables*. `expression`,`expressions`- Creates one or more expression holders; see
*Assignment and expression holders*. `maximise`,`maximize`- Specifies a maximization objective; see
*Objective functions*. `minimise`,`minimize`- Specifies a minimization objective; see
*Objective functions*. `variable`,`variables`- Creates one or more variables for use in the current CVX model; see
*Variables*.