reconfigure .gitpod.Dockerfile; add test_quasicvx2.py

Author: luk036

Diff for: .gitpod.Dockerfile

@@ -1,46 +1,79 @@
 FROM gitpod/workspace-full
 
 USER root
-
 # Install util tools.
+
 RUN apt-get update \
  && apt-get install -y \
   apt-utils \
   aria2 \
-  git \
-  less \
-  lcov \
-  neofetch \
-  neovim \
+# utilities (not ripgrep, gh) \
   asciinema \
-  tmux \
+  bat \
+  byobu \
+  curl \
+  elinks \
+  fd-find \
+  fish \
+  mdp \
+  ncdu \
+  neofetch \
+  patat \
+  pkg-config \
+  ranger \
   w3m \
-  wget
-
-RUN pip3 install --upgrade pip \
-    && pip3 install \
-	decorator \
-	networkx \
-	numpy \
-	numexpr \
-	scipy \
-	pre-commit \
-	mypy \
-	codecov \
-	coverage \
-	hypothesis \
-	pytest \
-	pytest-cov \
-	pytest-benchmark \
-	jupyter \
-	jupyterlab \
-	matplotlib
-
-USER gitpod
-
-# Install custom tools, runtime, etc. using apt-get
-# For example, the command below would install "bastet" - a command line tetris clone:
-#
-# RUN sudo apt-get -q update && #     sudo apt-get install -yq bastet && #     sudo rm -rf /var/lib/apt/lists/*
-#
-# More information: https://www.gitpod.io/docs/42_config_docker/
+# just for fun (not cmatrix) \
+  cowsay \
+  figlet \
+  fortune \
+  toilet \
+  tty-clock
+
+RUN mkdir -p /workspace/data \
+    && chown -R gitpod:gitpod /workspace/data
+
+RUN mkdir /home/gitpod/.conda
+# Install conda
+RUN wget --quiet https://repo.anaconda.com/miniconda/Miniconda3-latest-Linux-x86_64.sh -O ~/miniconda.sh && \
+    /bin/bash ~/miniconda.sh -b -p /opt/conda && \
+    rm ~/miniconda.sh && \
+    ln -s /opt/conda/etc/profile.d/conda.sh /etc/profile.d/conda.sh && \
+    echo ". /opt/conda/etc/profile.d/conda.sh" >> ~/.bashrc && \
+    echo "conda activate base" >> ~/.bashrc
+
+RUN /opt/conda/bin/conda config --set always_yes yes --set changeps1 no \
+    && /opt/conda/bin/conda update -q conda \
+    && /opt/conda/bin/conda info -a
+
+RUN /opt/conda/bin/conda install -y -c conda-forge \
+    pandoc-crossref \
+    pandoc
+
+RUN /opt/conda/bin/pip install \
+    codecov \
+    coverage \
+    coveralls \
+    cvxpy \
+    jupyter \
+    jupyterlab \
+    matplotlib \
+    networkx \
+    numexpr \
+    numpy \
+    scipy \
+    flake8 \
+    mypy \
+    pre-commit \
+    pyqt5 \
+    pytest-benchmark \
+    pytest-cov \
+    pytest \
+    yapf \
+    lolcat
+
+RUN chown -R gitpod:gitpod /opt/conda \
+    && chmod -R 777 /opt/conda \
+    && chown -R gitpod:gitpod /home/gitpod/.conda \
+    && chmod -R 777 /home/gitpod/.conda
+
+RUN apt-get clean && rm -rf /var/cache/apt/* && rm -rf /var/lib/apt/lists/* && rm -rf /tmp/*

Diff for: tests/test_quasicvx2.py

@@ -0,0 +1,81 @@
+# -*- coding: utf-8 -*-
+from __future__ import print_function
+
+import math
+
+from pytest import approx
+
+import numpy as np
+
+from ellpy.cutting_plane import CUTStatus, cutting_plane_dc
+from ellpy.ell import ell
+
+
+def my_quasicvx_oracle(z, t: float):
+    """[summary]
+
+    Arguments:
+        z ([type]): [description]
+        t (float): the best-so-far optimal value
+
+    Returns:
+        [type]: [description]
+    """
+    x, y = z
+
+    # constraint 1: exp(x) <= y
+    tmp = math.exp(x)
+    fj = tmp - y
+    if fj > 0:
+        return (np.array([tmp, -1.]), fj), None
+
+    # constraint 2: y > 0
+    if y <= 0.:
+        return (np.array([0., -1.]), -y), None
+
+    # constraint 3: x > 0
+    if x <= 0.:
+        return (np.array([-1., 0.]), -x), None
+
+    # objective: minimize -sqrt(x) / y
+    tmp2 = math.sqrt(x)
+    fj = -tmp2 - t * y
+    if fj < 0.:  # feasible
+        t = -tmp2 / y
+        return (np.array([-0.5/tmp2, -t]), 0), t
+
+    return (np.array([-0.5/tmp2, -t]), fj), None
+
+
+def test_case_feasible():
+    """[summary]
+    """
+    x0 = np.array([1., 1.])  # initial x0
+    E = ell(10., x0)
+    P = my_quasicvx_oracle
+    xb, fb, ell_info = cutting_plane_dc(P, E, 0.)
+    assert ell_info.feasible
+    assert fb == approx(-0.4288673396685956)
+    assert xb[0] == approx(0.501315956)
+    assert xb[1] == approx(1.650886769)
+
+
+def test_case_infeasible1():
+    """[summary]
+    """
+    x0 = np.array([100., 100.])  # wrong initial guess,
+    E = ell(10., x0)  # or ellipsoid is too small
+    P = my_quasicvx_oracle
+    _, _, ell_info = cutting_plane_dc(P, E, 0.)
+    assert not ell_info.feasible
+    assert ell_info.status == CUTStatus.nosoln  # no sol'n
+
+
+def test_case_infeasible2():
+    """[summary]
+    """
+    x0 = np.array([1., 1.])  # initial x0
+    E = ell(10., x0)
+    P = my_quasicvx_oracle
+    _, _, ell_info = cutting_plane_dc(P, E, -100)  # wrong initial best-so-far
+    assert not ell_info.feasible