{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import os\n", "\n", "import discretisedfield as df" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "test_filename = os.path.join(\n", " \"..\", \"discretisedfield\", \"tests\", \"test_sample/\", \"skyrmion.omf\"\n", ")" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "(20, 20, 1)" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" } ], "source": [ "f = df.Field.fromfile(test_filename)\n", "f.mesh.n" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAbQAAAF3CAYAAADNSL5cAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAAAnlElEQVR4nO3df3RU9Z3/8ddMGPKrRDbijxaBUAxDoV9p0eKmESESCocEqcQeEI4sq3UDAkXRLq5HVKhB0LCLSHNqVwVkVZDq2VXIoYaNmmIM2rPnLMdUpEnYQLOASoQkGxLCzP3+gU6N+WG4nzuZm8vzcc6cQ+7cz5137kx4z+fHvW+fZVmWAADo4/yxDgAAACeQ0AAAnkBCAwB4AgkNAOAJJDQAgCeQ0AAAntAv1gE47f8t/xej9m0DHArEgC9k2D5sHoO/zax9XKtpe/OrSeJazNoHzpidyH7N5m9Ev2azD0O//zN7I31nDD8IknytZ83anzWPQYbHsM4Z/lH6fWbtJSn1EqPmeyrXmMfQifDxkUbt/VcecigSDyY0AEDvCcvsi5uTw4QMOQIAPIEeGgDAtpBlODTvUBxOHwsAcJEJyz13TyShAQBsM51DcxJzaAAAT6CHBgCwLRTFgi2nTp3Sww8/rA8//FA+n0+zZs3S4sWLu9yfHhoAwLawLKNHdx588EGNHDlSpaWl2rlzp959913V1NR0uT89NACAbaEoLQo5ceKE/vjHP2rDhg2SpNTUVL300kvdtiGhAQBsi9Yqx48++khDhw7Vpk2btHfvXgUCAc2dO1ezZ8/usg1DjgAA1zl9+rQOHjyoUaNGqbi4WI8//rjWrl2r8vLyLtuQ0AAAtoUsy+jRlYEDB+ryyy/X9OnTJUmjR49WVlaW3nnnnS7bMOQIALAtWlehDR06VI2NjQqHw/L7z/e9/H6/4uLiumxDDw0AYFtIltGjK8OHD1cwGNQLL7wgSTp69Kjeeecd3XjjjV228VwPrX+D6QSlA2UeTCMw/MrjO2ceg/+c2Xk0Lj/TZj7RbBqD6ftgOVAyJBRv+J3TFzBq7u/f9bfhHh/jrFkMvjbD0i1OHMMF5WPODUwyPkY0hKJ456u1a9fqoYce0gsvvKBAIKB77rlHf/u3f9vl/p5LaAAAb7jqqqu0ZcuWHu9PQgMA2OaeOzmS0AAABkIumKb5EgkNAGBb2D3VY1jlCADwBnpoAADbGHIEAHgCCQ0A4Alhi4QGAPAAN/XQWBQCAPAEemgAANtCLuoXkdAAALa5aQ7NPan1Cw0NDZowYYIeeOABSVJ9fb0WL16sCRMmaOLEiXryySdldVNDBwDQe0LyGT2c5LoeWkFBgRISEiI/r1y5UikpKSorK1Nzc7PmzZunwYMHa+7cuTGMEgAgSSHLPf0i90Qi6a233tLRo0eVm5srSWpqatLbb7+t/Px8+Xw+JScna+7cuXr99ddjHCkAwG1ck9BOnz6tNWvWaM2aNZHqpLW1tbIsS0OHDo3sN2zYMFVVVcUqTADAV4TlN3o4yTVDjgUFBZo3b57S0tIi286cOaNAIBBJcJKUkJCg5ubm6AXixPSc6bAwU4T4gs/0s2Dcng+jJEcKdBpzQQid4Tq0ryktLVVdXZ3mz5/fbntSUpJaWlrU1vbX0sMNDQ1KTk7u7RABAJ0IWX6jh5Nc0UMrLi5WbW2tsrOzJZ1PWqFQSAcPHlQgENDhw4c1cuRISVJ1dbXS09NjGS4AwIVckdAKCwvb/fz000+rrq5Oa9eu1fLly/X8889rzZo1ampq0o4dO/Tzn/88RpECAL4qzJBjzz3yyCNqbm7WxIkTNWPGDE2fPl15eXmxDgsAoPN3CjF5OMkVPbSvW7p0aeTfl1xyiTZu3BjDaAAAXXHTdWiuTGgAgL7B6aX3JtwTCQAABuihAQBsC7no5sQkNACAbZSPAQB4QphFIQAAL3BTD809kQAAYIAeGgDANhaFAAA8wU3XoZHQAAC2caeQKDI9t1acM3HEknENLTlwHo3bOzCM4cSJMHp5B14/bHYMX8iwveHrS5LOhWPbXpLCpjGEzGMw5HPiPHic5xIaAKD3uOlu+yQ0AIBtDDkCADzBTdehkdAAALaFXbRs3z2pFQAAA/TQAAC2MeQIAPAEbk4MAPCEEMv2AQBe4KYemnsiAQDAAD00AIBtDDkCADwhWkOOf/nLXzR58mQNHz683faXXnpJqampnbYhoQEAbIv2ra/27NnT432ZQwMAeAI9NACAbdG+2/4vf/lLffjhh0pMTNSCBQt08803d7mv5xJa2PA3stxwRgxLLzlRNclv+CG1DOtomZavkiS/aU02079TB0qJ+QzPg2k9MydqcPlCZh9o0/aSpLZzhjHEvhaZry32MXQmWkOOSUlJmjVrlubMmaOxY8fqv/7rv3THHXfoyiuv1Pjx4ztt44b/vgEAfVS0bk6cmpqqxx9/PPLzuHHjNHnyZJWWlpLQAADOi9a9HE+dOqVTp04pLS3tr68VCikQCHTZhkUhAADX+e///m/Nnj1bR44ckSQdPHhQZWVlys7O7rINPTQAgG3RGnKcOHGiFi1apLvuukuSFB8fr4KCAo0dO7bLNiQ0AIBt4SgO9C1YsEALFizo8f4kNACAbSEXVawmoQEAbIvWkKMdLAoBAHgCPTQAgG1uqodGQgMA2Eb5GACAJzCHBgCAw+ihAQBsYw4NAOAJ0S4fcyFIaAAA27iwGgDgCQw5RlE4YPZtIRznUCAGfIafD58bvjAZfmvzWebVMcMhsxj8bYYxOFDg029Y1NF/1rC4Zpt5cU3TY/ha24xjkGmR0HMOFBk15D9rVqT0YuC5hAYA6D1uWrZPQgMA2MaiEACAJ7iph+ae2TwAAAzQQwMA2MYqRwCAJ7hpyJGEBgCwjUUhAABPcFMPzT2DnwAAGKCHBgCwzU09NBIaAMA2EhoAwBNIaJ2oqKjQ+vXr1djYKMuydNttt2nBggWqr6/XypUrdeDAAfn9fuXm5ur++++XzxV34AWAixurHL/ms88+06JFi1RUVKSMjAwdPXpUM2fO1Pe//31t3rxZKSkpKisrU3Nzs+bNm6fBgwdr7ty5sQ4bAOAirljl6Pf79cQTTygjI0OSNGTIEI0YMUKVlZV6++23lZ+fL5/Pp+TkZM2dO1evv/56jCMGAEjnhxxNHk5yRQ8tNTVVU6ZMifx85MgRHTp0SKNHj5ZlWRo6dGjkuWHDhqmqqqrLY4X7m8USDpi1d4LPsPSSI6OxsR5FcOCD7gsZFiQzPJE+B+qh+UJm9dB8hjW0TNtLDtQzO+tAPbRzZr+HZfg+OKL1bKwj6BRzaN04fvy4Fi5cqLvuuks+n0+BQEB+/187kgkJCWpubo5hhACAL7kpobliyPFLlZWVmjNnjvLy8rRkyRIlJSWppaVFbW1//YbW0NCg5OTkGEYJAHAj1/TQKisrlZ+fr0cffVTZ2dmSpLS0NAUCAR0+fFgjR46UJFVXVys9PT2WoQIAvkAP7WtaW1u1bNmydslMkpKSkvSTn/xEzz//vMLhsBoaGrRjxw7l5eXFMFoAwJcsy2f0cJIremglJSWqq6tTYWGhCgsLI9tzcnL0yCOPaOXKlZo4caL8fr9uvfVWEhoAuATXoX1Nbm6ucnNzu3x+48aNvRgNAKCnGHIEAMBhruihAQD6JqfnwUyQ0AAAtrlpyJGEBgCwjR4aAMAT3NRDY1EIAMAT6KEBAGyzHLgJt1NIaAAA27iwGgDgCSwKiaJQvFn7cMC8/2xaB8uKM6zDZVhPzREu+IyHDc+DZTjD7DtnXkPLf9bslzCuh9biQA0uw3pmlgvqoSkU+z8q3xkXFGuMgYaGBuXk5CgzM1Nr167tdl/PJTQAQO+J9irHgoICJSQk9GhfVjkCAGyzLLNHd9566y0dPXq023v9fhUJDQBgW7TKx5w+fVpr1qzRmjVr5Pf3LFUx5AgAsC1ai0IKCgo0b948paWl9bgNPTQAgKuUlpaqrq5O8+fPv6B29NAAALZFY1FIcXGxamtrlZ2dLen8SsdQKKTDhw9rx44dXbYjoQEAbIvGnUIKCwvb/fz000+rrq6OZfsAgOjhwmoAgCf0RkJbunRpj/ZjUQgAwBPooQEAbHPRzfZJaAAA+5hDAwB4g4u6aMyhAQA8gR4aAMA2hhwBAJ4QjQur7fJcQgslmJ1dy4Eaesbvb9jsCL440wDMi4xahvUUww78DqZVRvs1m716XJsDBT7/r9Wova/R7JewWh0o8GlYXNM660AMhgU6LRf8r+3W+SF6aAAAb3BRQnNr0gcA4ILQQwMA2OaC0dgIEhoAwD4SGgDAC1gUAgDwBhf10FgUAgDwBHpoAADbGHIEAHiDi4YcSWgAAAPu6aExhwYA8AR6aAAA+xhyBAB4Ql9LaE1NTSovL9fBgwd18uRJSdKll16qYDCozMxMfetb34pqkAAAl+orqxw//fRTFRUV6Xe/+51CoZC+/e1vKzU1VZJUX1+vY8eOKS4uTrNmzdLixYt1+eWX90rQAAB36BP3cnz33Xe1fPlyXX311SosLFRGRoZSUlLa7dPQ0KDy8nJt27ZNM2bM0D//8z8rMzMz6kF3J5Rg1t7yu+DdccE3Hp9hTbaQ4Wn0OfA2hBJNj2D2Plz6xybTAKRjnxg1D7eY1VOzQuY13SzDWmSyzGNw1f+6NoXaDIsMXgS6TGgrVqzQunXrNGnSpC4bp6SkaNq0aZo2bZpKS0u1YsUK7du3LxpxAgDcyEXfFbpMaK+++qquuOKKHh/opptu0pgxYxwJCgDQR7hgROlLXSa0ryYzy7JUWlqqw4cPq6WlpcO+S5Ys6dAGAOB9TkwPOKVHqxxXrFihPXv2KBgMKiGh/SSVz+ee7AwA6GV9LaGVlpbqtdde09VXXx3teAAAsKVHCS05OVlDhw6NdiwAgL7GRXNoPbqX48KFC/XUU0/p7Nmz0Y4HANCXWIYPB/WohzZmzBj967/+q7Zt26ZBgwZ1mDf7z//8T2ejAgD0DX1tDu2BBx7Q1VdfrZ///OdKTDS+WhUAAMf1KKH97//+r/793/9d/fv3j3Y8AIC+xEU9tB7NoY0bN07/8z//E+VQAAB9juUzezioRz20GTNm6L777tPkyZM1ePDgDnNot956q6NBAQD6hj53YfU//dM/SZL+/Oc/d3jO5/P1SkI7cOCAVq1apZMnTyoQCOgXv/iFZsyYEfXXBQB0o68ltMrKSsXFxUU7li61trbq7rvv1n333adbbrlF1dXVmjNnjoLBoEaOHBmzuAAA7tGjObRYJjNJqqiokN/v1y233CJJGjFihG688Ubt2rUrpnEBANyjRz20AwcO6Fe/+pWqqqo6vTnxRx995HhgX1VTU6Nhw4a125aWlhb11wUAdK/PzaGtXLlSgwYN0n333af4+Phox9RBc3Nzh5six8fHq7m5ucO+oUTDYoBhF9zGxedAQUNTPeq7d8PwNDpSaLWf2THixnT8fF2IhYt3G7WXpHHxZgU+B/jNRlcCMh+dCfjMjuE3/TBJivOZfqBjr80yLJQaLS669VWPElptba127twZs+vQkpOT1djY2G5bY2OjkpOTYxIPAOALLuqh9ehry5AhQ3TmzJlox9Kl9PR01dbWKhz+a8+lurpa6enpMYsJAOAuPUpoDz/8sFavXq2amhqdO3dO4XC43SPafvSjHykxMVGvvvqqpPNzdvv379fMmTOj/toAgG5E8ebE+/bt06xZszR16lRNnTpVW7Zs6Xb/Hg053nPPPWpoaFBxcXGnz0d7cUb//v1VVFSklStX6te//rXi4+O1du1aDR8+PKqvCwDoXrQWhXzyySdaunSpnn32WV177bU6evSoZs6cqe9973u6/vrrO23T44TWr1+Pdo2aUaNGaefOnTGNAQDwNVFKaD6fT08++aSuvfZaSeenvoYPH66ampoLT2gHDx7UqFGjJEk/+9nPehTAxx9/rGAweKFxAwDQzmWXXabs7OzIzxUVFaqrq1NmZmaXbbqcQ1uwYIFefPFFhULfvFQ0HA7rxRdf1IIFCy4sYgBA3xblAp/vvPOObrjhBi1btkyPPPKIhg4d2uW+XfbQtm7dqiVLluj555/XbbfdpszMTAWDQfn953NgOBzWoUOH9N577+nll1+WpG+csAMAeEu0L6yeOHGi9u3bp5qaGi1cuFChUEi5ubmd7ttlQgsGg3rjjTe0detWbd68WevXr5fP59OAAQMknb8OzLIsXXrppZo/f77mz5/f4eJnAIDHRenC6sOHD6u2tlaTJk2SJH33u9/V5MmTVVpaeuEJTZISEhKUn5+v/Px8VVZW6tChQ6qvr5ckpaamKhgMavTo0c7+FgCAviNKPbTGxkYtX75cL7/8soLBoBoaGrRv375uL9fq8dLFMWPGaMyYMY4ECgBAd6655ho9/PDDuv/++9XW1qZwOKzJkyd3u1YjtmvxAQB9WjTn0H7605/qpz/9aY/3J6EBAOxz0b0cSWgAANvcVD6mR/dyLCsrk2W5KGoAAL6mRz20RYsW6bLLLtMtt9yiWbNmaciQIdGOyzYrwfBmyW0uqO1jXIvMgS8fcWbH8JnWIguY3/Q6PuGsUfurBp42an/gTNcXgPaWwf1OGbUf4G8zjiHBZ/Y+JPnM/yYDhvXQ4hyoyWaqMXzOqP1gh+LowEV9nR730O6880699957mjZtmv7u7/5Or7/+us6eNfugAgD6uCjfKeRC9CihXXrppbr99tu1fft2/f73v1dmZqaeffZZ3XDDDVq9erWqqqqcjQoA0Cf4LLOHky64H37VVVcpKytLWVlZkqQ33nhDeXl5uvfeeztUlQYAoLf0OKE1NTXp5ZdfVl5enm6++WZ98MEHevDBB7Vv3z69+eabamxs1IoVK6IZKwAAXerRopD7779fe/fuVXx8vGbOnKl169bp6quvjjx/xRVXqLCwUBMnToxaoAAAF3LRopAeJbQTJ05o9erVmjZtmvr379/pPgMHDtTixYsdDQ4A4G5uug6tRwlt27ZtPTrYP/zDPxgFAwDoY/paQgMAoFMuSmiml/ACAOAK9NAAALb1uTk0AAA6RUIDAHiBm3pozKEBADyBHhoAwD4X9dBIaAAA+0ho0eNPMKsZFLYcOCWmpZNMB6UNa5lJkt+wHllcnFn7QH+z91GSkg3roSX2M6sF1hhKMGovScfPDTQ+homBcc3Gx0jxtRq1b/GbfxYSfCGj9gHjCMydCJlFEa16aG6aQ/NcQgMA9CIXJTQWhQAAPIEeGgDAPhf10EhoAADbmEMDAHgDCQ0A4AVu6qGxKAQA4An00AAA9rmoh0ZCAwDYR0IDAHiB6Y2RnMQcGgDAE+ihAQDsY8gRAOAFblq2T0IDANhHQgMAeIKLEhqLQgAAnuC5Hlog3qwYYOu52Od4n2GBTr8TBT7jDAsiBszaxwccKOrYz+wY/fxmv0Nr2PzPq8mwSGijP9E4BlNhv9nfVEgtxjG0mRb49JkVrHXCyXBSrEPoFHNoAABvIKEBALyAHhoAwBtclNBiP2EEAIAD6KEBAGxjyBEA4A0kNACAJ0QxoVVUVGj9+vVqbGyUZVm67bbbtGDBgi73J6EBAFzns88+06JFi1RUVKSMjAwdPXpUM2fO1Pe//31dd911nbZhUQgAwDafZfboit/v1xNPPKGMjAxJ0pAhQzRixAh9/PHHXbahhwYAsC9KQ46pqamaMmVK5OcjR47o0KFDuvbaa7tsQ0IDANjms6K/KuT48eNauHCh7rrrLo0aNarL/Vwx5GhZljZu3KipU6dq6tSpmjNnjj788MPI82VlZZo+fbqysrKUk5Oj8vLyGEYLAIiwDB/foLKyUnPmzFFeXp6WLFnS7b6u6KH927/9m0pKSrRz506lpKToueee07333quSkhJ98sknWrZsWWRi8P3339eiRYtUUlKi1NTUWIcOAIiSyspK5efn69FHH1V2dvY37u+KHtrYsWO1bt06paSkSJImTZqkI0eO6MyZM9q7d6+CwWBkYnD8+PFKT09XSUlJLEMGACh6i0JaW1u1bNmyHiczySU9tGuuuabdz2+++abGjBmjxMRE1dTUaPjw4e2eT0tLU1VVVW+GCADoTJSm0EpKSlRXV6fCwkIVFhZGtufk5Gjp0qWdtum1hLZ7926tXr26w/YBAwZo7969kZ+Li4u1detWbdmyRZLU3Nys+Pj4dm3i4+PV3Nzc6ev0729WA6vtbOxzvM9v9gmJizOv3dSvn1n9qP6G7Z2ohxZvWA/Nb3hPn7DlM2ovSS3hgFH7RsN6al7R5mszap/gM/88mjoVSo51CJ2K1q2vcnNzlZube0Fteu1/75ycHOXk5HS7zzPPPKPt27dry5YtkZUsSUlJ+vzzz9vt19jYqEGDBkUtVgBAD7no1leumEOTpA0bNuj3v/+9XnnllXbLMkeOHKnDhw+327e6ulrp6em9HSIAwMVckdD27dunN954Q5s3b9Zll13W7rns7GzV1tZGluqXlZXp2LFj7S64AwDERrQWhdgR+wkjSZs3b1ZjY6Nmz57dbvuGDRs0atQobdq0SQUFBWpqalJKSoqKioo0cODA2AQLAPgrFw05uiKhPffcc90+n5GRoV27dvVSNACAnnJTPTRXDDkCAGDKFT00AEAf1Qv3cuwpEhoAwDY3DTmS0AAA9pHQAABe4DO/MZFjWBQCAPAEemgAAPsYcgQAeAGLQgAA3sCyfQCAF9BDi6J4wzpcLQ7UEjPlN4zBkXpohscwrWcWH2f2Pp4/hlkM/QyXb4Ut8zVXzeH+Ru0DYbN6aKY14SQpTrH/mwr5zWrTnVWcQ5HY12j4Xl4MPJfQAAC9iB4aAMALGHIEAHiDixaFcGE1AMAT6KEBAGxjyBEA4A0kNACAF9BDAwB4Q9g9GY1FIQAAT6CHBgCwzz0dNBIaAMA+5tAAAN7gogurSWgAANvc1ENjUQgAwBPooQEA7HNRD42EBgCwzcccWvT072dW1DGuX+yLEfr9ZjGYFueUpP6GhVL7Gf4OAQcKfJoW6DQtbhmWWVFJSWqzzApLtoXN/sRbfQGj9pLU4jP7m/RbDvxNGh4iwddmHoOh/wvHxzqEzsX+v8wI5tAAAJ7guR4aAKD3MOQIAPAG9+QzEhoAwAA9NACAF3BhNQAADqOHBgCwjyFHAIAXGF7u6SgSGgDAPhf10JhDAwB4AgkNAGCfZfj4Bjt27NAPfvADPffcc9+4L0OOAADbonmnkFWrVqm+vl7f/e53e7Q/PTQAgH2WZfboRk5Ojp566iklJyf3KBR6aAAA+6K4yvG66667oP3poQEAPMFzPbQ4w/uwxDlQS8x0SNk0hjjDWmROHMO0HlqcAxe3+GN8gUzIMq+HFjY8hnE9NcP2ktQSNqup5negCxDnd8/ScrtaDc9jtHC3fQCAN5DQAACeQEIDAHhClEb2Q6GQcnJyJEnHjh1TVVWVdu7cqSlTpui+++7rtA0JDQDgOnFxcdqzZ88FtSGhAQBsY1EIAMAbSGgAAE9wUULjwmoAgCfQQwMA2OeiHhoJDQBgHxWrAQBe4KZVjq6bQysuLlYwGNT+/fsj28rKyjR9+nRlZWUpJydH5eXlMYwQABARxfIxF8pVPbSTJ0/qqaee0sCBAyPbPvnkEy1btkxFRUXKyMjQ+++/r0WLFqmkpESpqamxCxYA4Cqu6qGtWrVKd9xxR7tibnv37lUwGFRGRoYkafz48UpPT1dJSUmswgQAfClsmT0c5JqEVlxcrKamJs2ePbvd9pqaGg0fPrzdtrS0NFVVVfVmeACAzlyMQ467d+/W6tWrO2wfMGCAduzYofXr12vbtm0dnm9ublZ8fHy7bfHx8Wpubo5KnG6Y37QcqKMFbzCthxaWWfuQYXuvMK3PF7LM+w6xru/XJTf8p/mFXktoOTk5kTsnf93SpUuVn5+v73znOx2eS0pK0ueff95uW2NjowYNGhSVOAEAF+BiTGhdaWpqUkVFhSorK/Wb3/xGknT8+HEtX75cd955p0aOHKnt27e3a1NdXR2ZUwMAQHJBQvvWt76lDz74oN22m266SY8//riuv/561dfXa926dSovL9ePf/xjlZWV6dixY5oyZUqMIgYARDi8sMNEzBPaN0lNTdWmTZtUUFCgpqYmpaSkqKioqN3SfgBAjFjumdtzZUIrLS1t93NGRoZ27doVo2gAAF1y0Ryaa5btAwBgwpU9NABAH8EcGgDAE1w05EhCAwDYR0IDAHiCixIai0IAAJ5ADw0AYF+Y69AAAF7goiFHEhoAwD4SGgDAE1x0HRqLQgAAnuC5HtrJpiSj9mcaExyKxD5fnNkkqy/O/BtTv34ho/bx/c8ZtU+KjzNqL5kXSu0fZ3YOAg4UZEzp12LU/orAaaP2qf2ajNpL0kC/WTHeFL/ZOZCkJH+bUfsEn9lnIeThor0WNycGAHiCi4YcSWgAAPtctCiEOTQAgCfQQwMA2MeF1QAAT3DRkCMJDQBgm0UPDQDgCS7qobEoBADgCfTQAAD2cR0aAMATuFMIAMALLHpoAABPcFEPjUUhAABPoIcGALCNIUcAgDe4aMjRZ1kuuioOAACbmEMDAHgCCQ0A4AkkNACAJ5DQAACeQEKLguLiYgWDQe3fvz+yraysTNOnT1dWVpZycnJUXl4ewwijw7Isbdy4UVOnTtXUqVM1Z84cffjhh5HnL4ZzIEkHDhxQXl6eJk2apClTpuiNN96IdUi9oqKiQj/72c80bdo0TZ06VVu2bJEk1dfXa/HixZowYYImTpyoJ598Ul5fi9bQ0KAJEybogQcekHRxnoOYsOCozz77zPrJT35ijR8/3qqoqLAsy7JOnDhh/eAHP7DKy8sty7Ks/fv3W+PGjbNOnjwZy1Ad98ILL1i5ubnW6dOnLcuyrGeffdbKzs62LOviOQctLS1WZmam9dprr1mWZVlVVVXWddddZ3388ccxjiy6Pv3003bv75EjR6wf/vCH1gcffGDdfffd1gMPPGCFw2GrqanJmjlzpvXiiy/GOOLo+sd//EcrOzvbWrFihWVZ1kV5DmKBHprDVq1apTvuuEPJycmRbXv37lUwGFRGRoYkafz48UpPT1dJSUmswoyKsWPHat26dUpJSZEkTZo0SUeOHNGZM2cumnNQUVEhv9+vW265RZI0YsQI3Xjjjdq1a1eMI4suv9+vJ554IvL+DhkyRCNGjFBlZaXefvtt5efny+fzKTk5WXPnztXrr78e44ij56233tLRo0eVm5srSWpqarrozkGskNAcVFxcrKamJs2ePbvd9pqaGg0fPrzdtrS0NFVVVfVmeFF3zTXXaPTo0ZGf33zzTY0ZM0aJiYkXzTmoqanRsGHD2m3z4u/5dampqZoyZUrk5yNHjujQoUMaPXq0LMvS0KFDI88NGzbMs+fj9OnTWrNmjdasWSO///x/r7W1tRfVOYgl7hRyAXbv3q3Vq1d32D5gwADt2LFD69ev17Zt2zo839zcrPj4+Hbb4uPj1dzcHLVYo6W7c7B3797Iz8XFxdq6dWtkHsVL56A7zc3NSkhIaLfNi79nd44fP66FCxfqrrvuks/nUyAQiPznLkkJCQmePR8FBQWaN2+e0tLSItvOnDlzUZ2DWCKhXYCcnBzl5OR0+tzSpUuVn5+v73znOx2eS0pK0ueff95uW2NjowYNGhSVOKOpu3PwpWeeeUbbt2/Xli1bNGrUKEneOgfdSU5OVmNjY7ttjY2N7YagvayyslKLFy/W7bffrjvvvFN/+tOf1NLSora2NgUCAUnnF0x48XyUlpaqrq5Oa9eubbc9KSnpojkHsUZCc0BTU5MqKipUWVmp3/zmN5LOf0tdvny57rzzTo0cOVLbt29v16a6ujoy3+AlGzZsUFlZmV555RVddtllke0XyzlIT0/XM888o3A4HPlGXl1drZEjR8Y4suirrKxUfn6+Hn30UWVnZ0s6P9waCAR0+PDhyDmorq5Wenp6LEONiuLiYtXW1kZ+94aGBoVCIR08ePCiOQcxF+tVKV6VlZUVWeV48uRJa9y4cda7775rWZZlvfPOO9aPfvQj6/PPP49hhM77wx/+YN10003WqVOnOjx3sZyD1tZWKysry3rllVcsy7KsP/3pT9a4ceOsmpqaGEcWXS0tLdbkyZOtkpKSDs/de++91ooVK6xQKGSdPn3amjZtmvW73/0uBlH2ro0bN0ZWOV6s56C30UPrBampqdq0aZMKCgrU1NSklJQUFRUVaeDAgbEOzVGbN29WY2Njh0UxGzZs0KhRoy6Kc9C/f38VFRVp5cqV+vWvf634+HitXbu2w4IYrykpKVFdXZ0KCwtVWFgY2Z6Tk6NHHnlEK1eu1MSJE+X3+3XrrbcqLy8vhtH2Ps5B7+Bu+wAAT2DZPgDAE0hoAABPIKEBADyBhAYA8AQSGgDAE0hoAABPIKEBBj7++GNdf/31+uijjxw/9vHjx5WZmak//OEPjh8b8CISGmDTmTNntHjxYi1evFjf+973HD/+lVdeqccee0z333+/Tpw44fjxAa8hoQE2vfDCCwqHw7rtttui9hpZWVkaMWKEioqKovYagFeQ0ABJ//Ef/6Ef/vCH+stf/hLZ9uCDD2rWrFlqa2vrsL9lWdq6davmzJkTuYP6zp079eMf/1jvv/++ZsyYobFjx+rmm2/WgQMHJEnnzp1TMBjUa6+9pr//+7/X2LFjNW3aNB04cEAvvfSSbrzxRo0bN06//OUvFQqFIq81f/58vfrqqx3u4g+gPRIaIGnmzJmaMGGCVq1aJUnav3+/du/erSeffDKSsL6qsrJSJ0+eVGZmZmRbXFycGhoatH37dm3dulX79+/XoEGDIsfs1+/8rVOff/55rVixQu+//76GDBmiZcuW6ejRo3rzzTe1Y8cO7dmzR2+99VbkuBkZGQqFQnr33XejeQqAPo+EBnxh1apV+uijj/Taa6/p4Ycf1vLlyzVixIhO9z106JB8Pp+CwWC77W1tbbr77ruVmpqqhIQEZWVldahMfNNNN2nUqFGKj4/XpEmT9Omnn+qee+5RQkKC0tPTFQwGVVNTE9n/kksu0be//W0dPHjQ+V8a8BASGvCFv/mbv9Fjjz2mBx98UJdffrnmz5/f5b719fVKSUmJ9Lq+aujQoZF/JyQkqKWlpd3zgwcPjvw7MTFRgwYNalfNOzExUa2tre3apKamqr6+/oJ/J+BiQkIDvqK2tlZJSUk6duyYmpubbR3jy8KePX3+m/YH0DP8JQFfqKmp0YYNG/Tb3/5WV155pdasWdPlvqmpqWpoaNC5c+d6Jbb6+nqlpqb2ymsBfRUJDZAUCoW0YsUK5eXl6brrrtNjjz2mXbt2qbS0tNP9g8GgLMvSn//856jHdvr0aR07dqzDfB2A9khogKRnnnlGJ0+e1PLlyyVJaWlp+sUvfqGHHnqo07mr0aNHa9CgQSovL496bO+9957i4uJ0ww03RP21gL6MitWATb/97W/1yiuvaM+ePZ0uDnHK3LlzNWLECP3qV7+K2msAXkAPDbDp9ttvlyRt3749aq9RWlqqqqoqLVmyJGqvAXgFCQ2wKTExUZs2bdLTTz8dlWvEjh8/roceekiFhYW64oorHD8+4DUMOQIAPIEeGgDAE0hoAABPIKEBADyBhAYA8AQSGgDAE0hoAABP+P+ATqgE2jHSbwAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "f.plane(\"z\").angle.mpl_scalar(cmap=\"viridis\", filter_field=f.norm)" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAb8AAAF8CAYAAABfb7VoAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAAArJ0lEQVR4nO3dfXCU9b3//9cGk82NRBqjIgoJhmSj9GhLQSYgkFhSGALNxMgIODC2ticasJSbHqw1WhA4WKHiXU71qEA5VhSxViCDQFFToQn2255fxhDAJBgQQZFYsukSbpLr9wenW9eQEPeTZffK9XzM7Az57PW59s0G8sr7uq69Pi7LsiwBAOAgUeEuAACAi43wAwA4DuEHAHAcwg8A4DiEHwDAcQg/AIDjEH4AgIj1yiuv6Fvf+pZeeOGFDrd54403lJubq5ycHBUWFmrPnj0X3O8l3VkkAADdZeHChWpsbNR1113X4TZ79uzRI488opdfflkZGRn6wx/+oOLiYm3dulUxMTEdzqPzAwBEpLy8PD3xxBNKSEjocJvNmzcrOztbGRkZkqT8/Hy1tbVp9+7dne6b8AMARKShQ4decJu6ujoNHDgwYCwlJUW1tbWdziP8AAC2dfLkSbnd7oAxt9stn8/X6bwed87P5XKFuwQAiDghu42z62qz+dYRo+nx8fFqbm4OGPN6vZ0eKpV6YPgBAC6mji8quRgyMjJ04MAB/9etra366KOP/OcAO8JhTwCAbU2aNEk7d+7Uvn37JEnr169XYmKihg0b1uk8V09b0ojDngDQXugOe6abzbc+PO9wa2ur8vLyJElHjhxRfHy8LrvsMuXm5kqS4uLiVFxcLEkqKyvT448/rtOnT+vqq6/W4sWLNWjQoM7LJvwAoOcLXfgNNptvVXdPHV8T5/wAAAbCe84vWJzzAwA4Dp0fAMBAdLgLCArhBwAwYM/DnoQfAMAA4QcAcBx7hh8XvAAAHIfODwBgwJ6dH+EHADBA+AEAHIfwAwA4jj3DjwteAACOQ+cHADBgz86P8AMAGCD8AACOY8/w45wfAMBx6PwAAAZY1QEA4Dj2POxJ+AEADBB+AADHsWf4ccELAMBx6PwAAAbs2fkRfgAAA/YMv4g77NnU1KRRo0bp/vvvlyQ1NjZq5syZGjVqlMaMGaPHHntMlmWFuUoAwDkxho/wiLjOb8mSJYqNjfV/XVJSosTERJWXl8vn8+nOO+/UNddco2nTpoWxSgDAOXR+xt5++20dOnRIEydOlCQ1NzfrnXfeUVFRkVwulxISEjRt2jS9+eabYa4UAGBnERN+J06c0NKlS7V06VJFRZ0rq6GhQZZlacCAAf7tUlJSVFtbG64yAQABOOxpZMmSJbrzzjuVmprqHzt58qSio6P9YShJsbGx8vl8YagQF5Ur/Ptwmc43my5JMj69bTif0+u4MA57Bm3Hjh06fPiwZsyYETAeHx+vlpYWnTlzxj/W1NSkhISEi10iAOC8og0f4RERnV9ZWZkaGho0duxYSecCrrW1VXv37lV0dLQOHDigjIwMSVJdXZ3S09PDWS4AwOZcVgR+buCpp57S4cOHtWzZMs2dO1cxMTFaunSpmpubdccdd+hHP/qRCgsLzzvXZXqsCpGBw56SOOyJ7hOyH/WuLWbzrfHdU8fXFBGHPTvz8MMPy+fzacyYMZo0aZImTJjQYfABAC42e17wEpGdnwk6vx6Czk8SnR+6T+g6v/fM5lu3dE8dX1NEnPMDANgVV3sCAGALdH4AAAP27PwIPwCAAcIPAOA4hB8AwHHsGX5c8AIAcBw6PwCAAXt2foQfAMAA4QcAcBjLMEbCdU8uwg/nZ/ov0vBssqsbVjq5JNZsfrRhDb3MpkuSrDMX3qYzZ1rM5p81fH1Jamszm88t1hAKhB8AIGhthjHSHb8kBoPwAwAEjfADADiOafiFiz2rBgBEhLOGJ/jd3VTH18WH3AEAjkPnBwAI2mnD+QndUsXXR/gBAIJmGn7hQvgBAIJG+AEAHMeu4ccFLwAAx6HzAwAEjc4PAOA4pw0fnamqqlJhYaGys7OVm5urjRs3nne7l19+WRMmTND48eN12223qby8/IJ10/kBAILWDfc+P69Tp06puLhY8+bNU0FBgerq6jRlyhR5PB5lZGT4t/vLX/6iFStW6M0331S/fv1UUVGhoqIivfvuu+rTp0+H+6fzAwBEnIqKCkVFRamgoECSlJaWptGjR2vTpk0B2+3fv18DBw5Uv379JEnDhw/XmTNndOjQoU73T/gBAIIWqsOe9fX1SklJCRhLTU1VbW1twNjw4cPV0NCg/fv3S5L++Mc/Kjk5Wenp6Z3WzWHPnqg7Voc0XU/vUrPpl/QzfH1JfTr/t39B/a41m59sNl2S1Pqx2fzPPjSbf/Sw2XxJavaazT/Tal4DawKGTqguePH5fIqNDVyU0+12y+fzBYylpaVp9uzZKigo0GWXXaaWlhYtX7683dyvIvwAAEELVfglJCTI6w38zcnr9SohIfCGaO+++65eeOEFbdmyRf3791dNTY1+8IMf6Pnnn9c3v/nNDvfPYU8AQNBCddgzPT1dDQ0Namtr84/V1dW1O5z57rvvavjw4erfv78k6frrr1dGRoYqKys7rZvwAwBEnGHDhikuLk4bNmyQJNXU1KiyslL5+fkB22VmZur9999XY2OjJOnIkSPat2+frr/++k73z2FPAEDQQnXYMyYmRqWlpSopKdEzzzwjt9utZcuWaeDAgVqxYoXi4uJUXFysyZMn65NPPtHUqVPlcrkUFRWlWbNmacSIEZ3u32VZPetUsMvVHVd72FwPuOAlmgteJHHBi8QFL90lVD/q1xn+vJkSpu8NnR8AIGh2vb0Z4QcACJpdw48LXgAAjkPnBwAIml07P8IPABC0UN3YOtQIPwBA0Oza+XHODwDgOHR+AICg2bXzI/wAAEEj/AAAjkP4AQAch/BD9zG9N2d0N1zH9I1eZvOHmM1vnd5i9vqS0qb2Npo/VZ3fGPdCbjSafc512mU0f+cqsxtrrn3RaLok6a8fmM3/wvDeoJL5/UG5NWjPQ/gBAIJG5wcAcBzCDwDgOHYNPz7kDgBwHDo/AEDQ7Nr5EX4AgKARfgAAx2FVBwCA49i18+OCFwCA49D5AQCCZtfOj/ADAASN8DNUUVGhFStWyOv1yrIsTZ06VXfddZcaGxtVUlKiqqoqRUVFaeLEiZo/f75cLtMbYAIATBF+Bj7//HPde++9Ki0tVVZWlg4dOqT8/Hx985vf1KpVq5SYmKjy8nL5fD7deeeduuaaazRt2rRwlw0AjmfX8IuIC16ioqL0q1/9SllZWZKk/v37Ky0tTdXV1XrnnXdUVFQkl8ulhIQETZs2TW+++WaYKwYA2FlEdH5JSUnKzc31f33w4EHt379fN9xwgyzL0oABA/zPpaSkqLa2NhxlAgC+wq6dX0SE35cdPXpU99xzj3784x/L5XIpOjpaUVH/alBjY2Pl8/nCWOFFEGXYkF8aY17DkG+YzV90pdH0GUOzzV5f0lTdajTfo/5G8y8zmn3OpTpkNH/q+N8bzU8/sdpoviQ91Gg2v7LeuASdaDOb38qCfh2ya/hFxGHPf6qurtaUKVNUWFioWbNmKT4+Xi0tLTpz5l/3EGhqalJCQkIYqwQA/NNpw0e4REz4VVdXq6ioSA8++KDuvvtuSVJqaqqio6N14MAB/3Z1dXVKT08PV5kAgB4gIsLv1KlTmj17tn75y19q7Nix/vH4+Hh973vf04svvqi2tjY1NTXplVdeUWFhYRirBQD8k107v4g457dt2zYdPnxYy5cv1/Lly/3jeXl5evjhh1VSUqIxY8YoKipKt99+O+EHABHCruf8IiL8Jk6cqIkTJ3b4/JNPPnkRqwEAdBWrOgAAHMeunV9EnPMDAOBiovMDAAStzaYpYtOyAQCRgPADADgO4QcAcBy7hh8XvAAAHMemmQ0AiAjdcB/9cCD8AADBI/wAAI5D+MHP5TKbH93LbH6/JLP5kjT9ZqPpd30nz2h+kZVlNF+S0lxmawr2tuKN5kcb/jOQpF5t1xrNt2LN3oPvDDaaLkmaft1qo/mHDpvX4DO8DUmb4Xp+LAcYeQg/AEDw6PwAAI4THe4CgkP4AQCCR+cHAHAcm4YfH3IHADgOnR8AIHg27fwIPwBA8Ag/AIDjEH4AAMexafhxwQsAwHHo/AAAwbNp50f4AQCCZ9Pw47AnACB4MYaPTlRVVamwsFDZ2dnKzc3Vxo0bz7tdTU2NJk+erLFjx2rChAnavn37Bcum8wMARJxTp06puLhY8+bNU0FBgerq6jRlyhR5PB5lZGT4t/P5fCouLtbSpUuVlZWlv/zlL3riiSeUnZ2tSy7pOOLo/AAAwQtR51dRUaGoqCgVFBRIktLS0jR69Ght2rQpYLsdO3Zo0KBByso6twza0KFDtXbt2k6DT6LzAwCYCNGqDvX19UpJSQkYS01NVU1NTcBYdXW1rrzySs2ZM0d79uxRcnKy5syZo6FDh3a6f8IvFEwXs411m81PH2g2X1LvqdOM5hfoO0bz02S2CKskXSaz9/ESl9miwt2wlq1Mf7K4YtLMXj65wGy+pKnuDUbz17u8xjV8bDjfcC3cnr2YbYguePH5fIqNjQ0Yc7vd8vl8AWNNTU169913tWbNGl133XV67bXXdO+99+qtt95SUlLHC3tz2BMAELwQHfZMSEiQ1xv4i4/X61VCQkLA2GWXXaZbbrlFaWlpcrlcmjx5stxut/72t791WjbhBwCIOOnp6WpoaFBbW5t/rK6uTunp6QHbpaSk6O9//3vAWFRU1AXP+RF+AIDghajzGzZsmOLi4rRhw7nD5jU1NaqsrFR+fn7AdhMmTNDf/vY3VVVVSZK2b9+ulpYW3XTTTZ2WzTk/AEDwQnTOLyYmRqWlpSopKdEzzzwjt9utZcuWaeDAgVqxYoXi4uJUXFys3r176/HHH9cvfvELtbS0qHfv3nr66afVp0+fTvdP+AEAghfCO7xkZmZq/fr17cbnzZsX8PWIESM6/AB8Rwg/AEDwuL0ZAAD2QOcHAAieTTs/wg8AEDzCDwDgODYNP875AQAch84PABA8m3Z+XQq/5uZm7dq1S3v37tXx48clSZdffrk8Ho9GjhypSy+9NKRFAgAiVIhWdQi1TsPv2LFjKi0t1WuvvabW1lZdffXV/rtkNzY26siRI+rVq5duu+02zZw5U1deaX4nfgCAjfS0zm/nzp2aO3euBg0apOXLlysrK0uJiYkB2zQ1NWnXrl1au3atJk2apF//+tcaOXJkyIsGAESInhZ+CxYs0KOPPqrs7OwOJycmJmr8+PEaP368duzYoQULFui9994LRZ32YrqeX3TshbfpzLWZZvMljdAoo/n91cdofnw3HEu5xPB6LtP1+FzdsIibZXha3hUVb1ZAdH+z+ZKkEUazr/W+ZVxBdNuFt+lM96zNiEjS4f+sDRs26Kqrruryjm699VYNHjy4W4oCANhET+v8vhx8lmVpx44dOnDggFpaWtptO2vWrHZzAAAO0NPC78sWLFigLVu2yOPxtFtW3mV6iA8AYF89Ofx27Nih119/XYMGDQp1PQAAO7Fp+HXpioCEhAQNGDAg1LUAAHBRdCn87rnnHj3xxBM6ffp0qOsBANhJjOEjTLp02HPw4MH67//+b61du1bJycntzvP98Y9/DElxAIAIZ9PDnl0Kv/vvv1+DBg3Sj370I8XFxYW6JgCAXfTk8Pvkk0/0xhtvKCbGpn9LAAC+pEvn/IYMGaKPPvooxKUAAGwn2vARJl3q/CZNmqR58+bpu9/9rq655pp25/xuv/32kBQHAIhwNj0g2KXw+/nPfy5J+vDDD9s953K5Lkr4VVVVaeHChTp+/Liio6P1k5/8RJMmTQr56wIAOtGTw6+6ulq9evUKdS0dOnXqlIqLizVv3jwVFBSorq5OU6ZMkcfjUUZGRtjqAgDHs2n4demcXziDT5IqKioUFRWlgoICSVJaWppGjx6tTZs2hbUuAIA9danzq6qq0iOPPKLa2trz3ti6pqam2wv7svr6eqWkpASMpaamhvx1AQAX0KvVdAfdUsbX1aXwKykpUXJysubNmye32x3qmtrx+Xztbqjtdrvl8/kuei0AgC8zvfNXeD473qXwa2ho0Pr168P2Ob+EhAR5vd6AMa/Xq4SEhLDUE3KW4UoZPsMFTCXpaF+z+X0MXz+6G1aCjTJ7H80XLOmGv4OpNrMFfXWmG65Fb77UbH5rd3QGpt0JOmbP8OvS/4z+/fvr5MmToa6lQ+np6WpoaFBb27+WY66rq1N6enrYagIASOfCz+QRHl0Kv4ceekiLFi1SfX29zp49q7a2toBHqA0bNkxxcXHasGGDpHPnGCsrK5Wfnx/y1wYA9DxdOuz505/+VE1NTSorKzvv86G+8CQmJkalpaUqKSnRM888I7fbrWXLlmngwIEhfV0AwIXYc7WfLoffJZd0adOQyczM1Pr168NaAwDgq3pY+O3du1eZmZmSpMmTJ3dpZ/v27ZPH4+meygAANmDP8OvwnN9dd92ll156Sa2tF75Kqq2tTS+99JLuuuuu7qwNAICQ6LDzW7NmjWbNmqUXX3xRU6dO1ciRI+XxeBQVdS4v29ratH//fv35z3/Wyy+/LElavXr1RSkaABApzoS7gKB0GH4ej0cbN27UmjVrtGrVKq1YsUIul0u9e/eWdO5zdpZl6fLLL9eMGTM0Y8aMdh9EBwD0dPY87NnpVSyxsbEqKipSUVGRqqurtX//fjU2NkqSkpKS5PF4dMMNN1yUQgEAkagHht+XDR48WIMHDw5lLQAA27Fn+Bne+wgAAPsJ74f3AAA2Z8/Oj/ADABiwZ/h16bBneXm5LCsC7lAPAIgwPfjG1vfee69ycnL0xBNP6NChQ6GuCQCAkOrSYc/y8nKVlZVp8+bNeu655zR06FAVFhZq/PjxYVvjL6KZdskthh8a3fux2XxJu/5kNv/Qd8zmX3ul8WJ6ijZcd/kSw2XkXDL/O7jOGv5b8pmtumIdMv8As2uX2fyP1du4hjNqMppvKfSr19hXDz7sefnll2v69Olat26d3nrrLY0cOVLPP/+8brnlFi1atEi1tbWhrhMAEJF68GHPL7v22muVk5OjnJwcSdLGjRtVWFioOXPmtFttHQDQ0/Xw8GtubtbLL7+swsJCff/739f777+vBx54QO+99562bt0qr9erBQsWhLJWAEDEsWf4demc3/z587V9+3a53W7l5+fr0Ucf1aBBg/zPX3XVVVq+fLnGjBkTskIBAOguXQq/Tz/9VIsWLer0Apc+ffpo5syZ3VocACDS9bBVHb5s7dq1XdrZv//7vxsVAwCwG3te7ckdXgAABgg/AIDj2DP8WNUBAOA4dH4AAAP27PwIPwCAAXuGH4c9AQAGQvch96qqKhUWFio7O1u5ubnauHFjp9v/7//+r66//nq9/vrrF6yazg8AEHFOnTql4uJizZs3TwUFBaqrq9OUKVPk8XiUkZFx3u1LSkrUt2/fLu2fzg8AYCA0nV9FRYWioqJUUFAgSUpLS9Po0aO1adOm826/cuVK5eTk6Nprr+1S1YQfAMBAaMKvvr5eKSkpAWOpqannXUXor3/9q3bu3KlZs2Z1uWoOe4aC6Xp+p1rM5h/40Gy+JO/vXjaa//sTU4zmX5llvhZe2pVm83vHm82PNpsuSYrymc131ZldjOD6/edmBUh62XvKaP6HMvz/JMnwf1Q3VNCTheaCF5/Pp9jY2IAxt9stny/wP0VLS4tKSkq0bNmyr7W+LOEHADAQmnt7JiQktFsmz+v1KiEhIWBs5cqVGjt2rP7t3/7ta+2f8AMARJz09HQ9++yzamtrU1TUuTN0dXV17S522bp1q9ra2vxXgh47dkz79+/X3r179cADD3S4f8IPAGAgNIc9hw0bpri4OG3YsEGTJ09WTU2NKisr9bOf/Sxgux07dgR8PX36dBUUFOi2227rdP+EHwDAQGjCLyYmRqWlpSopKdEzzzwjt9utZcuWaeDAgVqxYoXi4uJUXFwc9P4JPwCAgdDd4SUzM1Pr169vNz5v3rwO53R1CT7CDwBggNubAQBgC3R+AAAD9uz8CD8AgAHCDwDgOPYMP875AQAch84PAGDAnp0f4QcAMED4AQAch/ADADhOaFZ1CDUueAEAOA6dXyiYLmbbavibVOMnZvMlqbJr98fryOoD6Wav/8FQs/mSpt5q9n3w9DdbUPcyo9nnXHrIrIZevzdbxvX/ra42mi9Ja1VvNP8TmS2GK0lnDJejZTHbznDYEwDgMFE6G+4SgkL4AQCCRvgBABzHruHHBS8AAMeh8wMABM2unR/hBwAIWky4CwgS4QcACJpdw49zfgAAx6HzAwAEza6dH+EHAAga4QcAcJzocBcQpIg452dZlp588kmNGzdO48aN05QpU/TBBx/4ny8vL9eECROUk5OjvLw87dq1K4zVAgD+KcbwES4REX7/8z//o23btmn9+vV66623lJubqzlz5kiSPvvsM82ePVslJSV6++239fDDD+u+++5TY2NjmKsGANhVRITfTTfdpEcffVSJiYmSpOzsbB08eFAnT57U9u3b5fF4lJWVJUm6+eablZ6erm3btoWzZACA7Nv5RcQ5vxtvvDHg661bt2rw4MGKi4tTfX29Bg4cGPB8amqqamtrL2aJAIDz4IKXC9i8ebMWLVrUbrx3797avn27/+uysjKtWbNGq1evliT5fD653e6AOW63Wz6fL6T1hlVbm9n8083mNXz+V7P5XzxkNP23/990s9eXtPfFKUbzp44we/0bL7zJBV23y2w9v53eLUbz12qD0XxJ+qvM1pds7oaVwttYkS9kCL8LyMvLU15eXqfbPPvss1q3bp1Wr16tzMxMSVJ8fLy++OKLgO28Xq+Sk5NDVisAoGvsGn4Rcc5PklauXKm33npLr776qj/4JCkjI0MHDhwI2Laurk7p6YYrhQMAHCsiwu+9997Txo0btWrVKl1xxRUBz40dO1YNDQ3+jzeUl5fryJEjys3NDUepAIAv4YIXA6tWrZLX69Udd9wRML5y5UplZmbq6aef1pIlS9Tc3KzExESVlpaqT58+4SkWAOBn18OeERF+L7zwQqfPZ2VladOmTRepGgBAV9k1/CLisCcAABdTRHR+AAB7smvnR/gBAIJG+AEAHMeuqzoQfgCAoNm18+OCFwCA49D5AQCCZtfOj/ADAASN8AMAOA7hBwBwHLuGn8uyrB610JXLZbb+Wc/QHe+B6QXMlxq+ej/D15f6yGzlj3661mh+dyy61aqPjeZ/pg+N5h/VYaP5ktQsr9H8M2o1rsFiPT+F6kf9EsOfN78I0/eGzg8AEDS7dn6EHwAgaIQfAMBx7Bp+fMgdAOA4dH4AgKDZtfMj/AAAQSP8AACOw6oOAADHsWvnxwUvAADHofMDAATNrp0f4QcACJpdw4/DngCAoMUYPjpTVVWlwsJCZWdnKzc3Vxs3bjzvdi+99JImTJigcePGqaCgQDt37rxg3XR+AICgharzO3XqlIqLizVv3jwVFBSorq5OU6ZMkcfjUUZGhn+7HTt26L/+67/02muvqW/fvtqyZYvuu+8+7dy5U3FxcR3un84PABBxKioqFBUVpYKCAklSWlqaRo8erU2bNgVsN2DAAK1cuVJ9+/aVJGVnZ+sf//iHDh/ufEUSOj8AQNBC1fnV19crJSUlYCw1NVU1NTUBY4MGDQr4euvWrbriiis0YMCATvdP+AEAghaq8PP5fIqNjQ0Yc7vd8vl8Hc6prKzU0qVLtXz5csXEdF4Z4dcjdcfikGcM558wmn1WHf8D76q/Gy4E+w/De1fUG80+xzL8PpxRi9H8s8b/DqQ2tRnNZyHayBaq8EtISJDXG7gQstfrVUJCwnm3f+ONN/TYY4/p8ccfV1ZW1gX3zzk/AEDESU9PV0NDg9ra/vXLU11dndLT09ttu379ej311FNau3Ztl4JPIvwAAAZC9VGHYcOGKS4uThs2bJAk1dTUqLKyUvn5+QHb1dbW6te//rXWrFmj6667rst1uyzL6lHHFFwuV7hL6CFM30ez36tc3XC73EsUe+GNOhFtWEMvo9nncNiTw57dJVQ/6g8Z/qzo38n3d+/evSopKdGxY8fkdrs1f/585ebmasWKFYqLi1NxcbEeeughbdy4UVdddVXA3J///OcaM2ZMh/sm/NABwo/wI/x6klD9qP/U8GfFVWH6/nLBCwAgaNzeDAAAm6DzAwAEza6dH+EHAAjaJWandMN2/JHwAwAELeqs4Q7C1Dpyzg8A4Dh0fgCAoNm18yP8AABBMw6/MCH8AABBc5nfByEsCD8AQPBOh7uA4HDBCwDAcej80AHT++21Gr666YeHpDOGv5KeNbxnYXfcZdb8vpZm87mvJi7Ipp0f4QcACB7n/AAAjkPnBwBwHJuGHxe8AAAch84PABA8m3Z+hB8AIHiEHwDAcWwafpzzAwA4Dp0fACB4Nu38CD8AQPAIPwCA49g0/CLunF9ZWZk8Ho8qKyv9Y+Xl5ZowYYJycnKUl5enXbt2hbFCAIDfacNHmERU53f8+HE98cQT6tOnj3/ss88+0+zZs1VaWqqsrCzt3r1b9957r7Zt26akpKTwFQsAsK2I6vwWLlyoH/7wh0pISPCPbd++XR6PR1lZWZKkm2++Wenp6dq2bVu4ygQA/JNNO7+ICb+ysjI1NzfrjjvuCBivr6/XwIEDA8ZSU1NVW1t7McsDAJzPGcNHmFy0w56bN2/WokWL2o337t1br7zyilasWKG1a9e2e97n88ntdgeMud1u+Xy+kNWKSNAd68iZrmUXzlcHbMKmF7xctPDLy8tTXl7eeZ+77777VFRUpH79+rV7Lj4+Xl988UXAmNfrVXJyckjqBAB8DYRfcJqbm1VRUaHq6mr95je/kSQdPXpUc+fO1d13362MjAytW7cuYE5dXZ3/HCAAAF9X2MPv0ksv1fvvvx8wduutt+o///M/NXz4cDU2NurRRx/Vrl27NGLECJWXl+vIkSPKzc0NU8UAAD86v9BISkrS008/rSVLlqi5uVmJiYkqLS0N+DgEACBMbBp+LsuyetR5eZfLFe4SACDihOxHfYnhz9xHwhNBEfNRBwAALpaIP+wJAIhgNj3sSfgBAIJH+AEAHIfwAwA4jk3DjwteAACOQ+cHAAheGG9ObYLwAwAEz6aHPQk/AEDwCD8AgOPYNPy44AUA4Dh0fgCA4Nm08yP8AADBI/wAAI5j0/DjnB8AwHHo/AAAwbNp50f4AQCCR/gBAByH8AMAOI5Nw48LXgAAjkPnBwAIHqs6AAAch8OeAADHOW346ERVVZUKCwuVnZ2t3Nxcbdy48bzbvfHGG8rNzVVOTo4KCwu1Z8+eC9dt9TCSePDgwYPHVx4hc6nMHh1oaWmxRo4cab3++uuWZVlWbW2tNXToUGvfvn0B21VXV1tDhgzxj7/xxhvWmDFjrFOnTnVaNp0fACDiVFRUKCoqSgUFBZKktLQ0jR49Wps2bQrYbvPmzcrOzlZGRoYkKT8/X21tbdq9e3en+yf8AADBC9Fhz/r6eqWkpASMpaamqra2NmCsrq5OAwcODBhLSUlpt91XccELACB4IbrgxefzKTY2NmDM7XbL5/MFjJ08eVJut/uC231Vjws/y7LCXQIAOEeIfuYmJCTI6/UGjHm9XiUkJASMxcfHq7m5+YLbfRWHPQEAESc9PV0NDQ1qa2vzj9XV1Sk9PT1gu4yMDB04cMD/dWtrqz766CP/OcCOEH4AgIgzbNgwxcXFacOGDZKkmpoaVVZWKj8/P2C7SZMmaefOndq3b58kaf369UpMTNSwYcM63b/L4jghACAC7d27VyUlJTp27Jjcbrfmz5+v3NxcrVixQnFxcSouLpYklZWV6fHHH9fp06d19dVXa/HixRo0aFCn+yb8AACOw2HPECgrK5PH41FlZaV/rLy8XBMmTFBOTo7y8vK0a9euMFYYGpZl6cknn9S4ceM0btw4TZkyRR988IH/eSe8B1LX70rR01RUVGjy5MkaP368xo0bp9WrV0uSGhsbNXPmTI0aNUpjxozRY4891uMvTGtqatKoUaN0//33S3LmexDxgvxMPzrw+eefW9/73vesm2++2aqoqLAsy7I+/fRT61vf+pa1a9cuy7Isq7Ky0hoyZIh1/PjxcJba7X77299aEydOtE6cOGFZlmU9//zz1tixYy3Lcs570NW7UvQ0x44dC/j+Hjx40Pr2t79tvf/++1ZxcbF1//33W21tbVZzc7OVn59vvfTSS2GuOLT+4z/+wxo7dqy1YMECy7IsR74HkY7Or5stXLhQP/zhDwMus92+fbs8Ho+ysrIkSTfffLPS09O1bdu2cJUZEjfddJMeffRRJSYmSpKys7N18OBBnTx50jHvQVfvStHTREVF6Ve/+pX/+9u/f3+lpaWpurpa77zzjoqKiuRyuZSQkKBp06bpzTffDHPFofP222/r0KFDmjhxoiSpubnZce+BHRB+3aisrEzNzc264447Asbr6+vb3YHgfHcqsLsbb7xRN9xwg//rrVu3avDgwYqLi3PMe9DVu1L0NElJScrNzfV/ffDgQe3fv1833HCDLMvSgAED/M915e4bdnXixAktXbpUS5cuVVTUuR+vDQ0NjnoP7KLHfcg9lDZv3qxFixa1G+/du7deeeUVrVixQmvXrm33vM/nC+oOBJGos/dg+/bt/q/Lysq0Zs0a/3mfnvQedKard6XoyY4ePap77rlHP/7xj+VyuRQdHe0PAkmKjY3tse/HkiVLdOeddyo1NdU/dvLkSUe9B3ZB+H0NeXl5ysvLO+9z9913n4qKitSvX792z8XHx+uLL74IGPN6vUpOTg5JnaHU2XvwT88++6zWrVun1atXKzMzU1LPeg8609W7UvRU1dXVmjlzpqZPn667775be/bsUUtLi86cOaPo6GhJ5y4G6Ynvx44dO3T48GEtW7YsYDw+Pt4x74GdEH7doLm5WRUVFaqurtZvfvMbSed++507d67uvvtuZWRkaN26dQFz6urq/OdHepKVK1eqvLxcr776qq644gr/uFPeg/T0dD377LNqa2vz/6ZfV1d3wbtN9ATV1dUqKirSL3/5S40dO1bSuUO+0dHROnDggP89ON9dOnqCsrIyNTQ0+P/uTU1Nam1t1d69ex3zHthKuK+46alycnL8V3seP37cGjJkiLVz507Lsizr3XfftYYNG2Z98cUXYayw+/3pT3+ybr31Vuvvf/97u+ec8h6cOnXKysnJsV599VXLsixrz5491pAhQ6z6+vowVxZaLS0t1ne/+11r27Zt7Z6bM2eOtWDBAqu1tdU6ceKENX78eOu1114LQ5UX15NPPum/2tOp70Eko/O7CJKSkvT0009ryZIlam5uVmJiokpLS9WnT59wl9atVq1aJa/X2+6Cn5UrVyozM9MR70FMTIxKS0tVUlKiZ555Rm63W8uWLWt3sU9Ps23bNh0+fFjLly/X8uXL/eN5eXl6+OGHVVJSojFjxigqKkq33367CgsLw1jtxcd7EHm4wwsAwHH4qAMAwHEIPwCA4xB+AADHIfwAAI5D+AEAHIfwAwA4DuEHGNi3b5+GDx+umpqabt/30aNHNXLkSP3pT3/q9n0DTkf4AUE6efKkZs6cqZkzZ+r666/v9v337dtXixcv1vz58/Xpp592+/4BJyP8gCD99re/VVtbm6ZOnRqy18jJyVFaWppKS0tD9hqAExF+gKQ//OEP+va3v62PP/7YP/bAAw/otttu05kzZ9ptb1mW1qxZoylTpvjv1L9+/XqNGDFCu3fv1qRJk3TTTTfp+9//vqqqqiRJZ8+elcfj0euvv64f/OAHuummmzR+/HhVVVXpd7/7nUaPHq0hQ4boZz/7mVpbW/2vNWPGDG3YsKHdahEAgkf4AZLy8/M1atQoLVy4UJJUWVmpzZs367HHHvOH25dVV1fr+PHjGjlypH+sV69eampq0rp167RmzRpVVlYqOTnZv89LLjl3K90XX3xRCxYs0O7du9W/f3/Nnj1bhw4d0tatW/XKK69oy5Ytevvtt/37zcrKUmtrq3bu3BnKtwBwFMIP+D8LFy5UTU2NXn/9dT300EOaO3eu0tLSzrvt/v375XK55PF4AsbPnDmj4uJiJSUlKTY2Vjk5Oe1W7L711luVmZkpt9ut7OxsHTt2TD/96U8VGxur9PR0eTwe1dfX+7e/7LLLdPXVV2vv3r3d/5cGHIrwA/7PN77xDS1evFgPPPCArrzySs2YMaPDbRsbG5WYmOjv5r5swIAB/j/HxsaqpaUl4PlrrrnG/+e4uDglJycHrHIfFxenU6dOBcxJSkpSY2Pj1/47ATg/wg/4koaGBsXHx+vIkSPy+XxB7eOfi9h29fkLbQ+g+/G/Dvg/9fX1WrlypZ577jn17dtXS5cu7XDbpKQkNTU16ezZsxeltsbGRiUlJV2U1wKcgPADJLW2tmrBggUqLCzU0KFDtXjxYm3atEk7duw47/Yej0eWZenDDz8MeW0nTpzQkSNH2p1fBBA8wg+Q9Oyzz+r48eOaO3euJCk1NVU/+clP9OCDD573XNsNN9yg5ORk7dq1K+S1/fnPf1avXr10yy23hPy1AKdgJXcgSM8995xeffVVbdmy5bwXvnSXadOmKS0tTY888kjIXgNwGjo/IEjTp0+XJK1bty5kr7Fjxw7V1tZq1qxZIXsNwIkIPyBIcXFxevrpp/XUU0+F5DN4R48e1YMPPqjly5frqquu6vb9A07GYU8AgOPQ+QEAHIfwAwA4DuEHAHAcwg8A4DiEHwDAcQg/AIDj/P+5ywZUvaP2pwAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "f.plane(\"z\").angle.mpl_scalar(lightness_field=f.z, filter_field=f.norm)" ] } ], "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.8.5" }, "widgets": { "application/vnd.jupyter.widget-state+json": { "state": {}, "version_major": 2, "version_minor": 0 } } }, "nbformat": 4, "nbformat_minor": 4 }