SympyTools.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "from sympy import *\n",
    "import numpy as np"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "IPython console for SymPy 1.0 (Python 3.5.2-64-bit) (ground types: python)\n",
      "\n",
      "These commands were executed:\n",
      ">>> from __future__ import division\n",
      ">>> from sympy import *\n",
      ">>> x, y, z, t = symbols('x y z t')\n",
      ">>> k, m, n = symbols('k m n', integer=True)\n",
      ">>> f, g, h = symbols('f g h', cls=Function)\n",
      ">>> init_printing()\n",
      "\n",
      "Documentation can be found at http://docs.sympy.org/1.0/\n"
     ]
    }
   ],
   "source": [
    "init_session()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "def lambdify2_0(Function):\n",
    "    Function_Num=lambdify(sorted(Function.free_symbols,key=default_sort_key) , Function)\n",
    "    \n",
    "    arguments=sorted(Function.free_symbols,key=default_sort_key)\n",
    "    setattr(Function_Num, 'arguements', arguments)\n",
    "    \n",
    "    def N(vals):\n",
    "        return Function_Num(*vals)\n",
    "    setattr(Function_Num, 'N', N)\n",
    "    \n",
    "    return Function_Num"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAASsAAAAVBAMAAADsjdXuAAAAMFBMVEX///8AAAAAAAAAAAAAAAAA\nAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAv3aB7AAAAD3RSTlMAIma7zZnddlTvRIky\nEKtZsEGBAAAACXBIWXMAAA7EAAAOxAGVKw4bAAADO0lEQVRIDcVWP0wTURj/3fXP0aMHTQguJqaw\nOCmVxITFUBfiJOjiIIYOxjhg2kXRkJAajW6mDi5ON7m2Jho2e1Hj4kAdjHHQ3mDCgAkgVQwK+L13\nf3q89zg2/HK9932/3/ev976+K6DlcZhiDR1YjXsMwCh3qsiWn6j941kYq50LysCF3fMtrHYckRyx\nRURsoC9HHmeBdpHWR5K3D8Sz6GfBKukQeI9V2Cvpwl6bWXtLWEVAnwIyNaJc+iglnkW9oYyC+YeS\nV2XO2JYxocQikM4ByTXAcGRvD4ln0bTVgb1bsJ4rKG1TBoUSn4Ej5NS3DhyVnX0knsXkPoFU6r2S\nmpBRoUTSxk1yoi+2/x4ewIKNkEr6i2w+FDLiSKDQQGoKo+Sjb8bsYZe9o5oik+/JjFyq3miWpPoM\naMtpggb8NHoR18nR/Bmzhz6bdJC2yVeU3r8MaYsw0JxP/ZZRQpZqEhw04KfpKeEp+VgbMXvos+0W\n6JKFTStwXyYm85pyf7OjNDKCBA34abJrvC380m3BMWpy9gHwLAoGen+BaS8Ds7tuAONO1wy1d62d\nUA8VvwE/DbXFNhHlQbrpiwvA4u3QNVQYO3h1Gq+Gv4dYqNQrwMe5c6EdKGzm6m5gdVfzMWa6VqDx\nBnga1oBZ4iOP8TPEL+A0UvabwLO7cpYKT7Sob1GaDqzZrHxEsplLyzDqNpaqYhLeAE/DG9CL3uOf\nsQFtHWNIjH2BJY4EY0FH2xYu++8JsgMpA6mpbE2KStPMGXTQC9nMH9RtQQABVoKnYQ0g1fAG5iQV\n0UvYRnaVvuS1VlDTWxlr1kC1acfTN6Kk9m13BZmq4YpRXyd3Kri1e0rEE7TnPevKEjwNb4D96nN+\nmV6X6t7FMoUFkM+wxXA/GAVtzQG5iNLOJU60lFHMU5VNCfI0vAF684RHcY/b52IW03Toi3XJTlUa\niYpeaAB0CZJxll6ro5ijKpsS5Gl4A/SqBj1nTy69yGP4uA18CpDIai47yZy20oJViaCeai4/vKiO\nYrwqmxLkaVgD7I8NBlgsl4zjrXlvUd91NQzsF6XElWCQOcGGSMvTjeSY6qj2qEO+D0XrvZ27EjX/\nr/4PCTzuDufrpwMAAAAASUVORK5CYII=\n",
      "text/latex": [
       "$$\\left ( V_{gs}, \\quad V_{th}, \\quad I_{d}, \\quad K, \\quad \\lambda, \\quad V_{ds}\\right )$$"
      ],
      "text/plain": [
       "(V_{gs}, V_{th}, I_d, K, \\lambda, V_{ds})"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "Vgs, Vth, Id, K, lam, Vds=symbols('V_{gs}, V_{th}, I_d, K, \\lambda, V_{ds}')\n",
    "Vgs, Vth, Id, K, lam, Vds"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 41,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAIUAAAAqBAMAAABvvoU0AAAAMFBMVEX///8AAAAAAAAAAAAAAAAA\nAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAv3aB7AAAAD3RSTlMAdqvNEDJUuyJEiWaZ\n3e/xv6KKAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAC2klEQVRIDaVUS2sTURg9ec9kJm1aELRUGRU3\nbhoF3Yhk0I0rG3QnPkIRuxKzqOiusxSNNH9AMnZV3KS4qlZxEKXoxlFBqiDtxp1oalrTihq/O7mZ\nmUwTBmc+mnvP+R7n9pv7ADDcOqxirFFAGGtQ8aQWRgHiL0AqhpJAegvxuXASiNVxLqQEBnOSEVaj\nqtfMsBq1EflnYI300gtWm58X2O4GsznsZ4VNYDrwAVvGJQ0Q14CqwsSC2EfMVID0BpBZD1Lfrrmm\nUn0diNBRDWpvgRP5H1nsbt3kEnf+Q2qB5aYNb4VQIr2jwIW/20KUGllpFJFY+cqrohqBJ5w4U0xD\nnG10UnF8LjSeIzLRccSJibpY6HA+M9Vv9It6/JwOrBJQOAHuA2evTKo2b4NlmvL02+Hxc0rXExFn\n3dPA91bLkyqY5BjTIOieAKfRTeC2E4pVHGyjiEGwVunXCntqXK1AZuleS82T51m2XyuQ1ngre3RW\nKuXY6LFMkRxVo18rEH+zVmIFZKwukma7vmUb8QEWGi9178rJz8yOsPR4k7UyrrI/skSdjR6zNAbN\nRx63Q/9ItMpV4J7l6qlh9ZLZNCjj8ZCC9FDRynWGlTJQfv8Ss4sHyCmaTsRG1jdNNVU6BUrKwKtE\nyQ61wfQUzXQXjqvDBHp+U+sCybSBmNFSRdy4rBJ02xtqBXRItvCagKy7Yxxb/5ykEHuKQRVnWgVP\n0kPi4ioSq/hAiDaoY+Wl6x3Izrpls7iIpBbNdrhrjiinIiWhTvV3bW9cYUezbfZNnsA+ZAqys5Cd\nDjmrp7JSSbfuHHdHVVj3kdGUxkYyeeQYpMVdbdI9insLMU04qILdfW4DJqIbHAuljhdfbNQP2AvS\nl1l3NMDfQkFJOov001joCmx71xOHRtWuDH9S1f1z/DLe+SX4x3u+Jf5lXRmjXSwQSRq4FajQVXQe\n2OmiQaDw6cFzM0ihqyZGr2EojX9B3a0eEbvlLAAAAABJRU5ErkJggg==\n",
      "text/latex": [
       "$$\\frac{K}{2} \\left(V_{gs} - V_{th}\\right)^{2}$$"
      ],
      "text/plain": [
       "                   2\n",
       "K⋅(V_{gs} - V_{th}) \n",
       "────────────────────\n",
       "         2          "
      ]
     },
     "execution_count": 41,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "Idrain_sat=(K/2)*(Vgs-Vth)**2; Idrain_sat# *(1-lam*Vds); simplify(Idrain_sat)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 42,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAJAAAAAVBAMAAABbFXrgAAAAMFBMVEX///8AAAAAAAAAAAAAAAAA\nAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAv3aB7AAAAD3RSTlMAdt3NMolEqxBUuyJm\nme9KAxjRAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAB30lEQVQ4Ea2TP0gbURzHv/lzyV2S4+KRKhRK\nM5SQLSm4iEoOFzfJFLvp5tZccdChkFs6J1sLXdJVl9scA4qKopg6dO2iq6BpVdLS9Jd3l9R77+Wm\nPg7y7vv5/T7v/SAHcxb/YR0UUAMijd42kgPeaA7mgKU/Xfk5HK0NRTDKgLbLN6g9SjSbj/13jnqi\nehU4ERtuKEqKsZ8EqSdaW8cHV+yoUJQTYz8JUk9UQZzuJKyGi7Qs9wqD1BP1cC5YKKCLTp6Mo0yU\n+KncyUQX2ZDJMKKR62EvEym31e8tianeDZkMHqX/TGEsiv1GpywRNdshk4FR9RH00GI3MtpI/ZCI\njPyRJB1FjGp5JPLDhInqWSQeRvzJb/RXl94Op2woU+LoQ5q+XLSU98eWL1qjTcmhHm6l+g59P3aq\ni2/xNscARmMtxIoR2xeVqOhlFWiWg9X6Pb133FQLO58cOTVcGOs69dJo6ubgHSKlvotoMSjK0EnY\nh+FgeWDJ6VtgBdGWf6Nx++fx7t9mlSo1N5kFZHRGdZ+jE+dFdEVhfcQWopZuATL6VcEXnCqcSKVj\nhaU/W0DmbBqQ0o0c5rH3hhNlBAsLrrx4AvV7ajBfydv9NG1r5dACD5rsOwkrjL9+4YTxMfsLtvJx\nI/WPjpEAAAAASUVORK5CYII=\n",
      "text/latex": [
       "$$\\left [ K, \\quad V_{gs}, \\quad V_{th}\\right ]$$"
      ],
      "text/plain": [
       "[K, V_{gs}, V_{th}]"
      ]
     },
     "execution_count": 42,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "sorted(Idrain_sat.free_symbols, key=lambda i: default_sort_key(i, order='rev-lex'))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 43,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAJAAAAAVBAMAAABbFXrgAAAAMFBMVEX///8AAAAAAAAAAAAAAAAA\nAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAv3aB7AAAAD3RSTlMAdt3NMolEqxBUuyJm\nme9KAxjRAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAB30lEQVQ4Ea2TP0gbURzHv/lzyV2S4+KRKhRK\nM5SQLSm4iEoOFzfJFLvp5tZccdChkFs6J1sLXdJVl9scA4qKopg6dO2iq6BpVdLS9Jd3l9R77+Wm\nPg7y7vv5/T7v/SAHcxb/YR0UUAMijd42kgPeaA7mgKU/Xfk5HK0NRTDKgLbLN6g9SjSbj/13jnqi\nehU4ERtuKEqKsZ8EqSdaW8cHV+yoUJQTYz8JUk9UQZzuJKyGi7Qs9wqD1BP1cC5YKKCLTp6Mo0yU\n+KncyUQX2ZDJMKKR62EvEym31e8tianeDZkMHqX/TGEsiv1GpywRNdshk4FR9RH00GI3MtpI/ZCI\njPyRJB1FjGp5JPLDhInqWSQeRvzJb/RXl94Op2woU+LoQ5q+XLSU98eWL1qjTcmhHm6l+g59P3aq\ni2/xNscARmMtxIoR2xeVqOhlFWiWg9X6Pb133FQLO58cOTVcGOs69dJo6ubgHSKlvotoMSjK0EnY\nh+FgeWDJ6VtgBdGWf6Nx++fx7t9mlSo1N5kFZHRGdZ+jE+dFdEVhfcQWopZuATL6VcEXnCqcSKVj\nhaU/W0DmbBqQ0o0c5rH3hhNlBAsLrrx4AvV7ajBfydv9NG1r5dACD5rsOwkrjL9+4YTxMfsLtvJx\nI/WPjpEAAAAASUVORK5CYII=\n",
      "text/latex": [
       "$$\\left [ K, \\quad V_{gs}, \\quad V_{th}\\right ]$$"
      ],
      "text/plain": [
       "[K, V_{gs}, V_{th}]"
      ]
     },
     "execution_count": 43,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "list(ordered(Idrain_sat.free_symbols))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 44,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "NIdrain_sat=lambdify2_0(Idrain_sat)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 45,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAJAAAAAVBAMAAABbFXrgAAAAMFBMVEX///8AAAAAAAAAAAAAAAAA\nAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAv3aB7AAAAD3RSTlMAdt3NMolEqxBUuyJm\nme9KAxjRAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAB30lEQVQ4Ea2TP0gbURzHv/lzyV2S4+KRKhRK\nM5SQLSm4iEoOFzfJFLvp5tZccdChkFs6J1sLXdJVl9scA4qKopg6dO2iq6BpVdLS9Jd3l9R77+Wm\nPg7y7vv5/T7v/SAHcxb/YR0UUAMijd42kgPeaA7mgKU/Xfk5HK0NRTDKgLbLN6g9SjSbj/13jnqi\nehU4ERtuKEqKsZ8EqSdaW8cHV+yoUJQTYz8JUk9UQZzuJKyGi7Qs9wqD1BP1cC5YKKCLTp6Mo0yU\n+KncyUQX2ZDJMKKR62EvEym31e8tianeDZkMHqX/TGEsiv1GpywRNdshk4FR9RH00GI3MtpI/ZCI\njPyRJB1FjGp5JPLDhInqWSQeRvzJb/RXl94Op2woU+LoQ5q+XLSU98eWL1qjTcmhHm6l+g59P3aq\ni2/xNscARmMtxIoR2xeVqOhlFWiWg9X6Pb133FQLO58cOTVcGOs69dJo6ubgHSKlvotoMSjK0EnY\nh+FgeWDJ6VtgBdGWf6Nx++fx7t9mlSo1N5kFZHRGdZ+jE+dFdEVhfcQWopZuATL6VcEXnCqcSKVj\nhaU/W0DmbBqQ0o0c5rH3hhNlBAsLrrx4AvV7ajBfydv9NG1r5dACD5rsOwkrjL9+4YTxMfsLtvJx\nI/WPjpEAAAAASUVORK5CYII=\n",
      "text/latex": [
       "$$\\left [ K, \\quad V_{gs}, \\quad V_{th}\\right ]$$"
      ],
      "text/plain": [
       "[K, V_{gs}, V_{th}]"
      ]
     },
     "execution_count": 45,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ArgList=NIdrain_sat.arguements; ArgList"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 46,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "def test(arg):\n",
    "    print(*arg)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 47,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "K V_{gs} V_{th}\n"
     ]
    }
   ],
   "source": [
    "test(ArgList)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 49,
   "metadata": {},
   "outputs": [
    {
     "ename": "SyntaxError",
     "evalue": "invalid syntax (<ipython-input-49-71a28903c861>, line 1)",
     "output_type": "error",
     "traceback": [
      "\u001b[1;36m  File \u001b[1;32m\"<ipython-input-49-71a28903c861>\"\u001b[1;36m, line \u001b[1;32m1\u001b[0m\n\u001b[1;33m    NIdrain_sat2=lambda(K, V_{gs}, V_{th}: K+ V_{ds}+ V_{gs}+ V_{th})\u001b[0m\n\u001b[1;37m                       ^\u001b[0m\n\u001b[1;31mSyntaxError\u001b[0m\u001b[1;31m:\u001b[0m invalid syntax\n"
     ]
    }
   ],
   "source": [
    "NIdrain_sat2=lambda(K, V_{gs}, V_{th}: K+ V_{ds}+ V_{gs}+ V_{th})"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 56,
   "metadata": {},
   "outputs": [],
   "source": [
    "import re\n",
    "#http://stackoverflow.com/questions/875968/how-to-remove-symbols-from-a-string-with-python"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 57,
   "metadata": {},
   "outputs": [
    {
     "ename": "AttributeError",
     "evalue": "module 're' has no attribute 'subs'",
     "output_type": "error",
     "traceback": [
      "\u001b[1;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[1;31mAttributeError\u001b[0m                            Traceback (most recent call last)",
      "\u001b[1;32m<ipython-input-57-bf44336fbd39>\u001b[0m in \u001b[0;36m<module>\u001b[1;34m()\u001b[0m\n\u001b[1;32m----> 1\u001b[1;33m \u001b[0mre\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0msubs\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;34mr'[^\\w]'\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;34m''\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mstr\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mArgList\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;36m1\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m",
      "\u001b[1;31mAttributeError\u001b[0m: module 're' has no attribute 'subs'"
     ]
    }
   ],
   "source": [
    "re.subs(r'[^\\w]', '', str(ArgList[1]))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python [conda root]",
   "language": "python",
   "name": "conda-root-py"
  },
  "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"
  },
  "latex_envs": {
   "LaTeX_envs_menu_present": true,
   "autocomplete": true,
   "bibliofile": "biblio.bib",
   "cite_by": "apalike",
   "current_citInitial": 1,
   "eqLabelWithNumbers": true,
   "eqNumInitial": 1,
   "hotkeys": {
    "equation": "Ctrl-E",
    "itemize": "Ctrl-I"
   },
   "labels_anchors": false,
   "latex_user_defs": false,
   "report_style_numbering": false,
   "user_envs_cfg": false
  },
  "toc": {
   "colors": {
    "hover_highlight": "#DAA520",
    "running_highlight": "#FF0000",
    "selected_highlight": "#FFD700"
   },
   "moveMenuLeft": true,
   "nav_menu": {
    "height": "0px",
    "width": "0px"
   },
   "navigate_menu": true,
   "number_sections": true,
   "sideBar": true,
   "threshold": 4,
   "toc_cell": false,
   "toc_section_display": "block",
   "toc_window_display": false,
   "widenNotebook": false
  },
  "varInspector": {
   "cols": {
    "lenName": 16,
    "lenType": 16,
    "lenVar": 40
   },
   "kernels_config": {
    "python": {
     "delete_cmd_postfix": "",
     "delete_cmd_prefix": "del ",
     "library": "var_list.py",
     "varRefreshCmd": "print(var_dic_list())"
    },
    "r": {
     "delete_cmd_postfix": ") ",
     "delete_cmd_prefix": "rm(",
     "library": "var_list.r",
     "varRefreshCmd": "cat(var_dic_list()) "
    }
   },
   "types_to_exclude": [
    "module",
    "function",
    "builtin_function_or_method",
    "instance",
    "_Feature"
   ],
   "window_display": false
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}