bump required numpy version to 2.0.0

needed for backward compatibility
groupeLIAMG · Dec 19, 2024 · 1317059 · 1317059
1 parent 21dd9ed
commit 1317059
Show file tree

Hide file tree

Showing 7 changed files with 203 additions and 8 deletions.
diff --git a/examples/example1.ipynb b/examples/example1.ipynb
@@ -2,7 +2,7 @@
  "cells": [
   {
    "cell_type": "code",
-   "execution_count": 1,
+   "execution_count": 4,
    "id": "secret-calculation",
    "metadata": {},
    "outputs": [],
@@ -12,7 +12,7 @@
     "import matplotlib.pyplot as plt\n",
     "\n",
     "plt.rcdefaults()\n",
-    "plt.style.use('seaborn-paper')\n",
+    "plt.style.use('seaborn-v0_8-paper')\n",
     "plt.rc('font', **{'family': 'sans-serif', 'sans-serif': ['Helvetica']})\n",
     "\n",
     "from ttcrpy.rgrid import Grid3d"
@@ -30,7 +30,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 2,
+   "execution_count": 5,
    "id": "available-dutch",
    "metadata": {},
    "outputs": [],
@@ -79,10 +79,19 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 3,
+   "execution_count": 6,
    "id": "choice-murder",
    "metadata": {},
-   "outputs": [],
+   "outputs": [
+    {
+     "name": "stderr",
+     "output_type": "stream",
+     "text": [
+      "/var/folders/kd/v1qqm8qd15506q5jp_2tz_bw0000gn/T/ipykernel_33100/3904946752.py:15: DeprecationWarning: Conversion of an array with ndim > 0 to a scalar is deprecated, and will error in future. Ensure you extract a single element from your array before performing this operation. (Deprecated NumPy 1.25.)\n",
+      "  tt, rays = grid.raytrace(src, rcv, return_rays=True)\n"
+     ]
+    }
+   ],
    "source": [
     "# Define the source location\n",
     "src = np.array([[0.5, 0.5, 0.5]])\n",
@@ -279,10 +288,59 @@
     "print('Compute time: ', compute_time3)"
    ]
   },
+  {
+   "cell_type": "code",
+   "execution_count": 3,
+   "id": "8bf80c5a-02c4-4509-9cad-b8d6883273d0",
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "text/plain": [
+       "['Solarize_Light2',\n",
+       " '_classic_test_patch',\n",
+       " '_mpl-gallery',\n",
+       " '_mpl-gallery-nogrid',\n",
+       " 'bmh',\n",
+       " 'classic',\n",
+       " 'dark_background',\n",
+       " 'fast',\n",
+       " 'fivethirtyeight',\n",
+       " 'ggplot',\n",
+       " 'grayscale',\n",
+       " 'petroff10',\n",
+       " 'seaborn-v0_8',\n",
+       " 'seaborn-v0_8-bright',\n",
+       " 'seaborn-v0_8-colorblind',\n",
+       " 'seaborn-v0_8-dark',\n",
+       " 'seaborn-v0_8-dark-palette',\n",
+       " 'seaborn-v0_8-darkgrid',\n",
+       " 'seaborn-v0_8-deep',\n",
+       " 'seaborn-v0_8-muted',\n",
+       " 'seaborn-v0_8-notebook',\n",
+       " 'seaborn-v0_8-paper',\n",
+       " 'seaborn-v0_8-pastel',\n",
+       " 'seaborn-v0_8-poster',\n",
+       " 'seaborn-v0_8-talk',\n",
+       " 'seaborn-v0_8-ticks',\n",
+       " 'seaborn-v0_8-white',\n",
+       " 'seaborn-v0_8-whitegrid',\n",
+       " 'tableau-colorblind10']"
+      ]
+     },
+     "execution_count": 3,
+     "metadata": {},
+     "output_type": "execute_result"
+    }
+   ],
+   "source": [
+    "plt.style.available"
+   ]
+  },
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "horizontal-naples",
+   "id": "d92e3c00-4aa3-4367-91cb-7ddfda2ba91c",
    "metadata": {},
    "outputs": [],
    "source": []
@@ -304,7 +362,7 @@
    "name": "python",
    "nbconvert_exporter": "python",
    "pygments_lexer": "ipython3",
-   "version": "3.9.7"
+   "version": "3.13.1"
   }
  },
  "nbformat": 4,

