{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import matplotlib.pyplot as plt"
   ],
   "metadata": {
    "collapsed": false,
    "pycharm": {
     "name": "#%%\n"
    }
   }
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "outputs": [],
   "source": [
    "data1 = np.loadtxt(\"../data/out/1Euler_msd_harmonic_konv.dat\")\n",
    "data2 = np.loadtxt(\"../data/out/2Heun_msd_harmonic_konv.dat\")\n",
    "data3 = np.loadtxt(\"../data/out/3Exact_msd_harmonic_konv.dat\")\n",
    "data4 = np.loadtxt(\"../data/out/4BDAS_msd_harmonic_konv.dat\")\n",
    "data5 = np.loadtxt(\"../data/out/5MBD_msd_harmonic_konv.dat\")"
   ],
   "metadata": {
    "collapsed": false,
    "pycharm": {
     "name": "#%%\n"
    }
   }
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "outputs": [
    {
     "data": {
      "text/plain": "<Figure size 432x288 with 1 Axes>",
      "image/png": "[encoded data, 23380 characters]\n"
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.loglog(data1[:,0], data1[:,1],label=\"Euler\")\n",
    "plt.loglog(data2[:,0], data2[:,1],label=\"Heun\")\n",
    "plt.loglog(data3[:,0], data3[:,1],label=\"Exact\")\n",
    "plt.loglog(data4[:,0], data4[:,1],label=\"BDAS\")\n",
    "plt.loglog(data5[:,0], data5[:,1],label=\"MBD\")\n",
    "plt.title(\"Mean square displacment\")\n",
    "plt.xlabel(\"Time-Step\")\n",
    "plt.ylabel(\"Error\")\n",
    "plt.legend()\n",
    "del data1,data2,data3,data4,data5"
   ],
   "metadata": {
    "collapsed": false,
    "pycharm": {
     "name": "#%%\n"
    }
   }
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "outputs": [],
   "source": [
    "data1 = np.loadtxt(\"../data/out/1Euler_oaf_harmonic_konv.dat\")\n",
    "data2 = np.loadtxt(\"../data/out/2Heun_oaf_harmonic_konv.dat\")\n",
    "#data3 = np.loadtxt(\"../data/3Exact_oaf_harmonic.dat\")\n",
    "#data4 = np.loadtxt(\"../data/4BDAS_oaf_harmonic.dat\")\n",
    "data5 = np.loadtxt(\"../data/out/5MBD_oaf_harmonic_konv.dat\")"
   ],
   "metadata": {
    "collapsed": false,
    "pycharm": {
     "name": "#%%\n"
    }
   }
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "outputs": [
    {
     "data": {
      "text/plain": "Text(0, 0.5, 'Error')"
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "text/plain": "<Figure size 432x288 with 1 Axes>",
      "image/png": "[encoded data, 20288 characters]\n"
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.loglog(data1[:,0], data1[:,1],label=\"Euler\")\n",
    "plt.loglog(data2[:,0], data2[:,1],label=\"Heun\")\n",
    "#plt.loglog(data3[:,0], data3[:,1])\n",
    "#plt.loglog(data4[:,0], data4[:,1])\n",
    "plt.loglog(data5[:,0], data5[:,1],label=\"MBD\")\n",
    "plt.legend()\n",
    "plt.title(\"Oaf\")\n",
    "plt.xlabel(\"Time-Step\")\n",
    "plt.ylabel(\"Error\")"
   ],
   "metadata": {
    "collapsed": false,
    "pycharm": {
     "name": "#%%\n"
    }
   }
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "outputs": [],
   "source": [],
   "metadata": {
    "collapsed": false,
    "pycharm": {
     "name": "#%%\n"
    }
   }
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "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.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}