Tuesday, December 3, 2019

macOS 10.15 (Catalina) MOSEK installation

MacOS 10.15 (Catalina) introduced stricter checks for unverified applications downloaded from the Internet. When you try to run MOSEK you may experience security exceptions such as:

"mosekopt.mexmaci64" cannot be opened because the developer cannot be verified. macOS cannot verify that this app is free from malware.


"mosek" cannot be opened because the developer cannot be verified. macOS cannot verify that this app is free from malware.


‘.../mosekopt.mexmaci64’:  dlopen(.../mosekopt.mexmaci64,6): no suitable image found. Did find:  mosekopt.mexmaci64:  code signature in (.../mosekopt.mexmaci64)  not valid for use in process using Library Validation:  library load disallowed by system policy.
or other errors of this kind.

Solution 1

Do not download the MOSEK package through a browser, use a command-line download tool such as curl instead:

curl https://download.mosek.com/stable/9.1.9/mosektoolsosx64x86.tar.bz2 -o mosektoolsosx64x86.tar.bz2

(replace 9.1.9 with the version you want to download). Continue with normal installation.

Solution 2

After you download the distribution remove the quarantine attribute from the downloaded file:

xattr -d com.apple.quarantine mosektoolsosx64x86.tar.bz2

Administrator permissions may be necessary. Continue with normal installation.

If you already unpacked the distribution remove the attribute recursively from the whole folder with MOSEK installation

xattr -dr com.apple.quarantine mosek

assuming you are in the folder which contains the unpacked mosek folder. Administrator permissions may be necessary. 

Monday, November 4, 2019

FlexLM update

In version 9.1.5 we upgraded the FlexLM license manager to from version 11.14 to 11.16 because of a potential security issue reported by a user.

This upgrade breaks client-server compatibility. That means that if you installed a MOSEK 9.1.5+ client, you are using a floating license, and you are getting an error such as

Version of vendor daemon is too old

then you should upgrade the license server to one that also comes from a 9.1.5+ distribution of MOSEK. It will be able to serve the new and all older client versions.

See also The licensing FAQ.

Tuesday, October 22, 2019

IPCO 2020

Just like this year, MOSEK is one of the sponsors of the upcoming IPCO 2020, The 21st Conference on Integer Programming and Combinatorial Optimization, taking place at LSE in London 8-10 June 2020.

Submission deadline is November 29, 2019.

Wednesday, October 2, 2019

MOSEK 9.1 and deprecation plans for version 10

Today we released MOSEK 9.1, which is a direct continuation of version 9.0 with no changes to the optimizer or the existing API.

It introduces one new feature: the possibility to use the remote OptServer from Fusion via Model.solve(server, port).

We take this opportunity to announce that the following will be dropped in the next major release:

  • support for Python 2.7 on all platforms,
  • support for Java on Windows 32bit,
  • support for all versions of Python on Windows 32bit.


Tuesday, August 20, 2019

Dropping 32bit Windows

Due to falling interest in our distribution for Windows 32bit platform and the fact that more and more software vendors are dropping it, we also plan to drop the Windows 32bit package in one of the minor or major releases in near future.

We will most likely continue to provide a minimal package consisting of the low-level C API and shared libraries.

If you rely on MOSEK for 32bit Windows please let us know at support@mosek.com so that we can best accommodate your needs.

Of course nothing changes on 64bit Windows.


Friday, July 19, 2019


Joachim Dahl speaks at ICCOPT 2019 in Berlin:

Title: A primal-dual algorithm for exponential-cone optimization

Session: Advances in Algorithms for Conic Optimization

Abstract: We discuss a primal-dual algorithm for exponential-cone optimization. The algorithm is based on primal-dual scaling proposed by Tuncel, with similar properties as Nesterov-Todd scalings for symmetric cones. We also discuss a new corrector for nonsymmetric cones, which considerably reduces the number of iterations required to solve a number of test problems. This algorithm has recently been developed for version 9 of the MOSEK software.

Time: Thursday 08/08/2019, 10:45, room H2032

Tuesday, June 4, 2019

Logarithmic mean temperature difference requires yet another cone?

The logarithmic mean temperature difference
$$LMTD(x,y) = \frac{x-y}{\ln(x/y)}$$

can be used to model heat exchange, e.g., when optimizing the reuse of excess process heat onsite at industrial facilities. In the heat exchanger network synthesis problem due to Yee and Grossmann (link goes to a recent paper by Mistry and Misener), we find that its reciprocal raised to some power, that is
$$RecLMTD^\beta(x,y) = \left(\frac{\ln(x/y)}{x-y}\right)^\beta,$$

can be extracted as a separately contributing term in the objective function. Capitalizing on the convexity of this term on $(x,y) \in \mathbb{R}_+^2$, for all considered $\beta \geq 0$, this leads to better performance when solving the otherwise nonconvex problem as argued in the paper.

A challenge to find the conic reformulation of this function was posed under the Oberwolfach Workshop on Mixed-Integer Nonlinear Optimization (2019) and we accepted. Of course, this is trivial if no restrictions are put on the set of cones as one may just define
$$\mathcal{K} = \mathrm{cl}\{(t,s,x,y) \in \mathbb{R}_{++}^4 : t \geq s \cdot RecLMTD^\beta(x/s, y/s) \}$$
and call it a day. This cone is nonempty, closed and convex and hence obeys $K = (K^{*})^{*}$ as well as all the usual properties of conic duality. Computationally, however, the cone is not particularly desirable and we can do better with a bit of reformulation:
$$\begin{array}{lll} t \geq \left(\frac{\ln(x/y)}{x-y}\right)^\beta,\\t \geq \left(\frac{\ln(u/y + 1)}{u}\right)^\beta,& u = x-y,\\ y \geq \frac{u}{\exp(ut^{1/\beta})-1},& u = x-y,\\ y \geq \frac{u}{\exp(u/s)-1},& u = x-y, & s \geq t^{-1/\beta},\\ \end{array}$$
where I substitute in the first step, rewrite assuming either $u > 0$ or $u <0$ (both leads to the same) in the second, and extract a power cone representable subexpression in the third. This means that the representation problem of $RecLMTD^\beta$ have been reduced to the representation problem of
$$\mathcal{K} = \mathrm{cl}\left\{(t,s,x) \in \mathbb{R}_{+}^2 \times \mathbb{R}_{++} : t \geq \frac{x}{\exp(x/s)-1} \right\},$$
which, just like the quadratic, power and exponential cones, is defined as the epigraph of the perspective of a univarite convex function; in this case $\frac{x}{\exp(x)-1}$. Whether this cone can be written in terms of the others, or has potential for computationally efficient implementations itself, remains open. We invite anyone interested in barrier functions and interior-point algorithms to take a crack at it.

What is possible right now?
The logarithmic mean is bounded from above and below by, respectively, the arithmetic and the geometric mean, both of which are representable in MOSEK. There is also Chen's approximation which corresponds to a power cone. Finally it can be observed that $\log(\exp(-x)+1)$ is a fairly good underestimator of $\frac{x}{\exp(x)-1}$. This means that the set described by the inequality from above,
$$y \geq \frac{u}{\exp(u/s)-1},$$

can be outer approximated by
$$y \geq s \log(\exp(-u/s)+1),$$

which is representable in MOSEK as the homogenization of the softplus function. Beware that the usefulness of any of these alternatives is unknown to the author.