diff --git a/examples/example1.vtr b/examples/example1.vtr
@@ -0,0 +1,25 @@
+<?xml version="1.0"?>
+<VTKFile type="RectilinearGrid" version="0.1" byte_order="LittleEndian" header_type="UInt32" compressor="vtkZLibDataCompressor">
+  <RectilinearGrid WholeExtent="0 40 0 40 0 40">
+  <Piece Extent="0 40 0 40 0 40">
+    <PointData>
+    </PointData>
+    <CellData>
+      <DataArray type="Float64" Name="Velocity" format="binary" RangeMin="1.1" RangeMax="2.9000000000000004">
+        EAAAAACAAAAAUAAATQAAAFYAAABOAAAAVQAAAFUAAABNAAAAVgAAAFUAAABMAAAAVAAAAFMAAABMAAAAUgAAAE0AAABUAAAAPQAAAA==eJztxTERAAAIBKBo5vRb2sAzh7CQPlOxbdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bfvxC5dldHV4nO3FMQ0AIAwAMGnoZJqGGBIELDOxcLRPY7e7wrZt27Zt27Zt27Zt27Zt27Zt27Zt27Zt27Ztj32yPdu2bdu2bdu2bdu2bdu2bdu2bdu2/dEF2Jfn7Xic7cVBDQAACASgaPbfvDYGcOYQPqTPVGzbtm3btm3btm3btm3btm3btm3btm3btm3btm3btm3btm3btm3btm3btm3btm3bfvwCWVDFknic7cUxDQAgDACwOcN/srnh4CSbAASQ9mll26ts27Zt2/a3xzi2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bfvxBRzg+6V4nO3FsQkAIAwEwN8p9u6/kWAj2SGNcNdc0s6Obdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bdsee1W7tm3btm3btm3btm3btm37ox+cGn94eJztxTERAAAIBKBQ/m7/TC6eOYSF1JmObdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bduPX9Ot0jt4nO3FsQkAIAwAsKPs7jf1/1mQ0hs6CcmSWO3usG3btm3btm3btm3b336yPdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2PbgA0A9PHHic7cUxDQAgDAAwnIEaZhYsLcsccJP2aeyWM2zbtm3btm3btm3btm3btm3btm3btm3btm3btm3btm3btm3btm3b9vP3tLFs27Zt27Zt2/Y/FzkIWKN4nO3FQQ0AAAgEIKPZf14nZw7hQ+ZUx7Zt27Zt27Zt27Zt27Zt27Zt27Zt27Zt27Zt27Zt27Zt27Zt27Zt27Zt27Zt27btxy+zu5SxeJztxTENACAMALBlyqYGtZsnAh5IONqn00fU2LZt27Zt27Zt27Zt27Zt27btb19X2rZt27Zt27Zt27Zt27Zt27Zt27Zt27Zt27Zt27bth2+SB7VJeJztxUENACAMBDCyIAw1518CYY+ZoP00eerEtm3btm3btm3btm3btm3btm3btm3btm3btm3btm3btm3btm3btm3btm3btm3Pq23btv3RF2rAAFl4nO3FMQ0AAAgDsB0Iw7+qBR20T5MzG9u2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2H1/wLUA9eJztxbERACAIBLAvdC/WxGk9pqBJmiTjVGzbtm3btm3btm3btm3btm3btm3btm3b9tqvx7Vt27Zt27Zt27Zt27Zt27Zt27Zt27Zt24t/9dwB8Hic7cUxEQAACASgyV7W9NN65hAWMqc6tm3btm3btm3btm3btm3btm3btm3btm3btm3btm3btm3btm3btm3btm3btm3btm0/fgHR38OUeJztxTENACAMALCLIIgJwCZTu+ybB9qn+dq6adu2vzlO27Zt27Zt27Zt27Zt27Zt27Zt27Zt27Zt27Zt27Zt27Zt27Zt27Zt27Zt27Zte1wZagOHeJztxTERAAAIBKDRQH4A+6fyDOEGC+lTE9u2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2bdu2/fYC9/TE0w==
+      </DataArray>
+    </CellData>
+    <Coordinates>
+      <DataArray type="Float64" Name="Array 0x1057c9970" format="binary" RangeMin="0" RangeMax="20">
+        AQAAAACAAABIAQAAYQAAAA==eJxNxykOwlAUQNEnEAgEIaQhDWmAzsyTb5fGUlkCEon4R3DNyY347z0kP/wyxsSEU84454JLZlwx55oFN9xyx5IVazZs2bHnngceeeKZF16TL8bNM+6e8fCM5/gD74oTzg==
+      </DataArray>
+      <DataArray type="Float64" Name="Array 0x1057acf60" format="binary" RangeMin="0" RangeMax="20">
+        AQAAAACAAABIAQAAYQAAAA==eJxNxykOwlAUQNEnEAgEIaQhDWmAzsyTb5fGUlkCEon4R3DNyY347z0kP/wyxsSEU84454JLZlwx55oFN9xyx5IVazZs2bHnngceeeKZF16TL8bNM+6e8fCM5/gD74oTzg==
+      </DataArray>
+      <DataArray type="Float64" Name="Array 0x1648ca490" format="binary" RangeMin="0" RangeMax="20">
+        AQAAAACAAABIAQAAYQAAAA==eJxNxykOwlAUQNEnEAgEIaQhDWmAzsyTb5fGUlkCEon4R3DNyY347z0kP/wyxsSEU84454JLZlwx55oFN9xyx5IVazZs2bHnngceeeKZF16TL8bNM+6e8fCM5/gD74oTzg==
+      </DataArray>
+    </Coordinates>
+  </Piece>
+  </RectilinearGrid>
+</VTKFile>
diff --git a/examples/example1_rays_raypaths for shot no 1.vtp b/examples/example1_rays_raypaths for shot no 1.vtp
diff --git a/examples/example1_rays_spm_raypaths for shot no 1.vtp b/examples/example1_rays_spm_raypaths for shot no 1.vtp
@@ -0,0 +1,56 @@
+<?xml version="1.0"?>
+<VTKFile type="PolyData" version="0.1" byte_order="LittleEndian" header_type="UInt32" compressor="vtkZLibDataCompressor">
+  <PolyData>
+    <Piece NumberOfPoints="2313" NumberOfVerts="0" NumberOfLines="38" NumberOfStrips="0" NumberOfPolys="0">
+      <PointData>
+      </PointData>
+      <CellData>
+      </CellData>
+      <Points>
+        <DataArray type="Float32" Name="Points" NumberOfComponents="3" format="binary" RangeMin="0.8660254037844386" RangeMax="33.77499074759311">
+          AQAAAACAAABsbAAAQg8AAA==eJztXH+oZVUZPUXFFBpP0XhJxi1UxiHjJRhTaLM/SblKxRTeGCXjNKRcpcIkZRo0juNVx/FHr+c4jabwCAJ5Q/AI+icIDkEgA8EQCBIIjyAIIZh/guivztprf/tbb8hCGNOc88dizuzz3ffe3Xuttb+9v31O0zR7moLZ7OIM/r8rQPuzGWzvC5ptePOf3Spo9pzYeD2j/KxUMFwvJSL//ASwfZKId8bfMJtdngC2ryQC7dckgO0pEWi/IQFs35sItM8SwPY2EefGd5zN7khA+V6JQPsjCWD7aiLGPvl3fTJ81wSwfT0RaP9FAti+mYixD89GH85mv05A6Z9EoP33CWD7qUSMff729PlWInB9OhHbxyL6v6t9PvyOoY9PZpQxyP3Mvv9QAthW+zZJvyXpkyT9MHzuUAaupU+GtiMZbH8qg5+tfTXcX8tgTO23oe35DLa/mMHP1v4c7v88gzG1b4e2jQy2/zKDn619Ptz/VQZjav8Pbb/JYPtvM/jZOi7D/d9lMKaO0dB2MoPtf8jgZ+vYDff/mMGYVzLY/qcMxnN8MS6z2d8SwPHiWGMchzgjML47jBg58G7hwGz2WobygTF/zmDMXzIYU31guP/PjDK2mScc3w8YwHEnZ8iHDxtAniwZcGLj7uF62YiRVyOvgkuMf68BjH+/AYypXoQxN4Ax5xvAGOcYuPERA8gZ8o1c+rgB5NjEAHJypxEjJ0dOnn1OMuZCAxhzsQGMqX44xFxiAGMuNaBwqXAVHLvCAHKPvCUnP20AubpiALm924iR2yO337ncZswnDWDM5QYwpvrzELPLAMZcZQBjnPPg6mcNIIfJf3L7CwaQ88kAamRqxKiRUSPvfo0w5moDGHONAYyp88UQ83kDGHOdAYxx7YDzNxpALVBH1MhXDKB29hpAre0zYrvWcI/3Q3MD7wdsZMQeT5/3d4YxzyjjVbCUYj9nUjH0Bb5LRmgxJe7XtALq8cTGQdEj2lxrnWinEy10wu1OuNoJ9zrhWCcc64RLnXCmE250woFOxhp/a533U3jgQRnHto4X96jMAPYBx4v9c7MB7DOOF/tzZgD718duaxiT2wzgeLT4PYbxms2+ZQDHcG7EONZnY6xDv52ML75DnddEmwdTzGttHVNy4KsGsG84puw3H1P0JceU/XyHAex3H1Nw4G5ounDgHgPIge8bwLE9YMTIgbeSAz7uvPZxx3erPpzChw+Klts61uSG6zfVsWZ/+lijjznW7P/7DOB4+FiDGwfxewo3Ovx+4/ji365wY2HEyI23gxvh853wAd/5mwaw3b39oGi/rRwgZ1zvqXKA/ewcQN+TAxiXoS8N4Dg5BzCOjxjAcT1sAMcd/x4unHnSiO2cGeaojO05+7rUiZQzyGf+vgeIHD3nU8KZJJy5OUUNqBXO3Jfv8X4n3Dmc7wHBla7k5c6f9Qrm2S8VbFYwp/b8ud+GyJ97yZ9x7+WCXnLpXnLpXjjZSy7dCz974Wcv/OwlH+4lH+6Ft73kw73kw73wuZd8uJd8uBee95IP95IP98L/XvLhXvLhXnTRSz7cSz7ci156yYd7yYd70VEv+XAv+XAv+kLM9QYw5gboo8S47nrJtRDzZQMY43rsZY7uqza5vvqGEZuizQ3R5rpo89nhd91rADnn/gx+PmhA4WbWJvkKXR4qGqc2we/QJrhPbUIvoU1opWpzaH/aAOpt1Qjo8ScGUJ9HjRi1PGr53avlExtfM4Axvp4KXTOmrqeGmNsNYIzn4b3kZn3VOz3hO0ZspsjDN0Tv66Jx7A89ZAC56/Pv4eFveNQActw1Dt4fMaKtGqcnPGUANUSN0xNWDaDmqsaH9mMGULfHjYCunzOAOn/BiNETRk8YPeG/ecKJjf0GMMbX6eEPbL/LALZ7bgDf+IERmyny9o0Uufp69QdyG97waPGNmp+nmPehgR8bQJ3QE4bvMPydawZQW/QEnqOCHxwrXkFPwB56eAK0657Qj54wesLoCWfBE2Tvboj5rgHqD4zxtUMv+zx9irwC7T80gO2eY8BbfB2xmWIdsTHEP2YAdeA5xoviIaspcgnUwOAZa0Y90TeGvxu6gicVn6FvYC/ixAb843jxGfeNrPfRN0bfGH3jHeIbsj8sexHhIbz2tQni4R+L0u6+AZ/xtclmirXJhuQh6+InL6ZYj6ymyDcOiW904hvtNt8Yxq3gzHMB6wXwBvcNeIGfP0/VN3AuPM4JtNU3yj5twar4w3r1BH5H94RePOGkeMIp8YRXxAe2qpaHv2n499WC01XLg8eJlmsNPGM2e09Gk8991pq4SU3c4ixoY1IfH+5/MIMxtVZuUiu3OC/amNTNh/sXZLC91sqHto9msL3Wx4e2T2TwZ9ZauUmt3OIMVWNSNx/uX5nBmFpDN6mhW9RbG5N6+nD/MxmMqbV1k9q6xbmUxqTOPtz/XAZjas3dpOZuUeNrLDwEMV7PbSz8pLHwE7R7jb4xOesy3L8pgzHuM42Fz6Dd60eNhec0FvuTjYX/NCb14qHt6xn8bK37D/dvzWCM+1Jj4Uto9xpEY+FRjcUeSGPhV41JjXJouzODn601aIt1UGPhY42Fj6H9exn8rHtaY+5d1MgDRpyufkVN1VqGeNQr4lGnxKNOptgX7cWXNpPsjYgXrY5eNHrR6EWjF1nUWxuLXKsxqcUPbfdnNOJdjPH1WmORgzUWORjaf5TBz3o+1ljsBUGbvhcEffpeELTs6zjo2PMuaN/XcdC9r+PgFe518AX3uny+cPS60etGrxu97m3xujinhHbfy0J7Pbticn7JYn2KdvdGtHtNDe2HC3DtOSGuPSfE83WeE8ILfJ8cZ5R9nxye4etWnHV2z4THjJ7ZjJ45eubomeeMZ/L66QJcu3/i2nNOXNdzCBZeimv3Unive2n2lOqlg7YyeA7Vz6b21Vd5/rS+jwJ78Xv4PD28kh6L85BxTjVJTSHvMWagrhF+Cw91v11N8d4J+Kh770tSU8jnwQr6FD58KoUPb6WoJZxO8iyPuffSH70uAE90H4YPbpWfsWTx3gDsxf8jEculTzGG2It3X4XfuU/i2r1xcgbc33DtfoVr95+JhZ9MxB8moveJ6HciepyIviail4nwfyJ8nghXJ8LDifBwUnnIeoNzb9lijXOh8G3JolZ1vgXf8twz8m3k2xvwzb1uIhybWHjaRDg2EY5dKhxbfkOORT00OEaOeK4YvCIPtvYQwSXWLL0euipcOp54FnzXGVxCzufvxumFSy8XXt6d4tw8fvZrKc7hbwnHwIn6bKM535gD+nOL4ITnmcjvvE66JDy8wKLWuSzcmwj3Lhfu7RTu7RLurVTuMZe6qGB35SFzJudgEh4m4WESHqYz4DxMwsMkPEzCwyQ8TMLDJDxMwkM8J+bnA5PM3Un4iZh6Hljm6yS8RbufDUa7e2ay8MzrhM+7LTzzauHzivD5KuHzzpHPI5//p3yezR7PYMyRDOUzY+qZWOE24v1MLGKeyWCMe3gSzifh/HXC+d0j588RznOdTM5zneyc/6LJs1cW5zZvFM5PZS2G3+Frsanwfyr8nwr/p8L/qfB/KvyfCv+nwn+0P5nBdl+XTUULU9HCVLQwFS0g5mgGY36awRivFyDGz2wg5mcZjHG9TN9yvcS7BFUv4K+fAdsUvfSil1Oil1dFL1spcunTSd7/YK4X7mcdT8QO0Qv2qPwdDktVL3wfkr+rYbnqhdx1vewUvaxY7K3trhohP6kRcu+8gmnVC/dlPN/2/HhS9lkuK9hX9cL9FNdLK3pxr3fvd720opfWYo5oz4DrpbWYI1qLOaIVvbSil1bmi1a004p2EPNwBmNcR63sabSiKcT4/kYr+kKMn4dsRWut5E6t6A4x9VkM0SBinshgjOuxFT22Mje1ok3E+LNaiKnnt0WzrWi2FZ3i2ues1mLOul00uG/U4DmiQe51+py1XzT4bYs56w6LOetO0eBctDYXrc1Fa3PR2ly0NhetzUVrc9HaXLQ2F63NRWtz0dpctDYXrc1Fa3PR2lx0NBcdzUVHc9HRXHQ0t5jX5qKd+ZvWTuwLoabjz8sn0VEr+0Kd7AutVk3xuVHfF0LN5pIMcpn6wtnLyAnzedfEZ+pPVa2hphP7RVuitdNy5hmcdq29z+LdRzss8kNw2t9ltGyxp/Qx0dpEtHaFaG2naG3F4mzz7qo71DPi3UGpahB5V7wrdWqRQ94ketxb9Yj6QdSP9lnkkLdbnGFuLXLI/RY55Lzqkfv7nyq4p2qT+/iuTV/nAA9Y1E0855uWa9dpZzFXdlWnob15ufa5srOYK7ttiPmxE512os1O9NiJBjvRXSda60RfnWiqkzmrE311oinE1GeKRVOd5JCd6OsB0deBUV+jvs66vjyP9Fwy9HWo6ouaob44t7m+FqKvhehrsQ2hr4XoayH6Woi+FqKvhehrIfpaiL4Woq+F6Gsh+sK1534L0Reuny9YiNa2v+fl/11rfKbxtpKfUXc4kxA5pOoO3PUcckl0d5HF+/pUdxPZ99hpse9xpcW79Vaq7rhf6jWU0F3TXGtxdiJ0xz2EvyYidNc0X7LYG9krGrzF4B3EPot9EryHwvPP1iL/3G+Rf3r9Y1J05/nnPRb5570W+eeBqkHWm68t6KoeY19jWnTiGjwkGjwkGvTc0fNH1+BCNPiIaNC1tCi6oqaYO7p2DleNMHf0szqPW9SUjwj/t7+zZuT/f+Z/vH8y+A+/jvdMBv/BlXh+LVnUEK+XeWdqMe/g/Mjr5WzR3sp/ng3x95zvq/zH+ijODrUW8xHOdPhZoLnFfHSX8N9rfvhO91uc+TlgMR89aHGGp6taILddCw9bnLdZVC34nh218JjF2ZjDVRfkpK/RnrA4x/KkaOSpqhGsW+LMyapFTrhqsXZbrXqJuQFYszgPvGZR63nGYk45egZ87jgqetn+Xpi3Qi9x5i30Qv75fkJopOQ3BaELcmpvQWih4XmbkoPtEC2cJzmYamFZcrCJaOEyycFUCysyF+wWLVxr8YynagH5iz+/GVogPzwH22eRg90qWmhFC+C8n52bW+xF3FX5T654DgbO+/k3r3lPCuf9zFsnWnjIYi5YVP7Th30u8JpiqvvRsQ9NzjOPuaVg1WJe8D2stvDW54U1i3lhzWJeAG/9eR9fS3RlXVHfjWCRdx2zOM95TDjv/MYc8ZzFHPG8cP4F+xc+WVnC
+          <InformationKey name="L2_NORM_RANGE" location="vtkDataArray" length="2">
+            <Value index="0">
+              0.86602540378
+            </Value>
+            <Value index="1">
+              33.774990748
+            </Value>
+          </InformationKey>
+        </DataArray>
+      </Points>
+      <Verts>
+        <DataArray type="Int64" Name="connectivity" format="binary" RangeMin="1e+299" RangeMax="-1e+299">
+          AAAAAACAAAAAAAAA
+        </DataArray>
+        <DataArray type="Int64" Name="offsets" format="binary" RangeMin="1e+299" RangeMax="-1e+299">
+          AAAAAACAAAAAAAAA
+        </DataArray>
+      </Verts>
+      <Lines>
+        <DataArray type="Int64" Name="connectivity" format="binary" RangeMin="0" RangeMax="2312">
+          AQAAAACAAABISAAAyg0AAA==eJw13MMWGIqSAMAX27Zt27Zt2/aNbdu2bdu2bXsWU+lNfUKf5v/+9/8RgAEZiIEZhEEZjMEZgiEZiqEZhmEZjuEZgREZiZEZhVEZjdEZgzEZi7EZh3EZj/GZgAmZiImZhEmZjMmZgimZiqmZhmmZjumZgRmZiZmZhVmZjdmZgzmZi7mZh3mZj/lZgAVZiIVZhEVZjMVZgiVZiqVZhmVZjuVZgRVZiZVZhVVZjdVZgzVZi7VZh3VZj/XZgA3ZiI3ZhE3ZjM3Zgi3Ziq3Zhm3Zju3ZgR3ZiZ3ZhV3Zjd3Zgz3Zi73Zh33Zj/+xPwdwIAdxMIdwKIdxOEdwJEdxNMdwLMdxPCdwIidxMqdwKqdxOmdwJmdxNudwLudxPhdwIRdxMZdwKZdxOVdwJVdxNddwLddxPTdwIzdxM7dwK7dxO3dwJ3dxN/dwL/dxPw/wIA/xMI/wKI/xOE/wJE/xNM/wLM/xPC/wIi/xMq/wKq/xOm/wJm/xNu/wLu/xPh/wIR/xMZ/wKZ/xOV/wJV/xNd/wLd/xPT/wIz/xM7/wK7/xO3/wJ3/xN//wL/8l/gAMyEAMzCAMymAMzhAMyVAMzTAMy3AMzwiMyEiMzCiMymiMzhiMyViMzTiMy3iMzwRMyERMzCRMymRMzhRMyVRMzTRMy3RMzwzMyEzMzCzMymzMzhzMyVzMzTzMy3zMzwIsyEIszCIsymIszhIsyVIszTIsy3IszwqsyEqszCqsymqszhqsyVqszTqsy3qszwZsyEZszCZsymZszhZsyVZszTZsy3Zszw7syE7szC7sym7szh7syV7szT7sy378j/05gAM5iIM5hEM5jMM5giM5iqM5hmM5juM5gRM5iZM5hVM5jdM5gzM5i7M5h3M5j/O5gAu5iIu5hEu5jMu5giu5iqu5hmu5juu5gRu5iZu5hVu5jdu5gzu5i7u5h3u5j/t5gAd5iId5hEd5jMd5gid5iqd5hmd5jud5gRd5iZd5hVd5jdd5gzd5i7d5h3d5j/f5gA/5iI/5hE/5jM/5gi/5iq/5hm/5ju/5gR/5iZ/5hV/5jd/5gz/5i7/5h3/5r+APwIAMxMAMwqAMxuAMwZAMxdAMw7AMx/CMwIiMxMiMwqiMxuiMwZiMxdiMw7iMx/hMwIRMxMRMwqRMxuRMwZRMxdRMw7RMx/TMwIzMxMzMwqzMxuzMwZzMxdzMw7zMx/wswIIsxMIswqIsxuIswZIsxdIsw7Isx/KswIqsxMqswqqsxuqswZqsxdqsw7qsx/pswIZsxMZswqZsxuZswZZsxdZsw7Zsx/bswI7sxM7swq7sxu7swZ7sxd7sw77sx//YnwM4kIM4mEM4lMM4nCM4kqM4mmM4luM4nhM4kZM4mVM4ldM4nTM4k7M4m3M4l/M4nwu4kIu4mEu4lMu4nCu4kqu4mmu4luu4nhu4kZu4mVu4ldu4nTu4k7u4m3u4l/u4nwd4kId4mEd4lMd4nCd4kqd4mmd4lud4nhd4kZd4mVd4ldd4nTd4k7d4m3d4l/d4nw/4kI/4mE/4lM/4nC/4kq/4mm/4lu/4nh/4kZ/4mV/4ld/4nT/4k7/4m3/4l/8a/QEYkIEYmEEYlMEYnCEYkqEYmmEYluEYnhEYkZEYmVEYldEYnTEYk7EYm3EYl/EYnwmYkImYmEmYlMmYnCmYkqmYmmmYlumYnhmYkZmYmVmYldmYnTmYk7mYm3mYl/mYnwVYkIVYmEVYlMVYnCVYkqVYmmVYluVYnhVYkZVYmVVYldVYnTVYk7VYm3VYl/VYnw3YkI3YmE3YlM3YnC3Ykq3Ymm3Ylu3Ynh3YkZ3YmV3Yld3YnT3Yk73Ym33Yl/34H/tzAAdyEAdzCIdyGIdzBEdyFEdzDMdyHMdzAidyEidzCqdyGqdzBmdyFmdzDudyHudzARdyERdzCZdyGZdzBVdyFVdzDddyHddzAzdyEzdzC7dyG7dzB3dyF3dzD/dyH/fzAA/yEA/zCI/yGI/zBE/yFE/zDM/yHM/zAi/yEi/zCq/yGq/zBm/yFm/zDu/yHu/zAR/yER/zCZ/yGZ/zBV/yFV/zDd/yHd/zAz/yEz/zC7/yG7/zB3/yF3/zD//y34A/AAMyEAMzCIMyGIMzBEMyFEMzDMMyHMMzAiMyEiMzCqMyGqMzBmMyFmMzDuMyHuMzARMyERMzCZMyGZMzBVMyFVMzDdMyHdMzAzMyEzMzC7MyG7MzB3MyF3MzD/MyH/OzAAuyEAuzCIuyGIuzBEuyFEuzDMuyHMuzAiuyEiuzCquyGquzBmuyFmuzDuuyHuuzARuyERuzCZuyGZuzBVuyFVuzDduyHduzAzuyEzuzC7uyG7uzB3uyF3uzD/uyH/9jfw7gQA7iYA7hUA7jcI7gSI7iaI7hWI7jeE7gRE7iZE7hVE7jdM7gTM7ibM7hXM7jfC7gQi7iYi7hUi7jcq7gSq7iaq7hWq7jem7gRm7iZm7hVm7jdu7gTu7ibu7hXu7jfh7gQR7iYR7hUR7jcZ7gSZ7iaZ7hWZ7jeV7gRV7iZV7hVV7jdd7gTd7ibd7hXd7jfT7gQz7iYz7hUz7jc77gS77ia77hW77je37gR37iZ37hV37jd/7gT/7ib/7hX/5b7AvAgAzEwAzCoAzG4AzBkAzF0AzDsAzH8IzAiIzEyIzCqIzG6IzBmIzF2IzDuIzH+EzAhEzExEzCpEzG5EzBlEzF1EzDtEzH9MzAjMzEzMzCrMzG7MzBnMzF3MzDvMzH/CzAgizEwizCoizG4izBkizF0izDsizH8qzAiqzEyqzCqqzG6qzBmqzF2qzDuqzH+mzAhmzExmzCpmzG5mzBlmzF1mzDtmzH9uzAjuzEzuzCruzG7uzBnuzF3uzDvuzH/9ifAziQgziYQziUwzicIziSoziaYziW4zieEziRkziZUziV0zidMziTszibcziX8zifC7iQi7iYS7iUy7icK7iSq7iaa7iW67ieG7iRm7iZW7iV27idO7iTu7ibe7iX+7ifB3iQh3iYR3iUx3icJ3iSp3iaZ3iW53ieF3iRl3iZV3iV13idN3iTt3ibd3iX93ifD/iQj/iYT/iUz/icL/iSr/iab/iW7/ieH/iRn/iZX/iV3/idP/iTv/ibf/iX/xb6AzAgAzEwgzAogzE4QzAkQzE0wzAswzE8IzAiIzEyozAqozE6YzAmYzE24zAu4zE+EzAhEzExkzApkzE5UzAlUzE10zAt0zE9MzAjMzEzszArszE7czAnczE38zAv8zE/C7AgC7Ewi7Aoi7E4S7AkS7E0y7Asy7E8K7AiK7Eyq7Aqq7E6a7Ama7E267Au67E+G7AhG7Exm7Apm7E5W7AlW7E127At27E9O7AjO7Ezu7Aru7E7e7Ane7E3+7Av+/E/9ucADuQgDuYQDuUwDucIjuQojuYYjuU4jucETuQkTuYUTuU0TucMzuQszuYczuU8zucCLuQiLuYSLuUyLucKruQqruYaruU6rucGbuQmbuYWbuU2bucO7uQu7uYe7uU+7ucBHuQhHuYRHuUxHucJnuQpnuYZnuU5nucFXuQlXuYVXuU1XucN3uQt3uYd3uU93ucDPuQjPuYTPuUzPucLvuQrvuYbvuU7vucHfuQnfuYXfuU3fucP/uQv/uYf/uW/Q74ADMhADMwgDMpgDM4QDMlQDM0wDMtwDM8IjMhIjMwojMpojM4YjMlYjM04jMt4jM8ETMhETMwkTMpkTM4UTMlUTM00TMt0TM8MzMhMzMwszMpszM4czMlczM08zMt8zM8CLMhCLMwiLMpiLM4SLMlSLM0yLMtyLM8KrMhKrMwqrMpqrM4arMlarM06rMt6rM8GbMhGbMwmbMpmbM4WbMlWbM02bMt2bM8O7MhO7Mwu7Mpu7M4e7Mle7M0+7Mt+/I/9OYADOYiDOYRDOYzDOYIjOYqjOYZjOY7jOYETOYmTOYVTOY3TOYMzOYuzOYdzOY/zuYALuYiLuYRLuYzLuYIruYqruYZruY7ruYEbuYmbuYVbuY3buYM7uYu7uYd7uY/7eYAHeYiHeYRHeYzHeYIneYqneYZneY7neYEXeYmXeYVXeY3XeYM3eYu3eYd3eY/3+YAP+YiP+YRP+YzP+YIv+Yqv+YZv+Y7v+YEf+Ymf+YVf+Y3f+YM/+Yu/+Yd/+e+APwADMhADMwiDMhiDMwRDMhRDMwzDMhzDMwIjMhIjMwqjMhqjMwZjMhZjMw7jMh7jMwETMhETMwmTMhmTMwVTMhVTMw3TMh3TMwMzMhMzMwuzMhuzMwdzMhdzMw/zMh/zswALshALswiLshiLswRLshRLswzLshzLswIrshIrswqrshqrswZrshZrsw7rsh7rswEbshEbswmbshmbswVbshVbsw3bsh3bswM7shM7swu7shu7swd7shd7sw/7sh//Y38O4EAO4mAO4VAO43CO4EiO4miO4ViO43hO4ERO4mRO4VRO43TO4EzO4mzO4VzO43wu4EIu4mIu4VIu43Ku4Equ4mqu4Vqu43pu4EZu4mZu4VZu43bu4E7u4m7u4V7u434e4EEe4mEe4VEe43Ge4Eme4mme4Vme43le4EVe4mVe4VVe43Xe4E3e4m3e4V3e430+4EM+4mM+4VM+43O+4Eu+4mu+4Vu+43t+4Ed+4md+4Vd+43f+4E/+4m/+4V/+e9wTgAEZiIEZhEEZjMH5f/PJoDI=
+        </DataArray>
+        <DataArray type="Int64" Name="offsets" format="binary" RangeMin="5" RangeMax="2313">
+          AQAAAACAAAAwAQAAbgAAAA==eJwtzqsRg1AABVHCH4ePQGQQaBwFIPCvBTQzpAYkjgLyXAoAjU4XtACSYYK4u+bYDRyV4gsrNPjGCVfcsXjIFr+449OVNfZo8cDGkx880fhywSyQI/6xDPlAixvmkexwxgvrWA74wySRN7jDEzg=
+        </DataArray>
+      </Lines>
+      <Strips>
+        <DataArray type="Int64" Name="connectivity" format="binary" RangeMin="1e+299" RangeMax="-1e+299">
+          AAAAAACAAAAAAAAA
+        </DataArray>
+        <DataArray type="Int64" Name="offsets" format="binary" RangeMin="1e+299" RangeMax="-1e+299">
+          AAAAAACAAAAAAAAA
+        </DataArray>
+      </Strips>
+      <Polys>
+        <DataArray type="Int64" Name="connectivity" format="binary" RangeMin="1e+299" RangeMax="-1e+299">
+          AAAAAACAAAAAAAAA
+        </DataArray>
+        <DataArray type="Int64" Name="offsets" format="binary" RangeMin="1e+299" RangeMax="-1e+299">
+          AAAAAACAAAAAAAAA
+        </DataArray>
+      </Polys>
+    </Piece>
+  </PolyData>
+</VTKFile>
diff --git a/examples/figs/example1_tt_fsm.pdf b/examples/figs/example1_tt_fsm.pdf
diff --git a/examples/figs/example1_tt_spm.pdf b/examples/figs/example1_tt_spm.pdf
diff --git a/requirements.txt b/requirements.txt
@@ -1,5 +1,5 @@
 setuptools
 cython
-numpy>=1.20.1
+numpy>=2.0.0
 scipy
 vtk