SymEngineTest.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": false,
    "scrolled": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "-(4.15725*B2_BCC2VA*T*Piecewise((-1.00000000000000, 0.5*B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*Piecewise((2.22000000000000, And(T < 6000.0, T >= 298.15)), (0, True)) <= 0), (1.00000000000000, True))*Piecewise((2.22000000000000, And(T < 6000.0, T >= 298.15)), (0, True))*Piecewise((B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(2*B2_BCC0AL - B2_BCC0FE + 2*B2_BCC1AL - B2_BCC1FE)*(0.0028689064946125*T**15/(B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True)))*Piecewise((-1.00000000000000, B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True))) <= 0), (1.00000000000000, True)) + 0.01)**16 + 0.00765041731896667*T**9/(B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True)))*Piecewise((-1.00000000000000, B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True))) <= 0), (1.00000000000000, True)) + 0.01)**10 + 0.05737812989225*T**3/(B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True)))*Piecewise((-1.00000000000000, B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True))) <= 0), (1.00000000000000, True)) + 0.01)**4 - 0.113162422843049/T)*Piecewise((-1.00000000000000, B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True))) <= 0), (1.00000000000000, True))*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)), T/(B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True)))*Piecewise((-1.00000000000000, B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True))) <= 0), (1.00000000000000, True)) + 0.01) < 1), (-B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True)))*Piecewise((-1.00000000000000, B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True))) <= 0), (1.00000000000000, True)) + 0.01)**4*(0.0401082005013337 + 0.00381982861917464*(B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True)))*Piecewise((-1.00000000000000, B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True))) <= 0), (1.00000000000000, True)) + 0.01)**10/T**10 + 0.00133694001671112*(B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True)))*Piecewise((-1.00000000000000, B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True))) <= 0), (1.00000000000000, True)) + 0.01)**20/T**20)*(2*B2_BCC0AL - B2_BCC0FE + 2*B2_BCC1AL - B2_BCC1FE)*Piecewise((-1.00000000000000, B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True))) <= 0), (1.00000000000000, True))*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True))/T**5, True))/(0.5*B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*Piecewise((-1.00000000000000, 0.5*B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*Piecewise((2.22000000000000, And(T < 6000.0, T >= 298.15)), (0, True)) <= 0), (1.00000000000000, True))*Piecewise((2.22000000000000, And(T < 6000.0, T >= 298.15)), (0, True)) + 1) + B2_BCC2VA*(0.125*(B2_BCC0AL + B2_BCC1AL)*Piecewise((2945.20000000000, And(T < 6000.0, T >= 298.15)), (0, True)) - 0.125*(B2_BCC0FE + B2_BCC1FE)*Piecewise((2945.20000000000, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.125*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((2945.20000000000, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.25*Piecewise((31.9888*T - 122960.0, And(T < 6000.0, T >= 298.15)), (0, True))) + 0.5*B2_BCC2VA*Piecewise((-14300.9572000000, And(T < 6000.0, T >= 298.15)), (0, True)) + 8.3145*T*Piecewise((B2_BCC2VA*(-0.0459025039138*B2_BCC2VA*T**15*(B2_BCC0FE + B2_BCC1FE)*(-0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0FE + B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True)))*(2*B2_BCC0AL - B2_BCC0FE + 2*B2_BCC1AL - B2_BCC1FE)/(B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True)))*Piecewise((-1.00000000000000, B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True))) <= 0), (1.00000000000000, True)) + 0.01)**17 - 0.0765041731896667*B2_BCC2VA*T**9*(B2_BCC0FE + B2_BCC1FE)*(-0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0FE + B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True)))*(2*B2_BCC0AL - B2_BCC0FE + 2*B2_BCC1AL - B2_BCC1FE)/(B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True)))*Piecewise((-1.00000000000000, B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True))) <= 0), (1.00000000000000, True)) + 0.01)**11 - 0.229512519569*B2_BCC2VA*T**3*(B2_BCC0FE + B2_BCC1FE)*(-0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0FE + B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True)))*(2*B2_BCC0AL - B2_BCC0FE + 2*B2_BCC1AL - B2_BCC1FE)/(B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True)))*Piecewise((-1.00000000000000, B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True))) <= 0), (1.00000000000000, True)) + 0.01)**5 + 0.005737812989225*T**15*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((-1.00000000000000, B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True))) <= 0), (1.00000000000000, True))/(B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True)))*Piecewise((-1.00000000000000, B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True))) <= 0), (1.00000000000000, True)) + 0.01)**16 + 0.0153008346379333*T**9*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((-1.00000000000000, B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True))) <= 0), (1.00000000000000, True))/(B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True)))*Piecewise((-1.00000000000000, B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True))) <= 0), (1.00000000000000, True)) + 0.01)**10 + 0.1147562597845*T**3*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((-1.00000000000000, B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True))) <= 0), (1.00000000000000, True))/(B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True)))*Piecewise((-1.00000000000000, B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True))) <= 0), (1.00000000000000, True)) + 0.01)**4 - 0.226324845686097*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((-1.00000000000000, B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True))) <= 0), (1.00000000000000, True))/T)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)), T/(B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True)))*Piecewise((-1.00000000000000, B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True))) <= 0), (1.00000000000000, True)) + 0.01) < 1), (-B2_BCC2VA*(B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True)))*Piecewise((-1.00000000000000, B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True))) <= 0), (1.00000000000000, True)) + 0.01)**3*(0.160432802005335*B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(-0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0FE + B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True)))*(2*B2_BCC0AL - B2_BCC0FE + 2*B2_BCC1AL - B2_BCC1FE) + 0.0534776006684449*B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True)))*Piecewise((-1.00000000000000, B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True))) <= 0), (1.00000000000000, True)) + 0.01)**10*(-0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0FE + B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True)))*(2*B2_BCC0AL - B2_BCC0FE + 2*B2_BCC1AL - B2_BCC1FE)/T**10 + 0.0320865604010669*B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True)))*Piecewise((-1.00000000000000, B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True))) <= 0), (1.00000000000000, True)) + 0.01)**20*(-0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0FE + B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True)))*(2*B2_BCC0AL - B2_BCC0FE + 2*B2_BCC1AL - B2_BCC1FE)/T**20 + 0.0802164010026674*(B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True)))*Piecewise((-1.00000000000000, B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True))) <= 0), (1.00000000000000, True)) + 0.01)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((-1.00000000000000, B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True))) <= 0), (1.00000000000000, True)) + 0.00763965723834927*(B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True)))*Piecewise((-1.00000000000000, B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True))) <= 0), (1.00000000000000, True)) + 0.01)**11*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((-1.00000000000000, B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True))) <= 0), (1.00000000000000, True))/T**10 + 0.00267388003342225*(B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True)))*Piecewise((-1.00000000000000, B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True))) <= 0), (1.00000000000000, True)) + 0.01)**21*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((-1.00000000000000, B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*(0.125*(B2_BCC0AL + B2_BCC1AL)*(B2_BCC0AL - B2_BCC0FE + B2_BCC1AL - B2_BCC1FE)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True)) + 0.5*Piecewise((1043, And(T < 6000.0, T >= 298.15)), (0, True))) <= 0), (1.00000000000000, True))/T**20)*Piecewise((504, And(T < 6000.0, T >= 298.15)), (0, True))/T**5, True))*log(0.5*B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*Piecewise((-1.00000000000000, 0.5*B2_BCC2VA*(B2_BCC0FE + B2_BCC1FE)*Piecewise((2.22000000000000, And(T < 6000.0, T >= 298.15)), (0, True)) <= 0), (1.00000000000000, True))*Piecewise((2.22000000000000, And(T < 6000.0, T >= 298.15)), (0, True)) + 1))/(3.0*B2_BCC2VA - 4.0)\n"
     ]
    },
    {
     "ename": "KeyboardInterrupt",
     "evalue": "",
     "output_type": "error",
     "traceback": [
      "\u001b[1;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[1;31mKeyboardInterrupt\u001b[0m                         Traceback (most recent call last)",
      "\u001b[1;32m<ipython-input-1-a6777b429aee>\u001b[0m in \u001b[0;36m<module>\u001b[1;34m()\u001b[0m\n\u001b[0;32m     19\u001b[0m \u001b[0mt\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mt\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mreplace\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;32mlambda\u001b[0m \u001b[0mexpr\u001b[0m\u001b[1;33m:\u001b[0m \u001b[0misinstance\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mexpr\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mPiecewise\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mfunc\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m     20\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0m_\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mrange\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;36m5\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m---> 21\u001b[1;33m     \u001b[0mt\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mt\u001b[0m \u001b[1;33m+\u001b[0m \u001b[0mt\u001b[0m\u001b[1;33m**\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mt\u001b[0m\u001b[1;33m+\u001b[0m\u001b[1;36m3\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m     22\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m     23\u001b[0m \u001b[1;32mfrom\u001b[0m \u001b[0mtimeit\u001b[0m \u001b[1;32mimport\u001b[0m \u001b[0mdefault_timer\u001b[0m \u001b[1;32mas\u001b[0m \u001b[0mclock\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m/home/rotis/anaconda/envs/calphadpy3/lib/python3.5/site-packages/sympy/core/decorators.py\u001b[0m in \u001b[0;36m__sympifyit_wrapper\u001b[1;34m(a, b)\u001b[0m\n\u001b[0;32m     75\u001b[0m                 \u001b[1;32mif\u001b[0m \u001b[1;32mnot\u001b[0m \u001b[0mhasattr\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mb\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;34m'_op_priority'\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m     76\u001b[0m                     \u001b[0mb\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0msympify\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mb\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mstrict\u001b[0m\u001b[1;33m=\u001b[0m\u001b[1;32mTrue\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m---> 77\u001b[1;33m                 \u001b[1;32mreturn\u001b[0m \u001b[0mfunc\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mb\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m     78\u001b[0m             \u001b[1;32mexcept\u001b[0m \u001b[0mSympifyError\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m     79\u001b[0m                 \u001b[1;32mreturn\u001b[0m \u001b[0mretval\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m/home/rotis/anaconda/envs/calphadpy3/lib/python3.5/site-packages/sympy/core/decorators.py\u001b[0m in \u001b[0;36mbinary_op_wrapper\u001b[1;34m(self, other)\u001b[0m\n\u001b[0;32m    116\u001b[0m                     \u001b[1;32melse\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m    117\u001b[0m                         \u001b[1;32mreturn\u001b[0m \u001b[0mf\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 118\u001b[1;33m             \u001b[1;32mreturn\u001b[0m \u001b[0mfunc\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mother\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m    119\u001b[0m         \u001b[1;32mreturn\u001b[0m \u001b[0mbinary_op_wrapper\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m    120\u001b[0m     \u001b[1;32mreturn\u001b[0m \u001b[0mpriority_decorator\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m/home/rotis/anaconda/envs/calphadpy3/lib/python3.5/site-packages/sympy/core/expr.py\u001b[0m in \u001b[0;36m__pow__\u001b[1;34m(self, other)\u001b[0m\n\u001b[0;32m    149\u001b[0m     \u001b[1;33m@\u001b[0m\u001b[0mcall_highest_priority\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;34m'__rpow__'\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m    150\u001b[0m     \u001b[1;32mdef\u001b[0m \u001b[0m__pow__\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mother\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 151\u001b[1;33m         \u001b[1;32mreturn\u001b[0m \u001b[0mPow\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mother\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m    152\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m    153\u001b[0m     \u001b[1;33m@\u001b[0m\u001b[0m_sympifyit\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;34m'other'\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mNotImplemented\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m/home/rotis/anaconda/envs/calphadpy3/lib/python3.5/site-packages/sympy/core/power.py\u001b[0m in \u001b[0;36m__new__\u001b[1;34m(cls, b, e, evaluate)\u001b[0m\n\u001b[0;32m    192\u001b[0m                 \u001b[1;32mif\u001b[0m \u001b[1;32mnot\u001b[0m \u001b[0me\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mis_Atom\u001b[0m \u001b[1;32mand\u001b[0m \u001b[0mb\u001b[0m \u001b[1;32mis\u001b[0m \u001b[1;32mnot\u001b[0m \u001b[0mS\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mExp1\u001b[0m \u001b[1;32mand\u001b[0m \u001b[0mb\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mfunc\u001b[0m \u001b[1;32mis\u001b[0m \u001b[1;32mnot\u001b[0m \u001b[0mexp_polar\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m    193\u001b[0m                     \u001b[1;32mfrom\u001b[0m \u001b[0msympy\u001b[0m \u001b[1;32mimport\u001b[0m \u001b[0mnumer\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mdenom\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mlog\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0msign\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mim\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mfactor_terms\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 194\u001b[1;33m                     \u001b[0mc\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mex\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mfactor_terms\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0me\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0msign\u001b[0m\u001b[1;33m=\u001b[0m\u001b[1;32mFalse\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mas_coeff_Mul\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m    195\u001b[0m                     \u001b[0mden\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mdenom\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mex\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m    196\u001b[0m                     \u001b[1;32mif\u001b[0m \u001b[0mden\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mfunc\u001b[0m \u001b[1;32mis\u001b[0m \u001b[0mlog\u001b[0m \u001b[1;32mand\u001b[0m \u001b[0mden\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0margs\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;36m0\u001b[0m\u001b[1;33m]\u001b[0m \u001b[1;33m==\u001b[0m \u001b[0mb\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m/home/rotis/anaconda/envs/calphadpy3/lib/python3.5/site-packages/sympy/core/exprtools.py\u001b[0m in \u001b[0;36mfactor_terms\u001b[1;34m(expr, radical, clear, fraction, sign)\u001b[0m\n\u001b[0;32m   1186\u001b[0m         \u001b[1;32mreturn\u001b[0m \u001b[0mrv\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1187\u001b[0m     \u001b[0mexpr\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0msympify\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mexpr\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m-> 1188\u001b[1;33m     \u001b[1;32mreturn\u001b[0m \u001b[0mdo\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mexpr\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m   1189\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1190\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m/home/rotis/anaconda/envs/calphadpy3/lib/python3.5/site-packages/sympy/core/exprtools.py\u001b[0m in \u001b[0;36mdo\u001b[1;34m(expr)\u001b[0m\n\u001b[0;32m   1162\u001b[0m         \u001b[0mcont\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mp\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mexpr\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mas_content_primitive\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mradical\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mradical\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mclear\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mclear\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1163\u001b[0m         \u001b[1;32mif\u001b[0m \u001b[0mp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mis_Add\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m-> 1164\u001b[1;33m             \u001b[0mlist_args\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;33m[\u001b[0m\u001b[0mdo\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0ma\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mAdd\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mmake_args\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mp\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m   1165\u001b[0m             \u001b[1;31m# get a common negative (if there) which gcd_terms does not remove\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1166\u001b[0m             \u001b[1;32mif\u001b[0m \u001b[0mall\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mas_coeff_Mul\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;36m0\u001b[0m\u001b[1;33m]\u001b[0m \u001b[1;33m<\u001b[0m \u001b[1;36m0\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0ma\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mlist_args\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m/home/rotis/anaconda/envs/calphadpy3/lib/python3.5/site-packages/sympy/core/exprtools.py\u001b[0m in \u001b[0;36m<listcomp>\u001b[1;34m(.0)\u001b[0m\n\u001b[0;32m   1162\u001b[0m         \u001b[0mcont\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mp\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mexpr\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mas_content_primitive\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mradical\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mradical\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mclear\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mclear\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1163\u001b[0m         \u001b[1;32mif\u001b[0m \u001b[0mp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mis_Add\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m-> 1164\u001b[1;33m             \u001b[0mlist_args\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;33m[\u001b[0m\u001b[0mdo\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0ma\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mAdd\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mmake_args\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mp\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m   1165\u001b[0m             \u001b[1;31m# get a common negative (if there) which gcd_terms does not remove\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1166\u001b[0m             \u001b[1;32mif\u001b[0m \u001b[0mall\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mas_coeff_Mul\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;36m0\u001b[0m\u001b[1;33m]\u001b[0m \u001b[1;33m<\u001b[0m \u001b[1;36m0\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0ma\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mlist_args\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m/home/rotis/anaconda/envs/calphadpy3/lib/python3.5/site-packages/sympy/core/exprtools.py\u001b[0m in \u001b[0;36mdo\u001b[1;34m(expr)\u001b[0m\n\u001b[0;32m   1155\u001b[0m                 \u001b[0mis_iterable\u001b[0m \u001b[1;32mor\u001b[0m \u001b[1;32mnot\u001b[0m \u001b[0mhasattr\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mexpr\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;34m'args_cnc'\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1156\u001b[0m             \u001b[0margs\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mexpr\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0margs\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m-> 1157\u001b[1;33m             \u001b[0mnewargs\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mtuple\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mdo\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mi\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0mi\u001b[0m \u001b[1;32min\u001b[0m \u001b[0margs\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m   1158\u001b[0m             \u001b[1;32mif\u001b[0m \u001b[0mnewargs\u001b[0m \u001b[1;33m==\u001b[0m \u001b[0margs\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1159\u001b[0m                 \u001b[1;32mreturn\u001b[0m \u001b[0mexpr\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m/home/rotis/anaconda/envs/calphadpy3/lib/python3.5/site-packages/sympy/core/exprtools.py\u001b[0m in \u001b[0;36m<listcomp>\u001b[1;34m(.0)\u001b[0m\n\u001b[0;32m   1155\u001b[0m                 \u001b[0mis_iterable\u001b[0m \u001b[1;32mor\u001b[0m \u001b[1;32mnot\u001b[0m \u001b[0mhasattr\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mexpr\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;34m'args_cnc'\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1156\u001b[0m             \u001b[0margs\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mexpr\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0margs\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m-> 1157\u001b[1;33m             \u001b[0mnewargs\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mtuple\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mdo\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mi\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0mi\u001b[0m \u001b[1;32min\u001b[0m \u001b[0margs\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m   1158\u001b[0m             \u001b[1;32mif\u001b[0m \u001b[0mnewargs\u001b[0m \u001b[1;33m==\u001b[0m \u001b[0margs\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1159\u001b[0m                 \u001b[1;32mreturn\u001b[0m \u001b[0mexpr\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m/home/rotis/anaconda/envs/calphadpy3/lib/python3.5/site-packages/sympy/core/exprtools.py\u001b[0m in \u001b[0;36mdo\u001b[1;34m(expr)\u001b[0m\n\u001b[0;32m   1162\u001b[0m         \u001b[0mcont\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mp\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mexpr\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mas_content_primitive\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mradical\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mradical\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mclear\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mclear\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1163\u001b[0m         \u001b[1;32mif\u001b[0m \u001b[0mp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mis_Add\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m-> 1164\u001b[1;33m             \u001b[0mlist_args\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;33m[\u001b[0m\u001b[0mdo\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0ma\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mAdd\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mmake_args\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mp\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m   1165\u001b[0m             \u001b[1;31m# get a common negative (if there) which gcd_terms does not remove\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1166\u001b[0m             \u001b[1;32mif\u001b[0m \u001b[0mall\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mas_coeff_Mul\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;36m0\u001b[0m\u001b[1;33m]\u001b[0m \u001b[1;33m<\u001b[0m \u001b[1;36m0\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0ma\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mlist_args\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m/home/rotis/anaconda/envs/calphadpy3/lib/python3.5/site-packages/sympy/core/exprtools.py\u001b[0m in \u001b[0;36m<listcomp>\u001b[1;34m(.0)\u001b[0m\n\u001b[0;32m   1162\u001b[0m         \u001b[0mcont\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mp\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mexpr\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mas_content_primitive\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mradical\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mradical\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mclear\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mclear\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1163\u001b[0m         \u001b[1;32mif\u001b[0m \u001b[0mp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mis_Add\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m-> 1164\u001b[1;33m             \u001b[0mlist_args\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;33m[\u001b[0m\u001b[0mdo\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0ma\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mAdd\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mmake_args\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mp\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m   1165\u001b[0m             \u001b[1;31m# get a common negative (if there) which gcd_terms does not remove\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1166\u001b[0m             \u001b[1;32mif\u001b[0m \u001b[0mall\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mas_coeff_Mul\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;36m0\u001b[0m\u001b[1;33m]\u001b[0m \u001b[1;33m<\u001b[0m \u001b[1;36m0\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0ma\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mlist_args\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m/home/rotis/anaconda/envs/calphadpy3/lib/python3.5/site-packages/sympy/core/exprtools.py\u001b[0m in \u001b[0;36mdo\u001b[1;34m(expr)\u001b[0m\n\u001b[0;32m   1155\u001b[0m                 \u001b[0mis_iterable\u001b[0m \u001b[1;32mor\u001b[0m \u001b[1;32mnot\u001b[0m \u001b[0mhasattr\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mexpr\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;34m'args_cnc'\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1156\u001b[0m             \u001b[0margs\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mexpr\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0margs\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m-> 1157\u001b[1;33m             \u001b[0mnewargs\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mtuple\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mdo\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mi\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0mi\u001b[0m \u001b[1;32min\u001b[0m \u001b[0margs\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m   1158\u001b[0m             \u001b[1;32mif\u001b[0m \u001b[0mnewargs\u001b[0m \u001b[1;33m==\u001b[0m \u001b[0margs\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1159\u001b[0m                 \u001b[1;32mreturn\u001b[0m \u001b[0mexpr\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m/home/rotis/anaconda/envs/calphadpy3/lib/python3.5/site-packages/sympy/core/exprtools.py\u001b[0m in \u001b[0;36m<listcomp>\u001b[1;34m(.0)\u001b[0m\n\u001b[0;32m   1155\u001b[0m                 \u001b[0mis_iterable\u001b[0m \u001b[1;32mor\u001b[0m \u001b[1;32mnot\u001b[0m \u001b[0mhasattr\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mexpr\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;34m'args_cnc'\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1156\u001b[0m             \u001b[0margs\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mexpr\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0margs\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m-> 1157\u001b[1;33m             \u001b[0mnewargs\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mtuple\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mdo\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mi\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0mi\u001b[0m \u001b[1;32min\u001b[0m \u001b[0margs\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m   1158\u001b[0m             \u001b[1;32mif\u001b[0m \u001b[0mnewargs\u001b[0m \u001b[1;33m==\u001b[0m \u001b[0margs\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1159\u001b[0m                 \u001b[1;32mreturn\u001b[0m \u001b[0mexpr\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m/home/rotis/anaconda/envs/calphadpy3/lib/python3.5/site-packages/sympy/core/exprtools.py\u001b[0m in \u001b[0;36mdo\u001b[1;34m(expr)\u001b[0m\n\u001b[0;32m   1162\u001b[0m         \u001b[0mcont\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mp\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mexpr\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mas_content_primitive\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mradical\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mradical\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mclear\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mclear\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1163\u001b[0m         \u001b[1;32mif\u001b[0m \u001b[0mp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mis_Add\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m-> 1164\u001b[1;33m             \u001b[0mlist_args\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;33m[\u001b[0m\u001b[0mdo\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0ma\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mAdd\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mmake_args\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mp\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m   1165\u001b[0m             \u001b[1;31m# get a common negative (if there) which gcd_terms does not remove\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1166\u001b[0m             \u001b[1;32mif\u001b[0m \u001b[0mall\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mas_coeff_Mul\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;36m0\u001b[0m\u001b[1;33m]\u001b[0m \u001b[1;33m<\u001b[0m \u001b[1;36m0\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0ma\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mlist_args\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m/home/rotis/anaconda/envs/calphadpy3/lib/python3.5/site-packages/sympy/core/exprtools.py\u001b[0m in \u001b[0;36m<listcomp>\u001b[1;34m(.0)\u001b[0m\n\u001b[0;32m   1162\u001b[0m         \u001b[0mcont\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mp\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mexpr\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mas_content_primitive\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mradical\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mradical\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mclear\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mclear\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1163\u001b[0m         \u001b[1;32mif\u001b[0m \u001b[0mp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mis_Add\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m-> 1164\u001b[1;33m             \u001b[0mlist_args\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;33m[\u001b[0m\u001b[0mdo\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0ma\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mAdd\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mmake_args\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mp\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m   1165\u001b[0m             \u001b[1;31m# get a common negative (if there) which gcd_terms does not remove\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1166\u001b[0m             \u001b[1;32mif\u001b[0m \u001b[0mall\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mas_coeff_Mul\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;36m0\u001b[0m\u001b[1;33m]\u001b[0m \u001b[1;33m<\u001b[0m \u001b[1;36m0\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0ma\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mlist_args\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m/home/rotis/anaconda/envs/calphadpy3/lib/python3.5/site-packages/sympy/core/exprtools.py\u001b[0m in \u001b[0;36mdo\u001b[1;34m(expr)\u001b[0m\n\u001b[0;32m   1155\u001b[0m                 \u001b[0mis_iterable\u001b[0m \u001b[1;32mor\u001b[0m \u001b[1;32mnot\u001b[0m \u001b[0mhasattr\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mexpr\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;34m'args_cnc'\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1156\u001b[0m             \u001b[0margs\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mexpr\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0margs\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m-> 1157\u001b[1;33m             \u001b[0mnewargs\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mtuple\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mdo\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mi\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0mi\u001b[0m \u001b[1;32min\u001b[0m \u001b[0margs\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m   1158\u001b[0m             \u001b[1;32mif\u001b[0m \u001b[0mnewargs\u001b[0m \u001b[1;33m==\u001b[0m \u001b[0margs\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1159\u001b[0m                 \u001b[1;32mreturn\u001b[0m \u001b[0mexpr\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m/home/rotis/anaconda/envs/calphadpy3/lib/python3.5/site-packages/sympy/core/exprtools.py\u001b[0m in \u001b[0;36m<listcomp>\u001b[1;34m(.0)\u001b[0m\n\u001b[0;32m   1155\u001b[0m                 \u001b[0mis_iterable\u001b[0m \u001b[1;32mor\u001b[0m \u001b[1;32mnot\u001b[0m \u001b[0mhasattr\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mexpr\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;34m'args_cnc'\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1156\u001b[0m             \u001b[0margs\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mexpr\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0margs\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m-> 1157\u001b[1;33m             \u001b[0mnewargs\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mtuple\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mdo\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mi\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0mi\u001b[0m \u001b[1;32min\u001b[0m \u001b[0margs\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m   1158\u001b[0m             \u001b[1;32mif\u001b[0m \u001b[0mnewargs\u001b[0m \u001b[1;33m==\u001b[0m \u001b[0margs\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1159\u001b[0m                 \u001b[1;32mreturn\u001b[0m \u001b[0mexpr\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m/home/rotis/anaconda/envs/calphadpy3/lib/python3.5/site-packages/sympy/core/exprtools.py\u001b[0m in \u001b[0;36mdo\u001b[1;34m(expr)\u001b[0m\n\u001b[0;32m   1162\u001b[0m         \u001b[0mcont\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mp\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mexpr\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mas_content_primitive\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mradical\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mradical\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mclear\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mclear\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1163\u001b[0m         \u001b[1;32mif\u001b[0m \u001b[0mp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mis_Add\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m-> 1164\u001b[1;33m             \u001b[0mlist_args\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;33m[\u001b[0m\u001b[0mdo\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0ma\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mAdd\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mmake_args\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mp\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m   1165\u001b[0m             \u001b[1;31m# get a common negative (if there) which gcd_terms does not remove\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1166\u001b[0m             \u001b[1;32mif\u001b[0m \u001b[0mall\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mas_coeff_Mul\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;36m0\u001b[0m\u001b[1;33m]\u001b[0m \u001b[1;33m<\u001b[0m \u001b[1;36m0\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0ma\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mlist_args\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m/home/rotis/anaconda/envs/calphadpy3/lib/python3.5/site-packages/sympy/core/exprtools.py\u001b[0m in \u001b[0;36m<listcomp>\u001b[1;34m(.0)\u001b[0m\n\u001b[0;32m   1162\u001b[0m         \u001b[0mcont\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mp\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mexpr\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mas_content_primitive\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mradical\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mradical\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mclear\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mclear\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1163\u001b[0m         \u001b[1;32mif\u001b[0m \u001b[0mp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mis_Add\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m-> 1164\u001b[1;33m             \u001b[0mlist_args\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;33m[\u001b[0m\u001b[0mdo\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0ma\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mAdd\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mmake_args\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mp\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m   1165\u001b[0m             \u001b[1;31m# get a common negative (if there) which gcd_terms does not remove\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1166\u001b[0m             \u001b[1;32mif\u001b[0m \u001b[0mall\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mas_coeff_Mul\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;36m0\u001b[0m\u001b[1;33m]\u001b[0m \u001b[1;33m<\u001b[0m \u001b[1;36m0\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0ma\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mlist_args\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m/home/rotis/anaconda/envs/calphadpy3/lib/python3.5/site-packages/sympy/core/exprtools.py\u001b[0m in \u001b[0;36mdo\u001b[1;34m(expr)\u001b[0m\n\u001b[0;32m   1182\u001b[0m         \u001b[1;32melif\u001b[0m \u001b[0mp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0margs\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1183\u001b[0m             p = p.func(\n\u001b[1;32m-> 1184\u001b[1;33m                 *[do(a) for a in p.args])\n\u001b[0m\u001b[0;32m   1185\u001b[0m         \u001b[0mrv\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0m_keep_coeff\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mcont\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mp\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mclear\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mclear\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0msign\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0msign\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1186\u001b[0m         \u001b[1;32mreturn\u001b[0m \u001b[0mrv\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m/home/rotis/anaconda/envs/calphadpy3/lib/python3.5/site-packages/sympy/core/exprtools.py\u001b[0m in \u001b[0;36m<listcomp>\u001b[1;34m(.0)\u001b[0m\n\u001b[0;32m   1182\u001b[0m         \u001b[1;32melif\u001b[0m \u001b[0mp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0margs\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1183\u001b[0m             p = p.func(\n\u001b[1;32m-> 1184\u001b[1;33m                 *[do(a) for a in p.args])\n\u001b[0m\u001b[0;32m   1185\u001b[0m         \u001b[0mrv\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0m_keep_coeff\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mcont\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mp\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mclear\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mclear\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0msign\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0msign\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1186\u001b[0m         \u001b[1;32mreturn\u001b[0m \u001b[0mrv\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m/home/rotis/anaconda/envs/calphadpy3/lib/python3.5/site-packages/sympy/core/exprtools.py\u001b[0m in \u001b[0;36mdo\u001b[1;34m(expr)\u001b[0m\n\u001b[0;32m   1162\u001b[0m         \u001b[0mcont\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mp\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mexpr\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mas_content_primitive\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mradical\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mradical\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mclear\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mclear\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1163\u001b[0m         \u001b[1;32mif\u001b[0m \u001b[0mp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mis_Add\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m-> 1164\u001b[1;33m             \u001b[0mlist_args\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;33m[\u001b[0m\u001b[0mdo\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0ma\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mAdd\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mmake_args\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mp\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m   1165\u001b[0m             \u001b[1;31m# get a common negative (if there) which gcd_terms does not remove\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1166\u001b[0m             \u001b[1;32mif\u001b[0m \u001b[0mall\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mas_coeff_Mul\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;36m0\u001b[0m\u001b[1;33m]\u001b[0m \u001b[1;33m<\u001b[0m \u001b[1;36m0\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0ma\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mlist_args\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m/home/rotis/anaconda/envs/calphadpy3/lib/python3.5/site-packages/sympy/core/exprtools.py\u001b[0m in \u001b[0;36m<listcomp>\u001b[1;34m(.0)\u001b[0m\n\u001b[0;32m   1162\u001b[0m         \u001b[0mcont\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mp\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mexpr\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mas_content_primitive\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mradical\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mradical\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mclear\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mclear\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1163\u001b[0m         \u001b[1;32mif\u001b[0m \u001b[0mp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mis_Add\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m-> 1164\u001b[1;33m             \u001b[0mlist_args\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;33m[\u001b[0m\u001b[0mdo\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0ma\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mAdd\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mmake_args\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mp\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m   1165\u001b[0m             \u001b[1;31m# get a common negative (if there) which gcd_terms does not remove\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1166\u001b[0m             \u001b[1;32mif\u001b[0m \u001b[0mall\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mas_coeff_Mul\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;36m0\u001b[0m\u001b[1;33m]\u001b[0m \u001b[1;33m<\u001b[0m \u001b[1;36m0\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0ma\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mlist_args\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m/home/rotis/anaconda/envs/calphadpy3/lib/python3.5/site-packages/sympy/core/exprtools.py\u001b[0m in \u001b[0;36mdo\u001b[1;34m(expr)\u001b[0m\n\u001b[0;32m   1182\u001b[0m         \u001b[1;32melif\u001b[0m \u001b[0mp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0margs\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1183\u001b[0m             p = p.func(\n\u001b[1;32m-> 1184\u001b[1;33m                 *[do(a) for a in p.args])\n\u001b[0m\u001b[0;32m   1185\u001b[0m         \u001b[0mrv\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0m_keep_coeff\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mcont\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mp\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mclear\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mclear\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0msign\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0msign\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1186\u001b[0m         \u001b[1;32mreturn\u001b[0m \u001b[0mrv\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m/home/rotis/anaconda/envs/calphadpy3/lib/python3.5/site-packages/sympy/core/exprtools.py\u001b[0m in \u001b[0;36m<listcomp>\u001b[1;34m(.0)\u001b[0m\n\u001b[0;32m   1182\u001b[0m         \u001b[1;32melif\u001b[0m \u001b[0mp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0margs\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1183\u001b[0m             p = p.func(\n\u001b[1;32m-> 1184\u001b[1;33m                 *[do(a) for a in p.args])\n\u001b[0m\u001b[0;32m   1185\u001b[0m         \u001b[0mrv\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0m_keep_coeff\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mcont\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mp\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mclear\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mclear\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0msign\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0msign\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1186\u001b[0m         \u001b[1;32mreturn\u001b[0m \u001b[0mrv\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m/home/rotis/anaconda/envs/calphadpy3/lib/python3.5/site-packages/sympy/core/exprtools.py\u001b[0m in \u001b[0;36mdo\u001b[1;34m(expr)\u001b[0m\n\u001b[0;32m   1155\u001b[0m                 \u001b[0mis_iterable\u001b[0m \u001b[1;32mor\u001b[0m \u001b[1;32mnot\u001b[0m \u001b[0mhasattr\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mexpr\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;34m'args_cnc'\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1156\u001b[0m             \u001b[0margs\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mexpr\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0margs\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m-> 1157\u001b[1;33m             \u001b[0mnewargs\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mtuple\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mdo\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mi\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0mi\u001b[0m \u001b[1;32min\u001b[0m \u001b[0margs\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m   1158\u001b[0m             \u001b[1;32mif\u001b[0m \u001b[0mnewargs\u001b[0m \u001b[1;33m==\u001b[0m \u001b[0margs\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1159\u001b[0m                 \u001b[1;32mreturn\u001b[0m \u001b[0mexpr\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m/home/rotis/anaconda/envs/calphadpy3/lib/python3.5/site-packages/sympy/core/exprtools.py\u001b[0m in \u001b[0;36m<listcomp>\u001b[1;34m(.0)\u001b[0m\n\u001b[0;32m   1155\u001b[0m                 \u001b[0mis_iterable\u001b[0m \u001b[1;32mor\u001b[0m \u001b[1;32mnot\u001b[0m \u001b[0mhasattr\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mexpr\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;34m'args_cnc'\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1156\u001b[0m             \u001b[0margs\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mexpr\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0margs\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m-> 1157\u001b[1;33m             \u001b[0mnewargs\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mtuple\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mdo\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mi\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0mi\u001b[0m \u001b[1;32min\u001b[0m \u001b[0margs\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m   1158\u001b[0m             \u001b[1;32mif\u001b[0m \u001b[0mnewargs\u001b[0m \u001b[1;33m==\u001b[0m \u001b[0margs\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1159\u001b[0m                 \u001b[1;32mreturn\u001b[0m \u001b[0mexpr\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m/home/rotis/anaconda/envs/calphadpy3/lib/python3.5/site-packages/sympy/core/exprtools.py\u001b[0m in \u001b[0;36mdo\u001b[1;34m(expr)\u001b[0m\n\u001b[0;32m   1182\u001b[0m         \u001b[1;32melif\u001b[0m \u001b[0mp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0margs\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1183\u001b[0m             p = p.func(\n\u001b[1;32m-> 1184\u001b[1;33m                 *[do(a) for a in p.args])\n\u001b[0m\u001b[0;32m   1185\u001b[0m         \u001b[0mrv\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0m_keep_coeff\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mcont\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mp\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mclear\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mclear\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0msign\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0msign\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1186\u001b[0m         \u001b[1;32mreturn\u001b[0m \u001b[0mrv\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m/home/rotis/anaconda/envs/calphadpy3/lib/python3.5/site-packages/sympy/core/exprtools.py\u001b[0m in \u001b[0;36m<listcomp>\u001b[1;34m(.0)\u001b[0m\n\u001b[0;32m   1182\u001b[0m         \u001b[1;32melif\u001b[0m \u001b[0mp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0margs\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1183\u001b[0m             p = p.func(\n\u001b[1;32m-> 1184\u001b[1;33m                 *[do(a) for a in p.args])\n\u001b[0m\u001b[0;32m   1185\u001b[0m         \u001b[0mrv\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0m_keep_coeff\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mcont\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mp\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mclear\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mclear\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0msign\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0msign\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1186\u001b[0m         \u001b[1;32mreturn\u001b[0m \u001b[0mrv\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m/home/rotis/anaconda/envs/calphadpy3/lib/python3.5/site-packages/sympy/core/exprtools.py\u001b[0m in \u001b[0;36mdo\u001b[1;34m(expr)\u001b[0m\n\u001b[0;32m   1162\u001b[0m         \u001b[0mcont\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mp\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mexpr\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mas_content_primitive\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mradical\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mradical\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mclear\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mclear\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1163\u001b[0m         \u001b[1;32mif\u001b[0m \u001b[0mp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mis_Add\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m-> 1164\u001b[1;33m             \u001b[0mlist_args\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;33m[\u001b[0m\u001b[0mdo\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0ma\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mAdd\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mmake_args\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mp\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m   1165\u001b[0m             \u001b[1;31m# get a common negative (if there) which gcd_terms does not remove\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1166\u001b[0m             \u001b[1;32mif\u001b[0m \u001b[0mall\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mas_coeff_Mul\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;36m0\u001b[0m\u001b[1;33m]\u001b[0m \u001b[1;33m<\u001b[0m \u001b[1;36m0\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0ma\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mlist_args\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m/home/rotis/anaconda/envs/calphadpy3/lib/python3.5/site-packages/sympy/core/exprtools.py\u001b[0m in \u001b[0;36m<listcomp>\u001b[1;34m(.0)\u001b[0m\n\u001b[0;32m   1162\u001b[0m         \u001b[0mcont\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mp\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mexpr\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mas_content_primitive\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mradical\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mradical\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mclear\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mclear\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1163\u001b[0m         \u001b[1;32mif\u001b[0m \u001b[0mp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mis_Add\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m-> 1164\u001b[1;33m             \u001b[0mlist_args\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;33m[\u001b[0m\u001b[0mdo\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0ma\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mAdd\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mmake_args\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mp\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m   1165\u001b[0m             \u001b[1;31m# get a common negative (if there) which gcd_terms does not remove\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1166\u001b[0m             \u001b[1;32mif\u001b[0m \u001b[0mall\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mas_coeff_Mul\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;36m0\u001b[0m\u001b[1;33m]\u001b[0m \u001b[1;33m<\u001b[0m \u001b[1;36m0\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0ma\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mlist_args\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m/home/rotis/anaconda/envs/calphadpy3/lib/python3.5/site-packages/sympy/core/exprtools.py\u001b[0m in \u001b[0;36mdo\u001b[1;34m(expr)\u001b[0m\n\u001b[0;32m   1182\u001b[0m         \u001b[1;32melif\u001b[0m \u001b[0mp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0margs\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1183\u001b[0m             p = p.func(\n\u001b[1;32m-> 1184\u001b[1;33m                 *[do(a) for a in p.args])\n\u001b[0m\u001b[0;32m   1185\u001b[0m         \u001b[0mrv\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0m_keep_coeff\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mcont\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mp\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mclear\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mclear\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0msign\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0msign\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1186\u001b[0m         \u001b[1;32mreturn\u001b[0m \u001b[0mrv\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m/home/rotis/anaconda/envs/calphadpy3/lib/python3.5/site-packages/sympy/core/exprtools.py\u001b[0m in \u001b[0;36m<listcomp>\u001b[1;34m(.0)\u001b[0m\n\u001b[0;32m   1182\u001b[0m         \u001b[1;32melif\u001b[0m \u001b[0mp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0margs\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1183\u001b[0m             p = p.func(\n\u001b[1;32m-> 1184\u001b[1;33m                 *[do(a) for a in p.args])\n\u001b[0m\u001b[0;32m   1185\u001b[0m         \u001b[0mrv\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0m_keep_coeff\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mcont\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mp\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mclear\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mclear\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0msign\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0msign\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1186\u001b[0m         \u001b[1;32mreturn\u001b[0m \u001b[0mrv\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m/home/rotis/anaconda/envs/calphadpy3/lib/python3.5/site-packages/sympy/core/exprtools.py\u001b[0m in \u001b[0;36mdo\u001b[1;34m(expr)\u001b[0m\n\u001b[0;32m   1164\u001b[0m             \u001b[0mlist_args\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;33m[\u001b[0m\u001b[0mdo\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0ma\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mAdd\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mmake_args\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mp\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1165\u001b[0m             \u001b[1;31m# get a common negative (if there) which gcd_terms does not remove\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m-> 1166\u001b[1;33m             \u001b[1;32mif\u001b[0m \u001b[0mall\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mas_coeff_Mul\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;36m0\u001b[0m\u001b[1;33m]\u001b[0m \u001b[1;33m<\u001b[0m \u001b[1;36m0\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0ma\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mlist_args\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m   1167\u001b[0m                 \u001b[0mcont\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;33m-\u001b[0m\u001b[0mcont\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1168\u001b[0m                 \u001b[0mlist_args\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;33m[\u001b[0m\u001b[1;33m-\u001b[0m\u001b[0ma\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0ma\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mlist_args\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m/home/rotis/anaconda/envs/calphadpy3/lib/python3.5/site-packages/sympy/core/exprtools.py\u001b[0m in \u001b[0;36m<genexpr>\u001b[1;34m(.0)\u001b[0m\n\u001b[0;32m   1164\u001b[0m             \u001b[0mlist_args\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;33m[\u001b[0m\u001b[0mdo\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0ma\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mAdd\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mmake_args\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mp\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1165\u001b[0m             \u001b[1;31m# get a common negative (if there) which gcd_terms does not remove\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m-> 1166\u001b[1;33m             \u001b[1;32mif\u001b[0m \u001b[0mall\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mas_coeff_Mul\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;36m0\u001b[0m\u001b[1;33m]\u001b[0m \u001b[1;33m<\u001b[0m \u001b[1;36m0\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0ma\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mlist_args\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m   1167\u001b[0m                 \u001b[0mcont\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;33m-\u001b[0m\u001b[0mcont\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1168\u001b[0m                 \u001b[0mlist_args\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;33m[\u001b[0m\u001b[1;33m-\u001b[0m\u001b[0ma\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0ma\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mlist_args\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m/home/rotis/anaconda/envs/calphadpy3/lib/python3.5/site-packages/sympy/core/numbers.py\u001b[0m in \u001b[0;36m__lt__\u001b[1;34m(self, other)\u001b[0m\n\u001b[0;32m   1853\u001b[0m     \u001b[1;32mdef\u001b[0m \u001b[0m__lt__\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mother\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1854\u001b[0m         \u001b[1;32mtry\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m-> 1855\u001b[1;33m             \u001b[0mother\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0m_sympify\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mother\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m   1856\u001b[0m         \u001b[1;32mexcept\u001b[0m \u001b[0mSympifyError\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   1857\u001b[0m             \u001b[1;32mraise\u001b[0m \u001b[0mTypeError\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;34m\"Invalid comparison %s < %s\"\u001b[0m \u001b[1;33m%\u001b[0m \u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mother\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m/home/rotis/anaconda/envs/calphadpy3/lib/python3.5/site-packages/sympy/core/sympify.py\u001b[0m in \u001b[0;36m_sympify\u001b[1;34m(a)\u001b[0m\n\u001b[0;32m    353\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m    354\u001b[0m     \"\"\"\n\u001b[1;32m--> 355\u001b[1;33m     \u001b[1;32mreturn\u001b[0m \u001b[0msympify\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mstrict\u001b[0m\u001b[1;33m=\u001b[0m\u001b[1;32mTrue\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m    356\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m    357\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;32m/home/rotis/anaconda/envs/calphadpy3/lib/python3.5/site-packages/sympy/core/sympify.py\u001b[0m in \u001b[0;36msympify\u001b[1;34m(a, locals, convert_xor, strict, rational, evaluate)\u001b[0m\n\u001b[0;32m    239\u001b[0m         \u001b[1;32mpass\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m    240\u001b[0m     \u001b[1;32mtry\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 241\u001b[1;33m         \u001b[0mcls\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0ma\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0m__class__\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m    242\u001b[0m     \u001b[1;32mexcept\u001b[0m \u001b[0mAttributeError\u001b[0m\u001b[1;33m:\u001b[0m  \u001b[1;31m# a is probably an old-style class object\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m    243\u001b[0m         \u001b[0mcls\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mtype\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;31mKeyboardInterrupt\u001b[0m: "
     ]
    }
   ],
   "source": [
    "import pycalphad\n",
    "from pycalphad.tests.datasets import ALFE_TDB\n",
    "from pycalphad import Database, Model\n",
    "import pycalphad.variables as v\n",
    "from sympy import Piecewise, Function\n",
    "\n",
    "dbf = Database(ALFE_TDB)\n",
    "mod = Model(dbf, ['AL','FE', 'VA'], 'B2_BCC')\n",
    "t = mod.ast.diff(v.Y('B2_BCC', 1, 'AL'), v.Y('B2_BCC', 0, 'FE'))\n",
    "#print(t)\n",
    "def func(x):\n",
    "    f = Function('f')\n",
    "    a = []\n",
    "    for t in x.args:\n",
    "        a.append(t[0])\n",
    "    #print(a)\n",
    "    return f(*a)\n",
    "\n",
    "t = t.replace(lambda expr: isinstance(expr, Piecewise), func)\n",
    "for _ in range(5):\n",
    "    t = t + t**(t+3)\n",
    "\n",
    "from timeit import default_timer as clock\n",
    "from symengine import sympify\n",
    "t1 = clock()\n",
    "p = str(t)\n",
    "t2 = clock()\n",
    "\n",
    "sympy_time = t2-t1\n",
    "print(t2-t1)\n",
    "\n",
    "t1 = clock()\n",
    "t = sympify(t)\n",
    "p = str(t)\n",
    "t2 = clock()\n",
    "\n",
    "print(t2-t1)\n",
    "symengine_time = t2-t1\n",
    "print(sympy_time / symengine_time)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "2.1349420790002114\n",
      "0.11995693399967422\n",
      "17.797571243409816\n"
     ]
    }
   ],
   "source": [
    "import pycalphad\n",
    "from pycalphad.tests.datasets import ALFE_TDB\n",
    "from pycalphad import Database, Model\n",
    "import pycalphad.variables as v\n",
    "from sympy import Piecewise, Function\n",
    "\n",
    "\n",
    "dbf = Database(ALFE_TDB)\n",
    "mod = Model(dbf, ['AL','FE', 'VA'], 'B2_BCC')\n",
    "t = mod.ast\n",
    "\n",
    "def func(x):\n",
    "    f = Function('f')\n",
    "    a = []\n",
    "    for t in x.args:\n",
    "        a.append(t[0])\n",
    "    return f(*a)\n",
    "\n",
    "t = t.replace(lambda expr: isinstance(expr, Piecewise), func)\n",
    "\n",
    "from timeit import default_timer as clock\n",
    "from symengine import sympify, diff\n",
    "t1 = clock()\n",
    "p = str(t.diff(v.Y('B2_BCC', 1, 'AL'), v.Y('B2_BCC', 0, 'FE')))\n",
    "t2 = clock()\n",
    "\n",
    "print(t2-t1)\n",
    "sympy_time = t2-t1\n",
    "\n",
    "t1 = clock()\n",
    "t = sympify(t)\n",
    "p = str(t.diff(v.Y('B2_BCC', 1, 'AL')).diff(v.Y('B2_BCC', 0, 'FE')))\n",
    "t2 = clock()\n",
    "\n",
    "print(t2-t1)\n",
    "\n",
    "symengine_time = t2-t1\n",
    "\n",
    "print(sympy_time/symengine_time)"
   ]
  },
  {
   "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
}