Skip to content

Instantly share code, notes, and snippets.

@jeanpat
Created October 6, 2015 22:16
Show Gist options
  • Select an option

  • Save jeanpat/addfefe5cae9bc9b9e9a to your computer and use it in GitHub Desktop.

Select an option

Save jeanpat/addfefe5cae9bc9b9e9a to your computer and use it in GitHub Desktop.
Getting the vertices of degree 1 and reading the weight of the bound edge using graph-tool (ipython notebook
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"import graph_tool.all as gt\n",
"#from graph_tool import topology"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### Make an undirectional weighted graph"
]
},
{
"cell_type": "code",
"execution_count": 73,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"T = gt.Graph(directed = False)\n",
"edge_weights = T.new_edge_property('double')\n",
"T.properties[(\"e\",\"weight\")] = edge_weights\n",
"\n",
"T.add_vertex(n=5)\n",
"\n",
"e_1 = T.add_edge(0,1)\n",
"e_2 = T.add_edge(1,2)\n",
"e_3 = T.add_edge(0,0)\n",
"e_4 = T.add_edge(1,3)\n",
"\n",
"edge_weights[e_1]= 2\n",
"edge_weights[e_2]= 6\n",
"edge_weights[e_3]= 1\n",
"edge_weights[e_4]= 3\n"
]
},
{
"cell_type": "code",
"execution_count": 74,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"weight (edge) (type: double)\n",
"None\n"
]
}
],
"source": [
"print T.list_properties()"
]
},
{
"cell_type": "code",
"execution_count": 75,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAMgAAADICAYAAACtWK6eAAAABmJLR0QA/wD/AP+gvaeTAAAgAElE\nQVR4nO3deXxd5Xkn8N/znu3ukrUvFraxMYuXGCTLGAx2CE0CxLTTJZluNGy+AruZhrbpMtOM0kkn\nnXSmbQigBU+apDSTkjSQ2IQmTVKHsgTbwuANbLCNba3W1Xr3s7zP/GGTppiALZ2rK+m+37/88ZWe\n85yPPz+fe855F0BRFEVRFEVRFEVRFEVRFEVRFEVRFEVRFEVRFEVRFEVRFEVRFEVRFEVRFEVRFEVR\nFEVRFEVRFEVRFEVRFGWuomI3MNts3759meu6GwEsBlADoAxAioiGmbmfmZ8dGhp6pb29XRa1UWVG\nqIAAePzxx7Xx8fFfklL+NhEthOOa5NhBllKHZI2IJGvCha7nYZpZAGNE9ITjOF/dunVrqtj9K4VT\n8gHp6OhoJaI/BvNikc5UIJcrg+eZIp93KZN1OZP1yDQEh8M6h4IGC+GyZSURCSdYiAQzP3zfffc9\nUezzUAqjpAPS2dn5MQAPUDZbSal0tUhnoB04NKa/cSwpBgfzb/95WV5muJddFnFXXLmAKysMBEOj\nHA0nQPQtAJ+Px+POzJ+FUkglG5DOzs77AdxNk8lakcmU63teSpi7d4/Cdt773oII7uqVZfb162s5\nFs3LsrI+CPFMRUXFJz/60Y96he9emSlasRsohs7OzluJ6AExPr5Qm5gMW9/e2WscODgBT/KF1hBD\nZ/LGseOTbmNjOQlRjkCgLJfPh3bu3PlCIXtXZlbJBWT79u2LpJQPUjJVT8lkLPD4t05q/f25KRXL\n5qT+2tFJb/myMmgiwoHAoo985CPHdu7cecLntpUiEcVuYKa5rruNHKdcZDKV1vd+0CfOnDnvXuNi\nUC7nWU/u6BXpTECk0wuI6BNdXV2GX/0qxVVSAeno6FhJRJuQTNZob745qR993ZdHtCIxYusvvTxC\nmWwVPO8SZv5lP+oqxVdSASGi28i2w8JxA8Yzzw37WVvfvXeUUmmJbLaciG7zs7ZSPKUUECKijcjm\nomJoKCsSI7avxR2HxRtvJJHPRwFc+eijj9b6WV8pjpIJyCOPPLIEQA3ZdkQcP1GQt9/6seNJ4XoB\neJ7ueV5rIY6hzKySCYgQ4uz/6FIaYnhkak+t3usYA0N5ZiZ4nklE6goyD5RMQABUwfN0ZiaRSruF\nOADlch55noSUOjNXF+IYyswqpYBYxHx25IDrXPALwYvmuoyzxwkU7BjKjCmZgBDRiBTCAwAOh/VC\nHIN1nWBZGmmaS0RnCnEMZWaVTEAAJEgISUJ4XB4ryIs8Li8zGAATuUQ0UohjKDOrlAJyjJltNvSM\nu2hRuBAH8JYuiZCmOTAMW0r5aiGOocyskglIPB7PENEeWFaSL2mKwDB8H8nsLr00Kk0zxcyjg4OD\nB/yur8y8kgkIABDRDzkQSMpAAHbL1RV+1paXLg7JuroQgsFxItqlpuTODyUVEGZ+mol6EQomvOZr\nKjkS8eVmnYUge8P1NQhYEzCNSV3Xv+xHXaX4Siog52b8PcLh8BiHQk7+F29rZF2f9lct5xc+UOvV\nVJsciQwz8zfvvvvufh/aVWaBkgoIALS1tX0PRD+U5WW9Xn29aW++pWGq9yMMwN5wfZW78qpyLov1\nQdePhEKhLp9bVoqo5CZMAcBNN930rGFZ62GaQQ6FquSSxVFx+lSacvkLv2+wLJH/yC0N7uqV5RyL\nDrBl9Wua1nbXXXepx7vzSMnOSe/q6qoC8NfwvNUYH18ocnbAOHR4VH9xzyglkz9/KIplCmfN6nKn\n+ZoqGQpxyjSkFQ6/rOv61ra2tsMzdwbKTCjZgADAgw8+aFmW9V8B3IpsNiZS6Wp4niGGh3PidF+a\nUimXslmXDUPjaEyX9bVBbmwISyE4bxp2GrBGx8dzz7/Y85Vn//Vffr/Y56P4r6QD8pZHH310ted5\nv0dEq2HbIcrmovDcAHlSZyl1EsKDEK4nhJcDjBzYzOTz7quvHUkceePYOLOUJOhj+3fvVnPR5xkV\nkH9H3d3dLcy8CcBGAHVv/wHP87KHXj2yqLevP3uqrzflut5PBz0y4XsH9+7+rzPYrzIDVEB+jr/7\nu78L2LZdw8wLpJRJy7KG77777uSq5tYHAPzGeb9AkMzerx/s6Tk2890qhaICcpGuaG2tND3+NoPO\nG85O4O/v79nzp8XoSymMknzMOx2Jvr5sTcPCcgCrzvuQaEl13SX/OjzQOzrznSmFUHIvCv0gTe3L\nBD5/2i5DaOTeVYSWlAJRV5ApGO7tzdY2NsUArD7vQ6IllTULdyUG+9RVZB5QV5ApEm7+ywAy533A\nELrGd894Q0pBqIBM0csvvzxO4H96xw8JH3hfa+vyGW6p4Nrb20WpLauqnmJNw5o1a8o9zfwOgNB5\nHzJ+eOCl3X808135o729XdTV1a0ioo0A1jJzNRFVABDMPApglIj2E9Gu/v7+ve3t7b4uxDdbqIBM\n08rmtZ8g0B3nfUBgIfCbr+zefbQIbU0HdXZ2foiI7gPQCNsOUt4Ow3NNMGvETCyEy0JzYZlpmGaG\ngQkAXxocHHx8vgVFBWSa3u0qwoQfHdy7+1NFaGtKOjo6GoUQnwPzCkqnF1A2V8GeZ2hj4zmMjtqU\nybgkmREO6TIWM2RNdRBCSBkIjCMcGoGmnSaiT2/ZsuXlYp+LX1RAfLCyufV3Cfid8z6YQ1eRrq6u\na5j588JxGjAx2SCyWaHv2z+iHT40KcYm3nFrOQ6FNPeK5VF3bXMVR6OCI+EhGQwmhBCf27Jly3dm\n+hwKQT3m9UFTQ90RD+LXALz9BpbAXDE00P8vxejrQnV3d68B8LDIZhtpYrJRP/xqMvDEjl7t+PF3\nnSNDjsPawGBOf+XAGEDEVZW1JKXJltWyefPmMzt37nxtBk+jINQVxCerWlq3gfHxn/07EoSaqqrg\n+2/c8O2qioowzu67XgEgycyjQoheKeULoVDo+TvuuCNdjL67urrqmfkrIptdQpPJeuPfnh8y9uwd\nm0otd9mlEeeWDzXIWGxCxqK9zHz/fffdt8/vnmeSCohPmpuby2xoOwCEhBC04orLFyxftrTSNE09\nrGuTASGykKyT9HSQ8KBpLuuazaaZJiFyAH4AoCMejw/MZN+dnZ3d5DgbxNj4Jcazz50xXtw7rRec\n3tKl4fwv3tokY7FBhEKvp9PpX37ggQeyfvU709RXLJ8MDAzk6xoagk2NDRvef8MNlzQ11EfKheZE\nPZcM2wlqiRFTO3NGaP2DTOPjushmg2Aup2yukhw7Ak1fBE37pc2bNwfuuOOOl77xjW8Ubv3gczo7\nO28koo/T2HiTcfT1tPmjH097uVQxNuYAJGVtTQ0HA45pWbmdO3e+5Ee/xaCuID76m7/5m3tyefuv\nw4ATsG2TUinPePmVUe3o68l3utGVpiHk0qVhb/XKcrexIYKANcHR6BAL8bxt23/yiU98YrKQ/XZ1\ndf0/ymTW0uhYbeBLXz0mUilfVr1nISh7529fKuvq0hyLHvc879atW7cWZE+WQlNXEJ90dXXdZ+h6\nPJjLB61cLqS/8JOEtfOf+7XTvdmfd6NLnmSRSNj6ocOT2tBQxquuXkCgcjasBZppXLt58+bv7ty5\n8x2fIE3Xo48+uoSZ22hissF4+ZVJ/cjrSb9qEzPI8zy56JJqDoWSuq4f3bFjx5ycJ6OGmvigu7v7\nQwDuovHxRm1iImj905MnzRd2j5DrXvDXJO34m5ngP/zjm9qp07YYH1skPO9KAJ9tb28vyL+R67ob\n4TgmXNfSXj3q+5VKf+1ICp7HwrbDUsqNftefKSog09TR0dHIzJ+myVQNpTNh65tPnhKnTk/pppRy\nOc/8xrd6tf7+HMYnmiDlprq6uvNnL/qAiNbAtsMik3W0gQH/d9yyHan19WeQy0eIaI3v9WeICsj0\n3Q/bKUMmXWV+/4f90953XUq2ntzZJyYnBSVT1QDufuyxx2I+9fpTQogaeNKgiYnCDQ0ZG7PBUgdQ\nVagrYaHNyaZni4ceemi5EOKDIpWs0U6dTuqvHfHlezxlMp7x/E/OUC63gFy3Mp1Of9yPuj/L87ga\n0jNEJlOQ7egAQKQzDjxpABAVFRWVhTpOIRVkp6VSYRjGh2E7QbadsPnMc8f9rK3tPzhJ11xdSYHA\nAi4zPgzgizi72ul7WrbsFiscHq+C4TR6LKoFUCXB1QSqAnEjM1UPnRlqrTVMF8C0rnjvhsFEABf8\neXUBqYBMAzNvErl8VBsdy033q9XbEQDt6NFJWV21AEBNZ2fnlW1tbYeXrVsXC0pZJZirPEkLSRNV\nxFzF4GoAVQCawCMRCQAsQDibKnrriT6f/VM+n3c9w4QXChZsfgdHIjprwgUgR0dH5+SSrCogU9TR\n0VFDRJfAzke0EycLs+/60TeSzrXratKTk037Dr/6tVXNrTm4rAMECQIRAMnnLisX90orl8/nZVnM\nQFm5wRf92xeGy8tNCM0BkJir+6WogExdDQDA80waHSnI1xRKjNjEzMKThmWalQCGLuwXYQM8DFAC\njGEG+gQhwcTDEkjotj68fNnyXxDS+aTM5ZfJhvqA1u/vkyxpGkI2NoTYMvuZec4Of1cBmboqllIQ\ns6DJd1nsehoIADIZVwQCHAwEjHN/aQM8zEx9BE4waBjMfcRIMGnDBucT+/btG37v6nfsgq5vg67n\nvSuviPkekCuWR6BpxKaZJqJdftaeSSogU0REBp27/SQpC3YfSp7HpqGnLr10yfd/9IN/2Xro0CFf\nHsved999b3Z2dh5FKFjurryq1ti9991Xtb8ILAQ5rWurORAYA1EyGAw+50fdYlCPeaeImRMQQoJI\nSp+2cnvH44TDOul6LmhZx/0Kx1uI6BGEQuNsWU7+/TfW+FXXuXZthSwr02UkPMLMXynWUH4/qIBM\nkZQyAQAkhIvyssLsux6J6KxpAkJzCrHvejwef5aZezgWHZCXLYs517ZOe2NTd+mSsHttazVHI0Ok\naf2BQOBrfvRaLCogUzQ8PNzHzKNsmml38aJIIY7hLrs0QkJ4bOh5Zi7IttKe5/13Now+jkYGnOuu\nrXFa1y6Ycq3Llkbyt354IQdDYxwMDhPRn955553+D2OZQSogU9Te3i6J6N+kZaZkfW3Irx1zf5a3\n/LIom0YKQIqI9vpdHwC2bt06qGnaH3AoNCxj0UHnhutq7Y/cWs+h0AWP9GbDIPuGDVX5zbctxNnZ\nhGcA/OV8WLxBDXefhltuucXVDONmyuXL2LR0/dhx375ry6aFQWf9uhqOxQZJ15+Ox+M/9Kv22+3Y\nsWPotttue4lMcy2bBnM0WuGtXlnJpkVaKuUgm3vHdxgcDmvu+1aV52+7pdG7dHEI0cggRyMDRPTn\nbW1tOwrV70xSE6amhzo7O79Eudz1YnyiwXrs6ye04eFpvxNhISj3Gx9dJJsWurK8/KSu6x+95557\nTvrR8Lv54he/2GAYxl8QsBrnlv2B5xlifDwvRsfyyGRdkpIRCuqyrMyU1VWBty37cwrAf4vH4wX5\nOlgMKiDT1NHRcbUQ4lEaG2sSwwkz8NjX36RczptOzfwHb671Vq8skwvK34RhfC0ej/+lX/1eAOrs\n7LyZiO4H0ATbCVA+H4HnGWCpg5kghAuhO2wZGTbNDIAxIvoSgG+c24t+3lAB8UF3d/d/kZ73O2Jk\ndIkYGPACT+w4Ten0RYeEATg3bKhyW5urZFmsF4HAQdd1P16M6art7e2ivr5+Bc5uR9cCoIaZK4lI\nABhl5hEi2u953i5N03rmWzDeogLig/b2diFBP4oGrKtjrsfa5CSs736vT5zuveCJUxwKafaHbq5z\nL10S4Vh0AMFgL4CPx+PxUwVs/aKcm9Mh2tvbCzZEfrZRAfHBqmvW3mmY+u9uuPbaumWXLklZ2WxU\n5PIx7cSbk/qLe0bebRgHRyK6/b5V5fKaNRUcCHiyrKyPDf0kEf1+PB6f8wuvzXUqINO0qrn11wD8\nEQhSSvzZ7265awmAO8lxokgma8hxgyKddsTp3jSlMg5l0h4buuBIROfa2oCsqQ6yEB6HgiMIh8dA\n9JLjOH+8bdu2OTk8fL5RAZmGVS3rbgW4HQCx5M8dfGnPt4CzK4ZIKbcC2ATHMZHPR8lxgiylTpJ1\nEElowoHQbAQDyXM3uqeZ+ZG2trYf4AInRimFpwIyRavWrt3IEp8nkEaEB/fv3f3Vt/9MR0fHYgCb\nzu2xcRnRv++My8weEQ0y8/NEtGtgYKCnlL7bzxUqIFOwoqWlVUD8LRgmBHUf2PNi94X83sMPPxzR\ndb0CQGZgYGB0rk4iKiUqIBdpVXPzKkB7GECIGP+4/6Xdf1XsnpTCUWOxLsKVa9ZdRtC+ACBEhJ37\nX9r9v4vdk1JYKiAXaEVz8yW6Jh9mIEbgf718yaL/AXUzPe+pr1gXYOW6dbXC5e0M1DP4Rc5lPun3\n5CVldlJXkPewYv36CnjyEQbqQbw/aup/oMJROkrmCtLe3m42NjYuyGazcnR0dORCniBdfv31UTPn\ndAFYDuCoHTDiR557zrdV0JXZb94GZPv27RWe593IzJsArCSi8p/5WAJIANjLzD/OZDLPv30XpNWr\nV4ehBzqYcBWAYya8LT09PRMzdwbKbDDvAtLV1VXGzHcT0a/C84KUz0dg22FIaUDy2QliZ99iO2yZ\naVhWioFxIcSX+vv7v9He3m4vW3aLFSxPPAimZgC9unTuvbCldJT5Zl4FpKOj4yYi+jPheRWcStWI\nvB1FPi+1vr40TUw6lE67JDR4kbDBCxaY3NgQlkIwB4Njb034yefzf9L9lb/fQqANID6je+Y9+/Y9\n11/sc1OKY76si0WdnZ13EVEbJVPVlMlUacOJnP7CT3q1YyfSb1+36qcnbZnCvfKKqL1ubTVy0XIZ\njYRztv29Sxcvzp44efIUCNtUOErbvLiCdHZ23k3A/ZiYaBDpTMTc9cyg9sqBiQs9OdZ1cjdcV+U2\nX12ZNY18CjASo6P3/Z/Pf/4fCtq4MuvN+YB0dXW9H8D/wvh4kzY2HjK/vfP0VJfRdFdcGbVv/kCj\nF4uMIhI5IaW88/777/d1WwNlbpnT70Eee+yxGDN/mpKpGpHORKYTDgDQD72aNHftGtAy2QrK52s1\nTfsM5sF/IsrUzfp7kGXLbrHCCxLLAFzuSdQCdLZn5rHX33hjXV1lZQ1lMlXmrmcG/FiAWX/l4IRX\nWxfk962qY8u6qru7+4Nbtmz53nTrKnPTbA2IWNHSsl5A/ApxolUyToBxhBi9RGfHP8Vi0UXhUPhX\nkUyaYjiR01454Ns7CuvZ54fl8mUxmU4vkOFwHIAKSImadV8fVrS0tAoWfwJggpj/KTNZ9f033nj6\nvLWmuru7b2fX/YxIjFxm7fhur3b0dV9X/rBvvKHKXddSJqurjjmO8xvbtm076md9ZW6YNVeQ9evX\nB5N575OCsZ4l/ueBfbtfeLefZ+ZNlM9HyLalOH7C99XD9deOJJ2Wq6vZdi1N1zYCUAEpQbMiIM3N\nzaGk7X2BiAcMyI/17OvJXMCvrYJth0VfX5pc1/dh5+LMmbxIpx2ORkIwQqv8rq/MDUUPyPr164NJ\n2/sCwMcP9Oz5S1zAHIuuri4DQDlJqdNEsmCrh9Nk0kZ1tYGzm2MqJajoj3nTjvdJAg8cvMBwAIAQ\nogJnVxLRKZ0u4D7faZekp0MFpGQV9QqyoqWlVQLXCif/67iI2Xn5fN41DAMEMAsq2IMGFoIAkkSk\nVhspUcW8ghCB/og8fG7//v0XdZOdSCTGmNljTbgcDhdw+7OQzkJzcXZovFKCihaQVVe3XktMufd6\nWvVOzk12SkBoDldUWgVoDyA6u8+3LhxmvrDtl5V5p3hXEIFfJZL/OI0Ke9kyU7KxLsSBgO8bAXmX\nNIWkZemwrJQQYo/f9ZW5oUgBaRcA1sK2p7xrEhHtgmWlmQR7V10R9bE5AIB75eUxMowsNM1xHOcZ\nv+src0NRAnLlmqeXAhi82HuPn5VOp3/CwBgHg2P22pYqNgzfbtZlZYXpXXVFuQyHRph539atWwf9\nqq3MLUUJiKbxcoCntbT/uTnk/xeRcIKjEeFuuM6fR7FEcG6+qRaGmUMgkATwsC91lTmpKAEhIASi\naY+dIqJvQoiTHIsOuNesqXRWXBWbbk37/TfWuAsbQ140NgDgR21tba9Mt6YydxUlICxYMPO0F26O\nx+MOM3+KLWtIhkIJ94MfaHDft6psSsWIYN+0sca9ek0Fl8X6YGjHwuHwZ6fbozK3FecrFuAIIl8e\nz7a1tb3OzJ/maOSMLIsN2jffVG9/+BfqLmafb1lVaeY/+itNztVrFsiyWB8CgUHTNB/4rd/6rUk/\nelTmrqLsk17V2GgRcMuZ/v5v+1HvqaeeevP2228/xIaxlg3d4/Lys/t8BwOC8rZHqdT5G2oSwVt0\nSci97tpK++ZN9V5lheSyBb1kmUc8z9u6ZcsWNdVWKc58kNWrV4fZCDx9oGf3JpxdxM0XHR0di4no\nMwSsRDq9gHL5criuJTIZB5OTDqXSDmmCOBzRZSxqIRgUUtdyCIdHEAgkmflHkUjks+rKobylaBOm\nVrWs/XsJ/uKhvXt3+1z6rX2+4wAWw3FMN5Ot1hgsIAEQs9Bc6MKBZaWgaQ4z79M07aF77713v8+9\nKHNc0QKyumXdrzDz2gM9u/+4UMfo6upamslkbj3V1/+HSxcv2m0YRhSAx8xDRDQspdyr6/qP7733\nXjWURHlHRQtIc3NzKA9tpwfvN1/t6Rko1HFWNq/9fQB0sGeP2uxGuWhFG4vV09OTIeav6qT9KQoU\n1NXNzdcQ4QO652wvRH1l/ivqhKkrli7+e2bEVjW3/qLftS+//vqoJNEuwX/x8ssvj/tdXykNRXnM\n+5bDhw9zXVPjK2D+8/qGphNDA32n/Ki7YsWmiI78QwCeO9Sz95t+1FRKU1EDAgBDfX3j1XVNPSD5\nudqFjSfP9PefnE69NWvWlLMpv0jA/oM9e77gV59KaZo162KtuHrdVULw5xn8nAX5hZ6eC1rZ5D9Y\nfU3rTSzwhyB64kL3LleUdzNrAgK89QLR+iSB1jHz17OG2PHGiy++50u7ldesW03gj4PQRCQ/s3/v\n3oMz0a8y/82qgLxlxbp1VwlXfoxBNwhwDzMOE2mHJeNcWNgA5DIiLAeoFQCY+euczzypNthU/DQr\nA/KWNWvWlHuG0cISVwnQlUwcAQCwcAF5ggWOsCsOHNr34uEit6ooiqIoiqIoiqIoiqIoiqIoyuwz\nq9+DKDPnkUceWaBpWgURRZl51DTNM3feeWfB9l6ZK1RASlRXV1eIiDYy8yZmXk9Eobf/DDP3A9hF\nRLvi8fhLRWiz6FRASkx7e7teV1d3OxHFWcpqyucjlM9H4HoWpNTBLEgIlzXhsm5kEbCSMM0sMx8i\noi+UWlBUQEpIR0dHIxH9NZgvo0xmAaUzVeS6EKd70+JMIidSSQeOIzkU1hEJ617TwrCsqgyyptkc\nCQ8jGJwE8C0AfxWPx51in89MUAEpEV1dXdcw818Jx6nH+ESjSKXZ2L0noR04NEGO83N39+JYTLev\nba3yVl5ZDtPMyLKyPhZiDxH9YTwe921v+tmq6BOmlMJ76KGHlgshHhHZbBNNTDbqBw9PBp54slf0\n9mdJvvuyZJTPS/3Y8ZT2+rFJb2FjOYgq2LBipImVLS0t/7xr167zF+WbR1RA5rnt27dXAOigfH6R\nNpls0J95bsh65tkEPHlRW2dTJusZh1+dkHV1IQSsCg5YoVhZWeWOHTt+XKDWZ4Wi73KrFJbneb9H\njreYJiYbtD0vJcw9e8emXCxvS/PJ7/RpwwlPTEw0MvPt3d3d1/nY7qyjriDzWHd393IAn6KJ8YXa\nqV7X/O7TA9O96SRPsjh5Mu2tuKqadc2DYSxsaWl5YteuXRd1RZor1BVkHpNS3oV8Pkq2EzJ//G9D\nfj2REWMTjrH/4KhIZ6rAvLy+vn6jT6VnHRWQeerBBx+0iOh6ZLPl4s2TKTE4mPezvv7Ci6Pkujry\n+QiA9/tZezZRAZmnTNNcB+YQ5e2IfuxY0u/6lMt5oq8/Tbl8BMCG9vb2gu1XX0wqIPMUEV0BxwmA\nWehvnJj2dnfvRLx5KkWOEwIQq6+vbyjEMYpNBWT+qiYpdXJdiXS6IO8qxMSEAykNABBCVBfiGMWm\nAjJ/VbLn6ZTNugU7QirlMjPB8zRm9meX4VlGBWSeYuYciBi6Ubh/Y/NcbSEkgHk5NF4FZP5KsBAu\nBwMaqDBD7jgS0UkID0QspRwuyEGKTAVknhJCDJGm2QyQrK31ZUfht5PVVQFomg0AruvOy126VEDm\nKWb+CQzDhq7n3cuWRQtxDO/SJRE2zRSA17dt2zZSiGMUmwrIPBWPx48xcy9MK+kuX1rGQvj6PUsu\nbAzKWMxiK5Bk5nk7YFEFZB4jop0yHBzj8nLdu3pNmZ+18zdcVwPLTLKhZXRdf9rP2rOJCsj89jUI\ncYaDoVFnXUs1x2K+vO1237eqjBsagjIcHSai79xzzz3T2vRoNlMBmcfi8XiGmTsQDSdkLGrn/tPt\nTT99NDtFsmlh0N60sV6GQ8Nk6iP5fH5eb1SkhrvPc0899dRrmzdvbmDTbIQQFXLxoqg48WaKbPvd\npxK+A/eypRH7I7cu5Eg4iVhsUEr5qW3btr1WiL5nCxWQEtDS0vJ8NBZbg4AVhmHEvMuXV1I+59KZ\n4fyF3LlzIKA5G2+odm7cUMuRyDiXlw0S0d+2tbU9VfDmi0wt2lAiHnzwQcuyrD8D8y00mayhXG6B\nSIzk9YOHx8TR11MimfwPQ1IYgGxsCMjll0Wdq65YwKGQ5Ej4DAcCCQCfbWtr+25xzmRmqYCUFurs\n7PzPALaQ51VQOlOJfD5KUuqUyTqUTrtsOxLhoI5QSJemqUHXcwhYExwOjzFwRAjxF1u2bDlU7BOZ\nKSogJeixxx6LpdPpu5h5MxGVwbaDcF2LPKkzsyBBLmuaC8PIQtcdZj4O4HCXIV4AAAB5SURBVMuD\ng4P/3N7eftH3LnOZCkgJe/zxx7XR0dGrieg6Zm4gompmDgMYAZBg5jeEED+Ox+Onit2roiiKoiiK\noiiKoiiKoiiKoiiKoiiKoiiKoiiKoiiKoiiKoiiKoiiKoiiKoiiKoiiKoiiKoiiKoiiKoiiKcrH+\nP0YGuQ11KIrQAAAAAElFTkSuQmCC\n",
"text/plain": [
"<IPython.core.display.Image object>"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"text/plain": [
"<PropertyMap object with key type 'Vertex' and value type 'vector<double>', for Graph 0x9b77bf4c, at 0x9b77b3cc>"
]
},
"execution_count": 75,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"gt.graph_draw(T,\n",
" edge_pen_width=edge_weights,\n",
" inline=True,\n",
" output_size=(200, 200))"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Read the property of an edge, here its weight:"
]
},
{
"cell_type": "code",
"execution_count": 76,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"6.0\n"
]
}
],
"source": [
"print edge_weights[e_2]"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### Now Collect the vertices of degree 1 in a dictionary:\n",
" keys : degree of a vertex\n",
" values: vertices"
]
},
{
"cell_type": "code",
"execution_count": 77,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"from collections import defaultdict\n",
"deg_dic = defaultdict(list)\n",
"\n",
"for v in T.vertices():\n",
" degree = v.out_degree()\n",
" deg_dic[degree].append(v)\n",
"#print deg_dic\n",
"#print deg_dic.keys()\n",
"#print deg_dic.values()"
]
},
{
"cell_type": "code",
"execution_count": 78,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"vertices of degree 1:\n",
"vertex 2 weight of out edge: [6.0]\n",
"vertex 3 weight of out edge: [3.0]\n"
]
}
],
"source": [
"print 'vertices of degree 1:'\n",
"for d1 in deg_dic[1]:\n",
" weight_of_bound_edge = [edge_weights[e] for e in d1.out_edges()]\n",
" print 'vertex',d1,' weight of out edge:', weight_of_bound_edge"
]
},
{
"cell_type": "code",
"execution_count": 79,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"# INVERT DICTIONARY\n",
"#inv_map = {}\n",
"#for k, v in map.iteritems():\n",
"# inv_map[v] = inv_map.get(v, [])\n",
"# inv_map[v].append(k)"
]
},
{
"cell_type": "code",
"execution_count": 82,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"{'weight': <PropertyMap object with key type 'Edge' and value type 'double', for Graph 0x9b77bf4c, at 0x9b77bf0c>}\n",
"{}\n",
"<bound method Graph.list_properties of <Graph object, undirected, with 5 vertices and 4 edges at 0x9b77bf4c>>\n",
"{}\n"
]
}
],
"source": [
"print T.edge_properties\n",
"print T.vertex_properties\n",
"print T.list_properties\n",
"print T.vp"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 2",
"language": "python",
"name": "python2"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 2
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython2",
"version": "2.7.9"
}
},
"nbformat": 4,
"nbformat_minor": 0
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment