MUSIC/tools/check_output.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "import h5py\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "%matplotlib widget"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 106,
   "metadata": {},
   "outputs": [],
   "source": [
    "class music_grid(object):\n",
    "\n",
    "  def __init__( self, fname, fieldname='rho' ):\n",
    "    with h5py.File(fname) as hf:\n",
    "      self.levelmin = hf['/header'].attrs['levelmin']\n",
    "      self.levelmax = hf['/header'].attrs['levelmax']\n",
    "      self.len = np.zeros((1+self.levelmax-self.levelmin,3), dtype=int)\n",
    "      self.len[:,0] = np.array( hf['/header/grid_len_x'])\n",
    "      self.len[:,1] = np.array( hf['/header/grid_len_y'])\n",
    "      self.len[:,2] = np.array( hf['/header/grid_len_z'])\n",
    "      self.off = np.zeros((1+self.levelmax-self.levelmin,3), dtype=int)\n",
    "      self.off[:,0] = np.array( hf['/header/grid_off_x'] )\n",
    "      self.off[:,1] = np.array( hf['/header/grid_off_y'] )\n",
    "      self.off[:,2] = np.array( hf['/header/grid_off_z'] )\n",
    "\n",
    "      self.absoff = np.zeros_like(self.off)\n",
    "      for i in range(self.off.shape[0] ):\n",
    "        self.absoff[i,:] = self.off[i,:] * 2**(1+self.levelmax-self.levelmin-i)\n",
    "      self.absoff = np.cumsum(self.absoff,axis=0)\n",
    "\n",
    "      self.density = {}\n",
    "      for i in range(self.levelmin, self.levelmax+1):\n",
    "        self.density[i] = np.array( hf[f'level_{i:03d}_DM_'+fieldname])[4:-4,4:-4,4:-4]\n",
    "      \n",
    "    \n",
    "  def get_slice( self, islz ):\n",
    "    scale = 2**(self.levelmax-self.levelmin)\n",
    "    iislz = islz//scale\n",
    "    img = self.density[self.levelmin][...,iislz].repeat( scale, axis=0 ).repeat( scale, axis=1 )# *0\n",
    "\n",
    "    for i,ilvl in enumerate(range(self.levelmin+1, self.levelmax+1)):\n",
    "      scale  = 2**(self.levelmax-ilvl)\n",
    "      if (islz >= self.absoff[i+1,2]) & (islz <  self.absoff[i+1,2]+scale*self.len[i+1,2]):\n",
    "        iislz = (islz - self.absoff[i+1,2])//scale\n",
    "        imgt = self.density[ilvl][...,iislz].repeat( scale, axis=0 ).repeat( scale, axis=1 )# *0+i+1\n",
    "        il = self.absoff[i+1,0]\n",
    "        ir = self.absoff[i+1,0]+scale*self.len[i+1,0]\n",
    "        jl = self.absoff[i+1,1]\n",
    "        jr = self.absoff[i+1,1]+scale*self.len[i+1,1]\n",
    "        \n",
    "        img[il:ir,jl:jr] = imgt\n",
    "\n",
    "    return img"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 220,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Linf = 0.05126058941694821, L2 = 0.008027430836950709\n"
     ]
    }
   ],
   "source": [
    "field = 'vz'\n",
    "\n",
    "mg1 = music_grid('../build/debug.hdf5',field)\n",
    "img1 = mg1.get_slice( 512 )\n",
    "\n",
    "mg2 = music_grid('../build/debug_full.hdf5',field)\n",
    "img2 = mg2.get_slice( 512 )\n",
    "\n",
    "diff = np.abs(img1-img2)[400:600,400:600]\n",
    "field = img2[400:600]\n",
    "print(f'Linf = {diff.flatten().max()/field.flatten().std()}, L2 = {diff.flatten().std()/field.flatten().std()}')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 221,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Text(0.5, 1.0, 'difference')"
      ]
     },
     "execution_count": 221,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
       "model_id": "44db17a691ef4d2bbddb1a9fdd8650e0",
       "version_major": 2,
       "version_minor": 0
      },
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABXgAAAK8CAYAAABV1dcbAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOzdd5gkV3n3/e85VR2nJ+3szkZJq4wQkkAJEBJJgABjLJIwYB5EtjEPNsnI2CTbJBsDxvACD7ZJjhgbjLEBg0gSiCyQBFpplbWrjZNDp6pz3j8qdPXszAZJgAf9PtclzU5Nd1fo6rv63HWfc4z33iMiIiIiIiIiIiIiq479ZW+AiIiIiIiIiIiIiNw9SvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiMivlE9+8pPc7373o1QqMTIy8sveHBEROUJvfvObMcbkv2/dupVLL7207zHbt2/ncY97HMPDwxhj+OxnPwvA97//fc477zwGBgYwxvDjH//4F7fhIiK/JOEvewNERERERO4t27Zt49JLL+Xxj388l112GfV6/Ze9SSIi8nPwvOc9j1tvvZW3vvWtjIyMcPbZZ9PtdnnGM55BtVrlPe95D/V6nWOOOeaXvakiIj93SvCKiIiIyK+Mr3/96zjn+Ku/+itOOOGEX/bmiIjIveCGG27A2l4H5GazyVVXXcUf/dEf8fKXvzxfvm3bNm6//XY+8pGP8KIXveiXsakiIr8UGqJBRERERP5XW1hYOOzH7t27F+BeHZrhSNYvIiL3vkqlQqlUyn/ft28fcGCs1zVARO6rlOAVkfukbFyvm266iUsvvZSRkRGGh4d5/vOfz+LiYt9j//7v/56zzjqLWq3GmjVr+M3f/E3uvPPO/O/ve9/7CIKA6enpfNlf/uVfYozhVa96Vb4sjmMGBwd53et
      "text/html": [
       "\n",
       "            <div style=\"display: inline-block;\">\n",
       "                <div class=\"jupyter-widgets widget-label\" style=\"text-align: center;\">\n",
       "                    Figure\n",
       "                </div>\n",
       "                <img src='data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAABXgAAAK8CAYAAABV1dcbAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOzdd5gkV3n3/e85VR2nJ+3szkZJq4wQkkAJEBJJgABjLJIwYB5EtjEPNsnI2CTbJBsDxvACD7ZJjhgbjLEBg0gSiCyQBFpplbWrjZNDp6pz3j8qdPXszAZJgAf9PtclzU5Nd1fo6rv63HWfc4z33iMiIiIiIiIiIiIiq479ZW+AiIiIiIiIiIiIiNw9SvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiIiIiIiIiKxSSvCKiIiIiIiIiIiIrFJK8IqIiIiIiIiIiIisUkrwioiIiMivlE9+8pPc7373o1QqMTIy8sveHBEROUJvfvObMcbkv2/dupVLL7207zHbt2/ncY97HMPDwxhj+OxnPwvA97//fc477zwGBgYwxvDjH//4F7fhIiK/JOEvewNERERERO4t27Zt49JLL+Xxj388l112GfV6/Ze9SSIi8nPwvOc9j1tvvZW3vvWtjIyMcPbZZ9PtdnnGM55BtVrlPe95D/V6nWOOOeaXvakiIj93SvCKiIiIyK+Mr3/96zjn+Ku/+itOOOGEX/bmiIjIveCGG27A2l4H5GazyVVXXcUf/dEf8fKXvzxfvm3bNm6//XY+8pGP8KIXveiXsakiIr8UGqJBRERERP5XW1hYOOzH7t27F+BeHZrhSNYvIiL3vkqlQqlUyn/ft28fcGCs1zVARO6rlOAVkfukbFyvm266iUsvvZSRkRGGh4d5/vOfz+LiYt9j//7v/56zzjqLWq3GmjVr+M3f/E3uvPPO/O/ve9/7CIK
       "            </div>\n",
       "        "
      ],
      "text/plain": [
       "Canvas(toolbar=Toolbar(toolitems=[('Home', 'Reset original view', 'home', 'home'), ('Back', 'Back to previous …"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig,ax  = plt.subplots(1,3,figsize=(14,7), sharex=True, sharey=True)\n",
    "ss = img1.std()\n",
    "\n",
    "cmap = 'RdBu'\n",
    "\n",
    "ax[0].imshow( img1, cmap=cmap, vmin=-3*ss,vmax=3*ss)\n",
    "ax[0].set_title('new')\n",
    "ax[1].imshow( img2, cmap=cmap, vmin=-3*ss,vmax=3*ss)\n",
    "ax[1].set_title('reference')\n",
    "ax[2].imshow( img2-img1, cmap=cmap, vmin=-1*ss/100,vmax=1*ss/100)\n",
    "ax[2].set_title('difference')\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "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.9.15"
  },
  "orig_nbformat": 4,
  "vscode": {
   "interpreter": {
    "hash": "397704579725e15f5c7cb49fe5f0341eb7531c82d19f2c29d197e8b64ab5776b"
   }
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}