EqPerformance-constraints.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "from pycalphad import Model, Database, calculate, equilibrium\n",
    "import pycalphad.variables as v\n",
    "\n",
    "#dbf = Database('2016-08-10-AlGdMgand18RLPSO-for 3d plot.tdb')\n",
    "dbf = Database('alfe_sei.TDB')\n",
    "models = {key: Model(dbf, ['AL', 'FE', 'VA'], key) for key in dbf.phases.keys()}"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/rotis/anaconda/envs/calphadpy3/lib/python3.5/site-packages/Cython/Distutils/old_build_ext.py:30: UserWarning: Cython.Distutils.old_build_ext does not properly handle dependencies and is deprecated.\n",
      "  \"Cython.Distutils.old_build_ext does not properly handle dependencies \"\n"
     ]
    }
   ],
   "source": [
    "#Set compiler directives (cf. http://docs.cython.org/src/reference/compilation.html)\n",
    "%load_ext cython\n",
    "from Cython.Compiler.Options import _directive_defaults\n",
    "\n",
    "_directive_defaults['linetrace'] = True\n",
    "_directive_defaults['binding'] = True"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "metadata": {
    "collapsed": false,
    "scrolled": false
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<!DOCTYPE html>\n",
       "<!-- Generated by Cython 0.25.1 -->\n",
       "<html>\n",
       "<head>\n",
       "    <meta http-equiv=\"Content-Type\" content=\"text/html; charset=utf-8\" />\n",
       "    <title>Cython: _cython_magic_07b8d211ecceae831b3e5546c95de0f8.pyx</title>\n",
       "    <style type=\"text/css\">\n",
       "    \n",
       "body.cython { font-family: courier; font-size: 12; }\n",
       "\n",
       ".cython.tag  {  }\n",
       ".cython.line { margin: 0em }\n",
       ".cython.code { font-size: 9; color: #444444; display: none; margin: 0px 0px 0px 8px; border-left: 8px none; }\n",
       "\n",
       ".cython.line .run { background-color: #B0FFB0; }\n",
       ".cython.line .mis { background-color: #FFB0B0; }\n",
       ".cython.code.run  { border-left: 8px solid #B0FFB0; }\n",
       ".cython.code.mis  { border-left: 8px solid #FFB0B0; }\n",
       "\n",
       ".cython.code .py_c_api  { color: red; }\n",
       ".cython.code .py_macro_api  { color: #FF7000; }\n",
       ".cython.code .pyx_c_api  { color: #FF3000; }\n",
       ".cython.code .pyx_macro_api  { color: #FF7000; }\n",
       ".cython.code .refnanny  { color: #FFA000; }\n",
       ".cython.code .trace  { color: #FFA000; }\n",
       ".cython.code .error_goto  { color: #FFA000; }\n",
       "\n",
       ".cython.code .coerce  { color: #008000; border: 1px dotted #008000 }\n",
       ".cython.code .py_attr { color: #FF0000; font-weight: bold; }\n",
       ".cython.code .c_attr  { color: #0000FF; }\n",
       ".cython.code .py_call { color: #FF0000; font-weight: bold; }\n",
       ".cython.code .c_call  { color: #0000FF; }\n",
       "\n",
       ".cython.score-0 {background-color: #FFFFff;}\n",
       ".cython.score-1 {background-color: #FFFFe7;}\n",
       ".cython.score-2 {background-color: #FFFFd4;}\n",
       ".cython.score-3 {background-color: #FFFFc4;}\n",
       ".cython.score-4 {background-color: #FFFFb6;}\n",
       ".cython.score-5 {background-color: #FFFFaa;}\n",
       ".cython.score-6 {background-color: #FFFF9f;}\n",
       ".cython.score-7 {background-color: #FFFF96;}\n",
       ".cython.score-8 {background-color: #FFFF8d;}\n",
       ".cython.score-9 {background-color: #FFFF86;}\n",
       ".cython.score-10 {background-color: #FFFF7f;}\n",
       ".cython.score-11 {background-color: #FFFF79;}\n",
       ".cython.score-12 {background-color: #FFFF73;}\n",
       ".cython.score-13 {background-color: #FFFF6e;}\n",
       ".cython.score-14 {background-color: #FFFF6a;}\n",
       ".cython.score-15 {background-color: #FFFF66;}\n",
       ".cython.score-16 {background-color: #FFFF62;}\n",
       ".cython.score-17 {background-color: #FFFF5e;}\n",
       ".cython.score-18 {background-color: #FFFF5b;}\n",
       ".cython.score-19 {background-color: #FFFF57;}\n",
       ".cython.score-20 {background-color: #FFFF55;}\n",
       ".cython.score-21 {background-color: #FFFF52;}\n",
       ".cython.score-22 {background-color: #FFFF4f;}\n",
       ".cython.score-23 {background-color: #FFFF4d;}\n",
       ".cython.score-24 {background-color: #FFFF4b;}\n",
       ".cython.score-25 {background-color: #FFFF48;}\n",
       ".cython.score-26 {background-color: #FFFF46;}\n",
       ".cython.score-27 {background-color: #FFFF44;}\n",
       ".cython.score-28 {background-color: #FFFF43;}\n",
       ".cython.score-29 {background-color: #FFFF41;}\n",
       ".cython.score-30 {background-color: #FFFF3f;}\n",
       ".cython.score-31 {background-color: #FFFF3e;}\n",
       ".cython.score-32 {background-color: #FFFF3c;}\n",
       ".cython.score-33 {background-color: #FFFF3b;}\n",
       ".cython.score-34 {background-color: #FFFF39;}\n",
       ".cython.score-35 {background-color: #FFFF38;}\n",
       ".cython.score-36 {background-color: #FFFF37;}\n",
       ".cython.score-37 {background-color: #FFFF36;}\n",
       ".cython.score-38 {background-color: #FFFF35;}\n",
       ".cython.score-39 {background-color: #FFFF34;}\n",
       ".cython.score-40 {background-color: #FFFF33;}\n",
       ".cython.score-41 {background-color: #FFFF32;}\n",
       ".cython.score-42 {background-color: #FFFF31;}\n",
       ".cython.score-43 {background-color: #FFFF30;}\n",
       ".cython.score-44 {background-color: #FFFF2f;}\n",
       ".cython.score-45 {background-color: #FFFF2e;}\n",
       ".cython.score-46 {background-color: #FFFF2d;}\n",
       ".cython.score-47 {background-color: #FFFF2c;}\n",
       ".cython.score-48 {background-color: #FFFF2b;}\n",
       ".cython.score-49 {background-color: #FFFF2b;}\n",
       ".cython.score-50 {background-color: #FFFF2a;}\n",
       ".cython.score-51 {background-color: #FFFF29;}\n",
       ".cython.score-52 {background-color: #FFFF29;}\n",
       ".cython.score-53 {background-color: #FFFF28;}\n",
       ".cython.score-54 {background-color: #FFFF27;}\n",
       ".cython.score-55 {background-color: #FFFF27;}\n",
       ".cython.score-56 {background-color: #FFFF26;}\n",
       ".cython.score-57 {background-color: #FFFF26;}\n",
       ".cython.score-58 {background-color: #FFFF25;}\n",
       ".cython.score-59 {background-color: #FFFF24;}\n",
       ".cython.score-60 {background-color: #FFFF24;}\n",
       ".cython.score-61 {background-color: #FFFF23;}\n",
       ".cython.score-62 {background-color: #FFFF23;}\n",
       ".cython.score-63 {background-color: #FFFF22;}\n",
       ".cython.score-64 {background-color: #FFFF22;}\n",
       ".cython.score-65 {background-color: #FFFF22;}\n",
       ".cython.score-66 {background-color: #FFFF21;}\n",
       ".cython.score-67 {background-color: #FFFF21;}\n",
       ".cython.score-68 {background-color: #FFFF20;}\n",
       ".cython.score-69 {background-color: #FFFF20;}\n",
       ".cython.score-70 {background-color: #FFFF1f;}\n",
       ".cython.score-71 {background-color: #FFFF1f;}\n",
       ".cython.score-72 {background-color: #FFFF1f;}\n",
       ".cython.score-73 {background-color: #FFFF1e;}\n",
       ".cython.score-74 {background-color: #FFFF1e;}\n",
       ".cython.score-75 {background-color: #FFFF1e;}\n",
       ".cython.score-76 {background-color: #FFFF1d;}\n",
       ".cython.score-77 {background-color: #FFFF1d;}\n",
       ".cython.score-78 {background-color: #FFFF1c;}\n",
       ".cython.score-79 {background-color: #FFFF1c;}\n",
       ".cython.score-80 {background-color: #FFFF1c;}\n",
       ".cython.score-81 {background-color: #FFFF1c;}\n",
       ".cython.score-82 {background-color: #FFFF1b;}\n",
       ".cython.score-83 {background-color: #FFFF1b;}\n",
       ".cython.score-84 {background-color: #FFFF1b;}\n",
       ".cython.score-85 {background-color: #FFFF1a;}\n",
       ".cython.score-86 {background-color: #FFFF1a;}\n",
       ".cython.score-87 {background-color: #FFFF1a;}\n",
       ".cython.score-88 {background-color: #FFFF1a;}\n",
       ".cython.score-89 {background-color: #FFFF19;}\n",
       ".cython.score-90 {background-color: #FFFF19;}\n",
       ".cython.score-91 {background-color: #FFFF19;}\n",
       ".cython.score-92 {background-color: #FFFF19;}\n",
       ".cython.score-93 {background-color: #FFFF18;}\n",
       ".cython.score-94 {background-color: #FFFF18;}\n",
       ".cython.score-95 {background-color: #FFFF18;}\n",
       ".cython.score-96 {background-color: #FFFF18;}\n",
       ".cython.score-97 {background-color: #FFFF17;}\n",
       ".cython.score-98 {background-color: #FFFF17;}\n",
       ".cython.score-99 {background-color: #FFFF17;}\n",
       ".cython.score-100 {background-color: #FFFF17;}\n",
       ".cython.score-101 {background-color: #FFFF16;}\n",
       ".cython.score-102 {background-color: #FFFF16;}\n",
       ".cython.score-103 {background-color: #FFFF16;}\n",
       ".cython.score-104 {background-color: #FFFF16;}\n",
       ".cython.score-105 {background-color: #FFFF16;}\n",
       ".cython.score-106 {background-color: #FFFF15;}\n",
       ".cython.score-107 {background-color: #FFFF15;}\n",
       ".cython.score-108 {background-color: #FFFF15;}\n",
       ".cython.score-109 {background-color: #FFFF15;}\n",
       ".cython.score-110 {background-color: #FFFF15;}\n",
       ".cython.score-111 {background-color: #FFFF15;}\n",
       ".cython.score-112 {background-color: #FFFF14;}\n",
       ".cython.score-113 {background-color: #FFFF14;}\n",
       ".cython.score-114 {background-color: #FFFF14;}\n",
       ".cython.score-115 {background-color: #FFFF14;}\n",
       ".cython.score-116 {background-color: #FFFF14;}\n",
       ".cython.score-117 {background-color: #FFFF14;}\n",
       ".cython.score-118 {background-color: #FFFF13;}\n",
       ".cython.score-119 {background-color: #FFFF13;}\n",
       ".cython.score-120 {background-color: #FFFF13;}\n",
       ".cython.score-121 {background-color: #FFFF13;}\n",
       ".cython.score-122 {background-color: #FFFF13;}\n",
       ".cython.score-123 {background-color: #FFFF13;}\n",
       ".cython.score-124 {background-color: #FFFF13;}\n",
       ".cython.score-125 {background-color: #FFFF12;}\n",
       ".cython.score-126 {background-color: #FFFF12;}\n",
       ".cython.score-127 {background-color: #FFFF12;}\n",
       ".cython.score-128 {background-color: #FFFF12;}\n",
       ".cython.score-129 {background-color: #FFFF12;}\n",
       ".cython.score-130 {background-color: #FFFF12;}\n",
       ".cython.score-131 {background-color: #FFFF12;}\n",
       ".cython.score-132 {background-color: #FFFF11;}\n",
       ".cython.score-133 {background-color: #FFFF11;}\n",
       ".cython.score-134 {background-color: #FFFF11;}\n",
       ".cython.score-135 {background-color: #FFFF11;}\n",
       ".cython.score-136 {background-color: #FFFF11;}\n",
       ".cython.score-137 {background-color: #FFFF11;}\n",
       ".cython.score-138 {background-color: #FFFF11;}\n",
       ".cython.score-139 {background-color: #FFFF11;}\n",
       ".cython.score-140 {background-color: #FFFF11;}\n",
       ".cython.score-141 {background-color: #FFFF10;}\n",
       ".cython.score-142 {background-color: #FFFF10;}\n",
       ".cython.score-143 {background-color: #FFFF10;}\n",
       ".cython.score-144 {background-color: #FFFF10;}\n",
       ".cython.score-145 {background-color: #FFFF10;}\n",
       ".cython.score-146 {background-color: #FFFF10;}\n",
       ".cython.score-147 {background-color: #FFFF10;}\n",
       ".cython.score-148 {background-color: #FFFF10;}\n",
       ".cython.score-149 {background-color: #FFFF10;}\n",
       ".cython.score-150 {background-color: #FFFF0f;}\n",
       ".cython.score-151 {background-color: #FFFF0f;}\n",
       ".cython.score-152 {background-color: #FFFF0f;}\n",
       ".cython.score-153 {background-color: #FFFF0f;}\n",
       ".cython.score-154 {background-color: #FFFF0f;}\n",
       ".cython.score-155 {background-color: #FFFF0f;}\n",
       ".cython.score-156 {background-color: #FFFF0f;}\n",
       ".cython.score-157 {background-color: #FFFF0f;}\n",
       ".cython.score-158 {background-color: #FFFF0f;}\n",
       ".cython.score-159 {background-color: #FFFF0f;}\n",
       ".cython.score-160 {background-color: #FFFF0f;}\n",
       ".cython.score-161 {background-color: #FFFF0e;}\n",
       ".cython.score-162 {background-color: #FFFF0e;}\n",
       ".cython.score-163 {background-color: #FFFF0e;}\n",
       ".cython.score-164 {background-color: #FFFF0e;}\n",
       ".cython.score-165 {background-color: #FFFF0e;}\n",
       ".cython.score-166 {background-color: #FFFF0e;}\n",
       ".cython.score-167 {background-color: #FFFF0e;}\n",
       ".cython.score-168 {background-color: #FFFF0e;}\n",
       ".cython.score-169 {background-color: #FFFF0e;}\n",
       ".cython.score-170 {background-color: #FFFF0e;}\n",
       ".cython.score-171 {background-color: #FFFF0e;}\n",
       ".cython.score-172 {background-color: #FFFF0e;}\n",
       ".cython.score-173 {background-color: #FFFF0d;}\n",
       ".cython.score-174 {background-color: #FFFF0d;}\n",
       ".cython.score-175 {background-color: #FFFF0d;}\n",
       ".cython.score-176 {background-color: #FFFF0d;}\n",
       ".cython.score-177 {background-color: #FFFF0d;}\n",
       ".cython.score-178 {background-color: #FFFF0d;}\n",
       ".cython.score-179 {background-color: #FFFF0d;}\n",
       ".cython.score-180 {background-color: #FFFF0d;}\n",
       ".cython.score-181 {background-color: #FFFF0d;}\n",
       ".cython.score-182 {background-color: #FFFF0d;}\n",
       ".cython.score-183 {background-color: #FFFF0d;}\n",
       ".cython.score-184 {background-color: #FFFF0d;}\n",
       ".cython.score-185 {background-color: #FFFF0d;}\n",
       ".cython.score-186 {background-color: #FFFF0d;}\n",
       ".cython.score-187 {background-color: #FFFF0c;}\n",
       ".cython.score-188 {background-color: #FFFF0c;}\n",
       ".cython.score-189 {background-color: #FFFF0c;}\n",
       ".cython.score-190 {background-color: #FFFF0c;}\n",
       ".cython.score-191 {background-color: #FFFF0c;}\n",
       ".cython.score-192 {background-color: #FFFF0c;}\n",
       ".cython.score-193 {background-color: #FFFF0c;}\n",
       ".cython.score-194 {background-color: #FFFF0c;}\n",
       ".cython.score-195 {background-color: #FFFF0c;}\n",
       ".cython.score-196 {background-color: #FFFF0c;}\n",
       ".cython.score-197 {background-color: #FFFF0c;}\n",
       ".cython.score-198 {background-color: #FFFF0c;}\n",
       ".cython.score-199 {background-color: #FFFF0c;}\n",
       ".cython.score-200 {background-color: #FFFF0c;}\n",
       ".cython.score-201 {background-color: #FFFF0c;}\n",
       ".cython.score-202 {background-color: #FFFF0c;}\n",
       ".cython.score-203 {background-color: #FFFF0b;}\n",
       ".cython.score-204 {background-color: #FFFF0b;}\n",
       ".cython.score-205 {background-color: #FFFF0b;}\n",
       ".cython.score-206 {background-color: #FFFF0b;}\n",
       ".cython.score-207 {background-color: #FFFF0b;}\n",
       ".cython.score-208 {background-color: #FFFF0b;}\n",
       ".cython.score-209 {background-color: #FFFF0b;}\n",
       ".cython.score-210 {background-color: #FFFF0b;}\n",
       ".cython.score-211 {background-color: #FFFF0b;}\n",
       ".cython.score-212 {background-color: #FFFF0b;}\n",
       ".cython.score-213 {background-color: #FFFF0b;}\n",
       ".cython.score-214 {background-color: #FFFF0b;}\n",
       ".cython.score-215 {background-color: #FFFF0b;}\n",
       ".cython.score-216 {background-color: #FFFF0b;}\n",
       ".cython.score-217 {background-color: #FFFF0b;}\n",
       ".cython.score-218 {background-color: #FFFF0b;}\n",
       ".cython.score-219 {background-color: #FFFF0b;}\n",
       ".cython.score-220 {background-color: #FFFF0b;}\n",
       ".cython.score-221 {background-color: #FFFF0b;}\n",
       ".cython.score-222 {background-color: #FFFF0a;}\n",
       ".cython.score-223 {background-color: #FFFF0a;}\n",
       ".cython.score-224 {background-color: #FFFF0a;}\n",
       ".cython.score-225 {background-color: #FFFF0a;}\n",
       ".cython.score-226 {background-color: #FFFF0a;}\n",
       ".cython.score-227 {background-color: #FFFF0a;}\n",
       ".cython.score-228 {background-color: #FFFF0a;}\n",
       ".cython.score-229 {background-color: #FFFF0a;}\n",
       ".cython.score-230 {background-color: #FFFF0a;}\n",
       ".cython.score-231 {background-color: #FFFF0a;}\n",
       ".cython.score-232 {background-color: #FFFF0a;}\n",
       ".cython.score-233 {background-color: #FFFF0a;}\n",
       ".cython.score-234 {background-color: #FFFF0a;}\n",
       ".cython.score-235 {background-color: #FFFF0a;}\n",
       ".cython.score-236 {background-color: #FFFF0a;}\n",
       ".cython.score-237 {background-color: #FFFF0a;}\n",
       ".cython.score-238 {background-color: #FFFF0a;}\n",
       ".cython.score-239 {background-color: #FFFF0a;}\n",
       ".cython.score-240 {background-color: #FFFF0a;}\n",
       ".cython.score-241 {background-color: #FFFF0a;}\n",
       ".cython.score-242 {background-color: #FFFF0a;}\n",
       ".cython.score-243 {background-color: #FFFF0a;}\n",
       ".cython.score-244 {background-color: #FFFF0a;}\n",
       ".cython.score-245 {background-color: #FFFF0a;}\n",
       ".cython.score-246 {background-color: #FFFF09;}\n",
       ".cython.score-247 {background-color: #FFFF09;}\n",
       ".cython.score-248 {background-color: #FFFF09;}\n",
       ".cython.score-249 {background-color: #FFFF09;}\n",
       ".cython.score-250 {background-color: #FFFF09;}\n",
       ".cython.score-251 {background-color: #FFFF09;}\n",
       ".cython.score-252 {background-color: #FFFF09;}\n",
       ".cython.score-253 {background-color: #FFFF09;}\n",
       ".cython.score-254 {background-color: #FFFF09;}\n",
       ".cython .hll { background-color: #ffffcc }\n",
       ".cython  { background: #f8f8f8; }\n",
       ".cython .c { color: #408080; font-style: italic } /* Comment */\n",
       ".cython .err { border: 1px solid #FF0000 } /* Error */\n",
       ".cython .k { color: #008000; font-weight: bold } /* Keyword */\n",
       ".cython .o { color: #666666 } /* Operator */\n",
       ".cython .ch { color: #408080; font-style: italic } /* Comment.Hashbang */\n",
       ".cython .cm { color: #408080; font-style: italic } /* Comment.Multiline */\n",
       ".cython .cp { color: #BC7A00 } /* Comment.Preproc */\n",
       ".cython .cpf { color: #408080; font-style: italic } /* Comment.PreprocFile */\n",
       ".cython .c1 { color: #408080; font-style: italic } /* Comment.Single */\n",
       ".cython .cs { color: #408080; font-style: italic } /* Comment.Special */\n",
       ".cython .gd { color: #A00000 } /* Generic.Deleted */\n",
       ".cython .ge { font-style: italic } /* Generic.Emph */\n",
       ".cython .gr { color: #FF0000 } /* Generic.Error */\n",
       ".cython .gh { color: #000080; font-weight: bold } /* Generic.Heading */\n",
       ".cython .gi { color: #00A000 } /* Generic.Inserted */\n",
       ".cython .go { color: #888888 } /* Generic.Output */\n",
       ".cython .gp { color: #000080; font-weight: bold } /* Generic.Prompt */\n",
       ".cython .gs { font-weight: bold } /* Generic.Strong */\n",
       ".cython .gu { color: #800080; font-weight: bold } /* Generic.Subheading */\n",
       ".cython .gt { color: #0044DD } /* Generic.Traceback */\n",
       ".cython .kc { color: #008000; font-weight: bold } /* Keyword.Constant */\n",
       ".cython .kd { color: #008000; font-weight: bold } /* Keyword.Declaration */\n",
       ".cython .kn { color: #008000; font-weight: bold } /* Keyword.Namespace */\n",
       ".cython .kp { color: #008000 } /* Keyword.Pseudo */\n",
       ".cython .kr { color: #008000; font-weight: bold } /* Keyword.Reserved */\n",
       ".cython .kt { color: #B00040 } /* Keyword.Type */\n",
       ".cython .m { color: #666666 } /* Literal.Number */\n",
       ".cython .s { color: #BA2121 } /* Literal.String */\n",
       ".cython .na { color: #7D9029 } /* Name.Attribute */\n",
       ".cython .nb { color: #008000 } /* Name.Builtin */\n",
       ".cython .nc { color: #0000FF; font-weight: bold } /* Name.Class */\n",
       ".cython .no { color: #880000 } /* Name.Constant */\n",
       ".cython .nd { color: #AA22FF } /* Name.Decorator */\n",
       ".cython .ni { color: #999999; font-weight: bold } /* Name.Entity */\n",
       ".cython .ne { color: #D2413A; font-weight: bold } /* Name.Exception */\n",
       ".cython .nf { color: #0000FF } /* Name.Function */\n",
       ".cython .nl { color: #A0A000 } /* Name.Label */\n",
       ".cython .nn { color: #0000FF; font-weight: bold } /* Name.Namespace */\n",
       ".cython .nt { color: #008000; font-weight: bold } /* Name.Tag */\n",
       ".cython .nv { color: #19177C } /* Name.Variable */\n",
       ".cython .ow { color: #AA22FF; font-weight: bold } /* Operator.Word */\n",
       ".cython .w { color: #bbbbbb } /* Text.Whitespace */\n",
       ".cython .mb { color: #666666 } /* Literal.Number.Bin */\n",
       ".cython .mf { color: #666666 } /* Literal.Number.Float */\n",
       ".cython .mh { color: #666666 } /* Literal.Number.Hex */\n",
       ".cython .mi { color: #666666 } /* Literal.Number.Integer */\n",
       ".cython .mo { color: #666666 } /* Literal.Number.Oct */\n",
       ".cython .sb { color: #BA2121 } /* Literal.String.Backtick */\n",
       ".cython .sc { color: #BA2121 } /* Literal.String.Char */\n",
       ".cython .sd { color: #BA2121; font-style: italic } /* Literal.String.Doc */\n",
       ".cython .s2 { color: #BA2121 } /* Literal.String.Double */\n",
       ".cython .se { color: #BB6622; font-weight: bold } /* Literal.String.Escape */\n",
       ".cython .sh { color: #BA2121 } /* Literal.String.Heredoc */\n",
       ".cython .si { color: #BB6688; font-weight: bold } /* Literal.String.Interpol */\n",
       ".cython .sx { color: #008000 } /* Literal.String.Other */\n",
       ".cython .sr { color: #BB6688 } /* Literal.String.Regex */\n",
       ".cython .s1 { color: #BA2121 } /* Literal.String.Single */\n",
       ".cython .ss { color: #19177C } /* Literal.String.Symbol */\n",
       ".cython .bp { color: #008000 } /* Name.Builtin.Pseudo */\n",
       ".cython .vc { color: #19177C } /* Name.Variable.Class */\n",
       ".cython .vg { color: #19177C } /* Name.Variable.Global */\n",
       ".cython .vi { color: #19177C } /* Name.Variable.Instance */\n",
       ".cython .il { color: #666666 } /* Literal.Number.Integer.Long */\n",
       "    </style>\n",
       "    <script>\n",
       "    function toggleDiv(id) {\n",
       "        theDiv = id.nextElementSibling\n",
       "        if (theDiv.style.display != 'block') theDiv.style.display = 'block';\n",
       "        else theDiv.style.display = 'none';\n",
       "    }\n",
       "    </script>\n",
       "</head>\n",
       "<body class=\"cython\">\n",
       "<p><span style=\"border-bottom: solid 1px grey;\">Generated by Cython 0.25.1</span></p>\n",
       "<p>\n",
       "    <span style=\"background-color: #FFFF00\">Yellow lines</span> hint at Python interaction.<br />\n",
       "    Click on a line that starts with a \"<code>+</code>\" to see the C code that Cython generated for it.\n",
       "</p>\n",
       "<div class=\"cython\"><pre class=\"cython line score-11\" onclick='toggleDiv(this)'>+<span class=\"\">01</span>: <span class=\"k\">cimport</span> <span class=\"nn\">numpy</span> <span class=\"k\">as</span> <span class=\"nn\">np</span></pre>\n",
       "<pre class='cython code score-11 '>  <span class='trace'>__Pyx_TraceLine</span>(1,0,__PYX_ERR(0, 1, __pyx_L1_error))\n",
       "  __pyx_t_2 = <span class='py_c_api'>PyDict_New</span>(); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 1, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "  if (<span class='py_c_api'>PyDict_SetItem</span>(__pyx_d, __pyx_n_s_test, __pyx_t_2) &lt; 0) __PYX_ERR(0, 1, __pyx_L1_error)\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "</pre><pre class=\"cython line score-8\" onclick='toggleDiv(this)'>+<span class=\"\">02</span>: <span class=\"k\">import</span> <span class=\"nn\">numpy</span> <span class=\"k\">as</span> <span class=\"nn\">np</span></pre>\n",
       "<pre class='cython code score-8 '>  <span class='trace'>__Pyx_TraceLine</span>(2,0,__PYX_ERR(0, 2, __pyx_L1_error))\n",
       "  __pyx_t_2 = <span class='pyx_c_api'>__Pyx_Import</span>(__pyx_n_s_numpy, 0, 0); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 2, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "  if (<span class='py_c_api'>PyDict_SetItem</span>(__pyx_d, __pyx_n_s_np, __pyx_t_2) &lt; 0) __PYX_ERR(0, 2, __pyx_L1_error)\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "</pre><pre class=\"cython line score-29\" onclick='toggleDiv(this)'>+<span class=\"\">03</span>: <span class=\"k\">from</span> <span class=\"nn\">collections</span> <span class=\"k\">import</span> <span class=\"n\">defaultdict</span><span class=\"p\">,</span> <span class=\"n\">OrderedDict</span></pre>\n",
       "<pre class='cython code score-29 '>  <span class='trace'>__Pyx_TraceLine</span>(3,0,__PYX_ERR(0, 3, __pyx_L1_error))\n",
       "  __pyx_t_2 = <span class='py_c_api'>PyList_New</span>(2); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 3, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "  <span class='pyx_macro_api'>__Pyx_INCREF</span>(__pyx_n_s_defaultdict);\n",
       "  <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_n_s_defaultdict);\n",
       "  <span class='py_macro_api'>PyList_SET_ITEM</span>(__pyx_t_2, 0, __pyx_n_s_defaultdict);\n",
       "  <span class='pyx_macro_api'>__Pyx_INCREF</span>(__pyx_n_s_OrderedDict);\n",
       "  <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_n_s_OrderedDict);\n",
       "  <span class='py_macro_api'>PyList_SET_ITEM</span>(__pyx_t_2, 1, __pyx_n_s_OrderedDict);\n",
       "  __pyx_t_3 = <span class='pyx_c_api'>__Pyx_Import</span>(__pyx_n_s_collections, __pyx_t_2, 0); if (unlikely(!__pyx_t_3)) __PYX_ERR(0, 3, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_3);\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "  __pyx_t_2 = <span class='pyx_c_api'>__Pyx_ImportFrom</span>(__pyx_t_3, __pyx_n_s_defaultdict); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 3, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "  if (<span class='py_c_api'>PyDict_SetItem</span>(__pyx_d, __pyx_n_s_defaultdict, __pyx_t_2) &lt; 0) __PYX_ERR(0, 3, __pyx_L1_error)\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "  __pyx_t_2 = <span class='pyx_c_api'>__Pyx_ImportFrom</span>(__pyx_t_3, __pyx_n_s_OrderedDict); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 3, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "  if (<span class='py_c_api'>PyDict_SetItem</span>(__pyx_d, __pyx_n_s_OrderedDict, __pyx_t_2) &lt; 0) __PYX_ERR(0, 3, __pyx_L1_error)\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_3); __pyx_t_3 = 0;\n",
       "</pre><pre class=\"cython line score-8\" onclick='toggleDiv(this)'>+<span class=\"\">04</span>: <span class=\"k\">import</span> <span class=\"nn\">copy</span></pre>\n",
       "<pre class='cython code score-8 '>  <span class='trace'>__Pyx_TraceLine</span>(4,0,__PYX_ERR(0, 4, __pyx_L1_error))\n",
       "  __pyx_t_3 = <span class='pyx_c_api'>__Pyx_Import</span>(__pyx_n_s_copy, 0, 0); if (unlikely(!__pyx_t_3)) __PYX_ERR(0, 4, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_3);\n",
       "  if (<span class='py_c_api'>PyDict_SetItem</span>(__pyx_d, __pyx_n_s_copy, __pyx_t_3) &lt; 0) __PYX_ERR(0, 4, __pyx_L1_error)\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_3); __pyx_t_3 = 0;\n",
       "</pre><pre class=\"cython line score-8\" onclick='toggleDiv(this)'>+<span class=\"\">05</span>: <span class=\"k\">import</span> <span class=\"nn\">itertools</span></pre>\n",
       "<pre class='cython code score-8 '>  <span class='trace'>__Pyx_TraceLine</span>(5,0,__PYX_ERR(0, 5, __pyx_L1_error))\n",
       "  __pyx_t_3 = <span class='pyx_c_api'>__Pyx_Import</span>(__pyx_n_s_itertools, 0, 0); if (unlikely(!__pyx_t_3)) __PYX_ERR(0, 5, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_3);\n",
       "  if (<span class='py_c_api'>PyDict_SetItem</span>(__pyx_d, __pyx_n_s_itertools, __pyx_t_3) &lt; 0) __PYX_ERR(0, 5, __pyx_L1_error)\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_3); __pyx_t_3 = 0;\n",
       "</pre><pre class=\"cython line score-8\" onclick='toggleDiv(this)'>+<span class=\"\">06</span>: <span class=\"k\">import</span> <span class=\"nn\">scipy.spatial</span></pre>\n",
       "<pre class='cython code score-8 '>  <span class='trace'>__Pyx_TraceLine</span>(6,0,__PYX_ERR(0, 6, __pyx_L1_error))\n",
       "  __pyx_t_3 = <span class='pyx_c_api'>__Pyx_Import</span>(__pyx_n_s_scipy_spatial, 0, 0); if (unlikely(!__pyx_t_3)) __PYX_ERR(0, 6, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_3);\n",
       "  if (<span class='py_c_api'>PyDict_SetItem</span>(__pyx_d, __pyx_n_s_scipy, __pyx_t_3) &lt; 0) __PYX_ERR(0, 6, __pyx_L1_error)\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_3); __pyx_t_3 = 0;\n",
       "</pre><pre class=\"cython line score-19\" onclick='toggleDiv(this)'>+<span class=\"\">07</span>: <span class=\"k\">from</span> <span class=\"nn\">pycalphad.core.sympydiff_utils</span> <span class=\"k\">import</span> <span class=\"n\">build_functions</span> <span class=\"k\">as</span> <span class=\"n\">compiled_build_functions</span></pre>\n",
       "<pre class='cython code score-19 '>  <span class='trace'>__Pyx_TraceLine</span>(7,0,__PYX_ERR(0, 7, __pyx_L1_error))\n",
       "  __pyx_t_3 = <span class='py_c_api'>PyList_New</span>(1); if (unlikely(!__pyx_t_3)) __PYX_ERR(0, 7, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_3);\n",
       "  <span class='pyx_macro_api'>__Pyx_INCREF</span>(__pyx_n_s_build_functions);\n",
       "  <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_n_s_build_functions);\n",
       "  <span class='py_macro_api'>PyList_SET_ITEM</span>(__pyx_t_3, 0, __pyx_n_s_build_functions);\n",
       "  __pyx_t_2 = <span class='pyx_c_api'>__Pyx_Import</span>(__pyx_n_s_pycalphad_core_sympydiff_utils, __pyx_t_3, 0); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 7, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_3); __pyx_t_3 = 0;\n",
       "  __pyx_t_3 = <span class='pyx_c_api'>__Pyx_ImportFrom</span>(__pyx_t_2, __pyx_n_s_build_functions); if (unlikely(!__pyx_t_3)) __PYX_ERR(0, 7, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_3);\n",
       "  if (<span class='py_c_api'>PyDict_SetItem</span>(__pyx_d, __pyx_n_s_compiled_build_functions, __pyx_t_3) &lt; 0) __PYX_ERR(0, 7, __pyx_L1_error)\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_3); __pyx_t_3 = 0;\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "</pre><pre class=\"cython line score-29\" onclick='toggleDiv(this)'>+<span class=\"\">08</span>: <span class=\"k\">from</span> <span class=\"nn\">pycalphad.core.constants</span> <span class=\"k\">import</span> <span class=\"n\">MIN_SITE_FRACTION</span><span class=\"p\">,</span> <span class=\"n\">COMP_DIFFERENCE_TOL</span></pre>\n",
       "<pre class='cython code score-29 '>  <span class='trace'>__Pyx_TraceLine</span>(8,0,__PYX_ERR(0, 8, __pyx_L1_error))\n",
       "  __pyx_t_2 = <span class='py_c_api'>PyList_New</span>(2); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 8, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "  <span class='pyx_macro_api'>__Pyx_INCREF</span>(__pyx_n_s_MIN_SITE_FRACTION);\n",
       "  <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_n_s_MIN_SITE_FRACTION);\n",
       "  <span class='py_macro_api'>PyList_SET_ITEM</span>(__pyx_t_2, 0, __pyx_n_s_MIN_SITE_FRACTION);\n",
       "  <span class='pyx_macro_api'>__Pyx_INCREF</span>(__pyx_n_s_COMP_DIFFERENCE_TOL);\n",
       "  <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_n_s_COMP_DIFFERENCE_TOL);\n",
       "  <span class='py_macro_api'>PyList_SET_ITEM</span>(__pyx_t_2, 1, __pyx_n_s_COMP_DIFFERENCE_TOL);\n",
       "  __pyx_t_3 = <span class='pyx_c_api'>__Pyx_Import</span>(__pyx_n_s_pycalphad_core_constants, __pyx_t_2, 0); if (unlikely(!__pyx_t_3)) __PYX_ERR(0, 8, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_3);\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "  __pyx_t_2 = <span class='pyx_c_api'>__Pyx_ImportFrom</span>(__pyx_t_3, __pyx_n_s_MIN_SITE_FRACTION); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 8, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "  if (<span class='py_c_api'>PyDict_SetItem</span>(__pyx_d, __pyx_n_s_MIN_SITE_FRACTION, __pyx_t_2) &lt; 0) __PYX_ERR(0, 8, __pyx_L1_error)\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "  __pyx_t_2 = <span class='pyx_c_api'>__Pyx_ImportFrom</span>(__pyx_t_3, __pyx_n_s_COMP_DIFFERENCE_TOL); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 8, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "  if (<span class='py_c_api'>PyDict_SetItem</span>(__pyx_d, __pyx_n_s_COMP_DIFFERENCE_TOL, __pyx_t_2) &lt; 0) __PYX_ERR(0, 8, __pyx_L1_error)\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_3); __pyx_t_3 = 0;\n",
       "</pre><pre class=\"cython line score-0\">&#xA0;<span class=\"\">09</span>: <span class=\"k\">from</span> <span class=\"nn\">pycalphad.core.phase_rec</span> <span class=\"k\">cimport</span> <span class=\"n\">PhaseRecord</span><span class=\"p\">,</span> <span class=\"n\">obj</span><span class=\"p\">,</span> <span class=\"n\">grad</span><span class=\"p\">,</span> <span class=\"n\">hess</span><span class=\"p\">,</span> <span class=\"n\">mass_obj</span><span class=\"p\">,</span> <span class=\"n\">mass_grad</span><span class=\"p\">,</span> <span class=\"n\">mass_hess</span></pre>\n",
       "<pre class=\"cython line score-19\" onclick='toggleDiv(this)'>+<span class=\"\">10</span>: <span class=\"k\">from</span> <span class=\"nn\">pycalphad.core.cache</span> <span class=\"k\">import</span> <span class=\"n\">cacheit</span></pre>\n",
       "<pre class='cython code score-19 '>  <span class='trace'>__Pyx_TraceLine</span>(10,0,__PYX_ERR(0, 10, __pyx_L1_error))\n",
       "  __pyx_t_3 = <span class='py_c_api'>PyList_New</span>(1); if (unlikely(!__pyx_t_3)) __PYX_ERR(0, 10, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_3);\n",
       "  <span class='pyx_macro_api'>__Pyx_INCREF</span>(__pyx_n_s_cacheit);\n",
       "  <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_n_s_cacheit);\n",
       "  <span class='py_macro_api'>PyList_SET_ITEM</span>(__pyx_t_3, 0, __pyx_n_s_cacheit);\n",
       "  __pyx_t_2 = <span class='pyx_c_api'>__Pyx_Import</span>(__pyx_n_s_pycalphad_core_cache, __pyx_t_3, 0); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 10, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_3); __pyx_t_3 = 0;\n",
       "  __pyx_t_3 = <span class='pyx_c_api'>__Pyx_ImportFrom</span>(__pyx_t_2, __pyx_n_s_cacheit); if (unlikely(!__pyx_t_3)) __PYX_ERR(0, 10, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_3);\n",
       "  if (<span class='py_c_api'>PyDict_SetItem</span>(__pyx_d, __pyx_n_s_cacheit, __pyx_t_3) &lt; 0) __PYX_ERR(0, 10, __pyx_L1_error)\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_3); __pyx_t_3 = 0;\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "</pre><pre class=\"cython line score-16\" onclick='toggleDiv(this)'>+<span class=\"\">11</span>: <span class=\"k\">import</span> <span class=\"nn\">pycalphad.variables</span> <span class=\"k\">as</span> <span class=\"nn\">v</span></pre>\n",
       "<pre class='cython code score-16 '>  <span class='trace'>__Pyx_TraceLine</span>(11,0,__PYX_ERR(0, 11, __pyx_L1_error))\n",
       "  __pyx_t_2 = <span class='py_c_api'>PyList_New</span>(1); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 11, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "  <span class='pyx_macro_api'>__Pyx_INCREF</span>(__pyx_n_s__33);\n",
       "  <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_n_s__33);\n",
       "  <span class='py_macro_api'>PyList_SET_ITEM</span>(__pyx_t_2, 0, __pyx_n_s__33);\n",
       "  __pyx_t_3 = <span class='pyx_c_api'>__Pyx_Import</span>(__pyx_n_s_pycalphad_variables, __pyx_t_2, 0); if (unlikely(!__pyx_t_3)) __PYX_ERR(0, 11, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_3);\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "  if (<span class='py_c_api'>PyDict_SetItem</span>(__pyx_d, __pyx_n_s_v, __pyx_t_3) &lt; 0) __PYX_ERR(0, 11, __pyx_L1_error)\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_3); __pyx_t_3 = 0;\n",
       "</pre><pre class=\"cython line score-19\" onclick='toggleDiv(this)'>+<span class=\"\">12</span>: <span class=\"k\">from</span> <span class=\"nn\">pycalphad.core.eqsolver</span> <span class=\"k\">import</span> <span class=\"n\">_compute_phase_dof</span></pre>\n",
       "<pre class='cython code score-19 '>  <span class='trace'>__Pyx_TraceLine</span>(12,0,__PYX_ERR(0, 12, __pyx_L1_error))\n",
       "  __pyx_t_3 = <span class='py_c_api'>PyList_New</span>(1); if (unlikely(!__pyx_t_3)) __PYX_ERR(0, 12, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_3);\n",
       "  <span class='pyx_macro_api'>__Pyx_INCREF</span>(__pyx_n_s_compute_phase_dof);\n",
       "  <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_n_s_compute_phase_dof);\n",
       "  <span class='py_macro_api'>PyList_SET_ITEM</span>(__pyx_t_3, 0, __pyx_n_s_compute_phase_dof);\n",
       "  __pyx_t_2 = <span class='pyx_c_api'>__Pyx_Import</span>(__pyx_n_s_pycalphad_core_eqsolver, __pyx_t_3, 0); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 12, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_3); __pyx_t_3 = 0;\n",
       "  __pyx_t_3 = <span class='pyx_c_api'>__Pyx_ImportFrom</span>(__pyx_t_2, __pyx_n_s_compute_phase_dof); if (unlikely(!__pyx_t_3)) __PYX_ERR(0, 12, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_3);\n",
       "  if (<span class='py_c_api'>PyDict_SetItem</span>(__pyx_d, __pyx_n_s_compute_phase_dof, __pyx_t_3) &lt; 0) __PYX_ERR(0, 12, __pyx_L1_error)\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_3); __pyx_t_3 = 0;\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "</pre><pre class=\"cython line score-0\">&#xA0;<span class=\"\">13</span>: <span class=\"k\">cimport</span> <span class=\"nn\">cython</span></pre>\n",
       "<pre class=\"cython line score-0\">&#xA0;<span class=\"\">14</span>: </pre>\n",
       "<pre class=\"cython line score-0\">&#xA0;<span class=\"\">15</span>: <span class=\"nd\">@cython</span><span class=\"o\">.</span><span class=\"n\">boundscheck</span><span class=\"p\">(</span><span class=\"bp\">False</span><span class=\"p\">)</span></pre>\n",
       "<pre class=\"cython line score-0\">&#xA0;<span class=\"\">16</span>: <span class=\"nd\">@cython</span><span class=\"o\">.</span><span class=\"n\">wraparound</span><span class=\"p\">(</span><span class=\"bp\">False</span><span class=\"p\">)</span></pre>\n",
       "<pre class=\"cython line score-144\" onclick='toggleDiv(this)'>+<span class=\"\">17</span>: <span class=\"k\">def</span> <span class=\"nf\">_compute_constraints</span><span class=\"p\">(</span><span class=\"nb\">object</span> <span class=\"n\">dbf</span><span class=\"p\">,</span> <span class=\"nb\">object</span> <span class=\"n\">comps</span><span class=\"p\">,</span> <span class=\"nb\">object</span> <span class=\"n\">phases</span><span class=\"p\">,</span></pre>\n",
       "<pre class='cython code score-144 '>/* Python wrapper */\n",
       "static PyObject *__pyx_pw_46_cython_magic_07b8d211ecceae831b3e5546c95de0f8_1_compute_constraints(PyObject *__pyx_self, PyObject *__pyx_args, PyObject *__pyx_kwds); /*proto*/\n",
       "static char __pyx_doc_46_cython_magic_07b8d211ecceae831b3e5546c95de0f8__compute_constraints[] = \"\\n    Compute the constraint vector and constraint Jacobian matrix.\\n    \";\n",
       "static PyMethodDef __pyx_mdef_46_cython_magic_07b8d211ecceae831b3e5546c95de0f8_1_compute_constraints = {\"_compute_constraints\", (PyCFunction)__pyx_pw_46_cython_magic_07b8d211ecceae831b3e5546c95de0f8_1_compute_constraints, METH_VARARGS|METH_KEYWORDS, __pyx_doc_46_cython_magic_07b8d211ecceae831b3e5546c95de0f8__compute_constraints};\n",
       "static PyObject *__pyx_pw_46_cython_magic_07b8d211ecceae831b3e5546c95de0f8_1_compute_constraints(PyObject *__pyx_self, PyObject *__pyx_args, PyObject *__pyx_kwds) {\n",
       "  PyObject *__pyx_v_dbf = 0;\n",
       "  PyObject *__pyx_v_comps = 0;\n",
       "  PyObject *__pyx_v_phases = 0;\n",
       "  PyObject *__pyx_v_cur_conds = 0;\n",
       "  PyArrayObject *__pyx_v_site_fracs = 0;\n",
       "  PyArrayObject *__pyx_v_phase_fracs = 0;\n",
       "  PyObject *__pyx_v_phase_records = 0;\n",
       "  PyObject *__pyx_r = 0;\n",
       "  <span class='refnanny'>__Pyx_RefNannyDeclarations</span>\n",
       "  <span class='refnanny'>__Pyx_RefNannySetupContext</span>(\"_compute_constraints (wrapper)\", 0);\n",
       "  {\n",
       "    static PyObject **__pyx_pyargnames[] = {&amp;__pyx_n_s_dbf,&amp;__pyx_n_s_comps,&amp;__pyx_n_s_phases,&amp;__pyx_n_s_cur_conds,&amp;__pyx_n_s_site_fracs,&amp;__pyx_n_s_phase_fracs,&amp;__pyx_n_s_phase_records,0};\n",
       "    PyObject* values[7] = {0,0,0,0,0,0,0};\n",
       "    if (unlikely(__pyx_kwds)) {\n",
       "      Py_ssize_t kw_args;\n",
       "      const Py_ssize_t pos_args = <span class='py_macro_api'>PyTuple_GET_SIZE</span>(__pyx_args);\n",
       "      switch (pos_args) {\n",
       "        case  7: values[6] = <span class='py_macro_api'>PyTuple_GET_ITEM</span>(__pyx_args, 6);\n",
       "        case  6: values[5] = <span class='py_macro_api'>PyTuple_GET_ITEM</span>(__pyx_args, 5);\n",
       "        case  5: values[4] = <span class='py_macro_api'>PyTuple_GET_ITEM</span>(__pyx_args, 4);\n",
       "        case  4: values[3] = <span class='py_macro_api'>PyTuple_GET_ITEM</span>(__pyx_args, 3);\n",
       "        case  3: values[2] = <span class='py_macro_api'>PyTuple_GET_ITEM</span>(__pyx_args, 2);\n",
       "        case  2: values[1] = <span class='py_macro_api'>PyTuple_GET_ITEM</span>(__pyx_args, 1);\n",
       "        case  1: values[0] = <span class='py_macro_api'>PyTuple_GET_ITEM</span>(__pyx_args, 0);\n",
       "        case  0: break;\n",
       "        default: goto __pyx_L5_argtuple_error;\n",
       "      }\n",
       "      kw_args = <span class='py_c_api'>PyDict_Size</span>(__pyx_kwds);\n",
       "      switch (pos_args) {\n",
       "        case  0:\n",
       "        if (likely((values[0] = <span class='py_c_api'>PyDict_GetItem</span>(__pyx_kwds, __pyx_n_s_dbf)) != 0)) kw_args--;\n",
       "        else goto __pyx_L5_argtuple_error;\n",
       "        case  1:\n",
       "        if (likely((values[1] = <span class='py_c_api'>PyDict_GetItem</span>(__pyx_kwds, __pyx_n_s_comps)) != 0)) kw_args--;\n",
       "        else {\n",
       "          <span class='pyx_c_api'>__Pyx_RaiseArgtupleInvalid</span>(\"_compute_constraints\", 1, 7, 7, 1); __PYX_ERR(0, 17, __pyx_L3_error)\n",
       "        }\n",
       "        case  2:\n",
       "        if (likely((values[2] = <span class='py_c_api'>PyDict_GetItem</span>(__pyx_kwds, __pyx_n_s_phases)) != 0)) kw_args--;\n",
       "        else {\n",
       "          <span class='pyx_c_api'>__Pyx_RaiseArgtupleInvalid</span>(\"_compute_constraints\", 1, 7, 7, 2); __PYX_ERR(0, 17, __pyx_L3_error)\n",
       "        }\n",
       "        case  3:\n",
       "        if (likely((values[3] = <span class='py_c_api'>PyDict_GetItem</span>(__pyx_kwds, __pyx_n_s_cur_conds)) != 0)) kw_args--;\n",
       "        else {\n",
       "          <span class='pyx_c_api'>__Pyx_RaiseArgtupleInvalid</span>(\"_compute_constraints\", 1, 7, 7, 3); __PYX_ERR(0, 17, __pyx_L3_error)\n",
       "        }\n",
       "        case  4:\n",
       "        if (likely((values[4] = <span class='py_c_api'>PyDict_GetItem</span>(__pyx_kwds, __pyx_n_s_site_fracs)) != 0)) kw_args--;\n",
       "        else {\n",
       "          <span class='pyx_c_api'>__Pyx_RaiseArgtupleInvalid</span>(\"_compute_constraints\", 1, 7, 7, 4); __PYX_ERR(0, 17, __pyx_L3_error)\n",
       "        }\n",
       "        case  5:\n",
       "        if (likely((values[5] = <span class='py_c_api'>PyDict_GetItem</span>(__pyx_kwds, __pyx_n_s_phase_fracs)) != 0)) kw_args--;\n",
       "        else {\n",
       "          <span class='pyx_c_api'>__Pyx_RaiseArgtupleInvalid</span>(\"_compute_constraints\", 1, 7, 7, 5); __PYX_ERR(0, 17, __pyx_L3_error)\n",
       "        }\n",
       "        case  6:\n",
       "        if (likely((values[6] = <span class='py_c_api'>PyDict_GetItem</span>(__pyx_kwds, __pyx_n_s_phase_records)) != 0)) kw_args--;\n",
       "        else {\n",
       "          <span class='pyx_c_api'>__Pyx_RaiseArgtupleInvalid</span>(\"_compute_constraints\", 1, 7, 7, 6); __PYX_ERR(0, 17, __pyx_L3_error)\n",
       "        }\n",
       "      }\n",
       "      if (unlikely(kw_args &gt; 0)) {\n",
       "        if (unlikely(<span class='pyx_c_api'>__Pyx_ParseOptionalKeywords</span>(__pyx_kwds, __pyx_pyargnames, 0, values, pos_args, \"_compute_constraints\") &lt; 0)) __PYX_ERR(0, 17, __pyx_L3_error)\n",
       "      }\n",
       "    } else if (<span class='py_macro_api'>PyTuple_GET_SIZE</span>(__pyx_args) != 7) {\n",
       "      goto __pyx_L5_argtuple_error;\n",
       "    } else {\n",
       "      values[0] = <span class='py_macro_api'>PyTuple_GET_ITEM</span>(__pyx_args, 0);\n",
       "      values[1] = <span class='py_macro_api'>PyTuple_GET_ITEM</span>(__pyx_args, 1);\n",
       "      values[2] = <span class='py_macro_api'>PyTuple_GET_ITEM</span>(__pyx_args, 2);\n",
       "      values[3] = <span class='py_macro_api'>PyTuple_GET_ITEM</span>(__pyx_args, 3);\n",
       "      values[4] = <span class='py_macro_api'>PyTuple_GET_ITEM</span>(__pyx_args, 4);\n",
       "      values[5] = <span class='py_macro_api'>PyTuple_GET_ITEM</span>(__pyx_args, 5);\n",
       "      values[6] = <span class='py_macro_api'>PyTuple_GET_ITEM</span>(__pyx_args, 6);\n",
       "    }\n",
       "    __pyx_v_dbf = values[0];\n",
       "    __pyx_v_comps = values[1];\n",
       "    __pyx_v_phases = values[2];\n",
       "    __pyx_v_cur_conds = values[3];\n",
       "    __pyx_v_site_fracs = ((PyArrayObject *)values[4]);\n",
       "    __pyx_v_phase_fracs = ((PyArrayObject *)values[5]);\n",
       "    __pyx_v_phase_records = values[6];\n",
       "  }\n",
       "  goto __pyx_L4_argument_unpacking_done;\n",
       "  __pyx_L5_argtuple_error:;\n",
       "  <span class='pyx_c_api'>__Pyx_RaiseArgtupleInvalid</span>(\"_compute_constraints\", 1, 7, 7, <span class='py_macro_api'>PyTuple_GET_SIZE</span>(__pyx_args)); __PYX_ERR(0, 17, __pyx_L3_error)\n",
       "  __pyx_L3_error:;\n",
       "  <span class='pyx_c_api'>__Pyx_AddTraceback</span>(\"_cython_magic_07b8d211ecceae831b3e5546c95de0f8._compute_constraints\", __pyx_clineno, __pyx_lineno, __pyx_filename);\n",
       "  <span class='refnanny'>__Pyx_RefNannyFinishContext</span>();\n",
       "  return NULL;\n",
       "  __pyx_L4_argument_unpacking_done:;\n",
       "  if (unlikely(!<span class='pyx_c_api'>__Pyx_ArgTypeTest</span>(((PyObject *)__pyx_v_site_fracs), __pyx_ptype_5numpy_ndarray, 1, \"site_fracs\", 0))) __PYX_ERR(0, 18, __pyx_L1_error)\n",
       "  if (unlikely(!<span class='pyx_c_api'>__Pyx_ArgTypeTest</span>(((PyObject *)__pyx_v_phase_fracs), __pyx_ptype_5numpy_ndarray, 1, \"phase_fracs\", 0))) __PYX_ERR(0, 19, __pyx_L1_error)\n",
       "  __pyx_r = __pyx_pf_46_cython_magic_07b8d211ecceae831b3e5546c95de0f8__compute_constraints(__pyx_self, __pyx_v_dbf, __pyx_v_comps, __pyx_v_phases, __pyx_v_cur_conds, __pyx_v_site_fracs, __pyx_v_phase_fracs, __pyx_v_phase_records);\n",
       "\n",
       "  /* function exit code */\n",
       "  goto __pyx_L0;\n",
       "  __pyx_L1_error:;\n",
       "  __pyx_r = NULL;\n",
       "  __pyx_L0:;\n",
       "  <span class='refnanny'>__Pyx_RefNannyFinishContext</span>();\n",
       "  return __pyx_r;\n",
       "}\n",
       "\n",
       "static PyObject *__pyx_pf_46_cython_magic_07b8d211ecceae831b3e5546c95de0f8__compute_constraints(CYTHON_UNUSED PyObject *__pyx_self, PyObject *__pyx_v_dbf, PyObject *__pyx_v_comps, PyObject *__pyx_v_phases, PyObject *__pyx_v_cur_conds, PyArrayObject *__pyx_v_site_fracs, PyArrayObject *__pyx_v_phase_fracs, PyObject *__pyx_v_phase_records) {\n",
       "  int __pyx_v_num_sitefrac_bals;\n",
       "  int __pyx_v_num_mass_bals;\n",
       "  double __pyx_v_indep_sum;\n",
       "  __Pyx_memviewslice __pyx_v_comp_obj_value = { 0, 0, { 0 }, { 0 }, { 0 } };\n",
       "  PyObject *__pyx_v_dependent_comp = 0;\n",
       "  int __pyx_v_num_constraints;\n",
       "  int __pyx_v_num_phases;\n",
       "  int __pyx_v_num_vars;\n",
       "  PyArrayObject *__pyx_v_phase_dof = 0;\n",
       "  PyArrayObject *__pyx_v_l_constraints = 0;\n",
       "  PyArrayObject *__pyx_v_constraint_jac = 0;\n",
       "  PyArrayObject *__pyx_v_constraint_hess = 0;\n",
       "  __Pyx_memviewslice __pyx_v_sfview = { 0, 0, { 0 }, { 0 }, { 0 } };\n",
       "  __Pyx_memviewslice __pyx_v_comp_grad_value = { 0, 0, { 0 }, { 0 }, { 0 } };\n",
       "  __Pyx_memviewslice __pyx_v_comp_hess_value = { 0, 0, { 0 }, { 0 }, { 0 } };\n",
       "  int __pyx_v_phase_idx;\n",
       "  int __pyx_v_var_offset;\n",
       "  int __pyx_v_constraint_offset;\n",
       "  int __pyx_v_var_idx;\n",
       "  int __pyx_v_grad_idx;\n",
       "  int __pyx_v_hess_idx;\n",
       "  int __pyx_v_comp_idx;\n",
       "  int __pyx_v_idx;\n",
       "  struct PhaseRecordObject *__pyx_v_prn = 0;\n",
       "  double __pyx_v_phase_frac;\n",
       "  PyObject *__pyx_v_name = NULL;\n",
       "  PyObject *__pyx_v_active_in_subl = NULL;\n",
       "  PyObject *__pyx_v_ais_len = NULL;\n",
       "  PyObject *__pyx_v_comp = NULL;\n",
       "  Py_ssize_t __pyx_v_spidx;\n",
       "  __Pyx_LocalBuf_ND __pyx_pybuffernd_constraint_hess;\n",
       "  __Pyx_Buffer __pyx_pybuffer_constraint_hess;\n",
       "  __Pyx_LocalBuf_ND __pyx_pybuffernd_constraint_jac;\n",
       "  __Pyx_Buffer __pyx_pybuffer_constraint_jac;\n",
       "  __Pyx_LocalBuf_ND __pyx_pybuffernd_l_constraints;\n",
       "  __Pyx_Buffer __pyx_pybuffer_l_constraints;\n",
       "  __Pyx_LocalBuf_ND __pyx_pybuffernd_phase_dof;\n",
       "  __Pyx_Buffer __pyx_pybuffer_phase_dof;\n",
       "  __Pyx_LocalBuf_ND __pyx_pybuffernd_phase_fracs;\n",
       "  __Pyx_Buffer __pyx_pybuffer_phase_fracs;\n",
       "  __Pyx_LocalBuf_ND __pyx_pybuffernd_site_fracs;\n",
       "  __Pyx_Buffer __pyx_pybuffer_site_fracs;\n",
       "  PyObject *__pyx_r = NULL;\n",
       "  <span class='trace'>__Pyx_TraceDeclarations</span>\n",
       "  <span class='trace'>__Pyx_TraceFrameInit</span>(__pyx_codeobj_)\n",
       "  <span class='refnanny'>__Pyx_RefNannyDeclarations</span>\n",
       "  <span class='refnanny'>__Pyx_RefNannySetupContext</span>(\"_compute_constraints\", 0);\n",
       "  <span class='trace'>__Pyx_TraceCall</span>(\"_compute_constraints\", __pyx_f[0], 17, 0, __PYX_ERR(0, 17, __pyx_L1_error));\n",
       "  __pyx_pybuffer_phase_dof.pybuffer.buf = NULL;\n",
       "  __pyx_pybuffer_phase_dof.refcount = 0;\n",
       "  __pyx_pybuffernd_phase_dof.data = NULL;\n",
       "  __pyx_pybuffernd_phase_dof.rcbuffer = &amp;__pyx_pybuffer_phase_dof;\n",
       "  __pyx_pybuffer_l_constraints.pybuffer.buf = NULL;\n",
       "  __pyx_pybuffer_l_constraints.refcount = 0;\n",
       "  __pyx_pybuffernd_l_constraints.data = NULL;\n",
       "  __pyx_pybuffernd_l_constraints.rcbuffer = &amp;__pyx_pybuffer_l_constraints;\n",
       "  __pyx_pybuffer_constraint_jac.pybuffer.buf = NULL;\n",
       "  __pyx_pybuffer_constraint_jac.refcount = 0;\n",
       "  __pyx_pybuffernd_constraint_jac.data = NULL;\n",
       "  __pyx_pybuffernd_constraint_jac.rcbuffer = &amp;__pyx_pybuffer_constraint_jac;\n",
       "  __pyx_pybuffer_constraint_hess.pybuffer.buf = NULL;\n",
       "  __pyx_pybuffer_constraint_hess.refcount = 0;\n",
       "  __pyx_pybuffernd_constraint_hess.data = NULL;\n",
       "  __pyx_pybuffernd_constraint_hess.rcbuffer = &amp;__pyx_pybuffer_constraint_hess;\n",
       "  __pyx_pybuffer_site_fracs.pybuffer.buf = NULL;\n",
       "  __pyx_pybuffer_site_fracs.refcount = 0;\n",
       "  __pyx_pybuffernd_site_fracs.data = NULL;\n",
       "  __pyx_pybuffernd_site_fracs.rcbuffer = &amp;__pyx_pybuffer_site_fracs;\n",
       "  __pyx_pybuffer_phase_fracs.pybuffer.buf = NULL;\n",
       "  __pyx_pybuffer_phase_fracs.refcount = 0;\n",
       "  __pyx_pybuffernd_phase_fracs.data = NULL;\n",
       "  __pyx_pybuffernd_phase_fracs.rcbuffer = &amp;__pyx_pybuffer_phase_fracs;\n",
       "  {\n",
       "    __Pyx_BufFmt_StackElem __pyx_stack[1];\n",
       "    if (unlikely(<span class='pyx_c_api'>__Pyx_GetBufferAndValidate</span>(&amp;__pyx_pybuffernd_site_fracs.rcbuffer-&gt;pybuffer, (PyObject*)__pyx_v_site_fracs, &amp;__Pyx_TypeInfo_nn___pyx_t_5numpy_float64_t, PyBUF_FORMAT| PyBUF_STRIDES, 1, 0, __pyx_stack) == -1)) __PYX_ERR(0, 17, __pyx_L1_error)\n",
       "  }\n",
       "  __pyx_pybuffernd_site_fracs.diminfo[0].strides = __pyx_pybuffernd_site_fracs.rcbuffer-&gt;pybuffer.strides[0]; __pyx_pybuffernd_site_fracs.diminfo[0].shape = __pyx_pybuffernd_site_fracs.rcbuffer-&gt;pybuffer.shape[0];\n",
       "  {\n",
       "    __Pyx_BufFmt_StackElem __pyx_stack[1];\n",
       "    if (unlikely(<span class='pyx_c_api'>__Pyx_GetBufferAndValidate</span>(&amp;__pyx_pybuffernd_phase_fracs.rcbuffer-&gt;pybuffer, (PyObject*)__pyx_v_phase_fracs, &amp;__Pyx_TypeInfo_nn___pyx_t_5numpy_float64_t, PyBUF_FORMAT| PyBUF_STRIDES, 1, 0, __pyx_stack) == -1)) __PYX_ERR(0, 17, __pyx_L1_error)\n",
       "  }\n",
       "  __pyx_pybuffernd_phase_fracs.diminfo[0].strides = __pyx_pybuffernd_phase_fracs.rcbuffer-&gt;pybuffer.strides[0]; __pyx_pybuffernd_phase_fracs.diminfo[0].shape = __pyx_pybuffernd_phase_fracs.rcbuffer-&gt;pybuffer.shape[0];\n",
       "/* … */\n",
       "  /* function exit code */\n",
       "  __pyx_L1_error:;\n",
       "  <span class='pyx_macro_api'>__Pyx_XDECREF</span>(__pyx_t_1);\n",
       "  <span class='pyx_macro_api'>__Pyx_XDECREF</span>(__pyx_t_2);\n",
       "  <span class='pyx_macro_api'>__Pyx_XDECREF</span>(__pyx_t_5);\n",
       "  <span class='pyx_macro_api'>__Pyx_XDECREF</span>(__pyx_t_6);\n",
       "  <span class='pyx_macro_api'>__Pyx_XDECREF</span>(__pyx_t_12);\n",
       "  __PYX_XDEC_MEMVIEW(&amp;__pyx_t_13, 1);\n",
       "  __PYX_XDEC_MEMVIEW(&amp;__pyx_t_26, 1);\n",
       "  <span class='pyx_macro_api'>__Pyx_XDECREF</span>(__pyx_t_28);\n",
       "  { PyObject *__pyx_type, *__pyx_value, *__pyx_tb;\n",
       "    __Pyx_PyThreadState_declare\n",
       "    __Pyx_PyThreadState_assign\n",
       "    <span class='pyx_c_api'>__Pyx_ErrFetch</span>(&amp;__pyx_type, &amp;__pyx_value, &amp;__pyx_tb);\n",
       "    <span class='pyx_c_api'>__Pyx_SafeReleaseBuffer</span>(&amp;__pyx_pybuffernd_constraint_hess.rcbuffer-&gt;pybuffer);\n",
       "    <span class='pyx_c_api'>__Pyx_SafeReleaseBuffer</span>(&amp;__pyx_pybuffernd_constraint_jac.rcbuffer-&gt;pybuffer);\n",
       "    <span class='pyx_c_api'>__Pyx_SafeReleaseBuffer</span>(&amp;__pyx_pybuffernd_l_constraints.rcbuffer-&gt;pybuffer);\n",
       "    <span class='pyx_c_api'>__Pyx_SafeReleaseBuffer</span>(&amp;__pyx_pybuffernd_phase_dof.rcbuffer-&gt;pybuffer);\n",
       "    <span class='pyx_c_api'>__Pyx_SafeReleaseBuffer</span>(&amp;__pyx_pybuffernd_phase_fracs.rcbuffer-&gt;pybuffer);\n",
       "    <span class='pyx_c_api'>__Pyx_SafeReleaseBuffer</span>(&amp;__pyx_pybuffernd_site_fracs.rcbuffer-&gt;pybuffer);\n",
       "  <span class='pyx_c_api'>__Pyx_ErrRestore</span>(__pyx_type, __pyx_value, __pyx_tb);}\n",
       "  <span class='pyx_c_api'>__Pyx_AddTraceback</span>(\"_cython_magic_07b8d211ecceae831b3e5546c95de0f8._compute_constraints\", __pyx_clineno, __pyx_lineno, __pyx_filename);\n",
       "  __pyx_r = NULL;\n",
       "  goto __pyx_L2;\n",
       "  __pyx_L0:;\n",
       "  <span class='pyx_c_api'>__Pyx_SafeReleaseBuffer</span>(&amp;__pyx_pybuffernd_constraint_hess.rcbuffer-&gt;pybuffer);\n",
       "  <span class='pyx_c_api'>__Pyx_SafeReleaseBuffer</span>(&amp;__pyx_pybuffernd_constraint_jac.rcbuffer-&gt;pybuffer);\n",
       "  <span class='pyx_c_api'>__Pyx_SafeReleaseBuffer</span>(&amp;__pyx_pybuffernd_l_constraints.rcbuffer-&gt;pybuffer);\n",
       "  <span class='pyx_c_api'>__Pyx_SafeReleaseBuffer</span>(&amp;__pyx_pybuffernd_phase_dof.rcbuffer-&gt;pybuffer);\n",
       "  <span class='pyx_c_api'>__Pyx_SafeReleaseBuffer</span>(&amp;__pyx_pybuffernd_phase_fracs.rcbuffer-&gt;pybuffer);\n",
       "  <span class='pyx_c_api'>__Pyx_SafeReleaseBuffer</span>(&amp;__pyx_pybuffernd_site_fracs.rcbuffer-&gt;pybuffer);\n",
       "  __pyx_L2:;\n",
       "  __PYX_XDEC_MEMVIEW(&amp;__pyx_v_comp_obj_value, 1);\n",
       "  <span class='pyx_macro_api'>__Pyx_XDECREF</span>(__pyx_v_dependent_comp);\n",
       "  <span class='pyx_macro_api'>__Pyx_XDECREF</span>((PyObject *)__pyx_v_phase_dof);\n",
       "  <span class='pyx_macro_api'>__Pyx_XDECREF</span>((PyObject *)__pyx_v_l_constraints);\n",
       "  <span class='pyx_macro_api'>__Pyx_XDECREF</span>((PyObject *)__pyx_v_constraint_jac);\n",
       "  <span class='pyx_macro_api'>__Pyx_XDECREF</span>((PyObject *)__pyx_v_constraint_hess);\n",
       "  __PYX_XDEC_MEMVIEW(&amp;__pyx_v_sfview, 1);\n",
       "  __PYX_XDEC_MEMVIEW(&amp;__pyx_v_comp_grad_value, 1);\n",
       "  __PYX_XDEC_MEMVIEW(&amp;__pyx_v_comp_hess_value, 1);\n",
       "  <span class='pyx_macro_api'>__Pyx_XDECREF</span>((PyObject *)__pyx_v_prn);\n",
       "  <span class='pyx_macro_api'>__Pyx_XDECREF</span>(__pyx_v_name);\n",
       "  <span class='pyx_macro_api'>__Pyx_XDECREF</span>(__pyx_v_active_in_subl);\n",
       "  <span class='pyx_macro_api'>__Pyx_XDECREF</span>(__pyx_v_ais_len);\n",
       "  <span class='pyx_macro_api'>__Pyx_XDECREF</span>(__pyx_v_comp);\n",
       "  <span class='refnanny'>__Pyx_XGIVEREF</span>(__pyx_r);\n",
       "  <span class='trace'>__Pyx_TraceReturn</span>(__pyx_r, 0);\n",
       "  <span class='refnanny'>__Pyx_RefNannyFinishContext</span>();\n",
       "  return __pyx_r;\n",
       "}\n",
       "/* … */\n",
       "  __pyx_tuple__34 = <span class='py_c_api'>PyTuple_Pack</span>(38, __pyx_n_s_dbf, __pyx_n_s_comps, __pyx_n_s_phases, __pyx_n_s_cur_conds, __pyx_n_s_site_fracs, __pyx_n_s_phase_fracs, __pyx_n_s_phase_records, __pyx_n_s_num_sitefrac_bals, __pyx_n_s_num_mass_bals, __pyx_n_s_indep_sum, __pyx_n_s_comp_obj_value, __pyx_n_s_dependent_comp, __pyx_n_s_num_constraints, __pyx_n_s_num_phases, __pyx_n_s_num_vars, __pyx_n_s_phase_dof, __pyx_n_s_l_constraints, __pyx_n_s_constraint_jac, __pyx_n_s_constraint_hess, __pyx_n_s_sfview, __pyx_n_s_comp_grad_value, __pyx_n_s_comp_hess_value, __pyx_n_s_phase_idx, __pyx_n_s_var_offset, __pyx_n_s_constraint_offset, __pyx_n_s_var_idx, __pyx_n_s_iter_idx, __pyx_n_s_grad_idx, __pyx_n_s_hess_idx, __pyx_n_s_comp_idx, __pyx_n_s_idx, __pyx_n_s_prn, __pyx_n_s_phase_frac, __pyx_n_s_name, __pyx_n_s_active_in_subl, __pyx_n_s_ais_len, __pyx_n_s_comp, __pyx_n_s_spidx); if (unlikely(!__pyx_tuple__34)) __PYX_ERR(0, 17, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_tuple__34);\n",
       "  <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_tuple__34);\n",
       "/* … */\n",
       "  <span class='trace'>__Pyx_TraceLine</span>(17,0,__PYX_ERR(0, 17, __pyx_L1_error))\n",
       "  __pyx_t_2 = <span class='pyx_c_api'>__Pyx_CyFunction_NewEx</span>(&amp;__pyx_mdef_46_cython_magic_07b8d211ecceae831b3e5546c95de0f8_1_compute_constraints, 0, __pyx_n_s_compute_constraints, NULL, __pyx_n_s_cython_magic_07b8d211ecceae831b, __pyx_d, ((PyObject *)__pyx_codeobj_)); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 17, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "  if (<span class='py_c_api'>PyDict_SetItem</span>(__pyx_d, __pyx_n_s_compute_constraints, __pyx_t_2) &lt; 0) __PYX_ERR(0, 17, __pyx_L1_error)\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "  __pyx_codeobj_ = (PyObject*)<span class='pyx_c_api'>__Pyx_PyCode_New</span>(7, 0, 38, 0, 0, __pyx_empty_bytes, __pyx_empty_tuple, __pyx_empty_tuple, __pyx_tuple__34, __pyx_empty_tuple, __pyx_empty_tuple, __pyx_kp_s_home_rotis_cache_ipython_cython, __pyx_n_s_compute_constraints, 17, __pyx_empty_bytes); if (unlikely(!__pyx_codeobj_)) __PYX_ERR(0, 17, __pyx_L1_error)\n",
       "</pre><pre class=\"cython line score-0\">&#xA0;<span class=\"\">18</span>:                          <span class=\"nb\">object</span> <span class=\"n\">cur_conds</span><span class=\"p\">,</span> <span class=\"n\">np</span><span class=\"o\">.</span><span class=\"n\">ndarray</span><span class=\"p\">[</span><span class=\"n\">dtype</span><span class=\"o\">=</span><span class=\"n\">np</span><span class=\"o\">.</span><span class=\"n\">float64_t</span><span class=\"p\">,</span> <span class=\"n\">ndim</span><span class=\"o\">=</span><span class=\"mf\">1</span><span class=\"p\">]</span> <span class=\"n\">site_fracs</span><span class=\"p\">,</span></pre>\n",
       "<pre class=\"cython line score-0\">&#xA0;<span class=\"\">19</span>:                          <span class=\"n\">np</span><span class=\"o\">.</span><span class=\"n\">ndarray</span><span class=\"p\">[</span><span class=\"n\">dtype</span><span class=\"o\">=</span><span class=\"n\">np</span><span class=\"o\">.</span><span class=\"n\">float64_t</span><span class=\"p\">,</span> <span class=\"n\">ndim</span><span class=\"o\">=</span><span class=\"mf\">1</span><span class=\"p\">]</span> <span class=\"n\">phase_fracs</span><span class=\"p\">,</span> <span class=\"nb\">object</span> <span class=\"n\">phase_records</span><span class=\"p\">):</span></pre>\n",
       "<pre class=\"cython line score-0\">&#xA0;<span class=\"\">20</span>:     <span class=\"sd\">&quot;&quot;&quot;</span></pre>\n",
       "<pre class=\"cython line score-0\">&#xA0;<span class=\"\">21</span>: <span class=\"sd\">    Compute the constraint vector and constraint Jacobian matrix.</span></pre>\n",
       "<pre class=\"cython line score-0\">&#xA0;<span class=\"\">22</span>: <span class=\"sd\">    &quot;&quot;&quot;</span></pre>\n",
       "<pre class=\"cython line score-95\" onclick='toggleDiv(this)'>+<span class=\"\">23</span>:     <span class=\"k\">cdef</span> <span class=\"kt\">int</span> <span class=\"nf\">num_sitefrac_bals</span> <span class=\"o\">=</span> <span class=\"nb\">sum</span><span class=\"p\">([</span><span class=\"nb\">len</span><span class=\"p\">(</span><span class=\"n\">dbf</span><span class=\"o\">.</span><span class=\"n\">phases</span><span class=\"p\">[</span><span class=\"n\">i</span><span class=\"p\">]</span><span class=\"o\">.</span><span class=\"n\">sublattices</span><span class=\"p\">)</span> <span class=\"k\">for</span> <span class=\"n\">i</span> <span class=\"ow\">in</span> <span class=\"n\">phases</span><span class=\"p\">])</span></pre>\n",
       "<pre class='cython code score-95 '>  <span class='trace'>__Pyx_TraceLine</span>(23,0,__PYX_ERR(0, 23, __pyx_L1_error))\n",
       "  { /* enter inner scope */\n",
       "    PyObject *__pyx_7genexpr__pyx_v_i = NULL;\n",
       "    __pyx_t_1 = <span class='py_c_api'>PyList_New</span>(0); if (unlikely(!__pyx_t_1)) __PYX_ERR(0, 23, __pyx_L5_error)\n",
       "    <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_1);\n",
       "    if (likely(<span class='py_c_api'>PyList_CheckExact</span>(__pyx_v_phases)) || <span class='py_c_api'>PyTuple_CheckExact</span>(__pyx_v_phases)) {\n",
       "      __pyx_t_2 = __pyx_v_phases; <span class='pyx_macro_api'>__Pyx_INCREF</span>(__pyx_t_2); __pyx_t_3 = 0;\n",
       "      __pyx_t_4 = NULL;\n",
       "    } else {\n",
       "      __pyx_t_3 = -1; __pyx_t_2 = <span class='py_c_api'>PyObject_GetIter</span>(__pyx_v_phases); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 23, __pyx_L5_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "      __pyx_t_4 = Py_TYPE(__pyx_t_2)-&gt;tp_iternext; if (unlikely(!__pyx_t_4)) __PYX_ERR(0, 23, __pyx_L5_error)\n",
       "    }\n",
       "    for (;;) {\n",
       "      if (likely(!__pyx_t_4)) {\n",
       "        if (likely(<span class='py_c_api'>PyList_CheckExact</span>(__pyx_t_2))) {\n",
       "          if (__pyx_t_3 &gt;= <span class='py_macro_api'>PyList_GET_SIZE</span>(__pyx_t_2)) break;\n",
       "          #if CYTHON_ASSUME_SAFE_MACROS &amp;&amp; !CYTHON_AVOID_BORROWED_REFS\n",
       "          __pyx_t_5 = <span class='py_macro_api'>PyList_GET_ITEM</span>(__pyx_t_2, __pyx_t_3); <span class='pyx_macro_api'>__Pyx_INCREF</span>(__pyx_t_5); __pyx_t_3++; if (unlikely(0 &lt; 0)) __PYX_ERR(0, 23, __pyx_L5_error)\n",
       "          #else\n",
       "          __pyx_t_5 = <span class='py_macro_api'>PySequence_ITEM</span>(__pyx_t_2, __pyx_t_3); __pyx_t_3++; if (unlikely(!__pyx_t_5)) __PYX_ERR(0, 23, __pyx_L5_error)\n",
       "          <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_5);\n",
       "          #endif\n",
       "        } else {\n",
       "          if (__pyx_t_3 &gt;= <span class='py_macro_api'>PyTuple_GET_SIZE</span>(__pyx_t_2)) break;\n",
       "          #if CYTHON_ASSUME_SAFE_MACROS &amp;&amp; !CYTHON_AVOID_BORROWED_REFS\n",
       "          __pyx_t_5 = <span class='py_macro_api'>PyTuple_GET_ITEM</span>(__pyx_t_2, __pyx_t_3); <span class='pyx_macro_api'>__Pyx_INCREF</span>(__pyx_t_5); __pyx_t_3++; if (unlikely(0 &lt; 0)) __PYX_ERR(0, 23, __pyx_L5_error)\n",
       "          #else\n",
       "          __pyx_t_5 = <span class='py_macro_api'>PySequence_ITEM</span>(__pyx_t_2, __pyx_t_3); __pyx_t_3++; if (unlikely(!__pyx_t_5)) __PYX_ERR(0, 23, __pyx_L5_error)\n",
       "          <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_5);\n",
       "          #endif\n",
       "        }\n",
       "      } else {\n",
       "        __pyx_t_5 = __pyx_t_4(__pyx_t_2);\n",
       "        if (unlikely(!__pyx_t_5)) {\n",
       "          PyObject* exc_type = <span class='py_c_api'>PyErr_Occurred</span>();\n",
       "          if (exc_type) {\n",
       "            if (likely(exc_type == PyExc_StopIteration || <span class='py_c_api'>PyErr_GivenExceptionMatches</span>(exc_type, PyExc_StopIteration))) <span class='py_c_api'>PyErr_Clear</span>();\n",
       "            else __PYX_ERR(0, 23, __pyx_L5_error)\n",
       "          }\n",
       "          break;\n",
       "        }\n",
       "        <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_5);\n",
       "      }\n",
       "      <span class='pyx_macro_api'>__Pyx_XDECREF_SET</span>(__pyx_7genexpr__pyx_v_i, __pyx_t_5);\n",
       "      __pyx_t_5 = 0;\n",
       "      __pyx_t_5 = <span class='pyx_c_api'>__Pyx_PyObject_GetAttrStr</span>(__pyx_v_dbf, __pyx_n_s_phases); if (unlikely(!__pyx_t_5)) __PYX_ERR(0, 23, __pyx_L5_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_5);\n",
       "      __pyx_t_6 = <span class='py_c_api'>PyObject_GetItem</span>(__pyx_t_5, __pyx_7genexpr__pyx_v_i); if (unlikely(!__pyx_t_6)) __PYX_ERR(0, 23, __pyx_L5_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_6);\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_5); __pyx_t_5 = 0;\n",
       "      __pyx_t_5 = <span class='pyx_c_api'>__Pyx_PyObject_GetAttrStr</span>(__pyx_t_6, __pyx_n_s_sublattices); if (unlikely(!__pyx_t_5)) __PYX_ERR(0, 23, __pyx_L5_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_5);\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_6); __pyx_t_6 = 0;\n",
       "      __pyx_t_7 = <span class='py_c_api'>PyObject_Length</span>(__pyx_t_5); if (unlikely(__pyx_t_7 == -1)) __PYX_ERR(0, 23, __pyx_L5_error)\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_5); __pyx_t_5 = 0;\n",
       "      __pyx_t_5 = <span class='py_c_api'>PyInt_FromSsize_t</span>(__pyx_t_7); if (unlikely(!__pyx_t_5)) __PYX_ERR(0, 23, __pyx_L5_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_5);\n",
       "      if (unlikely(<span class='pyx_c_api'>__Pyx_ListComp_Append</span>(__pyx_t_1, (PyObject*)__pyx_t_5))) __PYX_ERR(0, 23, __pyx_L5_error)\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_5); __pyx_t_5 = 0;\n",
       "    }\n",
       "    <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "    <span class='pyx_macro_api'>__Pyx_XDECREF</span>(__pyx_7genexpr__pyx_v_i);\n",
       "    goto __pyx_L8_exit_scope;\n",
       "    __pyx_L5_error:;\n",
       "    <span class='pyx_macro_api'>__Pyx_XDECREF</span>(__pyx_7genexpr__pyx_v_i);\n",
       "    goto __pyx_L1_error;\n",
       "    __pyx_L8_exit_scope:;\n",
       "  } /* exit inner scope */\n",
       "  __pyx_t_2 = <span class='py_c_api'>PyTuple_New</span>(1); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 23, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "  <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_t_1);\n",
       "  <span class='py_macro_api'>PyTuple_SET_ITEM</span>(__pyx_t_2, 0, __pyx_t_1);\n",
       "  __pyx_t_1 = 0;\n",
       "  __pyx_t_1 = <span class='pyx_c_api'>__Pyx_PyObject_Call</span>(__pyx_builtin_sum, __pyx_t_2, NULL); if (unlikely(!__pyx_t_1)) __PYX_ERR(0, 23, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_1);\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "  __pyx_t_8 = <span class='pyx_c_api'>__Pyx_PyInt_As_int</span>(__pyx_t_1); if (unlikely((__pyx_t_8 == (int)-1) &amp;&amp; <span class='py_c_api'>PyErr_Occurred</span>())) __PYX_ERR(0, 23, __pyx_L1_error)\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_1); __pyx_t_1 = 0;\n",
       "  __pyx_v_num_sitefrac_bals = __pyx_t_8;\n",
       "</pre><pre class=\"cython line score-32\" onclick='toggleDiv(this)'>+<span class=\"\">24</span>:     <span class=\"k\">cdef</span> <span class=\"kt\">int</span> <span class=\"nf\">num_mass_bals</span> <span class=\"o\">=</span> <span class=\"nb\">len</span><span class=\"p\">([</span><span class=\"n\">i</span> <span class=\"k\">for</span> <span class=\"n\">i</span> <span class=\"ow\">in</span> <span class=\"n\">cur_conds</span><span class=\"o\">.</span><span class=\"n\">keys</span><span class=\"p\">()</span> <span class=\"k\">if</span> <span class=\"n\">i</span><span class=\"o\">.</span><span class=\"n\">startswith</span><span class=\"p\">(</span><span class=\"s\">&#39;X_&#39;</span><span class=\"p\">)])</span> <span class=\"o\">+</span> <span class=\"mf\">1</span></pre>\n",
       "<pre class='cython code score-32 '>  <span class='trace'>__Pyx_TraceLine</span>(24,0,__PYX_ERR(0, 24, __pyx_L1_error))\n",
       "  { /* enter inner scope */\n",
       "    PyObject *__pyx_8genexpr1__pyx_v_i = NULL;\n",
       "    __pyx_t_1 = <span class='py_c_api'>PyList_New</span>(0); if (unlikely(!__pyx_t_1)) __PYX_ERR(0, 24, __pyx_L11_error)\n",
       "    <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_1);\n",
       "    __pyx_t_3 = 0;\n",
       "    if (unlikely(__pyx_v_cur_conds == Py_None)) {\n",
       "      <span class='py_c_api'>PyErr_Format</span>(PyExc_AttributeError, \"'NoneType' object has no attribute '%s'\", \"keys\");\n",
       "      __PYX_ERR(0, 24, __pyx_L11_error)\n",
       "    }\n",
       "    __pyx_t_5 = __Pyx_dict_iterator(__pyx_v_cur_conds, 0, __pyx_n_s_keys, (&amp;__pyx_t_7), (&amp;__pyx_t_8)); if (unlikely(!__pyx_t_5)) __PYX_ERR(0, 24, __pyx_L11_error)\n",
       "    <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_5);\n",
       "    <span class='pyx_macro_api'>__Pyx_XDECREF</span>(__pyx_t_2);\n",
       "    __pyx_t_2 = __pyx_t_5;\n",
       "    __pyx_t_5 = 0;\n",
       "    while (1) {\n",
       "      __pyx_t_9 = __Pyx_dict_iter_next(__pyx_t_2, __pyx_t_7, &amp;__pyx_t_3, &amp;__pyx_t_5, NULL, NULL, __pyx_t_8);\n",
       "      if (unlikely(__pyx_t_9 == 0)) break;\n",
       "      if (unlikely(__pyx_t_9 == -1)) __PYX_ERR(0, 24, __pyx_L11_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_5);\n",
       "      <span class='pyx_macro_api'>__Pyx_XDECREF_SET</span>(__pyx_8genexpr1__pyx_v_i, __pyx_t_5);\n",
       "      __pyx_t_5 = 0;\n",
       "      __pyx_t_5 = <span class='pyx_c_api'>__Pyx_PyObject_GetAttrStr</span>(__pyx_8genexpr1__pyx_v_i, __pyx_n_s_startswith); if (unlikely(!__pyx_t_5)) __PYX_ERR(0, 24, __pyx_L11_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_5);\n",
       "      __pyx_t_6 = <span class='pyx_c_api'>__Pyx_PyObject_Call</span>(__pyx_t_5, __pyx_tuple__2, NULL); if (unlikely(!__pyx_t_6)) __PYX_ERR(0, 24, __pyx_L11_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_6);\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_5); __pyx_t_5 = 0;\n",
       "      __pyx_t_10 = <span class='pyx_c_api'>__Pyx_PyObject_IsTrue</span>(__pyx_t_6); if (unlikely(__pyx_t_10 &lt; 0)) __PYX_ERR(0, 24, __pyx_L11_error)\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_6); __pyx_t_6 = 0;\n",
       "      if (__pyx_t_10) {\n",
       "        if (unlikely(<span class='pyx_c_api'>__Pyx_ListComp_Append</span>(__pyx_t_1, (PyObject*)__pyx_8genexpr1__pyx_v_i))) __PYX_ERR(0, 24, __pyx_L11_error)\n",
       "      }\n",
       "    }\n",
       "    <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "    <span class='pyx_macro_api'>__Pyx_XDECREF</span>(__pyx_8genexpr1__pyx_v_i);\n",
       "    goto __pyx_L15_exit_scope;\n",
       "    __pyx_L11_error:;\n",
       "    <span class='pyx_macro_api'>__Pyx_XDECREF</span>(__pyx_8genexpr1__pyx_v_i);\n",
       "    goto __pyx_L1_error;\n",
       "    __pyx_L15_exit_scope:;\n",
       "  } /* exit inner scope */\n",
       "  __pyx_t_7 = <span class='py_macro_api'>PyList_GET_SIZE</span>(__pyx_t_1); if (unlikely(__pyx_t_7 == -1)) __PYX_ERR(0, 24, __pyx_L1_error)\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_1); __pyx_t_1 = 0;\n",
       "  __pyx_v_num_mass_bals = (__pyx_t_7 + 1);\n",
       "/* … */\n",
       "  __pyx_tuple__2 = <span class='py_c_api'>PyTuple_Pack</span>(1, __pyx_n_u_X); if (unlikely(!__pyx_tuple__2)) __PYX_ERR(0, 24, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_tuple__2);\n",
       "  <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_tuple__2);\n",
       "</pre><pre class=\"cython line score-56\" onclick='toggleDiv(this)'>+<span class=\"\">25</span>:     <span class=\"k\">cdef</span> <span class=\"kt\">double</span> <span class=\"nf\">indep_sum</span> <span class=\"o\">=</span> <span class=\"nb\">sum</span><span class=\"p\">([</span><span class=\"nb\">float</span><span class=\"p\">(</span><span class=\"n\">val</span><span class=\"p\">)</span> <span class=\"k\">for</span> <span class=\"n\">i</span><span class=\"p\">,</span> <span class=\"n\">val</span> <span class=\"ow\">in</span> <span class=\"n\">cur_conds</span><span class=\"o\">.</span><span class=\"n\">items</span><span class=\"p\">()</span> <span class=\"k\">if</span> <span class=\"n\">i</span><span class=\"o\">.</span><span class=\"n\">startswith</span><span class=\"p\">(</span><span class=\"s\">&#39;X_&#39;</span><span class=\"p\">)])</span></pre>\n",
       "<pre class='cython code score-56 '>  <span class='trace'>__Pyx_TraceLine</span>(25,0,__PYX_ERR(0, 25, __pyx_L1_error))\n",
       "  { /* enter inner scope */\n",
       "    PyObject *__pyx_8genexpr2__pyx_v_i = NULL;\n",
       "    PyObject *__pyx_8genexpr2__pyx_v_val = NULL;\n",
       "    __pyx_t_1 = <span class='py_c_api'>PyList_New</span>(0); if (unlikely(!__pyx_t_1)) __PYX_ERR(0, 25, __pyx_L18_error)\n",
       "    <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_1);\n",
       "    __pyx_t_7 = 0;\n",
       "    if (unlikely(__pyx_v_cur_conds == Py_None)) {\n",
       "      <span class='py_c_api'>PyErr_Format</span>(PyExc_AttributeError, \"'NoneType' object has no attribute '%s'\", \"items\");\n",
       "      __PYX_ERR(0, 25, __pyx_L18_error)\n",
       "    }\n",
       "    __pyx_t_6 = __Pyx_dict_iterator(__pyx_v_cur_conds, 0, __pyx_n_s_items, (&amp;__pyx_t_3), (&amp;__pyx_t_8)); if (unlikely(!__pyx_t_6)) __PYX_ERR(0, 25, __pyx_L18_error)\n",
       "    <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_6);\n",
       "    <span class='pyx_macro_api'>__Pyx_XDECREF</span>(__pyx_t_2);\n",
       "    __pyx_t_2 = __pyx_t_6;\n",
       "    __pyx_t_6 = 0;\n",
       "    while (1) {\n",
       "      __pyx_t_9 = __Pyx_dict_iter_next(__pyx_t_2, __pyx_t_3, &amp;__pyx_t_7, &amp;__pyx_t_6, &amp;__pyx_t_5, NULL, __pyx_t_8);\n",
       "      if (unlikely(__pyx_t_9 == 0)) break;\n",
       "      if (unlikely(__pyx_t_9 == -1)) __PYX_ERR(0, 25, __pyx_L18_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_6);\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_5);\n",
       "      <span class='pyx_macro_api'>__Pyx_XDECREF_SET</span>(__pyx_8genexpr2__pyx_v_i, __pyx_t_6);\n",
       "      __pyx_t_6 = 0;\n",
       "      <span class='pyx_macro_api'>__Pyx_XDECREF_SET</span>(__pyx_8genexpr2__pyx_v_val, __pyx_t_5);\n",
       "      __pyx_t_5 = 0;\n",
       "      __pyx_t_5 = <span class='pyx_c_api'>__Pyx_PyObject_GetAttrStr</span>(__pyx_8genexpr2__pyx_v_i, __pyx_n_s_startswith); if (unlikely(!__pyx_t_5)) __PYX_ERR(0, 25, __pyx_L18_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_5);\n",
       "      __pyx_t_6 = <span class='pyx_c_api'>__Pyx_PyObject_Call</span>(__pyx_t_5, __pyx_tuple__3, NULL); if (unlikely(!__pyx_t_6)) __PYX_ERR(0, 25, __pyx_L18_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_6);\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_5); __pyx_t_5 = 0;\n",
       "      __pyx_t_10 = <span class='pyx_c_api'>__Pyx_PyObject_IsTrue</span>(__pyx_t_6); if (unlikely(__pyx_t_10 &lt; 0)) __PYX_ERR(0, 25, __pyx_L18_error)\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_6); __pyx_t_6 = 0;\n",
       "      if (__pyx_t_10) {\n",
       "        __pyx_t_6 = <span class='pyx_c_api'>__Pyx_PyNumber_Float</span>(__pyx_8genexpr2__pyx_v_val); if (unlikely(!__pyx_t_6)) __PYX_ERR(0, 25, __pyx_L18_error)\n",
       "        <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_6);\n",
       "        if (unlikely(<span class='pyx_c_api'>__Pyx_ListComp_Append</span>(__pyx_t_1, (PyObject*)__pyx_t_6))) __PYX_ERR(0, 25, __pyx_L18_error)\n",
       "        <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_6); __pyx_t_6 = 0;\n",
       "      }\n",
       "    }\n",
       "    <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "    <span class='pyx_macro_api'>__Pyx_XDECREF</span>(__pyx_8genexpr2__pyx_v_i);\n",
       "    <span class='pyx_macro_api'>__Pyx_XDECREF</span>(__pyx_8genexpr2__pyx_v_val);\n",
       "    goto __pyx_L22_exit_scope;\n",
       "    __pyx_L18_error:;\n",
       "    <span class='pyx_macro_api'>__Pyx_XDECREF</span>(__pyx_8genexpr2__pyx_v_i);\n",
       "    <span class='pyx_macro_api'>__Pyx_XDECREF</span>(__pyx_8genexpr2__pyx_v_val);\n",
       "    goto __pyx_L1_error;\n",
       "    __pyx_L22_exit_scope:;\n",
       "  } /* exit inner scope */\n",
       "  __pyx_t_2 = <span class='py_c_api'>PyTuple_New</span>(1); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 25, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "  <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_t_1);\n",
       "  <span class='py_macro_api'>PyTuple_SET_ITEM</span>(__pyx_t_2, 0, __pyx_t_1);\n",
       "  __pyx_t_1 = 0;\n",
       "  __pyx_t_1 = <span class='pyx_c_api'>__Pyx_PyObject_Call</span>(__pyx_builtin_sum, __pyx_t_2, NULL); if (unlikely(!__pyx_t_1)) __PYX_ERR(0, 25, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_1);\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "  __pyx_t_11 = __pyx_<span class='py_c_api'>PyFloat_AsDouble</span>(__pyx_t_1); if (unlikely((__pyx_t_11 == (double)-1) &amp;&amp; <span class='py_c_api'>PyErr_Occurred</span>())) __PYX_ERR(0, 25, __pyx_L1_error)\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_1); __pyx_t_1 = 0;\n",
       "  __pyx_v_indep_sum = __pyx_t_11;\n",
       "/* … */\n",
       "  __pyx_tuple__3 = <span class='py_c_api'>PyTuple_Pack</span>(1, __pyx_n_u_X); if (unlikely(!__pyx_tuple__3)) __PYX_ERR(0, 25, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_tuple__3);\n",
       "  <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_tuple__3);\n",
       "</pre><pre class=\"cython line score-60\" onclick='toggleDiv(this)'>+<span class=\"\">26</span>:     <span class=\"k\">cdef</span> <span class=\"kt\">double</span>[<span class=\"p\">::</span><span class=\"mf\">1</span><span class=\"p\">]</span> <span class=\"n\">comp_obj_value</span> <span class=\"o\">=</span> <span class=\"n\">np</span><span class=\"o\">.</span><span class=\"n\">atleast_1d</span><span class=\"p\">(</span><span class=\"n\">np</span><span class=\"o\">.</span><span class=\"n\">zeros</span><span class=\"p\">(</span><span class=\"mf\">1</span><span class=\"p\">))</span></pre>\n",
       "<pre class='cython code score-60 '>  <span class='trace'>__Pyx_TraceLine</span>(26,0,__PYX_ERR(0, 26, __pyx_L1_error))\n",
       "  __pyx_t_2 = <span class='pyx_c_api'>__Pyx_GetModuleGlobalName</span>(__pyx_n_s_np); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 26, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "  __pyx_t_6 = <span class='pyx_c_api'>__Pyx_PyObject_GetAttrStr</span>(__pyx_t_2, __pyx_n_s_atleast_1d); if (unlikely(!__pyx_t_6)) __PYX_ERR(0, 26, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_6);\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "  __pyx_t_2 = <span class='pyx_c_api'>__Pyx_GetModuleGlobalName</span>(__pyx_n_s_np); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 26, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "  __pyx_t_5 = <span class='pyx_c_api'>__Pyx_PyObject_GetAttrStr</span>(__pyx_t_2, __pyx_n_s_zeros); if (unlikely(!__pyx_t_5)) __PYX_ERR(0, 26, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_5);\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "  __pyx_t_2 = <span class='pyx_c_api'>__Pyx_PyObject_Call</span>(__pyx_t_5, __pyx_tuple__4, NULL); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 26, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_5); __pyx_t_5 = 0;\n",
       "  __pyx_t_5 = NULL;\n",
       "  if (CYTHON_UNPACK_METHODS &amp;&amp; unlikely(<span class='py_c_api'>PyMethod_Check</span>(__pyx_t_6))) {\n",
       "    __pyx_t_5 = <span class='py_macro_api'>PyMethod_GET_SELF</span>(__pyx_t_6);\n",
       "    if (likely(__pyx_t_5)) {\n",
       "      PyObject* function = <span class='py_macro_api'>PyMethod_GET_FUNCTION</span>(__pyx_t_6);\n",
       "      <span class='pyx_macro_api'>__Pyx_INCREF</span>(__pyx_t_5);\n",
       "      <span class='pyx_macro_api'>__Pyx_INCREF</span>(function);\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF_SET</span>(__pyx_t_6, function);\n",
       "    }\n",
       "  }\n",
       "  if (!__pyx_t_5) {\n",
       "    __pyx_t_1 = <span class='pyx_c_api'>__Pyx_PyObject_CallOneArg</span>(__pyx_t_6, __pyx_t_2); if (unlikely(!__pyx_t_1)) __PYX_ERR(0, 26, __pyx_L1_error)\n",
       "    <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "    <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_1);\n",
       "  } else {\n",
       "    #if CYTHON_FAST_PYCALL\n",
       "    if (<span class='py_c_api'>PyFunction_Check</span>(__pyx_t_6)) {\n",
       "      PyObject *__pyx_temp[2] = {__pyx_t_5, __pyx_t_2};\n",
       "      __pyx_t_1 = <span class='pyx_c_api'>__Pyx_PyFunction_FastCall</span>(__pyx_t_6, __pyx_temp+1-1, 1+1); if (unlikely(!__pyx_t_1)) __PYX_ERR(0, 26, __pyx_L1_error)\n",
       "      <span class='pyx_macro_api'>__Pyx_XDECREF</span>(__pyx_t_5); __pyx_t_5 = 0;\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_1);\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "    } else\n",
       "    #endif\n",
       "    #if CYTHON_FAST_PYCCALL\n",
       "    if (<span class='pyx_c_api'>__Pyx_PyFastCFunction_Check</span>(__pyx_t_6)) {\n",
       "      PyObject *__pyx_temp[2] = {__pyx_t_5, __pyx_t_2};\n",
       "      __pyx_t_1 = <span class='pyx_c_api'>__Pyx_PyCFunction_FastCall</span>(__pyx_t_6, __pyx_temp+1-1, 1+1); if (unlikely(!__pyx_t_1)) __PYX_ERR(0, 26, __pyx_L1_error)\n",
       "      <span class='pyx_macro_api'>__Pyx_XDECREF</span>(__pyx_t_5); __pyx_t_5 = 0;\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_1);\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "    } else\n",
       "    #endif\n",
       "    {\n",
       "      __pyx_t_12 = <span class='py_c_api'>PyTuple_New</span>(1+1); if (unlikely(!__pyx_t_12)) __PYX_ERR(0, 26, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_12);\n",
       "      <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_t_5); <span class='py_macro_api'>PyTuple_SET_ITEM</span>(__pyx_t_12, 0, __pyx_t_5); __pyx_t_5 = NULL;\n",
       "      <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_t_2);\n",
       "      <span class='py_macro_api'>PyTuple_SET_ITEM</span>(__pyx_t_12, 0+1, __pyx_t_2);\n",
       "      __pyx_t_2 = 0;\n",
       "      __pyx_t_1 = <span class='pyx_c_api'>__Pyx_PyObject_Call</span>(__pyx_t_6, __pyx_t_12, NULL); if (unlikely(!__pyx_t_1)) __PYX_ERR(0, 26, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_1);\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_12); __pyx_t_12 = 0;\n",
       "    }\n",
       "  }\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_6); __pyx_t_6 = 0;\n",
       "  __pyx_t_13 = <span class='pyx_c_api'>__Pyx_PyObject_to_MemoryviewSlice_dc_double</span>(__pyx_t_1);\n",
       "  if (unlikely(!__pyx_t_13.memview)) __PYX_ERR(0, 26, __pyx_L1_error)\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_1); __pyx_t_1 = 0;\n",
       "  __pyx_v_comp_obj_value = __pyx_t_13;\n",
       "  __pyx_t_13.memview = NULL;\n",
       "  __pyx_t_13.data = NULL;\n",
       "/* … */\n",
       "  __pyx_tuple__4 = <span class='py_c_api'>PyTuple_Pack</span>(1, __pyx_int_1); if (unlikely(!__pyx_tuple__4)) __PYX_ERR(0, 26, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_tuple__4);\n",
       "  <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_tuple__4);\n",
       "</pre><pre class=\"cython line score-73\" onclick='toggleDiv(this)'>+<span class=\"\">27</span>:     <span class=\"k\">cdef</span> <span class=\"kt\">object</span> <span class=\"nf\">dependent_comp</span> <span class=\"o\">=</span> <span class=\"nb\">set</span><span class=\"p\">(</span><span class=\"n\">comps</span><span class=\"p\">)</span> <span class=\"o\">-</span> <span class=\"nb\">set</span><span class=\"p\">([</span><span class=\"n\">i</span><span class=\"p\">[</span><span class=\"mf\">2</span><span class=\"p\">:]</span> <span class=\"k\">for</span> <span class=\"n\">i</span> <span class=\"ow\">in</span> <span class=\"n\">cur_conds</span><span class=\"o\">.</span><span class=\"n\">keys</span><span class=\"p\">()</span> <span class=\"k\">if</span> <span class=\"n\">i</span><span class=\"o\">.</span><span class=\"n\">startswith</span><span class=\"p\">(</span><span class=\"s\">&#39;X_&#39;</span><span class=\"p\">)])</span> <span class=\"o\">-</span> <span class=\"p\">{</span><span class=\"s\">&#39;VA&#39;</span><span class=\"p\">}</span></pre>\n",
       "<pre class='cython code score-73 '>  <span class='trace'>__Pyx_TraceLine</span>(27,0,__PYX_ERR(0, 27, __pyx_L1_error))\n",
       "  __pyx_t_1 = <span class='py_c_api'>PySet_New</span>(__pyx_v_comps); if (unlikely(!__pyx_t_1)) __PYX_ERR(0, 27, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_1);\n",
       "  { /* enter inner scope */\n",
       "    PyObject *__pyx_8genexpr3__pyx_v_i = NULL;\n",
       "    __pyx_t_6 = <span class='py_c_api'>PyList_New</span>(0); if (unlikely(!__pyx_t_6)) __PYX_ERR(0, 27, __pyx_L25_error)\n",
       "    <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_6);\n",
       "    __pyx_t_3 = 0;\n",
       "    if (unlikely(__pyx_v_cur_conds == Py_None)) {\n",
       "      <span class='py_c_api'>PyErr_Format</span>(PyExc_AttributeError, \"'NoneType' object has no attribute '%s'\", \"keys\");\n",
       "      __PYX_ERR(0, 27, __pyx_L25_error)\n",
       "    }\n",
       "    __pyx_t_2 = __Pyx_dict_iterator(__pyx_v_cur_conds, 0, __pyx_n_s_keys, (&amp;__pyx_t_7), (&amp;__pyx_t_8)); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 27, __pyx_L25_error)\n",
       "    <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "    <span class='pyx_macro_api'>__Pyx_XDECREF</span>(__pyx_t_12);\n",
       "    __pyx_t_12 = __pyx_t_2;\n",
       "    __pyx_t_2 = 0;\n",
       "    while (1) {\n",
       "      __pyx_t_9 = __Pyx_dict_iter_next(__pyx_t_12, __pyx_t_7, &amp;__pyx_t_3, &amp;__pyx_t_2, NULL, NULL, __pyx_t_8);\n",
       "      if (unlikely(__pyx_t_9 == 0)) break;\n",
       "      if (unlikely(__pyx_t_9 == -1)) __PYX_ERR(0, 27, __pyx_L25_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "      <span class='pyx_macro_api'>__Pyx_XDECREF_SET</span>(__pyx_8genexpr3__pyx_v_i, __pyx_t_2);\n",
       "      __pyx_t_2 = 0;\n",
       "      __pyx_t_2 = <span class='pyx_c_api'>__Pyx_PyObject_GetAttrStr</span>(__pyx_8genexpr3__pyx_v_i, __pyx_n_s_startswith); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 27, __pyx_L25_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "      __pyx_t_5 = <span class='pyx_c_api'>__Pyx_PyObject_Call</span>(__pyx_t_2, __pyx_tuple__5, NULL); if (unlikely(!__pyx_t_5)) __PYX_ERR(0, 27, __pyx_L25_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_5);\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "      __pyx_t_10 = <span class='pyx_c_api'>__Pyx_PyObject_IsTrue</span>(__pyx_t_5); if (unlikely(__pyx_t_10 &lt; 0)) __PYX_ERR(0, 27, __pyx_L25_error)\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_5); __pyx_t_5 = 0;\n",
       "      if (__pyx_t_10) {\n",
       "/* … */\n",
       "  __pyx_tuple__5 = <span class='py_c_api'>PyTuple_Pack</span>(1, __pyx_n_u_X); if (unlikely(!__pyx_tuple__5)) __PYX_ERR(0, 27, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_tuple__5);\n",
       "  <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_tuple__5);\n",
       "        __pyx_t_5 = <span class='pyx_c_api'>__Pyx_PyObject_GetSlice</span>(__pyx_8genexpr3__pyx_v_i, 2, 0, NULL, NULL, &amp;__pyx_slice__6, 1, 0, 0); if (unlikely(!__pyx_t_5)) __PYX_ERR(0, 27, __pyx_L25_error)\n",
       "        <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_5);\n",
       "        if (unlikely(<span class='pyx_c_api'>__Pyx_ListComp_Append</span>(__pyx_t_6, (PyObject*)__pyx_t_5))) __PYX_ERR(0, 27, __pyx_L25_error)\n",
       "        <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_5); __pyx_t_5 = 0;\n",
       "      }\n",
       "    }\n",
       "    <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_12); __pyx_t_12 = 0;\n",
       "    <span class='pyx_macro_api'>__Pyx_XDECREF</span>(__pyx_8genexpr3__pyx_v_i);\n",
       "    goto __pyx_L29_exit_scope;\n",
       "    __pyx_L25_error:;\n",
       "    <span class='pyx_macro_api'>__Pyx_XDECREF</span>(__pyx_8genexpr3__pyx_v_i);\n",
       "    goto __pyx_L1_error;\n",
       "    __pyx_L29_exit_scope:;\n",
       "  } /* exit inner scope */\n",
       "  __pyx_t_12 = <span class='py_c_api'>PySet_New</span>(__pyx_t_6); if (unlikely(!__pyx_t_12)) __PYX_ERR(0, 27, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_12);\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_6); __pyx_t_6 = 0;\n",
       "  __pyx_t_6 = <span class='py_c_api'>PyNumber_Subtract</span>(__pyx_t_1, __pyx_t_12); if (unlikely(!__pyx_t_6)) __PYX_ERR(0, 27, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_6);\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_1); __pyx_t_1 = 0;\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_12); __pyx_t_12 = 0;\n",
       "  __pyx_t_12 = <span class='py_c_api'>PySet_New</span>(0); if (unlikely(!__pyx_t_12)) __PYX_ERR(0, 27, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_12);\n",
       "  if (<span class='py_c_api'>PySet_Add</span>(__pyx_t_12, __pyx_n_u_VA) &lt; 0) __PYX_ERR(0, 27, __pyx_L1_error)\n",
       "  __pyx_t_1 = <span class='py_c_api'>PyNumber_Subtract</span>(__pyx_t_6, __pyx_t_12); if (unlikely(!__pyx_t_1)) __PYX_ERR(0, 27, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_1);\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_6); __pyx_t_6 = 0;\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_12); __pyx_t_12 = 0;\n",
       "  __pyx_v_dependent_comp = __pyx_t_1;\n",
       "  __pyx_t_1 = 0;\n",
       "  __pyx_slice__6 = <span class='py_c_api'>PySlice_New</span>(__pyx_int_2, Py_None, Py_None); if (unlikely(!__pyx_slice__6)) __PYX_ERR(0, 27, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_slice__6);\n",
       "  <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_slice__6);\n",
       "</pre><pre class=\"cython line score-9\" onclick='toggleDiv(this)'>+<span class=\"\">28</span>:     <span class=\"n\">dependent_comp</span> <span class=\"o\">=</span> <span class=\"nb\">list</span><span class=\"p\">(</span><span class=\"n\">dependent_comp</span><span class=\"p\">)[</span><span class=\"mf\">0</span><span class=\"p\">]</span></pre>\n",
       "<pre class='cython code score-9 '>  <span class='trace'>__Pyx_TraceLine</span>(28,0,__PYX_ERR(0, 28, __pyx_L1_error))\n",
       "  __pyx_t_1 = <span class='py_c_api'>PySequence_List</span>(__pyx_v_dependent_comp); if (unlikely(!__pyx_t_1)) __PYX_ERR(0, 28, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_1);\n",
       "  __pyx_t_12 = <span class='py_macro_api'>PyList_GET_ITEM</span>(__pyx_t_1, 0);\n",
       "  <span class='pyx_macro_api'>__Pyx_INCREF</span>(__pyx_t_12);\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_1); __pyx_t_1 = 0;\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF_SET</span>(__pyx_v_dependent_comp, __pyx_t_12);\n",
       "  __pyx_t_12 = 0;\n",
       "</pre><pre class=\"cython line score-0\" onclick='toggleDiv(this)'>+<span class=\"\">29</span>:     <span class=\"k\">cdef</span> <span class=\"kt\">int</span> <span class=\"nf\">num_constraints</span> <span class=\"o\">=</span> <span class=\"n\">num_sitefrac_bals</span> <span class=\"o\">+</span> <span class=\"n\">num_mass_bals</span></pre>\n",
       "<pre class='cython code score-0 '>  <span class='trace'>__Pyx_TraceLine</span>(29,0,__PYX_ERR(0, 29, __pyx_L1_error))\n",
       "  __pyx_v_num_constraints = (__pyx_v_num_sitefrac_bals + __pyx_v_num_mass_bals);\n",
       "</pre><pre class=\"cython line score-5\" onclick='toggleDiv(this)'>+<span class=\"\">30</span>:     <span class=\"k\">cdef</span> <span class=\"kt\">int</span> <span class=\"nf\">num_phases</span> <span class=\"o\">=</span> <span class=\"nb\">len</span><span class=\"p\">(</span><span class=\"n\">phases</span><span class=\"p\">)</span></pre>\n",
       "<pre class='cython code score-5 '>  <span class='trace'>__Pyx_TraceLine</span>(30,0,__PYX_ERR(0, 30, __pyx_L1_error))\n",
       "  __pyx_t_7 = <span class='py_c_api'>PyObject_Length</span>(__pyx_v_phases); if (unlikely(__pyx_t_7 == -1)) __PYX_ERR(0, 30, __pyx_L1_error)\n",
       "  __pyx_v_num_phases = __pyx_t_7;\n",
       "</pre><pre class=\"cython line score-5\" onclick='toggleDiv(this)'>+<span class=\"\">31</span>:     <span class=\"k\">cdef</span> <span class=\"kt\">int</span> <span class=\"nf\">num_vars</span> <span class=\"o\">=</span> <span class=\"nb\">len</span><span class=\"p\">(</span><span class=\"n\">site_fracs</span><span class=\"p\">)</span> <span class=\"o\">+</span> <span class=\"n\">num_phases</span></pre>\n",
       "<pre class='cython code score-5 '>  <span class='trace'>__Pyx_TraceLine</span>(31,0,__PYX_ERR(0, 31, __pyx_L1_error))\n",
       "  __pyx_t_7 = <span class='py_c_api'>PyObject_Length</span>(((PyObject *)__pyx_v_site_fracs)); if (unlikely(__pyx_t_7 == -1)) __PYX_ERR(0, 31, __pyx_L1_error)\n",
       "  __pyx_v_num_vars = (__pyx_t_7 + __pyx_v_num_phases);\n",
       "</pre><pre class=\"cython line score-46\" onclick='toggleDiv(this)'>+<span class=\"\">32</span>:     <span class=\"k\">cdef</span> <span class=\"kt\">np</span>.<span class=\"kt\">ndarray</span>[<span class=\"nf\">ndim</span><span class=\"o\">=</span><span class=\"mf\">1</span><span class=\"p\">,</span> <span class=\"n\">dtype</span><span class=\"o\">=</span><span class=\"n\">np</span><span class=\"o\">.</span><span class=\"n\">int_t</span><span class=\"p\">]</span> <span class=\"n\">phase_dof</span> <span class=\"o\">=</span> <span class=\"n\">_compute_phase_dof</span><span class=\"p\">(</span><span class=\"n\">dbf</span><span class=\"p\">,</span> <span class=\"n\">comps</span><span class=\"p\">,</span> <span class=\"n\">phases</span><span class=\"p\">)</span></pre>\n",
       "<pre class='cython code score-46 '>  <span class='trace'>__Pyx_TraceLine</span>(32,0,__PYX_ERR(0, 32, __pyx_L1_error))\n",
       "  __pyx_t_1 = <span class='pyx_c_api'>__Pyx_GetModuleGlobalName</span>(__pyx_n_s_compute_phase_dof); if (unlikely(!__pyx_t_1)) __PYX_ERR(0, 32, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_1);\n",
       "  __pyx_t_6 = NULL;\n",
       "  __pyx_t_8 = 0;\n",
       "  if (CYTHON_UNPACK_METHODS &amp;&amp; unlikely(<span class='py_c_api'>PyMethod_Check</span>(__pyx_t_1))) {\n",
       "    __pyx_t_6 = <span class='py_macro_api'>PyMethod_GET_SELF</span>(__pyx_t_1);\n",
       "    if (likely(__pyx_t_6)) {\n",
       "      PyObject* function = <span class='py_macro_api'>PyMethod_GET_FUNCTION</span>(__pyx_t_1);\n",
       "      <span class='pyx_macro_api'>__Pyx_INCREF</span>(__pyx_t_6);\n",
       "      <span class='pyx_macro_api'>__Pyx_INCREF</span>(function);\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF_SET</span>(__pyx_t_1, function);\n",
       "      __pyx_t_8 = 1;\n",
       "    }\n",
       "  }\n",
       "  #if CYTHON_FAST_PYCALL\n",
       "  if (<span class='py_c_api'>PyFunction_Check</span>(__pyx_t_1)) {\n",
       "    PyObject *__pyx_temp[4] = {__pyx_t_6, __pyx_v_dbf, __pyx_v_comps, __pyx_v_phases};\n",
       "    __pyx_t_12 = <span class='pyx_c_api'>__Pyx_PyFunction_FastCall</span>(__pyx_t_1, __pyx_temp+1-__pyx_t_8, 3+__pyx_t_8); if (unlikely(!__pyx_t_12)) __PYX_ERR(0, 32, __pyx_L1_error)\n",
       "    <span class='pyx_macro_api'>__Pyx_XDECREF</span>(__pyx_t_6); __pyx_t_6 = 0;\n",
       "    <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_12);\n",
       "  } else\n",
       "  #endif\n",
       "  #if CYTHON_FAST_PYCCALL\n",
       "  if (<span class='pyx_c_api'>__Pyx_PyFastCFunction_Check</span>(__pyx_t_1)) {\n",
       "    PyObject *__pyx_temp[4] = {__pyx_t_6, __pyx_v_dbf, __pyx_v_comps, __pyx_v_phases};\n",
       "    __pyx_t_12 = <span class='pyx_c_api'>__Pyx_PyCFunction_FastCall</span>(__pyx_t_1, __pyx_temp+1-__pyx_t_8, 3+__pyx_t_8); if (unlikely(!__pyx_t_12)) __PYX_ERR(0, 32, __pyx_L1_error)\n",
       "    <span class='pyx_macro_api'>__Pyx_XDECREF</span>(__pyx_t_6); __pyx_t_6 = 0;\n",
       "    <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_12);\n",
       "  } else\n",
       "  #endif\n",
       "  {\n",
       "    __pyx_t_5 = <span class='py_c_api'>PyTuple_New</span>(3+__pyx_t_8); if (unlikely(!__pyx_t_5)) __PYX_ERR(0, 32, __pyx_L1_error)\n",
       "    <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_5);\n",
       "    if (__pyx_t_6) {\n",
       "      <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_t_6); <span class='py_macro_api'>PyTuple_SET_ITEM</span>(__pyx_t_5, 0, __pyx_t_6); __pyx_t_6 = NULL;\n",
       "    }\n",
       "    <span class='pyx_macro_api'>__Pyx_INCREF</span>(__pyx_v_dbf);\n",
       "    <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_v_dbf);\n",
       "    <span class='py_macro_api'>PyTuple_SET_ITEM</span>(__pyx_t_5, 0+__pyx_t_8, __pyx_v_dbf);\n",
       "    <span class='pyx_macro_api'>__Pyx_INCREF</span>(__pyx_v_comps);\n",
       "    <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_v_comps);\n",
       "    <span class='py_macro_api'>PyTuple_SET_ITEM</span>(__pyx_t_5, 1+__pyx_t_8, __pyx_v_comps);\n",
       "    <span class='pyx_macro_api'>__Pyx_INCREF</span>(__pyx_v_phases);\n",
       "    <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_v_phases);\n",
       "    <span class='py_macro_api'>PyTuple_SET_ITEM</span>(__pyx_t_5, 2+__pyx_t_8, __pyx_v_phases);\n",
       "    __pyx_t_12 = <span class='pyx_c_api'>__Pyx_PyObject_Call</span>(__pyx_t_1, __pyx_t_5, NULL); if (unlikely(!__pyx_t_12)) __PYX_ERR(0, 32, __pyx_L1_error)\n",
       "    <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_12);\n",
       "    <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_5); __pyx_t_5 = 0;\n",
       "  }\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_1); __pyx_t_1 = 0;\n",
       "  if (!(likely(((__pyx_t_12) == Py_None) || likely(<span class='pyx_c_api'>__Pyx_TypeTest</span>(__pyx_t_12, __pyx_ptype_5numpy_ndarray))))) __PYX_ERR(0, 32, __pyx_L1_error)\n",
       "  __pyx_t_14 = ((PyArrayObject *)__pyx_t_12);\n",
       "  {\n",
       "    __Pyx_BufFmt_StackElem __pyx_stack[1];\n",
       "    if (unlikely(<span class='pyx_c_api'>__Pyx_GetBufferAndValidate</span>(&amp;__pyx_pybuffernd_phase_dof.rcbuffer-&gt;pybuffer, (PyObject*)__pyx_t_14, &amp;__Pyx_TypeInfo_nn___pyx_t_5numpy_int_t, PyBUF_FORMAT| PyBUF_STRIDES, 1, 0, __pyx_stack) == -1)) {\n",
       "      __pyx_v_phase_dof = ((PyArrayObject *)Py_None); <span class='pyx_macro_api'>__Pyx_INCREF</span>(Py_None); __pyx_pybuffernd_phase_dof.rcbuffer-&gt;pybuffer.buf = NULL;\n",
       "      __PYX_ERR(0, 32, __pyx_L1_error)\n",
       "    } else {__pyx_pybuffernd_phase_dof.diminfo[0].strides = __pyx_pybuffernd_phase_dof.rcbuffer-&gt;pybuffer.strides[0]; __pyx_pybuffernd_phase_dof.diminfo[0].shape = __pyx_pybuffernd_phase_dof.rcbuffer-&gt;pybuffer.shape[0];\n",
       "    }\n",
       "  }\n",
       "  __pyx_t_14 = 0;\n",
       "  __pyx_v_phase_dof = ((PyArrayObject *)__pyx_t_12);\n",
       "  __pyx_t_12 = 0;\n",
       "</pre><pre class=\"cython line score-51\" onclick='toggleDiv(this)'>+<span class=\"\">33</span>:     <span class=\"k\">cdef</span> <span class=\"kt\">np</span>.<span class=\"kt\">ndarray</span>[<span class=\"nf\">ndim</span><span class=\"o\">=</span><span class=\"mf\">1</span><span class=\"p\">,</span> <span class=\"n\">dtype</span><span class=\"o\">=</span><span class=\"n\">np</span><span class=\"o\">.</span><span class=\"n\">float64_t</span><span class=\"p\">]</span> <span class=\"n\">l_constraints</span> <span class=\"o\">=</span> <span class=\"n\">np</span><span class=\"o\">.</span><span class=\"n\">zeros</span><span class=\"p\">(</span><span class=\"n\">num_constraints</span><span class=\"p\">)</span></pre>\n",
       "<pre class='cython code score-51 '>  <span class='trace'>__Pyx_TraceLine</span>(33,0,__PYX_ERR(0, 33, __pyx_L1_error))\n",
       "  __pyx_t_1 = <span class='pyx_c_api'>__Pyx_GetModuleGlobalName</span>(__pyx_n_s_np); if (unlikely(!__pyx_t_1)) __PYX_ERR(0, 33, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_1);\n",
       "  __pyx_t_5 = <span class='pyx_c_api'>__Pyx_PyObject_GetAttrStr</span>(__pyx_t_1, __pyx_n_s_zeros); if (unlikely(!__pyx_t_5)) __PYX_ERR(0, 33, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_5);\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_1); __pyx_t_1 = 0;\n",
       "  __pyx_t_1 = <span class='pyx_c_api'>__Pyx_PyInt_From_int</span>(__pyx_v_num_constraints); if (unlikely(!__pyx_t_1)) __PYX_ERR(0, 33, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_1);\n",
       "  __pyx_t_6 = NULL;\n",
       "  if (CYTHON_UNPACK_METHODS &amp;&amp; unlikely(<span class='py_c_api'>PyMethod_Check</span>(__pyx_t_5))) {\n",
       "    __pyx_t_6 = <span class='py_macro_api'>PyMethod_GET_SELF</span>(__pyx_t_5);\n",
       "    if (likely(__pyx_t_6)) {\n",
       "      PyObject* function = <span class='py_macro_api'>PyMethod_GET_FUNCTION</span>(__pyx_t_5);\n",
       "      <span class='pyx_macro_api'>__Pyx_INCREF</span>(__pyx_t_6);\n",
       "      <span class='pyx_macro_api'>__Pyx_INCREF</span>(function);\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF_SET</span>(__pyx_t_5, function);\n",
       "    }\n",
       "  }\n",
       "  if (!__pyx_t_6) {\n",
       "    __pyx_t_12 = <span class='pyx_c_api'>__Pyx_PyObject_CallOneArg</span>(__pyx_t_5, __pyx_t_1); if (unlikely(!__pyx_t_12)) __PYX_ERR(0, 33, __pyx_L1_error)\n",
       "    <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_1); __pyx_t_1 = 0;\n",
       "    <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_12);\n",
       "  } else {\n",
       "    #if CYTHON_FAST_PYCALL\n",
       "    if (<span class='py_c_api'>PyFunction_Check</span>(__pyx_t_5)) {\n",
       "      PyObject *__pyx_temp[2] = {__pyx_t_6, __pyx_t_1};\n",
       "      __pyx_t_12 = <span class='pyx_c_api'>__Pyx_PyFunction_FastCall</span>(__pyx_t_5, __pyx_temp+1-1, 1+1); if (unlikely(!__pyx_t_12)) __PYX_ERR(0, 33, __pyx_L1_error)\n",
       "      <span class='pyx_macro_api'>__Pyx_XDECREF</span>(__pyx_t_6); __pyx_t_6 = 0;\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_12);\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_1); __pyx_t_1 = 0;\n",
       "    } else\n",
       "    #endif\n",
       "    #if CYTHON_FAST_PYCCALL\n",
       "    if (<span class='pyx_c_api'>__Pyx_PyFastCFunction_Check</span>(__pyx_t_5)) {\n",
       "      PyObject *__pyx_temp[2] = {__pyx_t_6, __pyx_t_1};\n",
       "      __pyx_t_12 = <span class='pyx_c_api'>__Pyx_PyCFunction_FastCall</span>(__pyx_t_5, __pyx_temp+1-1, 1+1); if (unlikely(!__pyx_t_12)) __PYX_ERR(0, 33, __pyx_L1_error)\n",
       "      <span class='pyx_macro_api'>__Pyx_XDECREF</span>(__pyx_t_6); __pyx_t_6 = 0;\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_12);\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_1); __pyx_t_1 = 0;\n",
       "    } else\n",
       "    #endif\n",
       "    {\n",
       "      __pyx_t_2 = <span class='py_c_api'>PyTuple_New</span>(1+1); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 33, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "      <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_t_6); <span class='py_macro_api'>PyTuple_SET_ITEM</span>(__pyx_t_2, 0, __pyx_t_6); __pyx_t_6 = NULL;\n",
       "      <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_t_1);\n",
       "      <span class='py_macro_api'>PyTuple_SET_ITEM</span>(__pyx_t_2, 0+1, __pyx_t_1);\n",
       "      __pyx_t_1 = 0;\n",
       "      __pyx_t_12 = <span class='pyx_c_api'>__Pyx_PyObject_Call</span>(__pyx_t_5, __pyx_t_2, NULL); if (unlikely(!__pyx_t_12)) __PYX_ERR(0, 33, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_12);\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "    }\n",
       "  }\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_5); __pyx_t_5 = 0;\n",
       "  if (!(likely(((__pyx_t_12) == Py_None) || likely(<span class='pyx_c_api'>__Pyx_TypeTest</span>(__pyx_t_12, __pyx_ptype_5numpy_ndarray))))) __PYX_ERR(0, 33, __pyx_L1_error)\n",
       "  __pyx_t_15 = ((PyArrayObject *)__pyx_t_12);\n",
       "  {\n",
       "    __Pyx_BufFmt_StackElem __pyx_stack[1];\n",
       "    if (unlikely(<span class='pyx_c_api'>__Pyx_GetBufferAndValidate</span>(&amp;__pyx_pybuffernd_l_constraints.rcbuffer-&gt;pybuffer, (PyObject*)__pyx_t_15, &amp;__Pyx_TypeInfo_nn___pyx_t_5numpy_float64_t, PyBUF_FORMAT| PyBUF_STRIDES| PyBUF_WRITABLE, 1, 0, __pyx_stack) == -1)) {\n",
       "      __pyx_v_l_constraints = ((PyArrayObject *)Py_None); <span class='pyx_macro_api'>__Pyx_INCREF</span>(Py_None); __pyx_pybuffernd_l_constraints.rcbuffer-&gt;pybuffer.buf = NULL;\n",
       "      __PYX_ERR(0, 33, __pyx_L1_error)\n",
       "    } else {__pyx_pybuffernd_l_constraints.diminfo[0].strides = __pyx_pybuffernd_l_constraints.rcbuffer-&gt;pybuffer.strides[0]; __pyx_pybuffernd_l_constraints.diminfo[0].shape = __pyx_pybuffernd_l_constraints.rcbuffer-&gt;pybuffer.shape[0];\n",
       "    }\n",
       "  }\n",
       "  __pyx_t_15 = 0;\n",
       "  __pyx_v_l_constraints = ((PyArrayObject *)__pyx_t_12);\n",
       "  __pyx_t_12 = 0;\n",
       "</pre><pre class=\"cython line score-60\" onclick='toggleDiv(this)'>+<span class=\"\">34</span>:     <span class=\"k\">cdef</span> <span class=\"kt\">np</span>.<span class=\"kt\">ndarray</span>[<span class=\"nf\">ndim</span><span class=\"o\">=</span><span class=\"mf\">2</span><span class=\"p\">,</span> <span class=\"n\">dtype</span><span class=\"o\">=</span><span class=\"n\">np</span><span class=\"o\">.</span><span class=\"n\">float64_t</span><span class=\"p\">]</span> <span class=\"n\">constraint_jac</span> <span class=\"o\">=</span> <span class=\"n\">np</span><span class=\"o\">.</span><span class=\"n\">zeros</span><span class=\"p\">((</span><span class=\"n\">num_constraints</span><span class=\"p\">,</span> <span class=\"n\">num_vars</span><span class=\"p\">))</span></pre>\n",
       "<pre class='cython code score-60 '>  <span class='trace'>__Pyx_TraceLine</span>(34,0,__PYX_ERR(0, 34, __pyx_L1_error))\n",
       "  __pyx_t_5 = <span class='pyx_c_api'>__Pyx_GetModuleGlobalName</span>(__pyx_n_s_np); if (unlikely(!__pyx_t_5)) __PYX_ERR(0, 34, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_5);\n",
       "  __pyx_t_2 = <span class='pyx_c_api'>__Pyx_PyObject_GetAttrStr</span>(__pyx_t_5, __pyx_n_s_zeros); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 34, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_5); __pyx_t_5 = 0;\n",
       "  __pyx_t_5 = <span class='pyx_c_api'>__Pyx_PyInt_From_int</span>(__pyx_v_num_constraints); if (unlikely(!__pyx_t_5)) __PYX_ERR(0, 34, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_5);\n",
       "  __pyx_t_1 = <span class='pyx_c_api'>__Pyx_PyInt_From_int</span>(__pyx_v_num_vars); if (unlikely(!__pyx_t_1)) __PYX_ERR(0, 34, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_1);\n",
       "  __pyx_t_6 = <span class='py_c_api'>PyTuple_New</span>(2); if (unlikely(!__pyx_t_6)) __PYX_ERR(0, 34, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_6);\n",
       "  <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_t_5);\n",
       "  <span class='py_macro_api'>PyTuple_SET_ITEM</span>(__pyx_t_6, 0, __pyx_t_5);\n",
       "  <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_t_1);\n",
       "  <span class='py_macro_api'>PyTuple_SET_ITEM</span>(__pyx_t_6, 1, __pyx_t_1);\n",
       "  __pyx_t_5 = 0;\n",
       "  __pyx_t_1 = 0;\n",
       "  __pyx_t_1 = NULL;\n",
       "  if (CYTHON_UNPACK_METHODS &amp;&amp; unlikely(<span class='py_c_api'>PyMethod_Check</span>(__pyx_t_2))) {\n",
       "    __pyx_t_1 = <span class='py_macro_api'>PyMethod_GET_SELF</span>(__pyx_t_2);\n",
       "    if (likely(__pyx_t_1)) {\n",
       "      PyObject* function = <span class='py_macro_api'>PyMethod_GET_FUNCTION</span>(__pyx_t_2);\n",
       "      <span class='pyx_macro_api'>__Pyx_INCREF</span>(__pyx_t_1);\n",
       "      <span class='pyx_macro_api'>__Pyx_INCREF</span>(function);\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF_SET</span>(__pyx_t_2, function);\n",
       "    }\n",
       "  }\n",
       "  if (!__pyx_t_1) {\n",
       "    __pyx_t_12 = <span class='pyx_c_api'>__Pyx_PyObject_CallOneArg</span>(__pyx_t_2, __pyx_t_6); if (unlikely(!__pyx_t_12)) __PYX_ERR(0, 34, __pyx_L1_error)\n",
       "    <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_6); __pyx_t_6 = 0;\n",
       "    <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_12);\n",
       "  } else {\n",
       "    #if CYTHON_FAST_PYCALL\n",
       "    if (<span class='py_c_api'>PyFunction_Check</span>(__pyx_t_2)) {\n",
       "      PyObject *__pyx_temp[2] = {__pyx_t_1, __pyx_t_6};\n",
       "      __pyx_t_12 = <span class='pyx_c_api'>__Pyx_PyFunction_FastCall</span>(__pyx_t_2, __pyx_temp+1-1, 1+1); if (unlikely(!__pyx_t_12)) __PYX_ERR(0, 34, __pyx_L1_error)\n",
       "      <span class='pyx_macro_api'>__Pyx_XDECREF</span>(__pyx_t_1); __pyx_t_1 = 0;\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_12);\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_6); __pyx_t_6 = 0;\n",
       "    } else\n",
       "    #endif\n",
       "    #if CYTHON_FAST_PYCCALL\n",
       "    if (<span class='pyx_c_api'>__Pyx_PyFastCFunction_Check</span>(__pyx_t_2)) {\n",
       "      PyObject *__pyx_temp[2] = {__pyx_t_1, __pyx_t_6};\n",
       "      __pyx_t_12 = <span class='pyx_c_api'>__Pyx_PyCFunction_FastCall</span>(__pyx_t_2, __pyx_temp+1-1, 1+1); if (unlikely(!__pyx_t_12)) __PYX_ERR(0, 34, __pyx_L1_error)\n",
       "      <span class='pyx_macro_api'>__Pyx_XDECREF</span>(__pyx_t_1); __pyx_t_1 = 0;\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_12);\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_6); __pyx_t_6 = 0;\n",
       "    } else\n",
       "    #endif\n",
       "    {\n",
       "      __pyx_t_5 = <span class='py_c_api'>PyTuple_New</span>(1+1); if (unlikely(!__pyx_t_5)) __PYX_ERR(0, 34, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_5);\n",
       "      <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_t_1); <span class='py_macro_api'>PyTuple_SET_ITEM</span>(__pyx_t_5, 0, __pyx_t_1); __pyx_t_1 = NULL;\n",
       "      <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_t_6);\n",
       "      <span class='py_macro_api'>PyTuple_SET_ITEM</span>(__pyx_t_5, 0+1, __pyx_t_6);\n",
       "      __pyx_t_6 = 0;\n",
       "      __pyx_t_12 = <span class='pyx_c_api'>__Pyx_PyObject_Call</span>(__pyx_t_2, __pyx_t_5, NULL); if (unlikely(!__pyx_t_12)) __PYX_ERR(0, 34, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_12);\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_5); __pyx_t_5 = 0;\n",
       "    }\n",
       "  }\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "  if (!(likely(((__pyx_t_12) == Py_None) || likely(<span class='pyx_c_api'>__Pyx_TypeTest</span>(__pyx_t_12, __pyx_ptype_5numpy_ndarray))))) __PYX_ERR(0, 34, __pyx_L1_error)\n",
       "  __pyx_t_16 = ((PyArrayObject *)__pyx_t_12);\n",
       "  {\n",
       "    __Pyx_BufFmt_StackElem __pyx_stack[1];\n",
       "    if (unlikely(<span class='pyx_c_api'>__Pyx_GetBufferAndValidate</span>(&amp;__pyx_pybuffernd_constraint_jac.rcbuffer-&gt;pybuffer, (PyObject*)__pyx_t_16, &amp;__Pyx_TypeInfo_nn___pyx_t_5numpy_float64_t, PyBUF_FORMAT| PyBUF_STRIDES| PyBUF_WRITABLE, 2, 0, __pyx_stack) == -1)) {\n",
       "      __pyx_v_constraint_jac = ((PyArrayObject *)Py_None); <span class='pyx_macro_api'>__Pyx_INCREF</span>(Py_None); __pyx_pybuffernd_constraint_jac.rcbuffer-&gt;pybuffer.buf = NULL;\n",
       "      __PYX_ERR(0, 34, __pyx_L1_error)\n",
       "    } else {__pyx_pybuffernd_constraint_jac.diminfo[0].strides = __pyx_pybuffernd_constraint_jac.rcbuffer-&gt;pybuffer.strides[0]; __pyx_pybuffernd_constraint_jac.diminfo[0].shape = __pyx_pybuffernd_constraint_jac.rcbuffer-&gt;pybuffer.shape[0]; __pyx_pybuffernd_constraint_jac.diminfo[1].strides = __pyx_pybuffernd_constraint_jac.rcbuffer-&gt;pybuffer.strides[1]; __pyx_pybuffernd_constraint_jac.diminfo[1].shape = __pyx_pybuffernd_constraint_jac.rcbuffer-&gt;pybuffer.shape[1];\n",
       "    }\n",
       "  }\n",
       "  __pyx_t_16 = 0;\n",
       "  __pyx_v_constraint_jac = ((PyArrayObject *)__pyx_t_12);\n",
       "  __pyx_t_12 = 0;\n",
       "</pre><pre class=\"cython line score-45\" onclick='toggleDiv(this)'>+<span class=\"\">35</span>:     <span class=\"k\">cdef</span> <span class=\"kt\">np</span>.<span class=\"kt\">ndarray</span>[<span class=\"nf\">ndim</span><span class=\"o\">=</span><span class=\"mf\">3</span><span class=\"p\">,</span> <span class=\"n\">dtype</span><span class=\"o\">=</span><span class=\"n\">np</span><span class=\"o\">.</span><span class=\"n\">float64_t</span><span class=\"p\">]</span> <span class=\"n\">constraint_hess</span> <span class=\"o\">=</span> <span class=\"n\">np</span><span class=\"o\">.</span><span class=\"n\">zeros</span><span class=\"p\">((</span><span class=\"n\">num_constraints</span><span class=\"p\">,</span> <span class=\"n\">num_vars</span><span class=\"p\">,</span> <span class=\"n\">num_vars</span><span class=\"p\">),</span> <span class=\"n\">order</span><span class=\"o\">=</span><span class=\"s\">&#39;F&#39;</span><span class=\"p\">)</span></pre>\n",
       "<pre class='cython code score-45 '>  <span class='trace'>__Pyx_TraceLine</span>(35,0,__PYX_ERR(0, 35, __pyx_L1_error))\n",
       "  __pyx_t_12 = <span class='pyx_c_api'>__Pyx_GetModuleGlobalName</span>(__pyx_n_s_np); if (unlikely(!__pyx_t_12)) __PYX_ERR(0, 35, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_12);\n",
       "  __pyx_t_2 = <span class='pyx_c_api'>__Pyx_PyObject_GetAttrStr</span>(__pyx_t_12, __pyx_n_s_zeros); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 35, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_12); __pyx_t_12 = 0;\n",
       "  __pyx_t_12 = <span class='pyx_c_api'>__Pyx_PyInt_From_int</span>(__pyx_v_num_constraints); if (unlikely(!__pyx_t_12)) __PYX_ERR(0, 35, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_12);\n",
       "  __pyx_t_5 = <span class='pyx_c_api'>__Pyx_PyInt_From_int</span>(__pyx_v_num_vars); if (unlikely(!__pyx_t_5)) __PYX_ERR(0, 35, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_5);\n",
       "  __pyx_t_6 = <span class='pyx_c_api'>__Pyx_PyInt_From_int</span>(__pyx_v_num_vars); if (unlikely(!__pyx_t_6)) __PYX_ERR(0, 35, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_6);\n",
       "  __pyx_t_1 = <span class='py_c_api'>PyTuple_New</span>(3); if (unlikely(!__pyx_t_1)) __PYX_ERR(0, 35, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_1);\n",
       "  <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_t_12);\n",
       "  <span class='py_macro_api'>PyTuple_SET_ITEM</span>(__pyx_t_1, 0, __pyx_t_12);\n",
       "  <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_t_5);\n",
       "  <span class='py_macro_api'>PyTuple_SET_ITEM</span>(__pyx_t_1, 1, __pyx_t_5);\n",
       "  <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_t_6);\n",
       "  <span class='py_macro_api'>PyTuple_SET_ITEM</span>(__pyx_t_1, 2, __pyx_t_6);\n",
       "  __pyx_t_12 = 0;\n",
       "  __pyx_t_5 = 0;\n",
       "  __pyx_t_6 = 0;\n",
       "  __pyx_t_6 = <span class='py_c_api'>PyTuple_New</span>(1); if (unlikely(!__pyx_t_6)) __PYX_ERR(0, 35, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_6);\n",
       "  <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_t_1);\n",
       "  <span class='py_macro_api'>PyTuple_SET_ITEM</span>(__pyx_t_6, 0, __pyx_t_1);\n",
       "  __pyx_t_1 = 0;\n",
       "  __pyx_t_1 = <span class='py_c_api'>PyDict_New</span>(); if (unlikely(!__pyx_t_1)) __PYX_ERR(0, 35, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_1);\n",
       "  if (<span class='py_c_api'>PyDict_SetItem</span>(__pyx_t_1, __pyx_n_s_order, __pyx_n_u_F) &lt; 0) __PYX_ERR(0, 35, __pyx_L1_error)\n",
       "  __pyx_t_5 = <span class='pyx_c_api'>__Pyx_PyObject_Call</span>(__pyx_t_2, __pyx_t_6, __pyx_t_1); if (unlikely(!__pyx_t_5)) __PYX_ERR(0, 35, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_5);\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_6); __pyx_t_6 = 0;\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_1); __pyx_t_1 = 0;\n",
       "  if (!(likely(((__pyx_t_5) == Py_None) || likely(<span class='pyx_c_api'>__Pyx_TypeTest</span>(__pyx_t_5, __pyx_ptype_5numpy_ndarray))))) __PYX_ERR(0, 35, __pyx_L1_error)\n",
       "  __pyx_t_17 = ((PyArrayObject *)__pyx_t_5);\n",
       "  {\n",
       "    __Pyx_BufFmt_StackElem __pyx_stack[1];\n",
       "    if (unlikely(<span class='pyx_c_api'>__Pyx_GetBufferAndValidate</span>(&amp;__pyx_pybuffernd_constraint_hess.rcbuffer-&gt;pybuffer, (PyObject*)__pyx_t_17, &amp;__Pyx_TypeInfo_nn___pyx_t_5numpy_float64_t, PyBUF_FORMAT| PyBUF_STRIDES| PyBUF_WRITABLE, 3, 0, __pyx_stack) == -1)) {\n",
       "      __pyx_v_constraint_hess = ((PyArrayObject *)Py_None); <span class='pyx_macro_api'>__Pyx_INCREF</span>(Py_None); __pyx_pybuffernd_constraint_hess.rcbuffer-&gt;pybuffer.buf = NULL;\n",
       "      __PYX_ERR(0, 35, __pyx_L1_error)\n",
       "    } else {__pyx_pybuffernd_constraint_hess.diminfo[0].strides = __pyx_pybuffernd_constraint_hess.rcbuffer-&gt;pybuffer.strides[0]; __pyx_pybuffernd_constraint_hess.diminfo[0].shape = __pyx_pybuffernd_constraint_hess.rcbuffer-&gt;pybuffer.shape[0]; __pyx_pybuffernd_constraint_hess.diminfo[1].strides = __pyx_pybuffernd_constraint_hess.rcbuffer-&gt;pybuffer.strides[1]; __pyx_pybuffernd_constraint_hess.diminfo[1].shape = __pyx_pybuffernd_constraint_hess.rcbuffer-&gt;pybuffer.shape[1]; __pyx_pybuffernd_constraint_hess.diminfo[2].strides = __pyx_pybuffernd_constraint_hess.rcbuffer-&gt;pybuffer.strides[2]; __pyx_pybuffernd_constraint_hess.diminfo[2].shape = __pyx_pybuffernd_constraint_hess.rcbuffer-&gt;pybuffer.shape[2];\n",
       "    }\n",
       "  }\n",
       "  __pyx_t_17 = 0;\n",
       "  __pyx_v_constraint_hess = ((PyArrayObject *)__pyx_t_5);\n",
       "  __pyx_t_5 = 0;\n",
       "</pre><pre class=\"cython line score-0\">&#xA0;<span class=\"\">36</span>:     <span class=\"k\">cdef</span> <span class=\"kt\">double</span>[<span class=\"p\">::</span><span class=\"mf\">1</span><span class=\"p\">]</span> <span class=\"n\">sfview</span></pre>\n",
       "<pre class=\"cython line score-0\">&#xA0;<span class=\"\">37</span>:     <span class=\"k\">cdef</span> <span class=\"kt\">double</span>[<span class=\"p\">::</span><span class=\"mf\">1</span><span class=\"p\">]</span> <span class=\"n\">comp_grad_value</span></pre>\n",
       "<pre class=\"cython line score-0\">&#xA0;<span class=\"\">38</span>:     <span class=\"k\">cdef</span> <span class=\"kt\">double</span>[<span class=\"p\">::</span><span class=\"mf\">1</span><span class=\"p\">,:]</span> <span class=\"n\">comp_hess_value</span></pre>\n",
       "<pre class=\"cython line score-0\">&#xA0;<span class=\"\">39</span>:     <span class=\"k\">cdef</span> <span class=\"kt\">int</span> <span class=\"nf\">phase_idx</span><span class=\"p\">,</span> <span class=\"nf\">var_offset</span><span class=\"p\">,</span> <span class=\"nf\">constraint_offset</span><span class=\"p\">,</span> <span class=\"nf\">var_idx</span><span class=\"p\">,</span> <span class=\"nf\">iter_idx</span><span class=\"p\">,</span> <span class=\"nf\">grad_idx</span><span class=\"p\">,</span> <span class=\"nf\">hess_idx</span><span class=\"p\">,</span> <span class=\"nf\">comp_idx</span><span class=\"p\">,</span> <span class=\"nf\">idx</span></pre>\n",
       "<pre class=\"cython line score-0\">&#xA0;<span class=\"\">40</span>:     <span class=\"k\">cdef</span> <span class=\"kt\">PhaseRecord</span> <span class=\"nf\">prn</span></pre>\n",
       "<pre class=\"cython line score-0\">&#xA0;<span class=\"\">41</span>:     <span class=\"k\">cdef</span> <span class=\"kt\">double</span> <span class=\"nf\">phase_frac</span></pre>\n",
       "<pre class=\"cython line score-0\">&#xA0;<span class=\"\">42</span>: </pre>\n",
       "<pre class=\"cython line score-0\">&#xA0;<span class=\"\">43</span>:     <span class=\"c\"># Ordering of constraints by row: sitefrac bal of each phase, then component mass balance</span></pre>\n",
       "<pre class=\"cython line score-0\">&#xA0;<span class=\"\">44</span>:     <span class=\"c\"># Ordering of constraints by column: site fractions of each phase, then phase fractions</span></pre>\n",
       "<pre class=\"cython line score-0\">&#xA0;<span class=\"\">45</span>:     <span class=\"c\"># First: Site fraction balance constraints</span></pre>\n",
       "<pre class=\"cython line score-0\" onclick='toggleDiv(this)'>+<span class=\"\">46</span>:     <span class=\"n\">var_idx</span> <span class=\"o\">=</span> <span class=\"mf\">0</span></pre>\n",
       "<pre class='cython code score-0 '>  <span class='trace'>__Pyx_TraceLine</span>(46,0,__PYX_ERR(0, 46, __pyx_L1_error))\n",
       "  __pyx_v_var_idx = 0;\n",
       "</pre><pre class=\"cython line score-0\" onclick='toggleDiv(this)'>+<span class=\"\">47</span>:     <span class=\"n\">constraint_offset</span> <span class=\"o\">=</span> <span class=\"mf\">0</span></pre>\n",
       "<pre class='cython code score-0 '>  <span class='trace'>__Pyx_TraceLine</span>(47,0,__PYX_ERR(0, 47, __pyx_L1_error))\n",
       "  __pyx_v_constraint_offset = 0;\n",
       "</pre><pre class=\"cython line score-0\" onclick='toggleDiv(this)'>+<span class=\"\">48</span>:     <span class=\"k\">for</span> <span class=\"n\">phase_idx</span> <span class=\"ow\">in</span> <span class=\"nb\">range</span><span class=\"p\">(</span><span class=\"n\">num_phases</span><span class=\"p\">):</span></pre>\n",
       "<pre class='cython code score-0 '>  <span class='trace'>__Pyx_TraceLine</span>(48,0,__PYX_ERR(0, 48, __pyx_L1_error))\n",
       "  __pyx_t_8 = __pyx_v_num_phases;\n",
       "  for (__pyx_t_9 = 0; __pyx_t_9 &lt; __pyx_t_8; __pyx_t_9+=1) {\n",
       "    __pyx_v_phase_idx = __pyx_t_9;\n",
       "</pre><pre class=\"cython line score-3\" onclick='toggleDiv(this)'>+<span class=\"\">49</span>:         <span class=\"n\">name</span> <span class=\"o\">=</span> <span class=\"n\">phases</span><span class=\"p\">[</span><span class=\"n\">phase_idx</span><span class=\"p\">]</span></pre>\n",
       "<pre class='cython code score-3 '>    <span class='trace'>__Pyx_TraceLine</span>(49,0,__PYX_ERR(0, 49, __pyx_L1_error))\n",
       "    __pyx_t_5 = <span class='pyx_c_api'>__Pyx_GetItemInt</span>(__pyx_v_phases, __pyx_v_phase_idx, int, 1, __Pyx_PyInt_From_int, 0, 0, 0); if (unlikely(!__pyx_t_5)) __PYX_ERR(0, 49, __pyx_L1_error)\n",
       "    <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_5);\n",
       "    <span class='pyx_macro_api'>__Pyx_XDECREF_SET</span>(__pyx_v_name, __pyx_t_5);\n",
       "    __pyx_t_5 = 0;\n",
       "</pre><pre class=\"cython line score-17\" onclick='toggleDiv(this)'>+<span class=\"\">50</span>:         <span class=\"k\">for</span> <span class=\"n\">idx</span> <span class=\"ow\">in</span> <span class=\"nb\">range</span><span class=\"p\">(</span><span class=\"nb\">len</span><span class=\"p\">(</span><span class=\"n\">dbf</span><span class=\"o\">.</span><span class=\"n\">phases</span><span class=\"p\">[</span><span class=\"n\">name</span><span class=\"p\">]</span><span class=\"o\">.</span><span class=\"n\">sublattices</span><span class=\"p\">)):</span></pre>\n",
       "<pre class='cython code score-17 '>    <span class='trace'>__Pyx_TraceLine</span>(50,0,__PYX_ERR(0, 50, __pyx_L1_error))\n",
       "    __pyx_t_5 = <span class='pyx_c_api'>__Pyx_PyObject_GetAttrStr</span>(__pyx_v_dbf, __pyx_n_s_phases); if (unlikely(!__pyx_t_5)) __PYX_ERR(0, 50, __pyx_L1_error)\n",
       "    <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_5);\n",
       "    __pyx_t_1 = <span class='py_c_api'>PyObject_GetItem</span>(__pyx_t_5, __pyx_v_name); if (unlikely(!__pyx_t_1)) __PYX_ERR(0, 50, __pyx_L1_error)\n",
       "    <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_1);\n",
       "    <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_5); __pyx_t_5 = 0;\n",
       "    __pyx_t_5 = <span class='pyx_c_api'>__Pyx_PyObject_GetAttrStr</span>(__pyx_t_1, __pyx_n_s_sublattices); if (unlikely(!__pyx_t_5)) __PYX_ERR(0, 50, __pyx_L1_error)\n",
       "    <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_5);\n",
       "    <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_1); __pyx_t_1 = 0;\n",
       "    __pyx_t_7 = <span class='py_c_api'>PyObject_Length</span>(__pyx_t_5); if (unlikely(__pyx_t_7 == -1)) __PYX_ERR(0, 50, __pyx_L1_error)\n",
       "    <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_5); __pyx_t_5 = 0;\n",
       "    for (__pyx_t_18 = 0; __pyx_t_18 &lt; __pyx_t_7; __pyx_t_18+=1) {\n",
       "      __pyx_v_idx = __pyx_t_18;\n",
       "</pre><pre class=\"cython line score-61\" onclick='toggleDiv(this)'>+<span class=\"\">51</span>:             <span class=\"n\">active_in_subl</span> <span class=\"o\">=</span> <span class=\"nb\">set</span><span class=\"p\">(</span><span class=\"n\">dbf</span><span class=\"o\">.</span><span class=\"n\">phases</span><span class=\"p\">[</span><span class=\"n\">name</span><span class=\"p\">]</span><span class=\"o\">.</span><span class=\"n\">constituents</span><span class=\"p\">[</span><span class=\"n\">idx</span><span class=\"p\">])</span><span class=\"o\">.</span><span class=\"n\">intersection</span><span class=\"p\">(</span><span class=\"n\">comps</span><span class=\"p\">)</span></pre>\n",
       "<pre class='cython code score-61 '>      <span class='trace'>__Pyx_TraceLine</span>(51,0,__PYX_ERR(0, 51, __pyx_L1_error))\n",
       "      __pyx_t_1 = <span class='pyx_c_api'>__Pyx_PyObject_GetAttrStr</span>(__pyx_v_dbf, __pyx_n_s_phases); if (unlikely(!__pyx_t_1)) __PYX_ERR(0, 51, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_1);\n",
       "      __pyx_t_6 = <span class='py_c_api'>PyObject_GetItem</span>(__pyx_t_1, __pyx_v_name); if (unlikely(!__pyx_t_6)) __PYX_ERR(0, 51, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_6);\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_1); __pyx_t_1 = 0;\n",
       "      __pyx_t_1 = <span class='pyx_c_api'>__Pyx_PyObject_GetAttrStr</span>(__pyx_t_6, __pyx_n_s_constituents); if (unlikely(!__pyx_t_1)) __PYX_ERR(0, 51, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_1);\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_6); __pyx_t_6 = 0;\n",
       "      __pyx_t_6 = <span class='pyx_c_api'>__Pyx_GetItemInt</span>(__pyx_t_1, __pyx_v_idx, int, 1, __Pyx_PyInt_From_int, 0, 0, 0); if (unlikely(!__pyx_t_6)) __PYX_ERR(0, 51, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_6);\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_1); __pyx_t_1 = 0;\n",
       "      __pyx_t_1 = <span class='py_c_api'>PySet_New</span>(__pyx_t_6); if (unlikely(!__pyx_t_1)) __PYX_ERR(0, 51, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_1);\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_6); __pyx_t_6 = 0;\n",
       "      __pyx_t_6 = <span class='pyx_c_api'>__Pyx_PyObject_GetAttrStr</span>(__pyx_t_1, __pyx_n_s_intersection); if (unlikely(!__pyx_t_6)) __PYX_ERR(0, 51, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_6);\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_1); __pyx_t_1 = 0;\n",
       "      __pyx_t_1 = NULL;\n",
       "      if (CYTHON_UNPACK_METHODS &amp;&amp; likely(<span class='py_c_api'>PyMethod_Check</span>(__pyx_t_6))) {\n",
       "        __pyx_t_1 = <span class='py_macro_api'>PyMethod_GET_SELF</span>(__pyx_t_6);\n",
       "        if (likely(__pyx_t_1)) {\n",
       "          PyObject* function = <span class='py_macro_api'>PyMethod_GET_FUNCTION</span>(__pyx_t_6);\n",
       "          <span class='pyx_macro_api'>__Pyx_INCREF</span>(__pyx_t_1);\n",
       "          <span class='pyx_macro_api'>__Pyx_INCREF</span>(function);\n",
       "          <span class='pyx_macro_api'>__Pyx_DECREF_SET</span>(__pyx_t_6, function);\n",
       "        }\n",
       "      }\n",
       "      if (!__pyx_t_1) {\n",
       "        __pyx_t_5 = <span class='pyx_c_api'>__Pyx_PyObject_CallOneArg</span>(__pyx_t_6, __pyx_v_comps); if (unlikely(!__pyx_t_5)) __PYX_ERR(0, 51, __pyx_L1_error)\n",
       "        <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_5);\n",
       "      } else {\n",
       "        #if CYTHON_FAST_PYCALL\n",
       "        if (<span class='py_c_api'>PyFunction_Check</span>(__pyx_t_6)) {\n",
       "          PyObject *__pyx_temp[2] = {__pyx_t_1, __pyx_v_comps};\n",
       "          __pyx_t_5 = <span class='pyx_c_api'>__Pyx_PyFunction_FastCall</span>(__pyx_t_6, __pyx_temp+1-1, 1+1); if (unlikely(!__pyx_t_5)) __PYX_ERR(0, 51, __pyx_L1_error)\n",
       "          <span class='pyx_macro_api'>__Pyx_XDECREF</span>(__pyx_t_1); __pyx_t_1 = 0;\n",
       "          <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_5);\n",
       "        } else\n",
       "        #endif\n",
       "        #if CYTHON_FAST_PYCCALL\n",
       "        if (<span class='pyx_c_api'>__Pyx_PyFastCFunction_Check</span>(__pyx_t_6)) {\n",
       "          PyObject *__pyx_temp[2] = {__pyx_t_1, __pyx_v_comps};\n",
       "          __pyx_t_5 = <span class='pyx_c_api'>__Pyx_PyCFunction_FastCall</span>(__pyx_t_6, __pyx_temp+1-1, 1+1); if (unlikely(!__pyx_t_5)) __PYX_ERR(0, 51, __pyx_L1_error)\n",
       "          <span class='pyx_macro_api'>__Pyx_XDECREF</span>(__pyx_t_1); __pyx_t_1 = 0;\n",
       "          <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_5);\n",
       "        } else\n",
       "        #endif\n",
       "        {\n",
       "          __pyx_t_2 = <span class='py_c_api'>PyTuple_New</span>(1+1); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 51, __pyx_L1_error)\n",
       "          <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "          <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_t_1); <span class='py_macro_api'>PyTuple_SET_ITEM</span>(__pyx_t_2, 0, __pyx_t_1); __pyx_t_1 = NULL;\n",
       "          <span class='pyx_macro_api'>__Pyx_INCREF</span>(__pyx_v_comps);\n",
       "          <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_v_comps);\n",
       "          <span class='py_macro_api'>PyTuple_SET_ITEM</span>(__pyx_t_2, 0+1, __pyx_v_comps);\n",
       "          __pyx_t_5 = <span class='pyx_c_api'>__Pyx_PyObject_Call</span>(__pyx_t_6, __pyx_t_2, NULL); if (unlikely(!__pyx_t_5)) __PYX_ERR(0, 51, __pyx_L1_error)\n",
       "          <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_5);\n",
       "          <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "        }\n",
       "      }\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_6); __pyx_t_6 = 0;\n",
       "      <span class='pyx_macro_api'>__Pyx_XDECREF_SET</span>(__pyx_v_active_in_subl, __pyx_t_5);\n",
       "      __pyx_t_5 = 0;\n",
       "</pre><pre class=\"cython line score-11\" onclick='toggleDiv(this)'>+<span class=\"\">52</span>:             <span class=\"n\">ais_len</span> <span class=\"o\">=</span> <span class=\"nb\">len</span><span class=\"p\">(</span><span class=\"n\">active_in_subl</span><span class=\"p\">)</span></pre>\n",
       "<pre class='cython code score-11 '>      <span class='trace'>__Pyx_TraceLine</span>(52,0,__PYX_ERR(0, 52, __pyx_L1_error))\n",
       "      __pyx_t_3 = <span class='py_c_api'>PyObject_Length</span>(__pyx_v_active_in_subl); if (unlikely(__pyx_t_3 == -1)) __PYX_ERR(0, 52, __pyx_L1_error)\n",
       "      __pyx_t_5 = <span class='py_c_api'>PyInt_FromSsize_t</span>(__pyx_t_3); if (unlikely(!__pyx_t_5)) __PYX_ERR(0, 52, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_5);\n",
       "      <span class='pyx_macro_api'>__Pyx_XDECREF_SET</span>(__pyx_v_ais_len, __pyx_t_5);\n",
       "      __pyx_t_5 = 0;\n",
       "</pre><pre class=\"cython line score-22\" onclick='toggleDiv(this)'>+<span class=\"\">53</span>:             <span class=\"n\">constraint_jac</span><span class=\"p\">[</span><span class=\"n\">constraint_offset</span> <span class=\"o\">+</span> <span class=\"n\">idx</span><span class=\"p\">,</span></pre>\n",
       "<pre class='cython code score-22 '>      <span class='trace'>__Pyx_TraceLine</span>(53,0,__PYX_ERR(0, 53, __pyx_L1_error))\n",
       "      __pyx_t_5 = <span class='pyx_c_api'>__Pyx_PyInt_From_int</span>((__pyx_v_constraint_offset + __pyx_v_idx)); if (unlikely(!__pyx_t_5)) __PYX_ERR(0, 53, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_5);\n",
       "/* … */\n",
       "      <span class='trace'>__Pyx_TraceLine</span>(53,0,__PYX_ERR(0, 53, __pyx_L1_error))\n",
       "      __pyx_t_2 = <span class='py_c_api'>PySlice_New</span>(__pyx_t_6, __pyx_t_1, Py_None); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 53, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_6); __pyx_t_6 = 0;\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_1); __pyx_t_1 = 0;\n",
       "      __pyx_t_1 = <span class='py_c_api'>PyTuple_New</span>(2); if (unlikely(!__pyx_t_1)) __PYX_ERR(0, 53, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_1);\n",
       "      <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_t_5);\n",
       "      <span class='py_macro_api'>PyTuple_SET_ITEM</span>(__pyx_t_1, 0, __pyx_t_5);\n",
       "      <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_t_2);\n",
       "      <span class='py_macro_api'>PyTuple_SET_ITEM</span>(__pyx_t_1, 1, __pyx_t_2);\n",
       "      __pyx_t_5 = 0;\n",
       "      __pyx_t_2 = 0;\n",
       "      if (unlikely(<span class='py_c_api'>PyObject_SetItem</span>(((PyObject *)__pyx_v_constraint_jac), __pyx_t_1, __pyx_int_1) &lt; 0)) __PYX_ERR(0, 53, __pyx_L1_error)\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_1); __pyx_t_1 = 0;\n",
       "</pre><pre class=\"cython line score-10\" onclick='toggleDiv(this)'>+<span class=\"\">54</span>:             <span class=\"n\">var_idx</span><span class=\"p\">:</span><span class=\"n\">var_idx</span> <span class=\"o\">+</span> <span class=\"n\">ais_len</span><span class=\"p\">]</span> <span class=\"o\">=</span> <span class=\"mf\">1</span></pre>\n",
       "<pre class='cython code score-10 '>      <span class='trace'>__Pyx_TraceLine</span>(54,0,__PYX_ERR(0, 54, __pyx_L1_error))\n",
       "      __pyx_t_6 = <span class='pyx_c_api'>__Pyx_PyInt_From_int</span>(__pyx_v_var_idx); if (unlikely(!__pyx_t_6)) __PYX_ERR(0, 54, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_6);\n",
       "      __pyx_t_2 = <span class='pyx_c_api'>__Pyx_PyInt_From_int</span>(__pyx_v_var_idx); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 54, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "      __pyx_t_1 = <span class='py_c_api'>PyNumber_Add</span>(__pyx_t_2, __pyx_v_ais_len); if (unlikely(!__pyx_t_1)) __PYX_ERR(0, 54, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_1);\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "</pre><pre class=\"cython line score-0\" onclick='toggleDiv(this)'>+<span class=\"\">55</span>:             <span class=\"n\">l_constraints</span><span class=\"p\">[</span><span class=\"n\">constraint_offset</span> <span class=\"o\">+</span> <span class=\"n\">idx</span><span class=\"p\">]</span> <span class=\"o\">=</span> \\</pre>\n",
       "<pre class='cython code score-0 '>      <span class='trace'>__Pyx_TraceLine</span>(55,0,__PYX_ERR(0, 55, __pyx_L1_error))\n",
       "      __pyx_t_20 = (__pyx_v_constraint_offset + __pyx_v_idx);\n",
       "      *__Pyx_BufPtrStrided1d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_l_constraints.rcbuffer-&gt;pybuffer.buf, __pyx_t_20, __pyx_pybuffernd_l_constraints.diminfo[0].strides) = __pyx_t_19;\n",
       "</pre><pre class=\"cython line score-34\" onclick='toggleDiv(this)'>+<span class=\"\">56</span>:                 <span class=\"p\">(</span><span class=\"nb\">sum</span><span class=\"p\">(</span><span class=\"n\">site_fracs</span><span class=\"p\">[</span><span class=\"n\">var_idx</span><span class=\"p\">:</span><span class=\"n\">var_idx</span> <span class=\"o\">+</span> <span class=\"n\">ais_len</span><span class=\"p\">])</span> <span class=\"o\">-</span> <span class=\"mf\">1</span><span class=\"p\">)</span></pre>\n",
       "<pre class='cython code score-34 '>      <span class='trace'>__Pyx_TraceLine</span>(56,0,__PYX_ERR(0, 56, __pyx_L1_error))\n",
       "      __pyx_t_1 = <span class='pyx_c_api'>__Pyx_PyInt_From_int</span>(__pyx_v_var_idx); if (unlikely(!__pyx_t_1)) __PYX_ERR(0, 56, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_1);\n",
       "      __pyx_t_2 = <span class='py_c_api'>PyNumber_Add</span>(__pyx_t_1, __pyx_v_ais_len); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 56, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_1); __pyx_t_1 = 0;\n",
       "      __pyx_t_1 = <span class='pyx_c_api'>__Pyx_PyObject_GetSlice</span>(((PyObject *)__pyx_v_site_fracs), __pyx_v_var_idx, 0, NULL, &amp;__pyx_t_2, NULL, 1, 0, 0); if (unlikely(!__pyx_t_1)) __PYX_ERR(0, 56, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_1);\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "      __pyx_t_2 = <span class='py_c_api'>PyTuple_New</span>(1); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 56, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "      <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_t_1);\n",
       "      <span class='py_macro_api'>PyTuple_SET_ITEM</span>(__pyx_t_2, 0, __pyx_t_1);\n",
       "      __pyx_t_1 = 0;\n",
       "      __pyx_t_1 = <span class='pyx_c_api'>__Pyx_PyObject_Call</span>(__pyx_builtin_sum, __pyx_t_2, NULL); if (unlikely(!__pyx_t_1)) __PYX_ERR(0, 56, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_1);\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "      __pyx_t_2 = <span class='pyx_c_api'>__Pyx_PyInt_SubtractObjC</span>(__pyx_t_1, __pyx_int_1, 1, 0); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 56, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_1); __pyx_t_1 = 0;\n",
       "      __pyx_t_19 = __pyx_<span class='py_c_api'>PyFloat_AsDouble</span>(__pyx_t_2); if (unlikely((__pyx_t_19 == ((npy_float64)-1)) &amp;&amp; <span class='py_c_api'>PyErr_Occurred</span>())) __PYX_ERR(0, 56, __pyx_L1_error)\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "</pre><pre class=\"cython line score-16\" onclick='toggleDiv(this)'>+<span class=\"\">57</span>:             <span class=\"n\">var_idx</span> <span class=\"o\">+=</span> <span class=\"n\">ais_len</span></pre>\n",
       "<pre class='cython code score-16 '>      <span class='trace'>__Pyx_TraceLine</span>(57,0,__PYX_ERR(0, 57, __pyx_L1_error))\n",
       "      __pyx_t_2 = <span class='pyx_c_api'>__Pyx_PyInt_From_int</span>(__pyx_v_var_idx); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 57, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "      __pyx_t_1 = <span class='py_c_api'>PyNumber_InPlaceAdd</span>(__pyx_t_2, __pyx_v_ais_len); if (unlikely(!__pyx_t_1)) __PYX_ERR(0, 57, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_1);\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "      __pyx_t_21 = <span class='pyx_c_api'>__Pyx_PyInt_As_int</span>(__pyx_t_1); if (unlikely((__pyx_t_21 == (int)-1) &amp;&amp; <span class='py_c_api'>PyErr_Occurred</span>())) __PYX_ERR(0, 57, __pyx_L1_error)\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_1); __pyx_t_1 = 0;\n",
       "      __pyx_v_var_idx = __pyx_t_21;\n",
       "    }\n",
       "</pre><pre class=\"cython line score-17\" onclick='toggleDiv(this)'>+<span class=\"\">58</span>:         <span class=\"n\">constraint_offset</span> <span class=\"o\">+=</span> <span class=\"nb\">len</span><span class=\"p\">(</span><span class=\"n\">dbf</span><span class=\"o\">.</span><span class=\"n\">phases</span><span class=\"p\">[</span><span class=\"n\">name</span><span class=\"p\">]</span><span class=\"o\">.</span><span class=\"n\">sublattices</span><span class=\"p\">)</span></pre>\n",
       "<pre class='cython code score-17 '>    <span class='trace'>__Pyx_TraceLine</span>(58,0,__PYX_ERR(0, 58, __pyx_L1_error))\n",
       "    __pyx_t_1 = <span class='pyx_c_api'>__Pyx_PyObject_GetAttrStr</span>(__pyx_v_dbf, __pyx_n_s_phases); if (unlikely(!__pyx_t_1)) __PYX_ERR(0, 58, __pyx_L1_error)\n",
       "    <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_1);\n",
       "    __pyx_t_2 = <span class='py_c_api'>PyObject_GetItem</span>(__pyx_t_1, __pyx_v_name); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 58, __pyx_L1_error)\n",
       "    <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "    <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_1); __pyx_t_1 = 0;\n",
       "    __pyx_t_1 = <span class='pyx_c_api'>__Pyx_PyObject_GetAttrStr</span>(__pyx_t_2, __pyx_n_s_sublattices); if (unlikely(!__pyx_t_1)) __PYX_ERR(0, 58, __pyx_L1_error)\n",
       "    <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_1);\n",
       "    <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "    __pyx_t_7 = <span class='py_c_api'>PyObject_Length</span>(__pyx_t_1); if (unlikely(__pyx_t_7 == -1)) __PYX_ERR(0, 58, __pyx_L1_error)\n",
       "    <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_1); __pyx_t_1 = 0;\n",
       "    __pyx_v_constraint_offset = (__pyx_v_constraint_offset + __pyx_t_7);\n",
       "  }\n",
       "</pre><pre class=\"cython line score-0\">&#xA0;<span class=\"\">59</span>:     <span class=\"c\"># Second: Mass balance of each component</span></pre>\n",
       "<pre class=\"cython line score-46\" onclick='toggleDiv(this)'>+<span class=\"\">60</span>:     <span class=\"k\">for</span> <span class=\"n\">comp_idx</span><span class=\"p\">,</span> <span class=\"n\">comp</span> <span class=\"ow\">in</span> <span class=\"nb\">enumerate</span><span class=\"p\">(</span><span class=\"n\">comps</span><span class=\"p\">):</span></pre>\n",
       "<pre class='cython code score-46 '>  <span class='trace'>__Pyx_TraceLine</span>(60,0,__PYX_ERR(0, 60, __pyx_L1_error))\n",
       "  __pyx_t_8 = 0;\n",
       "  if (likely(<span class='py_c_api'>PyList_CheckExact</span>(__pyx_v_comps)) || <span class='py_c_api'>PyTuple_CheckExact</span>(__pyx_v_comps)) {\n",
       "    __pyx_t_1 = __pyx_v_comps; <span class='pyx_macro_api'>__Pyx_INCREF</span>(__pyx_t_1); __pyx_t_7 = 0;\n",
       "    __pyx_t_4 = NULL;\n",
       "  } else {\n",
       "    __pyx_t_7 = -1; __pyx_t_1 = <span class='py_c_api'>PyObject_GetIter</span>(__pyx_v_comps); if (unlikely(!__pyx_t_1)) __PYX_ERR(0, 60, __pyx_L1_error)\n",
       "    <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_1);\n",
       "    __pyx_t_4 = Py_TYPE(__pyx_t_1)-&gt;tp_iternext; if (unlikely(!__pyx_t_4)) __PYX_ERR(0, 60, __pyx_L1_error)\n",
       "  }\n",
       "  for (;;) {\n",
       "    if (likely(!__pyx_t_4)) {\n",
       "      if (likely(<span class='py_c_api'>PyList_CheckExact</span>(__pyx_t_1))) {\n",
       "        if (__pyx_t_7 &gt;= <span class='py_macro_api'>PyList_GET_SIZE</span>(__pyx_t_1)) break;\n",
       "        #if CYTHON_ASSUME_SAFE_MACROS &amp;&amp; !CYTHON_AVOID_BORROWED_REFS\n",
       "        __pyx_t_2 = <span class='py_macro_api'>PyList_GET_ITEM</span>(__pyx_t_1, __pyx_t_7); <span class='pyx_macro_api'>__Pyx_INCREF</span>(__pyx_t_2); __pyx_t_7++; if (unlikely(0 &lt; 0)) __PYX_ERR(0, 60, __pyx_L1_error)\n",
       "        #else\n",
       "        __pyx_t_2 = <span class='py_macro_api'>PySequence_ITEM</span>(__pyx_t_1, __pyx_t_7); __pyx_t_7++; if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 60, __pyx_L1_error)\n",
       "        <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "        #endif\n",
       "      } else {\n",
       "        if (__pyx_t_7 &gt;= <span class='py_macro_api'>PyTuple_GET_SIZE</span>(__pyx_t_1)) break;\n",
       "        #if CYTHON_ASSUME_SAFE_MACROS &amp;&amp; !CYTHON_AVOID_BORROWED_REFS\n",
       "        __pyx_t_2 = <span class='py_macro_api'>PyTuple_GET_ITEM</span>(__pyx_t_1, __pyx_t_7); <span class='pyx_macro_api'>__Pyx_INCREF</span>(__pyx_t_2); __pyx_t_7++; if (unlikely(0 &lt; 0)) __PYX_ERR(0, 60, __pyx_L1_error)\n",
       "        #else\n",
       "        __pyx_t_2 = <span class='py_macro_api'>PySequence_ITEM</span>(__pyx_t_1, __pyx_t_7); __pyx_t_7++; if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 60, __pyx_L1_error)\n",
       "        <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "        #endif\n",
       "      }\n",
       "    } else {\n",
       "      __pyx_t_2 = __pyx_t_4(__pyx_t_1);\n",
       "      if (unlikely(!__pyx_t_2)) {\n",
       "        PyObject* exc_type = <span class='py_c_api'>PyErr_Occurred</span>();\n",
       "        if (exc_type) {\n",
       "          if (likely(exc_type == PyExc_StopIteration || <span class='py_c_api'>PyErr_GivenExceptionMatches</span>(exc_type, PyExc_StopIteration))) <span class='py_c_api'>PyErr_Clear</span>();\n",
       "          else __PYX_ERR(0, 60, __pyx_L1_error)\n",
       "        }\n",
       "        break;\n",
       "      }\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "    }\n",
       "    <span class='pyx_macro_api'>__Pyx_XDECREF_SET</span>(__pyx_v_comp, __pyx_t_2);\n",
       "    __pyx_t_2 = 0;\n",
       "    __pyx_v_comp_idx = __pyx_t_8;\n",
       "    __pyx_t_8 = (__pyx_t_8 + 1);\n",
       "/* … */\n",
       "    <span class='trace'>__Pyx_TraceLine</span>(60,0,__PYX_ERR(0, 60, __pyx_L1_error))\n",
       "    __pyx_L34_continue:;\n",
       "  }\n",
       "  <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_1); __pyx_t_1 = 0;\n",
       "</pre><pre class=\"cython line score-2\" onclick='toggleDiv(this)'>+<span class=\"\">61</span>:         <span class=\"k\">if</span> <span class=\"n\">comp</span> <span class=\"o\">==</span> <span class=\"s\">&#39;VA&#39;</span><span class=\"p\">:</span></pre>\n",
       "<pre class='cython code score-2 '>    <span class='trace'>__Pyx_TraceLine</span>(61,0,__PYX_ERR(0, 61, __pyx_L1_error))\n",
       "    __pyx_t_10 = (<span class='pyx_c_api'>__Pyx_PyUnicode_Equals</span>(__pyx_v_comp, __pyx_n_u_VA, Py_EQ)); if (unlikely(__pyx_t_10 &lt; 0)) __PYX_ERR(0, 61, __pyx_L1_error)\n",
       "    if (__pyx_t_10) {\n",
       "/* … */\n",
       "    }\n",
       "</pre><pre class=\"cython line score-0\" onclick='toggleDiv(this)'>+<span class=\"\">62</span>:             <span class=\"k\">continue</span></pre>\n",
       "<pre class='cython code score-0 '>      <span class='trace'>__Pyx_TraceLine</span>(62,0,__PYX_ERR(0, 62, __pyx_L1_error))\n",
       "      goto __pyx_L34_continue;\n",
       "</pre><pre class=\"cython line score-0\" onclick='toggleDiv(this)'>+<span class=\"\">63</span>:         <span class=\"n\">var_offset</span> <span class=\"o\">=</span> <span class=\"mf\">0</span></pre>\n",
       "<pre class='cython code score-0 '>    <span class='trace'>__Pyx_TraceLine</span>(63,0,__PYX_ERR(0, 63, __pyx_L1_error))\n",
       "    __pyx_v_var_offset = 0;\n",
       "</pre><pre class=\"cython line score-0\" onclick='toggleDiv(this)'>+<span class=\"\">64</span>:         <span class=\"n\">phase_idx</span> <span class=\"o\">=</span> <span class=\"mf\">0</span></pre>\n",
       "<pre class='cython code score-0 '>    <span class='trace'>__Pyx_TraceLine</span>(64,0,__PYX_ERR(0, 64, __pyx_L1_error))\n",
       "    __pyx_v_phase_idx = 0;\n",
       "</pre><pre class=\"cython line score-0\" onclick='toggleDiv(this)'>+<span class=\"\">65</span>:         <span class=\"k\">for</span> <span class=\"n\">phase_idx</span> <span class=\"ow\">in</span> <span class=\"nb\">range</span><span class=\"p\">(</span><span class=\"n\">num_phases</span><span class=\"p\">):</span></pre>\n",
       "<pre class='cython code score-0 '>    <span class='trace'>__Pyx_TraceLine</span>(65,0,__PYX_ERR(0, 65, __pyx_L1_error))\n",
       "    __pyx_t_9 = __pyx_v_num_phases;\n",
       "    for (__pyx_t_18 = 0; __pyx_t_18 &lt; __pyx_t_9; __pyx_t_18+=1) {\n",
       "      __pyx_v_phase_idx = __pyx_t_18;\n",
       "</pre><pre class=\"cython line score-3\" onclick='toggleDiv(this)'>+<span class=\"\">66</span>:             <span class=\"n\">name</span> <span class=\"o\">=</span> <span class=\"n\">phases</span><span class=\"p\">[</span><span class=\"n\">phase_idx</span><span class=\"p\">]</span></pre>\n",
       "<pre class='cython code score-3 '>      <span class='trace'>__Pyx_TraceLine</span>(66,0,__PYX_ERR(0, 66, __pyx_L1_error))\n",
       "      __pyx_t_2 = <span class='pyx_c_api'>__Pyx_GetItemInt</span>(__pyx_v_phases, __pyx_v_phase_idx, int, 1, __Pyx_PyInt_From_int, 0, 0, 0); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 66, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "      <span class='pyx_macro_api'>__Pyx_XDECREF_SET</span>(__pyx_v_name, __pyx_t_2);\n",
       "      __pyx_t_2 = 0;\n",
       "</pre><pre class=\"cython line score-8\" onclick='toggleDiv(this)'>+<span class=\"\">67</span>:             <span class=\"n\">prn</span> <span class=\"o\">=</span> <span class=\"n\">phase_records</span><span class=\"p\">[</span><span class=\"n\">name</span><span class=\"p\">]</span></pre>\n",
       "<pre class='cython code score-8 '>      <span class='trace'>__Pyx_TraceLine</span>(67,0,__PYX_ERR(0, 67, __pyx_L1_error))\n",
       "      __pyx_t_2 = <span class='py_c_api'>PyObject_GetItem</span>(__pyx_v_phase_records, __pyx_v_name); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 67, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "      if (!(likely(((__pyx_t_2) == Py_None) || likely(<span class='pyx_c_api'>__Pyx_TypeTest</span>(__pyx_t_2, __pyx_ptype_9pycalphad_4core_9phase_rec_PhaseRecord))))) __PYX_ERR(0, 67, __pyx_L1_error)\n",
       "      <span class='pyx_macro_api'>__Pyx_XDECREF_SET</span>(__pyx_v_prn, ((struct PhaseRecordObject *)__pyx_t_2));\n",
       "      __pyx_t_2 = 0;\n",
       "</pre><pre class=\"cython line score-0\" onclick='toggleDiv(this)'>+<span class=\"\">68</span>:             <span class=\"n\">phase_frac</span> <span class=\"o\">=</span> <span class=\"n\">phase_fracs</span><span class=\"p\">[</span><span class=\"n\">phase_idx</span><span class=\"p\">]</span></pre>\n",
       "<pre class='cython code score-0 '>      <span class='trace'>__Pyx_TraceLine</span>(68,0,__PYX_ERR(0, 68, __pyx_L1_error))\n",
       "      __pyx_t_22 = __pyx_v_phase_idx;\n",
       "      __pyx_v_phase_frac = (*__Pyx_BufPtrStrided1d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_phase_fracs.rcbuffer-&gt;pybuffer.buf, __pyx_t_22, __pyx_pybuffernd_phase_fracs.diminfo[0].strides));\n",
       "</pre><pre class=\"cython line score-5\" onclick='toggleDiv(this)'>+<span class=\"\">69</span>:             <span class=\"n\">spidx</span> <span class=\"o\">=</span> <span class=\"nb\">len</span><span class=\"p\">(</span><span class=\"n\">site_fracs</span><span class=\"p\">)</span> <span class=\"o\">+</span> <span class=\"n\">phase_idx</span></pre>\n",
       "<pre class='cython code score-5 '>      <span class='trace'>__Pyx_TraceLine</span>(69,0,__PYX_ERR(0, 69, __pyx_L1_error))\n",
       "      __pyx_t_3 = <span class='py_c_api'>PyObject_Length</span>(((PyObject *)__pyx_v_site_fracs)); if (unlikely(__pyx_t_3 == -1)) __PYX_ERR(0, 69, __pyx_L1_error)\n",
       "      __pyx_v_spidx = (__pyx_t_3 + __pyx_v_phase_idx);\n",
       "</pre><pre class=\"cython line score-5\" onclick='toggleDiv(this)'>+<span class=\"\">70</span>:             <span class=\"n\">sfview</span> <span class=\"o\">=</span> <span class=\"n\">site_fracs</span><span class=\"p\">[</span><span class=\"n\">var_offset</span><span class=\"p\">:</span><span class=\"n\">var_offset</span> <span class=\"o\">+</span> <span class=\"n\">phase_dof</span><span class=\"p\">[</span><span class=\"n\">phase_idx</span><span class=\"p\">]]</span></pre>\n",
       "<pre class='cython code score-5 '>      <span class='trace'>__Pyx_TraceLine</span>(70,0,__PYX_ERR(0, 70, __pyx_L1_error))\n",
       "      __pyx_t_23 = __pyx_v_phase_idx;\n",
       "      __pyx_t_2 = <span class='pyx_c_api'>__Pyx_PyObject_GetSlice</span>(((PyObject *)__pyx_v_site_fracs), __pyx_v_var_offset, (__pyx_v_var_offset + (*__Pyx_BufPtrStrided1d(__pyx_t_5numpy_int_t *, __pyx_pybuffernd_phase_dof.rcbuffer-&gt;pybuffer.buf, __pyx_t_23, __pyx_pybuffernd_phase_dof.diminfo[0].strides))), NULL, NULL, NULL, 1, 1, 0); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 70, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "      __pyx_t_13 = <span class='pyx_c_api'>__Pyx_PyObject_to_MemoryviewSlice_dc_double</span>(__pyx_t_2);\n",
       "      if (unlikely(!__pyx_t_13.memview)) __PYX_ERR(0, 70, __pyx_L1_error)\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "      __PYX_XDEC_MEMVIEW(&amp;__pyx_v_sfview, 1);\n",
       "      __pyx_v_sfview = __pyx_t_13;\n",
       "      __pyx_t_13.memview = NULL;\n",
       "      __pyx_t_13.data = NULL;\n",
       "</pre><pre class=\"cython line score-40\" onclick='toggleDiv(this)'>+<span class=\"\">71</span>:             <span class=\"n\">comp_hess_value</span> <span class=\"o\">=</span> <span class=\"n\">np</span><span class=\"o\">.</span><span class=\"n\">zeros</span><span class=\"p\">((</span><span class=\"n\">phase_dof</span><span class=\"p\">[</span><span class=\"n\">phase_idx</span><span class=\"p\">],</span> <span class=\"n\">phase_dof</span><span class=\"p\">[</span><span class=\"n\">phase_idx</span><span class=\"p\">]),</span> <span class=\"n\">order</span><span class=\"o\">=</span><span class=\"s\">&#39;F&#39;</span><span class=\"p\">)</span></pre>\n",
       "<pre class='cython code score-40 '>      <span class='trace'>__Pyx_TraceLine</span>(71,0,__PYX_ERR(0, 71, __pyx_L1_error))\n",
       "      __pyx_t_2 = <span class='pyx_c_api'>__Pyx_GetModuleGlobalName</span>(__pyx_n_s_np); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 71, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "      __pyx_t_5 = <span class='pyx_c_api'>__Pyx_PyObject_GetAttrStr</span>(__pyx_t_2, __pyx_n_s_zeros); if (unlikely(!__pyx_t_5)) __PYX_ERR(0, 71, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_5);\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "      __pyx_t_24 = __pyx_v_phase_idx;\n",
       "      __pyx_t_2 = <span class='pyx_c_api'>__Pyx_PyInt_From_npy_long</span>((*__Pyx_BufPtrStrided1d(__pyx_t_5numpy_int_t *, __pyx_pybuffernd_phase_dof.rcbuffer-&gt;pybuffer.buf, __pyx_t_24, __pyx_pybuffernd_phase_dof.diminfo[0].strides))); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 71, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "      __pyx_t_25 = __pyx_v_phase_idx;\n",
       "      __pyx_t_6 = <span class='pyx_c_api'>__Pyx_PyInt_From_npy_long</span>((*__Pyx_BufPtrStrided1d(__pyx_t_5numpy_int_t *, __pyx_pybuffernd_phase_dof.rcbuffer-&gt;pybuffer.buf, __pyx_t_25, __pyx_pybuffernd_phase_dof.diminfo[0].strides))); if (unlikely(!__pyx_t_6)) __PYX_ERR(0, 71, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_6);\n",
       "      __pyx_t_12 = <span class='py_c_api'>PyTuple_New</span>(2); if (unlikely(!__pyx_t_12)) __PYX_ERR(0, 71, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_12);\n",
       "      <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_t_2);\n",
       "      <span class='py_macro_api'>PyTuple_SET_ITEM</span>(__pyx_t_12, 0, __pyx_t_2);\n",
       "      <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_t_6);\n",
       "      <span class='py_macro_api'>PyTuple_SET_ITEM</span>(__pyx_t_12, 1, __pyx_t_6);\n",
       "      __pyx_t_2 = 0;\n",
       "      __pyx_t_6 = 0;\n",
       "      __pyx_t_6 = <span class='py_c_api'>PyTuple_New</span>(1); if (unlikely(!__pyx_t_6)) __PYX_ERR(0, 71, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_6);\n",
       "      <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_t_12);\n",
       "      <span class='py_macro_api'>PyTuple_SET_ITEM</span>(__pyx_t_6, 0, __pyx_t_12);\n",
       "      __pyx_t_12 = 0;\n",
       "      __pyx_t_12 = <span class='py_c_api'>PyDict_New</span>(); if (unlikely(!__pyx_t_12)) __PYX_ERR(0, 71, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_12);\n",
       "      if (<span class='py_c_api'>PyDict_SetItem</span>(__pyx_t_12, __pyx_n_s_order, __pyx_n_u_F) &lt; 0) __PYX_ERR(0, 71, __pyx_L1_error)\n",
       "      __pyx_t_2 = <span class='pyx_c_api'>__Pyx_PyObject_Call</span>(__pyx_t_5, __pyx_t_6, __pyx_t_12); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 71, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_5); __pyx_t_5 = 0;\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_6); __pyx_t_6 = 0;\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_12); __pyx_t_12 = 0;\n",
       "      __pyx_t_26 = <span class='pyx_c_api'>__Pyx_PyObject_to_MemoryviewSlice_dcd__double</span>(__pyx_t_2);\n",
       "      if (unlikely(!__pyx_t_26.memview)) __PYX_ERR(0, 71, __pyx_L1_error)\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "      __PYX_XDEC_MEMVIEW(&amp;__pyx_v_comp_hess_value, 1);\n",
       "      __pyx_v_comp_hess_value = __pyx_t_26;\n",
       "      __pyx_t_26.memview = NULL;\n",
       "      __pyx_t_26.data = NULL;\n",
       "</pre><pre class=\"cython line score-49\" onclick='toggleDiv(this)'>+<span class=\"\">72</span>:             <span class=\"n\">comp_grad_value</span> <span class=\"o\">=</span> <span class=\"n\">np</span><span class=\"o\">.</span><span class=\"n\">zeros</span><span class=\"p\">(</span><span class=\"n\">phase_dof</span><span class=\"p\">[</span><span class=\"n\">phase_idx</span><span class=\"p\">])</span></pre>\n",
       "<pre class='cython code score-49 '>      <span class='trace'>__Pyx_TraceLine</span>(72,0,__PYX_ERR(0, 72, __pyx_L1_error))\n",
       "      __pyx_t_12 = <span class='pyx_c_api'>__Pyx_GetModuleGlobalName</span>(__pyx_n_s_np); if (unlikely(!__pyx_t_12)) __PYX_ERR(0, 72, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_12);\n",
       "      __pyx_t_6 = <span class='pyx_c_api'>__Pyx_PyObject_GetAttrStr</span>(__pyx_t_12, __pyx_n_s_zeros); if (unlikely(!__pyx_t_6)) __PYX_ERR(0, 72, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_6);\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_12); __pyx_t_12 = 0;\n",
       "      __pyx_t_27 = __pyx_v_phase_idx;\n",
       "      __pyx_t_12 = <span class='pyx_c_api'>__Pyx_PyInt_From_npy_long</span>((*__Pyx_BufPtrStrided1d(__pyx_t_5numpy_int_t *, __pyx_pybuffernd_phase_dof.rcbuffer-&gt;pybuffer.buf, __pyx_t_27, __pyx_pybuffernd_phase_dof.diminfo[0].strides))); if (unlikely(!__pyx_t_12)) __PYX_ERR(0, 72, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_12);\n",
       "      __pyx_t_5 = NULL;\n",
       "      if (CYTHON_UNPACK_METHODS &amp;&amp; unlikely(<span class='py_c_api'>PyMethod_Check</span>(__pyx_t_6))) {\n",
       "        __pyx_t_5 = <span class='py_macro_api'>PyMethod_GET_SELF</span>(__pyx_t_6);\n",
       "        if (likely(__pyx_t_5)) {\n",
       "          PyObject* function = <span class='py_macro_api'>PyMethod_GET_FUNCTION</span>(__pyx_t_6);\n",
       "          <span class='pyx_macro_api'>__Pyx_INCREF</span>(__pyx_t_5);\n",
       "          <span class='pyx_macro_api'>__Pyx_INCREF</span>(function);\n",
       "          <span class='pyx_macro_api'>__Pyx_DECREF_SET</span>(__pyx_t_6, function);\n",
       "        }\n",
       "      }\n",
       "      if (!__pyx_t_5) {\n",
       "        __pyx_t_2 = <span class='pyx_c_api'>__Pyx_PyObject_CallOneArg</span>(__pyx_t_6, __pyx_t_12); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 72, __pyx_L1_error)\n",
       "        <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_12); __pyx_t_12 = 0;\n",
       "        <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "      } else {\n",
       "        #if CYTHON_FAST_PYCALL\n",
       "        if (<span class='py_c_api'>PyFunction_Check</span>(__pyx_t_6)) {\n",
       "          PyObject *__pyx_temp[2] = {__pyx_t_5, __pyx_t_12};\n",
       "          __pyx_t_2 = <span class='pyx_c_api'>__Pyx_PyFunction_FastCall</span>(__pyx_t_6, __pyx_temp+1-1, 1+1); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 72, __pyx_L1_error)\n",
       "          <span class='pyx_macro_api'>__Pyx_XDECREF</span>(__pyx_t_5); __pyx_t_5 = 0;\n",
       "          <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "          <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_12); __pyx_t_12 = 0;\n",
       "        } else\n",
       "        #endif\n",
       "        #if CYTHON_FAST_PYCCALL\n",
       "        if (<span class='pyx_c_api'>__Pyx_PyFastCFunction_Check</span>(__pyx_t_6)) {\n",
       "          PyObject *__pyx_temp[2] = {__pyx_t_5, __pyx_t_12};\n",
       "          __pyx_t_2 = <span class='pyx_c_api'>__Pyx_PyCFunction_FastCall</span>(__pyx_t_6, __pyx_temp+1-1, 1+1); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 72, __pyx_L1_error)\n",
       "          <span class='pyx_macro_api'>__Pyx_XDECREF</span>(__pyx_t_5); __pyx_t_5 = 0;\n",
       "          <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "          <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_12); __pyx_t_12 = 0;\n",
       "        } else\n",
       "        #endif\n",
       "        {\n",
       "          __pyx_t_28 = <span class='py_c_api'>PyTuple_New</span>(1+1); if (unlikely(!__pyx_t_28)) __PYX_ERR(0, 72, __pyx_L1_error)\n",
       "          <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_28);\n",
       "          <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_t_5); <span class='py_macro_api'>PyTuple_SET_ITEM</span>(__pyx_t_28, 0, __pyx_t_5); __pyx_t_5 = NULL;\n",
       "          <span class='refnanny'>__Pyx_GIVEREF</span>(__pyx_t_12);\n",
       "          <span class='py_macro_api'>PyTuple_SET_ITEM</span>(__pyx_t_28, 0+1, __pyx_t_12);\n",
       "          __pyx_t_12 = 0;\n",
       "          __pyx_t_2 = <span class='pyx_c_api'>__Pyx_PyObject_Call</span>(__pyx_t_6, __pyx_t_28, NULL); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 72, __pyx_L1_error)\n",
       "          <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "          <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_28); __pyx_t_28 = 0;\n",
       "        }\n",
       "      }\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_6); __pyx_t_6 = 0;\n",
       "      __pyx_t_13 = <span class='pyx_c_api'>__Pyx_PyObject_to_MemoryviewSlice_dc_double</span>(__pyx_t_2);\n",
       "      if (unlikely(!__pyx_t_13.memview)) __PYX_ERR(0, 72, __pyx_L1_error)\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "      __PYX_XDEC_MEMVIEW(&amp;__pyx_v_comp_grad_value, 1);\n",
       "      __pyx_v_comp_grad_value = __pyx_t_13;\n",
       "      __pyx_t_13.memview = NULL;\n",
       "      __pyx_t_13.data = NULL;\n",
       "</pre><pre class=\"cython line score-0\" onclick='toggleDiv(this)'>+<span class=\"\">73</span>:             <span class=\"k\">with</span> <span class=\"k\">nogil</span><span class=\"p\">:</span></pre>\n",
       "<pre class='cython code score-0 '>      <span class='trace'>__Pyx_TraceLine</span>(73,0,__PYX_ERR(0, 73, __pyx_L1_error))\n",
       "      {\n",
       "          #ifdef WITH_THREAD\n",
       "          PyThreadState *_save;\n",
       "          Py_UNBLOCK_THREADS\n",
       "          #endif\n",
       "          /*try:*/ {\n",
       "/* … */\n",
       "          <span class='trace'>__Pyx_TraceLine</span>(73,1,__PYX_ERR(0, 73, __pyx_L1_error))\n",
       "          /*finally:*/ {\n",
       "            /*normal exit:*/{\n",
       "              #ifdef WITH_THREAD\n",
       "              Py_BLOCK_THREADS\n",
       "              #endif\n",
       "              goto __pyx_L43;\n",
       "            }\n",
       "            __pyx_L42_error: {\n",
       "              #ifdef WITH_THREAD\n",
       "              Py_BLOCK_THREADS\n",
       "              #endif\n",
       "              goto __pyx_L1_error;\n",
       "            }\n",
       "            __pyx_L43:;\n",
       "          }\n",
       "      }\n",
       "    }\n",
       "</pre><pre class=\"cython line score-0\" onclick='toggleDiv(this)'>+<span class=\"\">74</span>:                 <span class=\"n\">mass_obj</span><span class=\"p\">(</span><span class=\"n\">prn</span><span class=\"p\">,</span> <span class=\"n\">comp_obj_value</span><span class=\"p\">,</span> <span class=\"n\">sfview</span><span class=\"p\">,</span> <span class=\"n\">comp_idx</span><span class=\"p\">)</span></pre>\n",
       "<pre class='cython code score-0 '>            <span class='trace'>__Pyx_TraceLine</span>(74,1,__PYX_ERR(0, 74, __pyx_L42_error))\n",
       "            __pyx_f_9pycalphad_4core_9phase_rec_mass_obj(__pyx_v_prn, __pyx_v_comp_obj_value, __pyx_v_sfview, __pyx_v_comp_idx);\n",
       "</pre><pre class=\"cython line score-0\" onclick='toggleDiv(this)'>+<span class=\"\">75</span>:                 <span class=\"n\">mass_grad</span><span class=\"p\">(</span><span class=\"n\">prn</span><span class=\"p\">,</span> <span class=\"n\">comp_grad_value</span><span class=\"p\">,</span> <span class=\"n\">sfview</span><span class=\"p\">,</span> <span class=\"n\">comp_idx</span><span class=\"p\">)</span></pre>\n",
       "<pre class='cython code score-0 '>            <span class='trace'>__Pyx_TraceLine</span>(75,1,__PYX_ERR(0, 75, __pyx_L42_error))\n",
       "            __pyx_f_9pycalphad_4core_9phase_rec_mass_grad(__pyx_v_prn, __pyx_v_comp_grad_value, __pyx_v_sfview, __pyx_v_comp_idx);\n",
       "</pre><pre class=\"cython line score-0\" onclick='toggleDiv(this)'>+<span class=\"\">76</span>:                 <span class=\"n\">mass_hess</span><span class=\"p\">(</span><span class=\"n\">prn</span><span class=\"p\">,</span> <span class=\"n\">comp_hess_value</span><span class=\"p\">,</span> <span class=\"n\">sfview</span><span class=\"p\">,</span> <span class=\"n\">comp_idx</span><span class=\"p\">)</span></pre>\n",
       "<pre class='cython code score-0 '>            <span class='trace'>__Pyx_TraceLine</span>(76,1,__PYX_ERR(0, 76, __pyx_L42_error))\n",
       "            __pyx_f_9pycalphad_4core_9phase_rec_mass_hess(__pyx_v_prn, __pyx_v_comp_hess_value, __pyx_v_sfview, __pyx_v_comp_idx);\n",
       "</pre><pre class=\"cython line score-0\">&#xA0;<span class=\"\">77</span>:                 <span class=\"c\"># current phase frac times the comp_grad</span></pre>\n",
       "<pre class=\"cython line score-0\" onclick='toggleDiv(this)'>+<span class=\"\">78</span>:                 <span class=\"k\">for</span> <span class=\"n\">grad_idx</span> <span class=\"ow\">in</span> <span class=\"nb\">range</span><span class=\"p\">(</span><span class=\"n\">var_offset</span><span class=\"p\">,</span> <span class=\"n\">var_offset</span> <span class=\"o\">+</span> <span class=\"n\">phase_dof</span><span class=\"p\">[</span><span class=\"n\">phase_idx</span><span class=\"p\">]):</span></pre>\n",
       "<pre class='cython code score-0 '>            <span class='trace'>__Pyx_TraceLine</span>(78,1,__PYX_ERR(0, 78, __pyx_L42_error))\n",
       "            __pyx_t_29 = __pyx_v_phase_idx;\n",
       "            __pyx_t_30 = (__pyx_v_var_offset + (*__Pyx_BufPtrStrided1d(__pyx_t_5numpy_int_t *, __pyx_pybuffernd_phase_dof.rcbuffer-&gt;pybuffer.buf, __pyx_t_29, __pyx_pybuffernd_phase_dof.diminfo[0].strides)));\n",
       "            for (__pyx_t_21 = __pyx_v_var_offset; __pyx_t_21 &lt; __pyx_t_30; __pyx_t_21+=1) {\n",
       "              __pyx_v_grad_idx = __pyx_t_21;\n",
       "</pre><pre class=\"cython line score-0\" onclick='toggleDiv(this)'>+<span class=\"\">79</span>:                     <span class=\"n\">constraint_jac</span><span class=\"p\">[</span><span class=\"n\">constraint_offset</span><span class=\"p\">,</span> <span class=\"n\">grad_idx</span><span class=\"p\">]</span> <span class=\"o\">=</span> \\</pre>\n",
       "<pre class='cython code score-0 '>              <span class='trace'>__Pyx_TraceLine</span>(79,1,__PYX_ERR(0, 79, __pyx_L42_error))\n",
       "              __pyx_t_32 = __pyx_v_constraint_offset;\n",
       "              __pyx_t_33 = __pyx_v_grad_idx;\n",
       "              *__Pyx_BufPtrStrided2d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_constraint_jac.rcbuffer-&gt;pybuffer.buf, __pyx_t_32, __pyx_pybuffernd_constraint_jac.diminfo[0].strides, __pyx_t_33, __pyx_pybuffernd_constraint_jac.diminfo[1].strides) = (__pyx_v_phase_frac * (*((double *) ( /* dim=0 */ ((char *) (((double *) __pyx_v_comp_grad_value.data) + __pyx_t_31)) ))));\n",
       "</pre><pre class=\"cython line score-0\" onclick='toggleDiv(this)'>+<span class=\"\">80</span>:                         <span class=\"n\">phase_frac</span> <span class=\"o\">*</span> <span class=\"n\">comp_grad_value</span><span class=\"p\">[</span><span class=\"n\">grad_idx</span> <span class=\"o\">-</span> <span class=\"n\">var_offset</span><span class=\"p\">]</span></pre>\n",
       "<pre class='cython code score-0 '>              <span class='trace'>__Pyx_TraceLine</span>(80,1,__PYX_ERR(0, 80, __pyx_L42_error))\n",
       "              __pyx_t_31 = (__pyx_v_grad_idx - __pyx_v_var_offset);\n",
       "</pre><pre class=\"cython line score-0\" onclick='toggleDiv(this)'>+<span class=\"\">81</span>:                     <span class=\"n\">constraint_hess</span><span class=\"p\">[</span><span class=\"n\">constraint_offset</span><span class=\"p\">,</span> <span class=\"n\">spidx</span><span class=\"p\">,</span> <span class=\"n\">grad_idx</span><span class=\"p\">]</span> <span class=\"o\">=</span> <span class=\"n\">comp_grad_value</span><span class=\"p\">[</span><span class=\"n\">grad_idx</span> <span class=\"o\">-</span> <span class=\"n\">var_offset</span><span class=\"p\">]</span></pre>\n",
       "<pre class='cython code score-0 '>              <span class='trace'>__Pyx_TraceLine</span>(81,1,__PYX_ERR(0, 81, __pyx_L42_error))\n",
       "              __pyx_t_34 = (__pyx_v_grad_idx - __pyx_v_var_offset);\n",
       "              __pyx_t_35 = __pyx_v_constraint_offset;\n",
       "              __pyx_t_36 = __pyx_v_spidx;\n",
       "              __pyx_t_37 = __pyx_v_grad_idx;\n",
       "              *__Pyx_BufPtrStrided3d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_constraint_hess.rcbuffer-&gt;pybuffer.buf, __pyx_t_35, __pyx_pybuffernd_constraint_hess.diminfo[0].strides, __pyx_t_36, __pyx_pybuffernd_constraint_hess.diminfo[1].strides, __pyx_t_37, __pyx_pybuffernd_constraint_hess.diminfo[2].strides) = (*((double *) ( /* dim=0 */ ((char *) (((double *) __pyx_v_comp_grad_value.data) + __pyx_t_34)) )));\n",
       "</pre><pre class=\"cython line score-0\" onclick='toggleDiv(this)'>+<span class=\"\">82</span>:                     <span class=\"n\">constraint_hess</span><span class=\"p\">[</span><span class=\"n\">constraint_offset</span><span class=\"p\">,</span> <span class=\"n\">grad_idx</span><span class=\"p\">,</span> <span class=\"n\">spidx</span><span class=\"p\">]</span> <span class=\"o\">=</span> <span class=\"n\">comp_grad_value</span><span class=\"p\">[</span><span class=\"n\">grad_idx</span> <span class=\"o\">-</span> <span class=\"n\">var_offset</span><span class=\"p\">]</span></pre>\n",
       "<pre class='cython code score-0 '>              <span class='trace'>__Pyx_TraceLine</span>(82,1,__PYX_ERR(0, 82, __pyx_L42_error))\n",
       "              __pyx_t_38 = (__pyx_v_grad_idx - __pyx_v_var_offset);\n",
       "              __pyx_t_39 = __pyx_v_constraint_offset;\n",
       "              __pyx_t_40 = __pyx_v_grad_idx;\n",
       "              __pyx_t_41 = __pyx_v_spidx;\n",
       "              *__Pyx_BufPtrStrided3d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_constraint_hess.rcbuffer-&gt;pybuffer.buf, __pyx_t_39, __pyx_pybuffernd_constraint_hess.diminfo[0].strides, __pyx_t_40, __pyx_pybuffernd_constraint_hess.diminfo[1].strides, __pyx_t_41, __pyx_pybuffernd_constraint_hess.diminfo[2].strides) = (*((double *) ( /* dim=0 */ ((char *) (((double *) __pyx_v_comp_grad_value.data) + __pyx_t_38)) )));\n",
       "</pre><pre class=\"cython line score-0\" onclick='toggleDiv(this)'>+<span class=\"\">83</span>:                     <span class=\"k\">for</span> <span class=\"n\">hess_idx</span> <span class=\"ow\">in</span> <span class=\"nb\">range</span><span class=\"p\">(</span><span class=\"n\">var_offset</span><span class=\"p\">,</span> <span class=\"n\">var_offset</span> <span class=\"o\">+</span> <span class=\"n\">phase_dof</span><span class=\"p\">[</span><span class=\"n\">phase_idx</span><span class=\"p\">]):</span></pre>\n",
       "<pre class='cython code score-0 '>              <span class='trace'>__Pyx_TraceLine</span>(83,1,__PYX_ERR(0, 83, __pyx_L42_error))\n",
       "              __pyx_t_42 = __pyx_v_phase_idx;\n",
       "              __pyx_t_43 = (__pyx_v_var_offset + (*__Pyx_BufPtrStrided1d(__pyx_t_5numpy_int_t *, __pyx_pybuffernd_phase_dof.rcbuffer-&gt;pybuffer.buf, __pyx_t_42, __pyx_pybuffernd_phase_dof.diminfo[0].strides)));\n",
       "              for (__pyx_t_44 = __pyx_v_var_offset; __pyx_t_44 &lt; __pyx_t_43; __pyx_t_44+=1) {\n",
       "                __pyx_v_hess_idx = __pyx_t_44;\n",
       "</pre><pre class=\"cython line score-0\" onclick='toggleDiv(this)'>+<span class=\"\">84</span>:                         <span class=\"n\">constraint_hess</span><span class=\"p\">[</span><span class=\"n\">constraint_offset</span><span class=\"p\">,</span> <span class=\"n\">grad_idx</span><span class=\"p\">,</span> <span class=\"n\">hess_idx</span><span class=\"p\">]</span> <span class=\"o\">=</span> <span class=\"n\">phase_frac</span> <span class=\"o\">*</span> <span class=\"n\">comp_hess_value</span><span class=\"p\">[</span><span class=\"n\">grad_idx</span> <span class=\"o\">-</span> <span class=\"n\">var_offset</span><span class=\"p\">,</span> <span class=\"n\">hess_idx</span> <span class=\"o\">-</span> <span class=\"n\">var_offset</span><span class=\"p\">]</span></pre>\n",
       "<pre class='cython code score-0 '>                <span class='trace'>__Pyx_TraceLine</span>(84,1,__PYX_ERR(0, 84, __pyx_L42_error))\n",
       "                __pyx_t_45 = (__pyx_v_grad_idx - __pyx_v_var_offset);\n",
       "                __pyx_t_46 = (__pyx_v_hess_idx - __pyx_v_var_offset);\n",
       "                __pyx_t_47 = __pyx_v_constraint_offset;\n",
       "                __pyx_t_48 = __pyx_v_grad_idx;\n",
       "                __pyx_t_49 = __pyx_v_hess_idx;\n",
       "                *__Pyx_BufPtrStrided3d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_constraint_hess.rcbuffer-&gt;pybuffer.buf, __pyx_t_47, __pyx_pybuffernd_constraint_hess.diminfo[0].strides, __pyx_t_48, __pyx_pybuffernd_constraint_hess.diminfo[1].strides, __pyx_t_49, __pyx_pybuffernd_constraint_hess.diminfo[2].strides) = (__pyx_v_phase_frac * (*((double *) ( /* dim=1 */ (( /* dim=0 */ ((char *) (((double *) __pyx_v_comp_hess_value.data) + __pyx_t_45)) ) + __pyx_t_46 * __pyx_v_comp_hess_value.strides[1]) ))));\n",
       "              }\n",
       "            }\n",
       "</pre><pre class=\"cython line score-0\" onclick='toggleDiv(this)'>+<span class=\"\">85</span>:                 <span class=\"n\">l_constraints</span><span class=\"p\">[</span><span class=\"n\">constraint_offset</span><span class=\"p\">]</span> <span class=\"o\">+=</span> <span class=\"n\">phase_frac</span> <span class=\"o\">*</span> <span class=\"n\">comp_obj_value</span><span class=\"p\">[</span><span class=\"mf\">0</span><span class=\"p\">]</span></pre>\n",
       "<pre class='cython code score-0 '>            <span class='trace'>__Pyx_TraceLine</span>(85,1,__PYX_ERR(0, 85, __pyx_L42_error))\n",
       "            __pyx_t_50 = 0;\n",
       "            __pyx_t_51 = __pyx_v_constraint_offset;\n",
       "            *__Pyx_BufPtrStrided1d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_l_constraints.rcbuffer-&gt;pybuffer.buf, __pyx_t_51, __pyx_pybuffernd_l_constraints.diminfo[0].strides) += (__pyx_v_phase_frac * (*((double *) ( /* dim=0 */ ((char *) (((double *) __pyx_v_comp_obj_value.data) + __pyx_t_50)) ))));\n",
       "</pre><pre class=\"cython line score-0\" onclick='toggleDiv(this)'>+<span class=\"\">86</span>:                 <span class=\"n\">constraint_jac</span><span class=\"p\">[</span><span class=\"n\">constraint_offset</span><span class=\"p\">,</span> <span class=\"n\">spidx</span><span class=\"p\">]</span> <span class=\"o\">+=</span> <span class=\"n\">comp_obj_value</span><span class=\"p\">[</span><span class=\"mf\">0</span><span class=\"p\">]</span></pre>\n",
       "<pre class='cython code score-0 '>            <span class='trace'>__Pyx_TraceLine</span>(86,1,__PYX_ERR(0, 86, __pyx_L42_error))\n",
       "            __pyx_t_52 = 0;\n",
       "            __pyx_t_53 = __pyx_v_constraint_offset;\n",
       "            __pyx_t_54 = __pyx_v_spidx;\n",
       "            *__Pyx_BufPtrStrided2d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_constraint_jac.rcbuffer-&gt;pybuffer.buf, __pyx_t_53, __pyx_pybuffernd_constraint_jac.diminfo[0].strides, __pyx_t_54, __pyx_pybuffernd_constraint_jac.diminfo[1].strides) += (*((double *) ( /* dim=0 */ ((char *) (((double *) __pyx_v_comp_obj_value.data) + __pyx_t_52)) )));\n",
       "</pre><pre class=\"cython line score-0\" onclick='toggleDiv(this)'>+<span class=\"\">87</span>:                 <span class=\"n\">var_offset</span> <span class=\"o\">+=</span> <span class=\"n\">phase_dof</span><span class=\"p\">[</span><span class=\"n\">phase_idx</span><span class=\"p\">]</span></pre>\n",
       "<pre class='cython code score-0 '>            <span class='trace'>__Pyx_TraceLine</span>(87,1,__PYX_ERR(0, 87, __pyx_L42_error))\n",
       "            __pyx_t_55 = __pyx_v_phase_idx;\n",
       "            __pyx_v_var_offset = (__pyx_v_var_offset + (*__Pyx_BufPtrStrided1d(__pyx_t_5numpy_int_t *, __pyx_pybuffernd_phase_dof.rcbuffer-&gt;pybuffer.buf, __pyx_t_55, __pyx_pybuffernd_phase_dof.diminfo[0].strides)));\n",
       "          }\n",
       "</pre><pre class=\"cython line score-8\" onclick='toggleDiv(this)'>+<span class=\"\">88</span>:         <span class=\"k\">if</span> <span class=\"n\">comp</span> <span class=\"o\">!=</span> <span class=\"n\">dependent_comp</span><span class=\"p\">:</span></pre>\n",
       "<pre class='cython code score-8 '>    <span class='trace'>__Pyx_TraceLine</span>(88,0,__PYX_ERR(0, 88, __pyx_L1_error))\n",
       "    __pyx_t_2 = <span class='py_c_api'>PyObject_RichCompare</span>(__pyx_v_comp, __pyx_v_dependent_comp, Py_NE); <span class='refnanny'>__Pyx_XGOTREF</span>(__pyx_t_2); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 88, __pyx_L1_error)\n",
       "    __pyx_t_10 = <span class='pyx_c_api'>__Pyx_PyObject_IsTrue</span>(__pyx_t_2); if (unlikely(__pyx_t_10 &lt; 0)) __PYX_ERR(0, 88, __pyx_L1_error)\n",
       "    <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "    if (__pyx_t_10) {\n",
       "/* … */\n",
       "      goto __pyx_L48;\n",
       "    }\n",
       "</pre><pre class=\"cython line score-19\" onclick='toggleDiv(this)'>+<span class=\"\">89</span>:             <span class=\"n\">l_constraints</span><span class=\"p\">[</span><span class=\"n\">constraint_offset</span><span class=\"p\">]</span> <span class=\"o\">-=</span> <span class=\"nb\">float</span><span class=\"p\">(</span><span class=\"n\">cur_conds</span><span class=\"p\">[</span><span class=\"s\">&#39;X_&#39;</span> <span class=\"o\">+</span> <span class=\"n\">comp</span><span class=\"p\">])</span></pre>\n",
       "<pre class='cython code score-19 '>      <span class='trace'>__Pyx_TraceLine</span>(89,0,__PYX_ERR(0, 89, __pyx_L1_error))\n",
       "      __pyx_t_2 = <span class='py_c_api'>PyNumber_Add</span>(__pyx_n_u_X, __pyx_v_comp); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 89, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_2);\n",
       "      __pyx_t_6 = <span class='py_c_api'>PyObject_GetItem</span>(__pyx_v_cur_conds, __pyx_t_2); if (unlikely(!__pyx_t_6)) __PYX_ERR(0, 89, __pyx_L1_error)\n",
       "      <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_6);\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_2); __pyx_t_2 = 0;\n",
       "      __pyx_t_11 = <span class='pyx_c_api'>__Pyx_PyObject_AsDouble</span>(__pyx_t_6); if (unlikely(__pyx_t_11 == ((double)-1) &amp;&amp; <span class='py_c_api'>PyErr_Occurred</span>())) __PYX_ERR(0, 89, __pyx_L1_error)\n",
       "      <span class='pyx_macro_api'>__Pyx_DECREF</span>(__pyx_t_6); __pyx_t_6 = 0;\n",
       "      __pyx_t_56 = __pyx_v_constraint_offset;\n",
       "      *__Pyx_BufPtrStrided1d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_l_constraints.rcbuffer-&gt;pybuffer.buf, __pyx_t_56, __pyx_pybuffernd_l_constraints.diminfo[0].strides) -= __pyx_t_11;\n",
       "</pre><pre class=\"cython line score-0\">&#xA0;<span class=\"\">90</span>:         <span class=\"k\">else</span><span class=\"p\">:</span></pre>\n",
       "<pre class=\"cython line score-0\">&#xA0;<span class=\"\">91</span>:             <span class=\"c\"># TODO: Assuming N=1 (fixed for dependent component)</span></pre>\n",
       "<pre class=\"cython line score-0\" onclick='toggleDiv(this)'>+<span class=\"\">92</span>:             <span class=\"n\">l_constraints</span><span class=\"p\">[</span><span class=\"n\">constraint_offset</span><span class=\"p\">]</span> <span class=\"o\">-=</span> <span class=\"p\">(</span><span class=\"mf\">1</span> <span class=\"o\">-</span> <span class=\"n\">indep_sum</span><span class=\"p\">)</span></pre>\n",
       "<pre class='cython code score-0 '>    <span class='trace'>__Pyx_TraceLine</span>(92,0,__PYX_ERR(0, 92, __pyx_L1_error))\n",
       "    /*else*/ {\n",
       "      __pyx_t_57 = __pyx_v_constraint_offset;\n",
       "      *__Pyx_BufPtrStrided1d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_l_constraints.rcbuffer-&gt;pybuffer.buf, __pyx_t_57, __pyx_pybuffernd_l_constraints.diminfo[0].strides) -= (1.0 - __pyx_v_indep_sum);\n",
       "    }\n",
       "    __pyx_L48:;\n",
       "</pre><pre class=\"cython line score-0\" onclick='toggleDiv(this)'>+<span class=\"\">93</span>:         <span class=\"n\">constraint_offset</span> <span class=\"o\">+=</span> <span class=\"mf\">1</span></pre>\n",
       "<pre class='cython code score-0 '>    <span class='trace'>__Pyx_TraceLine</span>(93,0,__PYX_ERR(0, 93, __pyx_L1_error))\n",
       "    __pyx_v_constraint_offset = (__pyx_v_constraint_offset + 1);\n",
       "</pre><pre class=\"cython line score-12\" onclick='toggleDiv(this)'>+<span class=\"\">94</span>:     <span class=\"k\">return</span> <span class=\"n\">l_constraints</span><span class=\"p\">,</span> <span class=\"n\">constraint_jac</span><span class=\"p\">,</span> <span class=\"n\">constraint_hess</span></pre>\n",
       "<pre class='cython code score-12 '>  <span class='trace'>__Pyx_TraceLine</span>(94,0,__PYX_ERR(0, 94, __pyx_L1_error))\n",
       "  <span class='pyx_macro_api'>__Pyx_XDECREF</span>(__pyx_r);\n",
       "  __pyx_t_1 = <span class='py_c_api'>PyTuple_New</span>(3); if (unlikely(!__pyx_t_1)) __PYX_ERR(0, 94, __pyx_L1_error)\n",
       "  <span class='refnanny'>__Pyx_GOTREF</span>(__pyx_t_1);\n",
       "  <span class='pyx_macro_api'>__Pyx_INCREF</span>(((PyObject *)__pyx_v_l_constraints));\n",
       "  <span class='refnanny'>__Pyx_GIVEREF</span>(((PyObject *)__pyx_v_l_constraints));\n",
       "  <span class='py_macro_api'>PyTuple_SET_ITEM</span>(__pyx_t_1, 0, ((PyObject *)__pyx_v_l_constraints));\n",
       "  <span class='pyx_macro_api'>__Pyx_INCREF</span>(((PyObject *)__pyx_v_constraint_jac));\n",
       "  <span class='refnanny'>__Pyx_GIVEREF</span>(((PyObject *)__pyx_v_constraint_jac));\n",
       "  <span class='py_macro_api'>PyTuple_SET_ITEM</span>(__pyx_t_1, 1, ((PyObject *)__pyx_v_constraint_jac));\n",
       "  <span class='pyx_macro_api'>__Pyx_INCREF</span>(((PyObject *)__pyx_v_constraint_hess));\n",
       "  <span class='refnanny'>__Pyx_GIVEREF</span>(((PyObject *)__pyx_v_constraint_hess));\n",
       "  <span class='py_macro_api'>PyTuple_SET_ITEM</span>(__pyx_t_1, 2, ((PyObject *)__pyx_v_constraint_hess));\n",
       "  __pyx_r = __pyx_t_1;\n",
       "  __pyx_t_1 = 0;\n",
       "  goto __pyx_L0;\n",
       "</pre></div></body></html>"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "execution_count": 29,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "%%cython -a -f --compile-args=-DCYTHON_TRACE=1\n",
    "cimport numpy as np\n",
    "import numpy as np\n",
    "from collections import defaultdict, OrderedDict\n",
    "import copy\n",
    "import itertools\n",
    "import scipy.spatial\n",
    "from pycalphad.core.sympydiff_utils import build_functions as compiled_build_functions\n",
    "from pycalphad.core.constants import MIN_SITE_FRACTION, COMP_DIFFERENCE_TOL\n",
    "from pycalphad.core.phase_rec cimport PhaseRecord, obj, grad, hess, mass_obj, mass_grad, mass_hess\n",
    "from pycalphad.core.cache import cacheit\n",
    "import pycalphad.variables as v\n",
    "from pycalphad.core.eqsolver import _compute_phase_dof\n",
    "cimport cython\n",
    "\n",
    "@cython.boundscheck(False)\n",
    "@cython.wraparound(False)\n",
    "def _compute_constraints(object dbf, object comps, object phases,\n",
    "                         object cur_conds, np.ndarray[dtype=np.float64_t, ndim=1] site_fracs,\n",
    "                         np.ndarray[dtype=np.float64_t, ndim=1] phase_fracs, object phase_records):\n",
    "    \"\"\"\n",
    "    Compute the constraint vector and constraint Jacobian matrix.\n",
    "    \"\"\"\n",
    "    cdef int num_sitefrac_bals = sum([len(dbf.phases[i].sublattices) for i in phases])\n",
    "    cdef int num_mass_bals = len([i for i in cur_conds.keys() if i.startswith('X_')]) + 1\n",
    "    cdef double indep_sum = sum([float(val) for i, val in cur_conds.items() if i.startswith('X_')])\n",
    "    cdef double[::1] comp_obj_value = np.atleast_1d(np.zeros(1))\n",
    "    cdef object dependent_comp = set(comps) - set([i[2:] for i in cur_conds.keys() if i.startswith('X_')]) - {'VA'}\n",
    "    dependent_comp = list(dependent_comp)[0]\n",
    "    cdef int num_constraints = num_sitefrac_bals + num_mass_bals\n",
    "    cdef int num_phases = len(phases)\n",
    "    cdef int num_vars = len(site_fracs) + num_phases\n",
    "    cdef np.ndarray[ndim=1, dtype=np.int_t] phase_dof = _compute_phase_dof(dbf, comps, phases)\n",
    "    cdef np.ndarray[ndim=1, dtype=np.float64_t] l_constraints = np.zeros(num_constraints)\n",
    "    cdef np.ndarray[ndim=2, dtype=np.float64_t] constraint_jac = np.zeros((num_constraints, num_vars))\n",
    "    cdef np.ndarray[ndim=3, dtype=np.float64_t] constraint_hess = np.zeros((num_constraints, num_vars, num_vars), order='F')\n",
    "    cdef double[::1] sfview\n",
    "    cdef double[::1] comp_grad_value\n",
    "    cdef double[::1,:] comp_hess_value\n",
    "    cdef int phase_idx, var_offset, constraint_offset, var_idx, iter_idx, grad_idx, hess_idx, comp_idx, idx\n",
    "    cdef PhaseRecord prn\n",
    "    cdef double phase_frac\n",
    "\n",
    "    # Ordering of constraints by row: sitefrac bal of each phase, then component mass balance\n",
    "    # Ordering of constraints by column: site fractions of each phase, then phase fractions\n",
    "    # First: Site fraction balance constraints\n",
    "    var_idx = 0\n",
    "    constraint_offset = 0\n",
    "    for phase_idx in range(num_phases):\n",
    "        name = phases[phase_idx]\n",
    "        for idx in range(len(dbf.phases[name].sublattices)):\n",
    "            active_in_subl = set(dbf.phases[name].constituents[idx]).intersection(comps)\n",
    "            ais_len = len(active_in_subl)\n",
    "            constraint_jac[constraint_offset + idx,\n",
    "            var_idx:var_idx + ais_len] = 1\n",
    "            l_constraints[constraint_offset + idx] = \\\n",
    "                (sum(site_fracs[var_idx:var_idx + ais_len]) - 1)\n",
    "            var_idx += ais_len\n",
    "        constraint_offset += len(dbf.phases[name].sublattices)\n",
    "    # Second: Mass balance of each component\n",
    "    for comp_idx, comp in enumerate(comps):\n",
    "        if comp == 'VA':\n",
    "            continue\n",
    "        var_offset = 0\n",
    "        phase_idx = 0\n",
    "        for phase_idx in range(num_phases):\n",
    "            name = phases[phase_idx]\n",
    "            prn = phase_records[name]\n",
    "            phase_frac = phase_fracs[phase_idx]\n",
    "            spidx = len(site_fracs) + phase_idx\n",
    "            sfview = site_fracs[var_offset:var_offset + phase_dof[phase_idx]]\n",
    "            comp_hess_value = np.zeros((phase_dof[phase_idx], phase_dof[phase_idx]), order='F')\n",
    "            comp_grad_value = np.zeros(phase_dof[phase_idx])\n",
    "            with nogil:\n",
    "                mass_obj(prn, comp_obj_value, sfview, comp_idx)\n",
    "                mass_grad(prn, comp_grad_value, sfview, comp_idx)\n",
    "                mass_hess(prn, comp_hess_value, sfview, comp_idx)\n",
    "                # current phase frac times the comp_grad\n",
    "                for grad_idx in range(var_offset, var_offset + phase_dof[phase_idx]):\n",
    "                    constraint_jac[constraint_offset, grad_idx] = \\\n",
    "                        phase_frac * comp_grad_value[grad_idx - var_offset]\n",
    "                    constraint_hess[constraint_offset, spidx, grad_idx] = comp_grad_value[grad_idx - var_offset]\n",
    "                    constraint_hess[constraint_offset, grad_idx, spidx] = comp_grad_value[grad_idx - var_offset]\n",
    "                    for hess_idx in range(var_offset, var_offset + phase_dof[phase_idx]):\n",
    "                        constraint_hess[constraint_offset, grad_idx, hess_idx] = phase_frac * comp_hess_value[grad_idx - var_offset, hess_idx - var_offset]\n",
    "                l_constraints[constraint_offset] += phase_frac * comp_obj_value[0]\n",
    "                constraint_jac[constraint_offset, spidx] += comp_obj_value[0]\n",
    "                var_offset += phase_dof[phase_idx]\n",
    "        if comp != dependent_comp:\n",
    "            l_constraints[constraint_offset] -= float(cur_conds['X_' + comp])\n",
    "        else:\n",
    "            # TODO: Assuming N=1 (fixed for dependent component)\n",
    "            l_constraints[constraint_offset] -= (1 - indep_sum)\n",
    "        constraint_offset += 1\n",
    "    return l_constraints, constraint_jac, constraint_hess"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {
    "collapsed": false,
    "scrolled": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "The line_profiler extension is already loaded. To reload it, use:\n",
      "  %reload_ext line_profiler\n"
     ]
    }
   ],
   "source": [
    "%load_ext line_profiler\n",
    "%lprun -f _compute_constraints equilibrium(dbf, ['AL', 'FE', 'VA'], ['B2_BCC', 'AL13FE4'], {v.T: 700, v.X('AL'): (0,1,0.02), v.P: 101325}, model=models, compute_constraints=_compute_constraints)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "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.5.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}