Test run_simulation.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "from run_sweep_expander import *\n",
    "\n",
    "%load_ext autoreload"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [],
   "source": [
    "n = 1000\n",
    "t = 4\n",
    "r = 3\n",
    "sigma = 0.0\n",
    "r_fit = 32\n",
    "#const = 12\n",
    "d = 12"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Help on function run_simulation in module run_sweep_expander:\n",
      "\n",
      "run_simulation(n, t, r, sigma, r_fit, rep, d=10, max_iter=None, inner_max_iter=10, tol=1e-10, alg='fro', verbosity=0, kappa=100, beta=1, epsilon=0.01, delta=None)\n",
      "\n"
     ]
    }
   ],
   "source": [
    "help(run_simulation)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "1.728e-06"
      ]
     },
     "execution_count": 16,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "(d/n)**(t-1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Running a simulation: n = 1000, t = 4, r = 3, sigma = 0.000000, r_fit = 32, alg=fro\n",
      "\n",
      "max_qnorm_ub_true = 3.941e+01\n",
      "expander degree = 12, sampling 1.73e-04%\n",
      "Initial cost: 7.697e+07\n",
      "Initial qnorm_ub: 8.482e-01\n",
      "|| r || = 9.938e-01, delta = 5.000e-02\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 1.062e+00\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.858e+00\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 6.283e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 7.928e+01\n",
      "\n",
      "=============================\n",
      "Iteration 0 complete\n",
      "\n",
      "\n",
      "scaled || r || = 9.928e-01, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 3.093e+01\n",
      "Cost function:         7.680e+07\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 1.944e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.252e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.479e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 6.667e+01\n",
      "\n",
      "=============================\n",
      "Iteration 1 complete\n",
      "\n",
      "\n",
      "scaled || r || = 9.871e-01, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.504e+02\n",
      "Cost function:         7.587e+07\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.639e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.813e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.539e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 5.736e+01\n",
      "\n",
      "=============================\n",
      "Iteration 2 complete\n",
      "\n",
      "\n",
      "scaled || r || = 9.807e-01, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.178e+02\n",
      "Cost function:         7.484e+07\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.419e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.103e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.924e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 4.617e+01\n",
      "\n",
      "=============================\n",
      "Iteration 3 complete\n",
      "\n",
      "\n",
      "scaled || r || = 9.718e-01, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 9.368e+01\n",
      "Cost function:         7.341e+07\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.159e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.150e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.126e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 4.312e+01\n",
      "\n",
      "=============================\n",
      "Iteration 4 complete\n",
      "\n",
      "\n",
      "scaled || r || = 9.472e-01, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 8.756e+01\n",
      "Cost function:         6.955e+07\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.778e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.836e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.500e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 5.394e+01\n",
      "\n",
      "=============================\n",
      "Iteration 5 complete\n",
      "\n",
      "\n",
      "scaled || r || = 8.425e-01, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.026e+02\n",
      "Cost function:         5.427e+07\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.022e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.626e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.589e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 5.871e+01\n",
      "\n",
      "=============================\n",
      "Iteration 6 complete\n",
      "\n",
      "\n",
      "scaled || r || = 6.013e-01, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.125e+02\n",
      "Cost function:         2.627e+07\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.570e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.705e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.694e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 7.090e+01\n",
      "\n",
      "=============================\n",
      "Iteration 7 complete\n",
      "\n",
      "\n",
      "scaled || r || = 5.221e-01, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.373e+02\n",
      "Cost function:         1.926e+07\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.196e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.403e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.955e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 9.790e+01\n",
      "\n",
      "=============================\n",
      "Iteration 8 complete\n",
      "\n",
      "\n",
      "scaled || r || = 4.604e-01, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.891e+02\n",
      "Cost function:         1.456e+07\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 5.427e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 4.472e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.065e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 1.300e+02\n",
      "\n",
      "=============================\n",
      "Iteration 9 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.813e-01, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.454e+02\n",
      "Cost function:         4.634e+06\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 6.799e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 5.211e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.793e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 1.385e+02\n",
      "\n",
      "=============================\n",
      "Iteration 10 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.404e-01, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.614e+02\n",
      "Cost function:         7.166e+05\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 7.759e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 5.979e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.807e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 1.554e+02\n",
      "\n",
      "=============================\n",
      "Iteration 11 complete\n",
      "\n",
      "\n",
      "scaled || r || = 9.997e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.944e+02\n",
      "Cost function:         2.266e+05\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 8.699e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 6.864e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.926e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 1.760e+02\n",
      "\n",
      "=============================\n",
      "Iteration 12 complete\n",
      "\n",
      "\n",
      "scaled || r || = 8.369e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 3.342e+02\n",
      "Cost function:         1.090e+05\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 9.644e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 7.405e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.902e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 1.850e+02\n",
      "\n",
      "=============================\n",
      "Iteration 13 complete\n",
      "\n",
      "\n",
      "scaled || r || = 7.450e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 3.518e+02\n",
      "Cost function:         6.271e+04\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 1.028e+02\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 7.814e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.855e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 1.909e+02\n",
      "\n",
      "=============================\n",
      "Iteration 14 complete\n",
      "\n",
      "\n",
      "scaled || r || = 6.809e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 3.634e+02\n",
      "Cost function:         3.908e+04\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 1.070e+02\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 8.225e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.824e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 1.966e+02\n",
      "\n",
      "=============================\n",
      "Iteration 15 complete\n",
      "\n",
      "\n",
      "scaled || r || = 6.321e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 3.748e+02\n",
      "Cost function:         2.586e+04\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 1.103e+02\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 8.587e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.822e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.020e+02\n",
      "\n",
      "=============================\n",
      "Iteration 16 complete\n",
      "\n",
      "\n",
      "scaled || r || = 5.957e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 3.855e+02\n",
      "Cost function:         1.868e+04\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 1.134e+02\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 8.819e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.811e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.050e+02\n",
      "\n",
      "=============================\n",
      "Iteration 17 complete\n",
      "\n",
      "\n",
      "scaled || r || = 5.631e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 3.927e+02\n",
      "Cost function:         1.420e+04\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 1.160e+02\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 8.975e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.797e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.070e+02\n",
      "\n",
      "=============================\n",
      "Iteration 18 complete\n",
      "\n",
      "\n",
      "scaled || r || = 5.369e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 3.978e+02\n",
      "Cost function:         1.191e+04\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 1.182e+02\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 9.104e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.780e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.076e+02\n",
      "\n",
      "=============================\n",
      "Iteration 19 complete\n",
      "\n",
      "\n",
      "scaled || r || = 5.134e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 4.004e+02\n",
      "Cost function:         1.087e+04\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 1.198e+02\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 9.130e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.723e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.008e+02\n",
      "\n",
      "=============================\n",
      "Iteration 20 complete\n",
      "\n",
      "\n",
      "scaled || r || = 4.903e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 3.891e+02\n",
      "Cost function:         1.054e+04\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 1.180e+02\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 8.998e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.677e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 1.945e+02\n",
      "\n",
      "=============================\n",
      "Iteration 21 complete\n",
      "\n",
      "\n",
      "scaled || r || = 4.704e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 3.786e+02\n",
      "Cost function:         1.032e+04\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 1.160e+02\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 8.822e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.628e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 1.871e+02\n",
      "\n",
      "=============================\n",
      "Iteration 22 complete\n",
      "\n",
      "\n",
      "scaled || r || = 4.510e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 3.658e+02\n",
      "Cost function:         1.011e+04\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 1.136e+02\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 8.621e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.596e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 1.810e+02\n",
      "\n",
      "=============================\n",
      "Iteration 23 complete\n",
      "\n",
      "\n",
      "scaled || r || = 4.335e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 3.552e+02\n",
      "Cost function:         9.929e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 1.117e+02\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 8.456e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.593e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 1.761e+02\n",
      "\n",
      "=============================\n",
      "Iteration 24 complete\n",
      "\n",
      "\n",
      "scaled || r || = 4.157e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 3.470e+02\n",
      "Cost function:         9.750e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 1.097e+02\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 8.258e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.594e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 1.712e+02\n",
      "\n",
      "=============================\n",
      "Iteration 25 complete\n",
      "\n",
      "\n",
      "scaled || r || = 3.999e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 3.386e+02\n",
      "Cost function:         9.593e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 1.077e+02\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 8.022e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.590e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 1.656e+02\n",
      "\n",
      "=============================\n",
      "Iteration 26 complete\n",
      "\n",
      "\n",
      "scaled || r || = 3.839e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 3.287e+02\n",
      "Cost function:         9.439e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 1.054e+02\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 7.787e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.594e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 1.604e+02\n",
      "\n",
      "=============================\n",
      "Iteration 27 complete\n",
      "\n",
      "\n",
      "scaled || r || = 3.698e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 3.194e+02\n",
      "Cost function:         9.305e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 1.030e+02\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 7.522e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.591e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 1.546e+02\n",
      "\n",
      "=============================\n",
      "Iteration 28 complete\n",
      "\n",
      "\n",
      "scaled || r || = 3.516e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 3.098e+02\n",
      "Cost function:         9.148e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 1.007e+02\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 7.316e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.612e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 1.499e+02\n",
      "\n",
      "=============================\n",
      "Iteration 29 complete\n",
      "\n",
      "\n",
      "scaled || r || = 3.342e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 3.019e+02\n",
      "Cost function:         9.000e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 9.828e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 7.107e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.635e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 1.452e+02\n",
      "\n",
      "=============================\n",
      "Iteration 30 complete\n",
      "\n",
      "\n",
      "scaled || r || = 3.174e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.943e+02\n",
      "Cost function:         8.860e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 9.590e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 6.898e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.658e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 1.404e+02\n",
      "\n",
      "=============================\n",
      "Iteration 31 complete\n",
      "\n",
      "\n",
      "scaled || r || = 3.021e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.886e+02\n",
      "Cost function:         8.735e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 9.415e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 6.740e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.702e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 1.356e+02\n",
      "\n",
      "=============================\n",
      "Iteration 32 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.867e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.827e+02\n",
      "Cost function:         8.613e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 9.230e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 6.579e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.744e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 1.325e+02\n",
      "\n",
      "=============================\n",
      "Iteration 33 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.727e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.798e+02\n",
      "Cost function:         8.508e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 9.146e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 6.492e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.782e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 1.294e+02\n",
      "\n",
      "=============================\n",
      "Iteration 34 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.598e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.765e+02\n",
      "Cost function:         8.413e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 9.053e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 6.401e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.817e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 1.262e+02\n",
      "\n",
      "=============================\n",
      "Iteration 35 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.473e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.731e+02\n",
      "Cost function:         8.322e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 8.954e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 6.309e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.849e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 1.232e+02\n",
      "\n",
      "=============================\n",
      "Iteration 36 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.358e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.694e+02\n",
      "Cost function:         8.239e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 8.850e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 6.216e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.879e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 1.202e+02\n",
      "\n",
      "=============================\n",
      "Iteration 37 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.249e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.658e+02\n",
      "Cost function:         8.161e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 8.745e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 6.128e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.909e+01\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 1.175e+02\n",
      "\n",
      "=============================\n",
      "Iteration 38 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.113e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.645e+02\n",
      "Cost function:         8.076e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 8.715e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 6.108e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.969e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 1.151e+02\n",
      "\n",
      "=============================\n",
      "Iteration 39 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.994e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.634e+02\n",
      "Cost function:         8.000e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 8.668e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 6.080e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.022e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 1.129e+02\n",
      "\n",
      "=============================\n",
      "Iteration 40 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.886e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.621e+02\n",
      "Cost function:         7.930e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 8.613e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 6.046e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.070e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 1.108e+02\n",
      "\n",
      "=============================\n",
      "Iteration 41 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.790e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.606e+02\n",
      "Cost function:         7.866e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 8.553e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 6.010e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.114e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 1.099e+02\n",
      "\n",
      "=============================\n",
      "Iteration 42 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.712e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.613e+02\n",
      "Cost function:         7.815e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 8.565e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 6.021e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.150e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 1.090e+02\n",
      "\n",
      "=============================\n",
      "Iteration 43 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.640e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.615e+02\n",
      "Cost function:         7.767e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 8.567e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 6.025e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.181e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 1.080e+02\n",
      "\n",
      "=============================\n",
      "Iteration 44 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.574e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.614e+02\n",
      "Cost function:         7.721e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 8.563e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 6.024e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.209e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 1.071e+02\n",
      "\n",
      "=============================\n",
      "Iteration 45 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.512e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.612e+02\n",
      "Cost function:         7.678e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 8.552e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 6.018e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.233e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 1.062e+02\n",
      "\n",
      "=============================\n",
      "Iteration 46 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.454e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.608e+02\n",
      "Cost function:         7.638e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 8.535e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 6.009e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.255e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 1.053e+02\n",
      "\n",
      "=============================\n",
      "Iteration 47 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.400e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.602e+02\n",
      "Cost function:         7.598e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 8.514e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 5.995e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.275e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 1.044e+02\n",
      "\n",
      "=============================\n",
      "Iteration 48 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.349e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.593e+02\n",
      "Cost function:         7.561e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 8.488e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 5.977e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.291e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 1.035e+02\n",
      "\n",
      "=============================\n",
      "Iteration 49 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.301e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.584e+02\n",
      "Cost function:         7.525e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 8.458e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 5.956e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.306e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 1.025e+02\n",
      "\n",
      "=============================\n",
      "Iteration 50 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.256e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.572e+02\n",
      "Cost function:         7.491e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 8.424e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 5.931e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.318e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 1.015e+02\n",
      "\n",
      "=============================\n",
      "Iteration 51 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.213e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.559e+02\n",
      "Cost function:         7.457e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 8.385e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 5.904e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.329e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 1.005e+02\n",
      "\n",
      "=============================\n",
      "Iteration 52 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.173e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.546e+02\n",
      "Cost function:         7.426e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 8.346e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 5.876e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.338e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 9.954e+01\n",
      "\n",
      "=============================\n",
      "Iteration 53 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.134e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.531e+02\n",
      "Cost function:         7.395e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 8.303e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 5.845e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.346e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 9.853e+01\n",
      "\n",
      "=============================\n",
      "Iteration 54 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.098e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.514e+02\n",
      "Cost function:         7.365e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 8.255e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 5.812e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.352e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 9.754e+01\n",
      "\n",
      "=============================\n",
      "Iteration 55 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.063e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.498e+02\n",
      "Cost function:         7.336e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 8.209e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 5.777e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.358e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 9.654e+01\n",
      "\n",
      "=============================\n",
      "Iteration 56 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.030e-02, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.480e+02\n",
      "Cost function:         7.309e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 8.159e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 5.739e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.362e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 9.554e+01\n",
      "\n",
      "=============================\n",
      "Iteration 57 complete\n",
      "\n",
      "\n",
      "scaled || r || = 9.989e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.461e+02\n",
      "Cost function:         7.282e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 8.107e+01\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 5.699e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.364e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 9.453e+01\n",
      "\n",
      "=============================\n",
      "Iteration 58 complete\n",
      "\n",
      "\n",
      "scaled || r || = 9.687e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.442e+02\n",
      "Cost function:         7.255e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 8.055e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 5.659e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.366e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 9.352e+01\n",
      "\n",
      "=============================\n",
      "Iteration 59 complete\n",
      "\n",
      "\n",
      "scaled || r || = 9.402e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.422e+02\n",
      "Cost function:         7.229e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 8.001e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 5.617e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.367e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 9.250e+01\n",
      "\n",
      "=============================\n",
      "Iteration 60 complete\n",
      "\n",
      "\n",
      "scaled || r || = 9.124e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.401e+02\n",
      "Cost function:         7.203e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 7.946e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 5.574e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.366e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 9.149e+01\n",
      "\n",
      "=============================\n",
      "Iteration 61 complete\n",
      "\n",
      "\n",
      "scaled || r || = 8.863e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.380e+02\n",
      "Cost function:         7.178e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 7.889e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 5.530e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.365e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 9.047e+01\n",
      "\n",
      "=============================\n",
      "Iteration 62 complete\n",
      "\n",
      "\n",
      "scaled || r || = 8.609e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.359e+02\n",
      "Cost function:         7.154e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 7.831e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 5.486e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.363e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 8.946e+01\n",
      "\n",
      "=============================\n",
      "Iteration 63 complete\n",
      "\n",
      "\n",
      "scaled || r || = 8.365e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.338e+02\n",
      "Cost function:         7.129e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 7.774e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 5.442e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.361e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 8.844e+01\n",
      "\n",
      "=============================\n",
      "Iteration 64 complete\n",
      "\n",
      "\n",
      "scaled || r || = 8.130e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.316e+02\n",
      "Cost function:         7.105e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 7.715e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 5.397e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.358e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 8.744e+01\n",
      "\n",
      "=============================\n",
      "Iteration 65 complete\n",
      "\n",
      "\n",
      "scaled || r || = 7.905e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.294e+02\n",
      "Cost function:         7.082e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 7.655e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 5.352e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.354e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 8.643e+01\n",
      "\n",
      "=============================\n",
      "Iteration 66 complete\n",
      "\n",
      "\n",
      "scaled || r || = 7.717e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.271e+02\n",
      "Cost function:         7.059e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 7.591e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 5.302e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.348e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 8.544e+01\n",
      "\n",
      "=============================\n",
      "Iteration 67 complete\n",
      "\n",
      "\n",
      "scaled || r || = 7.507e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.249e+02\n",
      "Cost function:         7.035e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 7.532e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 5.257e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.344e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 8.446e+01\n",
      "\n",
      "=============================\n",
      "Iteration 68 complete\n",
      "\n",
      "\n",
      "scaled || r || = 7.310e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.227e+02\n",
      "Cost function:         7.013e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 7.473e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 5.212e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.340e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 8.347e+01\n",
      "\n",
      "=============================\n",
      "Iteration 69 complete\n",
      "\n",
      "\n",
      "scaled || r || = 7.121e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.205e+02\n",
      "Cost function:         6.990e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 7.414e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 5.168e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.335e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 8.250e+01\n",
      "\n",
      "=============================\n",
      "Iteration 70 complete\n",
      "\n",
      "\n",
      "scaled || r || = 6.939e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.183e+02\n",
      "Cost function:         6.968e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 7.354e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 5.123e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.330e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 8.153e+01\n",
      "\n",
      "=============================\n",
      "Iteration 71 complete\n",
      "\n",
      "\n",
      "scaled || r || = 6.766e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.161e+02\n",
      "Cost function:         6.946e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 7.295e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 5.079e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.325e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 8.058e+01\n",
      "\n",
      "=============================\n",
      "Iteration 72 complete\n",
      "\n",
      "\n",
      "scaled || r || = 6.601e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.140e+02\n",
      "Cost function:         6.924e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 7.236e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 5.035e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.320e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 7.963e+01\n",
      "\n",
      "=============================\n",
      "Iteration 73 complete\n",
      "\n",
      "\n",
      "scaled || r || = 6.444e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.119e+02\n",
      "Cost function:         6.903e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 7.178e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 4.991e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.315e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 7.871e+01\n",
      "\n",
      "=============================\n",
      "Iteration 74 complete\n",
      "\n",
      "\n",
      "scaled || r || = 6.283e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.096e+02\n",
      "Cost function:         6.881e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 7.118e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 4.946e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.307e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 7.777e+01\n",
      "\n",
      "=============================\n",
      "Iteration 75 complete\n",
      "\n",
      "\n",
      "scaled || r || = 6.159e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.074e+02\n",
      "Cost function:         6.860e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 7.056e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 4.901e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.301e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 7.692e+01\n",
      "\n",
      "=============================\n",
      "Iteration 76 complete\n",
      "\n",
      "\n",
      "scaled || r || = 6.003e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.054e+02\n",
      "Cost function:         6.839e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 7.001e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 4.859e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.293e+01\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 7.602e+01\n",
      "\n",
      "=============================\n",
      "Iteration 77 complete\n",
      "\n",
      "\n",
      "scaled || r || = 5.863e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.032e+02\n",
      "Cost function:         6.818e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 6.940e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 4.813e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.284e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 7.516e+01\n",
      "\n",
      "=============================\n",
      "Iteration 78 complete\n",
      "\n",
      "\n",
      "scaled || r || = 5.732e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 2.011e+02\n",
      "Cost function:         6.797e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 6.884e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 4.770e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.275e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 7.435e+01\n",
      "\n",
      "=============================\n",
      "Iteration 79 complete\n",
      "\n",
      "\n",
      "scaled || r || = 5.605e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.991e+02\n",
      "Cost function:         6.777e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 6.830e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 4.730e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.266e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 7.355e+01\n",
      "\n",
      "=============================\n",
      "Iteration 80 complete\n",
      "\n",
      "\n",
      "scaled || r || = 5.485e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.971e+02\n",
      "Cost function:         6.757e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 6.777e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 4.690e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.257e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 7.277e+01\n",
      "\n",
      "=============================\n",
      "Iteration 81 complete\n",
      "\n",
      "\n",
      "scaled || r || = 5.369e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.951e+02\n",
      "Cost function:         6.738e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 6.723e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 4.650e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.248e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 7.195e+01\n",
      "\n",
      "=============================\n",
      "Iteration 82 complete\n",
      "\n",
      "\n",
      "scaled || r || = 5.254e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.931e+02\n",
      "Cost function:         6.718e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 6.667e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 4.608e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.239e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 7.115e+01\n",
      "\n",
      "=============================\n",
      "Iteration 83 complete\n",
      "\n",
      "\n",
      "scaled || r || = 5.144e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.911e+02\n",
      "Cost function:         6.699e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 6.611e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 4.566e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.230e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 7.036e+01\n",
      "\n",
      "=============================\n",
      "Iteration 84 complete\n",
      "\n",
      "\n",
      "scaled || r || = 5.037e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.891e+02\n",
      "Cost function:         6.680e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 6.556e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 4.525e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.221e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 6.958e+01\n",
      "\n",
      "=============================\n",
      "Iteration 85 complete\n",
      "\n",
      "\n",
      "scaled || r || = 4.933e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.871e+02\n",
      "Cost function:         6.661e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 6.500e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 4.485e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.212e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 6.882e+01\n",
      "\n",
      "=============================\n",
      "Iteration 86 complete\n",
      "\n",
      "\n",
      "scaled || r || = 4.834e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.852e+02\n",
      "Cost function:         6.642e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 6.445e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 4.444e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.203e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 6.807e+01\n",
      "\n",
      "=============================\n",
      "Iteration 87 complete\n",
      "\n",
      "\n",
      "scaled || r || = 4.739e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.832e+02\n",
      "Cost function:         6.624e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 6.391e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 4.405e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.193e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 6.732e+01\n",
      "\n",
      "=============================\n",
      "Iteration 88 complete\n",
      "\n",
      "\n",
      "scaled || r || = 4.647e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.813e+02\n",
      "Cost function:         6.606e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 6.337e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 4.365e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.184e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 6.659e+01\n",
      "\n",
      "=============================\n",
      "Iteration 89 complete\n",
      "\n",
      "\n",
      "scaled || r || = 4.558e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.794e+02\n",
      "Cost function:         6.587e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 6.283e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 4.325e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.175e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 6.584e+01\n",
      "\n",
      "=============================\n",
      "Iteration 90 complete\n",
      "\n",
      "\n",
      "scaled || r || = 4.468e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.774e+02\n",
      "Cost function:         6.569e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 6.226e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 4.283e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.164e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 6.510e+01\n",
      "\n",
      "=============================\n",
      "Iteration 91 complete\n",
      "\n",
      "\n",
      "scaled || r || = 4.380e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.755e+02\n",
      "Cost function:         6.551e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 6.169e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 4.241e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.153e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 6.438e+01\n",
      "\n",
      "=============================\n",
      "Iteration 92 complete\n",
      "\n",
      "\n",
      "scaled || r || = 4.295e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.735e+02\n",
      "Cost function:         6.533e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 6.113e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 4.199e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.142e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 6.366e+01\n",
      "\n",
      "=============================\n",
      "Iteration 93 complete\n",
      "\n",
      "\n",
      "scaled || r || = 4.213e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.715e+02\n",
      "Cost function:         6.515e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 6.057e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 4.157e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.131e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 6.293e+01\n",
      "\n",
      "=============================\n",
      "Iteration 94 complete\n",
      "\n",
      "\n",
      "scaled || r || = 4.140e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.698e+02\n",
      "Cost function:         6.497e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 6.008e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 4.121e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.125e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 6.221e+01\n",
      "\n",
      "=============================\n",
      "Iteration 95 complete\n",
      "\n",
      "\n",
      "scaled || r || = 4.069e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.684e+02\n",
      "Cost function:         6.479e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 5.970e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 4.092e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.124e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 6.150e+01\n",
      "\n",
      "=============================\n",
      "Iteration 96 complete\n",
      "\n",
      "\n",
      "scaled || r || = 3.999e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.670e+02\n",
      "Cost function:         6.461e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 5.932e+01\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 4.063e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.123e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 6.080e+01\n",
      "\n",
      "=============================\n",
      "Iteration 97 complete\n",
      "\n",
      "\n",
      "scaled || r || = 3.931e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.655e+02\n",
      "Cost function:         6.443e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 5.893e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 4.034e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.122e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 6.014e+01\n",
      "\n",
      "=============================\n",
      "Iteration 98 complete\n",
      "\n",
      "\n",
      "scaled || r || = 3.857e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.640e+02\n",
      "Cost function:         6.426e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 5.849e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 4.000e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.115e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 5.944e+01\n",
      "\n",
      "=============================\n",
      "Iteration 99 complete\n",
      "\n",
      "\n",
      "scaled || r || = 3.796e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.625e+02\n",
      "Cost function:         6.408e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 5.810e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.970e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.112e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 5.876e+01\n",
      "\n",
      "=============================\n",
      "Iteration 100 complete\n",
      "\n",
      "\n",
      "scaled || r || = 3.733e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.611e+02\n",
      "Cost function:         6.390e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 5.771e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.941e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.110e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 5.809e+01\n",
      "\n",
      "=============================\n",
      "Iteration 101 complete\n",
      "\n",
      "\n",
      "scaled || r || = 3.672e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.597e+02\n",
      "Cost function:         6.373e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 5.732e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.913e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.108e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 5.746e+01\n",
      "\n",
      "=============================\n",
      "Iteration 102 complete\n",
      "\n",
      "\n",
      "scaled || r || = 3.608e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.583e+02\n",
      "Cost function:         6.356e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 5.694e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.883e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.105e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 5.680e+01\n",
      "\n",
      "=============================\n",
      "Iteration 103 complete\n",
      "\n",
      "\n",
      "scaled || r || = 3.554e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.569e+02\n",
      "Cost function:         6.339e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 5.655e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.855e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.103e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 5.616e+01\n",
      "\n",
      "=============================\n",
      "Iteration 104 complete\n",
      "\n",
      "\n",
      "scaled || r || = 3.496e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.555e+02\n",
      "Cost function:         6.322e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 5.617e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.826e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.100e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 5.553e+01\n",
      "\n",
      "=============================\n",
      "Iteration 105 complete\n",
      "\n",
      "\n",
      "scaled || r || = 3.441e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.542e+02\n",
      "Cost function:         6.305e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 5.579e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.797e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.097e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 5.493e+01\n",
      "\n",
      "=============================\n",
      "Iteration 106 complete\n",
      "\n",
      "\n",
      "scaled || r || = 3.380e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.526e+02\n",
      "Cost function:         6.288e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 5.535e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.764e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.088e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 5.431e+01\n",
      "\n",
      "=============================\n",
      "Iteration 107 complete\n",
      "\n",
      "\n",
      "scaled || r || = 3.331e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.512e+02\n",
      "Cost function:         6.271e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 5.496e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.735e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.085e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 5.371e+01\n",
      "\n",
      "=============================\n",
      "Iteration 108 complete\n",
      "\n",
      "\n",
      "scaled || r || = 3.280e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.499e+02\n",
      "Cost function:         6.255e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 5.458e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.707e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.081e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 5.312e+01\n",
      "\n",
      "=============================\n",
      "Iteration 109 complete\n",
      "\n",
      "\n",
      "scaled || r || = 3.229e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.485e+02\n",
      "Cost function:         6.238e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 5.420e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.679e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.077e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 5.255e+01\n",
      "\n",
      "=============================\n",
      "Iteration 110 complete\n",
      "\n",
      "\n",
      "scaled || r || = 3.180e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.472e+02\n",
      "Cost function:         6.222e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 5.383e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.650e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.073e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 5.197e+01\n",
      "\n",
      "=============================\n",
      "Iteration 111 complete\n",
      "\n",
      "\n",
      "scaled || r || = 3.133e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.459e+02\n",
      "Cost function:         6.206e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 5.344e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.621e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.069e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 5.140e+01\n",
      "\n",
      "=============================\n",
      "Iteration 112 complete\n",
      "\n",
      "\n",
      "scaled || r || = 3.089e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.445e+02\n",
      "Cost function:         6.190e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 5.307e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.593e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.064e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 5.084e+01\n",
      "\n",
      "=============================\n",
      "Iteration 113 complete\n",
      "\n",
      "\n",
      "scaled || r || = 3.044e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.432e+02\n",
      "Cost function:         6.174e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 5.270e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.565e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.060e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 5.029e+01\n",
      "\n",
      "=============================\n",
      "Iteration 114 complete\n",
      "\n",
      "\n",
      "scaled || r || = 3.000e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.419e+02\n",
      "Cost function:         6.158e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 5.233e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.537e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.055e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 4.974e+01\n",
      "\n",
      "=============================\n",
      "Iteration 115 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.957e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.406e+02\n",
      "Cost function:         6.142e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 5.196e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.509e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.050e+01\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 4.920e+01\n",
      "\n",
      "=============================\n",
      "Iteration 116 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.917e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.393e+02\n",
      "Cost function:         6.126e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 5.159e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.481e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.045e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 4.868e+01\n",
      "\n",
      "=============================\n",
      "Iteration 117 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.876e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.381e+02\n",
      "Cost function:         6.111e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 5.123e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.454e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.040e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 4.816e+01\n",
      "\n",
      "=============================\n",
      "Iteration 118 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.837e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.368e+02\n",
      "Cost function:         6.095e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 5.086e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.426e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.035e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 4.764e+01\n",
      "\n",
      "=============================\n",
      "Iteration 119 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.800e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.355e+02\n",
      "Cost function:         6.080e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 5.050e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.399e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.029e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 4.714e+01\n",
      "\n",
      "=============================\n",
      "Iteration 120 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.763e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.343e+02\n",
      "Cost function:         6.065e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 5.013e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.371e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.024e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 4.664e+01\n",
      "\n",
      "=============================\n",
      "Iteration 121 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.727e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.331e+02\n",
      "Cost function:         6.050e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.977e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.344e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.018e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 4.615e+01\n",
      "\n",
      "=============================\n",
      "Iteration 122 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.693e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.318e+02\n",
      "Cost function:         6.035e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.941e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.317e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.012e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 4.567e+01\n",
      "\n",
      "=============================\n",
      "Iteration 123 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.658e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.306e+02\n",
      "Cost function:         6.020e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.905e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.290e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 3.006e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 4.519e+01\n",
      "\n",
      "=============================\n",
      "Iteration 124 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.626e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.294e+02\n",
      "Cost function:         6.006e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.869e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.263e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.999e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 4.472e+01\n",
      "\n",
      "=============================\n",
      "Iteration 125 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.595e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.282e+02\n",
      "Cost function:         5.991e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.833e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.236e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.993e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 4.426e+01\n",
      "\n",
      "=============================\n",
      "Iteration 126 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.565e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.269e+02\n",
      "Cost function:         5.977e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.798e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.209e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.986e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 4.380e+01\n",
      "\n",
      "=============================\n",
      "Iteration 127 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.534e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.257e+02\n",
      "Cost function:         5.962e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.762e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.183e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.979e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 4.335e+01\n",
      "\n",
      "=============================\n",
      "Iteration 128 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.506e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.246e+02\n",
      "Cost function:         5.948e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.727e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.157e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.972e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 4.291e+01\n",
      "\n",
      "=============================\n",
      "Iteration 129 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.477e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.234e+02\n",
      "Cost function:         5.934e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.692e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.131e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.965e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 4.248e+01\n",
      "\n",
      "=============================\n",
      "Iteration 130 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.449e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.222e+02\n",
      "Cost function:         5.921e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.656e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.105e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.958e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 4.220e+01\n",
      "\n",
      "=============================\n",
      "Iteration 131 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.423e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.215e+02\n",
      "Cost function:         5.907e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.637e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.089e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.951e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 4.194e+01\n",
      "\n",
      "=============================\n",
      "Iteration 132 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.397e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.208e+02\n",
      "Cost function:         5.894e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.620e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.076e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.943e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 4.169e+01\n",
      "\n",
      "=============================\n",
      "Iteration 133 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.372e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.201e+02\n",
      "Cost function:         5.881e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.602e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.062e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.936e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 4.144e+01\n",
      "\n",
      "=============================\n",
      "Iteration 134 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.348e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.194e+02\n",
      "Cost function:         5.868e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.585e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.048e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.928e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 4.118e+01\n",
      "\n",
      "=============================\n",
      "Iteration 135 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.325e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.187e+02\n",
      "Cost function:         5.855e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.567e+01\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.034e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.920e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 4.094e+01\n",
      "\n",
      "=============================\n",
      "Iteration 136 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.303e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.181e+02\n",
      "Cost function:         5.842e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.548e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.020e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.912e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 4.069e+01\n",
      "\n",
      "=============================\n",
      "Iteration 137 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.281e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.174e+02\n",
      "Cost function:         5.830e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.529e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 3.005e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.904e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 4.044e+01\n",
      "\n",
      "=============================\n",
      "Iteration 138 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.260e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.167e+02\n",
      "Cost function:         5.818e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.510e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.991e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.896e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 4.020e+01\n",
      "\n",
      "=============================\n",
      "Iteration 139 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.240e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.159e+02\n",
      "Cost function:         5.805e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.491e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.976e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.888e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.996e+01\n",
      "\n",
      "=============================\n",
      "Iteration 140 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.220e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.152e+02\n",
      "Cost function:         5.793e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.472e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.961e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.879e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.972e+01\n",
      "\n",
      "=============================\n",
      "Iteration 141 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.202e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.145e+02\n",
      "Cost function:         5.781e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.452e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.946e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.871e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.948e+01\n",
      "\n",
      "=============================\n",
      "Iteration 142 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.184e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.138e+02\n",
      "Cost function:         5.770e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.432e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.931e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.862e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.925e+01\n",
      "\n",
      "=============================\n",
      "Iteration 143 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.166e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.131e+02\n",
      "Cost function:         5.758e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.412e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.916e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.853e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.902e+01\n",
      "\n",
      "=============================\n",
      "Iteration 144 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.149e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.124e+02\n",
      "Cost function:         5.747e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.391e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.901e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.844e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.878e+01\n",
      "\n",
      "=============================\n",
      "Iteration 145 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.133e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.117e+02\n",
      "Cost function:         5.735e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.371e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.886e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.835e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.855e+01\n",
      "\n",
      "=============================\n",
      "Iteration 146 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.118e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.110e+02\n",
      "Cost function:         5.724e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.350e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.871e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.826e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.833e+01\n",
      "\n",
      "=============================\n",
      "Iteration 147 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.103e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.103e+02\n",
      "Cost function:         5.713e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.330e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.855e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.817e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.810e+01\n",
      "\n",
      "=============================\n",
      "Iteration 148 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.088e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.095e+02\n",
      "Cost function:         5.702e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.309e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.840e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.808e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.788e+01\n",
      "\n",
      "=============================\n",
      "Iteration 149 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.074e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.088e+02\n",
      "Cost function:         5.692e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.288e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.824e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.798e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.766e+01\n",
      "\n",
      "=============================\n",
      "Iteration 150 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.061e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.081e+02\n",
      "Cost function:         5.681e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.267e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.809e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.789e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.744e+01\n",
      "\n",
      "=============================\n",
      "Iteration 151 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.048e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.074e+02\n",
      "Cost function:         5.670e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.246e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.794e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.780e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.722e+01\n",
      "\n",
      "=============================\n",
      "Iteration 152 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.036e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.067e+02\n",
      "Cost function:         5.660e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.224e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.778e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.770e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.701e+01\n",
      "\n",
      "=============================\n",
      "Iteration 153 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.024e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.060e+02\n",
      "Cost function:         5.650e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.203e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.763e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.761e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.679e+01\n",
      "\n",
      "=============================\n",
      "Iteration 154 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.013e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.053e+02\n",
      "Cost function:         5.640e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.182e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.748e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.751e+01\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.658e+01\n",
      "\n",
      "=============================\n",
      "Iteration 155 complete\n",
      "\n",
      "\n",
      "scaled || r || = 2.002e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.046e+02\n",
      "Cost function:         5.630e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.161e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.732e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.741e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.637e+01\n",
      "\n",
      "=============================\n",
      "Iteration 156 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.991e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.039e+02\n",
      "Cost function:         5.620e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.140e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.717e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.732e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.617e+01\n",
      "\n",
      "=============================\n",
      "Iteration 157 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.981e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.032e+02\n",
      "Cost function:         5.610e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.118e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.702e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.722e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.596e+01\n",
      "\n",
      "=============================\n",
      "Iteration 158 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.972e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.025e+02\n",
      "Cost function:         5.601e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.097e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.687e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.713e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.576e+01\n",
      "\n",
      "=============================\n",
      "Iteration 159 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.962e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.019e+02\n",
      "Cost function:         5.591e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.076e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.672e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.703e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.556e+01\n",
      "\n",
      "=============================\n",
      "Iteration 160 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.953e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.012e+02\n",
      "Cost function:         5.582e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.055e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.657e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.693e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.536e+01\n",
      "\n",
      "=============================\n",
      "Iteration 161 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.945e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 1.005e+02\n",
      "Cost function:         5.572e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.034e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.642e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.684e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.516e+01\n",
      "\n",
      "=============================\n",
      "Iteration 162 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.936e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 9.985e+01\n",
      "Cost function:         5.563e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 4.013e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.628e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.674e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.497e+01\n",
      "\n",
      "=============================\n",
      "Iteration 163 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.928e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 9.918e+01\n",
      "Cost function:         5.554e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.992e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.613e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.664e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.478e+01\n",
      "\n",
      "=============================\n",
      "Iteration 164 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.920e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 9.852e+01\n",
      "Cost function:         5.546e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.971e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.598e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.655e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.459e+01\n",
      "\n",
      "=============================\n",
      "Iteration 165 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.913e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 9.787e+01\n",
      "Cost function:         5.537e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.950e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.584e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.645e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.440e+01\n",
      "\n",
      "=============================\n",
      "Iteration 166 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.911e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 9.720e+01\n",
      "Cost function:         5.528e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.929e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.569e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.635e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.421e+01\n",
      "\n",
      "=============================\n",
      "Iteration 167 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.903e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 9.655e+01\n",
      "Cost function:         5.519e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.908e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.555e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.625e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.403e+01\n",
      "\n",
      "=============================\n",
      "Iteration 168 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.897e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 9.590e+01\n",
      "Cost function:         5.511e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.887e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.540e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.615e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.385e+01\n",
      "\n",
      "=============================\n",
      "Iteration 169 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.891e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 9.524e+01\n",
      "Cost function:         5.502e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.865e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.525e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.605e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.367e+01\n",
      "\n",
      "=============================\n",
      "Iteration 170 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.885e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 9.459e+01\n",
      "Cost function:         5.493e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.844e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.511e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.595e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.349e+01\n",
      "\n",
      "=============================\n",
      "Iteration 171 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.879e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 9.395e+01\n",
      "Cost function:         5.485e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.823e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.497e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.585e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.332e+01\n",
      "\n",
      "=============================\n",
      "Iteration 172 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.873e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 9.331e+01\n",
      "Cost function:         5.477e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.802e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.482e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.575e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.314e+01\n",
      "\n",
      "=============================\n",
      "Iteration 173 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.868e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 9.268e+01\n",
      "Cost function:         5.468e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.781e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.468e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.565e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.297e+01\n",
      "\n",
      "=============================\n",
      "Iteration 174 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.862e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 9.205e+01\n",
      "Cost function:         5.460e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.760e+01\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.454e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.554e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.280e+01\n",
      "\n",
      "=============================\n",
      "Iteration 175 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.858e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 9.141e+01\n",
      "Cost function:         5.452e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.739e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.440e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.544e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.263e+01\n",
      "\n",
      "=============================\n",
      "Iteration 176 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.853e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 9.079e+01\n",
      "Cost function:         5.444e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.718e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.426e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.534e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.246e+01\n",
      "\n",
      "=============================\n",
      "Iteration 177 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.848e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 9.017e+01\n",
      "Cost function:         5.436e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.697e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.412e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.523e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.230e+01\n",
      "\n",
      "=============================\n",
      "Iteration 178 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.843e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 8.955e+01\n",
      "Cost function:         5.428e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.676e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.398e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.513e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.213e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.384e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.503e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.197e+01\n",
      "\n",
      "=============================\n",
      "Iteration 180 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.834e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 8.833e+01\n",
      "Cost function:         5.413e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.635e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.371e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.492e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.181e+01\n",
      "\n",
      "=============================\n",
      "Iteration 181 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.830e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 8.772e+01\n",
      "Cost function:         5.405e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.614e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.357e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.482e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.165e+01\n",
      "\n",
      "=============================\n",
      "Iteration 182 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.827e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 8.711e+01\n",
      "Cost function:         5.397e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.594e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.343e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.472e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.150e+01\n",
      "\n",
      "=============================\n",
      "Iteration 183 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.823e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 8.651e+01\n",
      "Cost function:         5.390e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.573e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.330e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.461e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.134e+01\n",
      "\n",
      "=============================\n",
      "Iteration 184 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.819e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 8.592e+01\n",
      "Cost function:         5.382e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.553e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.316e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.451e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.119e+01\n",
      "\n",
      "=============================\n",
      "Iteration 185 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.813e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 8.532e+01\n",
      "Cost function:         5.375e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.532e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.303e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.440e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.103e+01\n",
      "\n",
      "=============================\n",
      "Iteration 186 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.809e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 8.473e+01\n",
      "Cost function:         5.368e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.512e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.290e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.430e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.088e+01\n",
      "\n",
      "=============================\n",
      "Iteration 187 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.808e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 8.414e+01\n",
      "Cost function:         5.361e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.491e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.276e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.419e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.073e+01\n",
      "\n",
      "=============================\n",
      "Iteration 188 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.804e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 8.356e+01\n",
      "Cost function:         5.353e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.471e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.263e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.408e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.059e+01\n",
      "\n",
      "=============================\n",
      "Iteration 189 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.801e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 8.297e+01\n",
      "Cost function:         5.346e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.451e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.250e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.397e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.044e+01\n",
      "\n",
      "=============================\n",
      "Iteration 190 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.798e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 8.240e+01\n",
      "Cost function:         5.339e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.430e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.237e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.387e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.029e+01\n",
      "\n",
      "=============================\n",
      "Iteration 191 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.795e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 8.182e+01\n",
      "Cost function:         5.332e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.410e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.224e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.376e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.015e+01\n",
      "\n",
      "=============================\n",
      "Iteration 192 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.792e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 8.125e+01\n",
      "Cost function:         5.325e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.390e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.211e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.365e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 3.001e+01\n",
      "\n",
      "=============================\n",
      "Iteration 193 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.788e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 8.068e+01\n",
      "Cost function:         5.319e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.370e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.198e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.355e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.987e+01\n",
      "\n",
      "=============================\n",
      "Iteration 194 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.785e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 8.012e+01\n",
      "Cost function:         5.312e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.350e+01\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.185e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.344e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.973e+01\n",
      "\n",
      "=============================\n",
      "Iteration 195 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.783e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 7.956e+01\n",
      "Cost function:         5.305e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.330e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.173e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.333e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.959e+01\n",
      "\n",
      "=============================\n",
      "Iteration 196 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.779e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 7.900e+01\n",
      "Cost function:         5.298e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.311e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.160e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.322e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.945e+01\n",
      "\n",
      "=============================\n",
      "Iteration 197 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.776e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 7.844e+01\n",
      "Cost function:         5.292e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.291e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.147e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.311e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.931e+01\n",
      "\n",
      "=============================\n",
      "Iteration 198 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.776e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 7.789e+01\n",
      "Cost function:         5.285e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.271e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.135e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.301e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.918e+01\n",
      "\n",
      "=============================\n",
      "Iteration 199 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.770e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 7.735e+01\n",
      "Cost function:         5.279e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.252e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.123e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.290e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.905e+01\n",
      "\n",
      "=============================\n",
      "Iteration 200 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.768e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 7.680e+01\n",
      "Cost function:         5.273e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.232e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.110e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.279e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.891e+01\n",
      "\n",
      "=============================\n",
      "Iteration 201 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.766e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 7.626e+01\n",
      "Cost function:         5.266e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.213e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.098e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.268e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.878e+01\n",
      "\n",
      "=============================\n",
      "Iteration 202 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.762e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 7.573e+01\n",
      "Cost function:         5.260e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.193e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.086e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.257e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.865e+01\n",
      "\n",
      "=============================\n",
      "Iteration 203 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.764e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 7.519e+01\n",
      "Cost function:         5.254e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.174e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.074e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.247e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.853e+01\n",
      "\n",
      "=============================\n",
      "Iteration 204 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.757e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 7.469e+01\n",
      "Cost function:         5.248e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.156e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.063e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.236e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.845e+01\n",
      "\n",
      "=============================\n",
      "Iteration 205 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.758e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 7.428e+01\n",
      "Cost function:         5.242e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.142e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.054e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.225e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.836e+01\n",
      "\n",
      "=============================\n",
      "Iteration 206 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.752e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 7.387e+01\n",
      "Cost function:         5.236e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.128e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.045e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.214e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.828e+01\n",
      "\n",
      "=============================\n",
      "Iteration 207 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.750e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 7.346e+01\n",
      "Cost function:         5.230e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.113e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.037e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.203e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.820e+01\n",
      "\n",
      "=============================\n",
      "Iteration 208 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.749e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 7.305e+01\n",
      "Cost function:         5.224e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.099e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.028e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.192e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.811e+01\n",
      "\n",
      "=============================\n",
      "Iteration 209 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.749e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 7.264e+01\n",
      "Cost function:         5.218e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.085e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.019e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.181e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.803e+01\n",
      "\n",
      "=============================\n",
      "Iteration 210 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.743e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 7.223e+01\n",
      "Cost function:         5.213e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.070e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.011e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.171e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.795e+01\n",
      "\n",
      "=============================\n",
      "Iteration 211 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.741e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 7.182e+01\n",
      "Cost function:         5.207e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.056e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 2.002e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.160e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.786e+01\n",
      "\n",
      "=============================\n",
      "Iteration 212 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.739e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 7.141e+01\n",
      "Cost function:         5.202e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.041e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.993e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.149e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.778e+01\n",
      "\n",
      "=============================\n",
      "Iteration 213 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.738e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 7.101e+01\n",
      "Cost function:         5.196e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.027e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.985e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.138e+01\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.770e+01\n",
      "\n",
      "=============================\n",
      "Iteration 214 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.736e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 7.060e+01\n",
      "Cost function:         5.191e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 3.012e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.976e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.127e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.761e+01\n",
      "\n",
      "=============================\n",
      "Iteration 215 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.734e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 7.020e+01\n",
      "Cost function:         5.185e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.998e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.968e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.116e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.753e+01\n",
      "\n",
      "=============================\n",
      "Iteration 216 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.732e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 6.979e+01\n",
      "Cost function:         5.180e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.984e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.959e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.106e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.745e+01\n",
      "\n",
      "=============================\n",
      "Iteration 217 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.729e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 6.939e+01\n",
      "Cost function:         5.175e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.969e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.951e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.095e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.736e+01\n",
      "\n",
      "=============================\n",
      "Iteration 218 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.726e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 6.899e+01\n",
      "Cost function:         5.169e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.955e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.942e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.084e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.728e+01\n",
      "\n",
      "=============================\n",
      "Iteration 219 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.725e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 6.859e+01\n",
      "Cost function:         5.164e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.940e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.934e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.073e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.720e+01\n",
      "\n",
      "=============================\n",
      "Iteration 220 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.722e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 6.819e+01\n",
      "Cost function:         5.159e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.926e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.925e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.063e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.712e+01\n",
      "\n",
      "=============================\n",
      "Iteration 221 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.721e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 6.779e+01\n",
      "Cost function:         5.154e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.911e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.917e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.052e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.703e+01\n",
      "\n",
      "=============================\n",
      "Iteration 222 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.718e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 6.739e+01\n",
      "Cost function:         5.149e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.897e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.908e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.041e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.695e+01\n",
      "\n",
      "=============================\n",
      "Iteration 223 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.717e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 6.699e+01\n",
      "Cost function:         5.144e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.882e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.900e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.030e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.687e+01\n",
      "\n",
      "=============================\n",
      "Iteration 224 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.714e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 6.659e+01\n",
      "Cost function:         5.139e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.868e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.892e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.019e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.678e+01\n",
      "\n",
      "=============================\n",
      "Iteration 225 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.713e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 6.620e+01\n",
      "Cost function:         5.134e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.853e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.883e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 2.009e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.670e+01\n",
      "\n",
      "=============================\n",
      "Iteration 226 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.712e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 6.580e+01\n",
      "Cost function:         5.129e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.839e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.875e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.998e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.662e+01\n",
      "\n",
      "=============================\n",
      "Iteration 227 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.709e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 6.541e+01\n",
      "Cost function:         5.125e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.824e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.867e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.987e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.654e+01\n",
      "\n",
      "=============================\n",
      "Iteration 228 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.708e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 6.501e+01\n",
      "Cost function:         5.120e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.810e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.858e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.977e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.645e+01\n",
      "\n",
      "=============================\n",
      "Iteration 229 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.707e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 6.462e+01\n",
      "Cost function:         5.115e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.796e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.850e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.966e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.637e+01\n",
      "\n",
      "=============================\n",
      "Iteration 230 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.705e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 6.423e+01\n",
      "Cost function:         5.111e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.781e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.842e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.955e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.629e+01\n",
      "\n",
      "=============================\n",
      "Iteration 231 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.704e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 6.384e+01\n",
      "Cost function:         5.106e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.767e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.834e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.944e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.621e+01\n",
      "\n",
      "=============================\n",
      "Iteration 232 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.702e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 6.345e+01\n",
      "Cost function:         5.101e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.752e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.825e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.934e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.612e+01\n",
      "\n",
      "=============================\n",
      "Iteration 233 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.700e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 6.306e+01\n",
      "Cost function:         5.097e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.738e+01\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.817e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.923e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.604e+01\n",
      "\n",
      "=============================\n",
      "Iteration 234 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.699e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 6.267e+01\n",
      "Cost function:         5.092e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.723e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.809e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.912e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.595e+01\n",
      "\n",
      "=============================\n",
      "Iteration 235 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.704e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 6.228e+01\n",
      "Cost function:         5.088e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.709e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.801e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.902e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.588e+01\n",
      "\n",
      "=============================\n",
      "Iteration 236 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.703e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 6.191e+01\n",
      "Cost function:         5.083e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.694e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.794e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.893e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.581e+01\n",
      "\n",
      "=============================\n",
      "Iteration 237 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.701e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 6.155e+01\n",
      "Cost function:         5.078e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.679e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.787e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.883e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.573e+01\n",
      "\n",
      "=============================\n",
      "Iteration 238 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.699e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 6.118e+01\n",
      "Cost function:         5.073e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.664e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.779e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.873e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.566e+01\n",
      "\n",
      "=============================\n",
      "Iteration 239 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.698e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 6.081e+01\n",
      "Cost function:         5.068e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.649e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.772e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.863e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.559e+01\n",
      "\n",
      "=============================\n",
      "Iteration 240 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.696e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 6.045e+01\n",
      "Cost function:         5.064e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.634e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.764e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.853e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.552e+01\n",
      "\n",
      "=============================\n",
      "Iteration 241 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.694e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 6.008e+01\n",
      "Cost function:         5.059e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.619e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.757e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.843e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.544e+01\n",
      "\n",
      "=============================\n",
      "Iteration 242 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.693e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 5.971e+01\n",
      "Cost function:         5.055e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.604e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.749e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.833e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.537e+01\n",
      "\n",
      "=============================\n",
      "Iteration 243 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.691e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 5.934e+01\n",
      "Cost function:         5.050e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.588e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.741e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.823e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.529e+01\n",
      "\n",
      "=============================\n",
      "Iteration 244 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.689e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 5.897e+01\n",
      "Cost function:         5.046e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.573e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.734e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.812e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.522e+01\n",
      "\n",
      "=============================\n",
      "Iteration 245 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.688e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 5.859e+01\n",
      "Cost function:         5.042e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.557e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.726e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.802e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.514e+01\n",
      "\n",
      "=============================\n",
      "Iteration 246 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.686e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 5.822e+01\n",
      "Cost function:         5.037e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.542e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.719e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.792e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.507e+01\n",
      "\n",
      "=============================\n",
      "Iteration 247 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.685e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 5.784e+01\n",
      "Cost function:         5.033e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.526e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.711e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.781e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.499e+01\n",
      "\n",
      "=============================\n",
      "Iteration 248 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.684e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 5.747e+01\n",
      "Cost function:         5.029e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.510e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.704e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.771e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.491e+01\n",
      "\n",
      "=============================\n",
      "Iteration 249 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.682e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 5.709e+01\n",
      "Cost function:         5.025e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.495e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.696e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.760e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.484e+01\n",
      "\n",
      "=============================\n",
      "Iteration 250 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.681e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 5.671e+01\n",
      "Cost function:         5.021e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.479e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.689e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.750e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.476e+01\n",
      "\n",
      "=============================\n",
      "Iteration 251 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.680e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 5.634e+01\n",
      "Cost function:         5.017e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.463e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.681e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.739e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.468e+01\n",
      "\n",
      "=============================\n",
      "Iteration 252 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.678e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 5.596e+01\n",
      "Cost function:         5.013e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.447e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.674e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.728e+01\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.460e+01\n",
      "\n",
      "=============================\n",
      "Iteration 253 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.677e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 5.558e+01\n",
      "Cost function:         5.009e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.431e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.666e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.718e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.452e+01\n",
      "\n",
      "=============================\n",
      "Iteration 254 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.676e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 5.520e+01\n",
      "Cost function:         5.005e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.415e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.659e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.707e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.444e+01\n",
      "\n",
      "=============================\n",
      "Iteration 255 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.675e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 5.482e+01\n",
      "Cost function:         5.001e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.399e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.651e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.696e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.436e+01\n",
      "\n",
      "=============================\n",
      "Iteration 256 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.674e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 5.444e+01\n",
      "Cost function:         4.997e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.383e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.644e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.686e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.428e+01\n",
      "\n",
      "=============================\n",
      "Iteration 257 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.673e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 5.406e+01\n",
      "Cost function:         4.993e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.367e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.637e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.675e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.420e+01\n",
      "\n",
      "=============================\n",
      "Iteration 258 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.672e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 5.368e+01\n",
      "Cost function:         4.990e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.351e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.629e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.664e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.411e+01\n",
      "\n",
      "=============================\n",
      "Iteration 259 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.671e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 5.331e+01\n",
      "Cost function:         4.986e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.335e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.622e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.653e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.403e+01\n",
      "\n",
      "=============================\n",
      "Iteration 260 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.670e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 5.293e+01\n",
      "Cost function:         4.982e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.319e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.615e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.643e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.395e+01\n",
      "\n",
      "=============================\n",
      "Iteration 261 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.669e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 5.255e+01\n",
      "Cost function:         4.979e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.303e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.607e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.632e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.386e+01\n",
      "\n",
      "=============================\n",
      "Iteration 262 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.668e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 5.217e+01\n",
      "Cost function:         4.975e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.287e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.600e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.621e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.378e+01\n",
      "\n",
      "=============================\n",
      "Iteration 263 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.667e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 5.179e+01\n",
      "Cost function:         4.971e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.271e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.593e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.610e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.369e+01\n",
      "\n",
      "=============================\n",
      "Iteration 264 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.666e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 5.141e+01\n",
      "Cost function:         4.968e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.255e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.585e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.599e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.361e+01\n",
      "\n",
      "=============================\n",
      "Iteration 265 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.665e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 5.104e+01\n",
      "Cost function:         4.964e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.239e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.578e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.588e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.352e+01\n",
      "\n",
      "=============================\n",
      "Iteration 266 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.664e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 5.066e+01\n",
      "Cost function:         4.961e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.223e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.571e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.578e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.343e+01\n",
      "\n",
      "=============================\n",
      "Iteration 267 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.663e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 5.029e+01\n",
      "Cost function:         4.958e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.207e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.564e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.567e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.335e+01\n",
      "\n",
      "=============================\n",
      "Iteration 268 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.662e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 4.991e+01\n",
      "Cost function:         4.954e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.191e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.557e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.556e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.326e+01\n",
      "\n",
      "=============================\n",
      "Iteration 269 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.662e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 4.954e+01\n",
      "Cost function:         4.951e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "Entering inner loop for factor 0\n",
      "qnorm_factr = 2.175e+01\n",
      "Entering inner loop for factor 1\n",
      "qnorm_factr = 1.550e+01\n",
      "Entering inner loop for factor 2\n",
      "qnorm_factr = 1.545e+01\n",
      "Entering inner loop for factor 3\n",
      "qnorm_factr = 2.317e+01\n",
      "\n",
      "=============================\n",
      "Iteration 270 complete\n",
      "\n",
      "\n",
      "scaled || r || = 1.661e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 4.917e+01\n",
      "Cost function:         4.947e+03\n",
      "\n",
      "=============================\n",
      "\n",
      "\n",
      "finished in 271 iterations\n",
      "\n",
      "\n",
      "scaled || r || = 1.661e-03, delta = 5.000e-02\n",
      "Max-qnorm upper bound: 4.917e+01\n",
      "Cost function:         4.947e+03\n",
      "(0.0, 39.414819264217435, 0.001660774324614873, 49.168516501680054467, 0.015133934711390585)\n"
     ]
    }
   ],
   "source": [
    "%autoreload\n",
    "result = run_simulation(n, t, r, 0., r_fit, 4, d=d, verbosity=2, delta=0.05, kappa=100, beta=1, alg='fro', max_iter=271)\n",
    "print(result)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(0.0,\n",
       " 39.414819264217435,\n",
       " 0.001660774324614873,\n",
       " 49.168516501680054467,\n",
       " 0.015133934711390585)"
      ]
     },
     "execution_count": 19,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "result"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "3 * 5 *3 * np.log10(3) * 20"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "20**-2."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "20**1.25"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "5/(5-1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.8.13"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}