From b5c6a6c6db3d12ff2a33663d426af384a6931c13 Mon Sep 17 00:00:00 2001
From: utensil <utensilcandel@gmail.com>
Date: Fri, 10 May 2024 13:52:52 +0800
Subject: [PATCH] Fixes #516

---
 examples/ipython/issue-516.ipynb | 59 ++++++++++++++++++++++++++++++++
 galgebra/mv.py                   |  1 +
 2 files changed, 60 insertions(+)
 create mode 100644 examples/ipython/issue-516.ipynb

diff --git a/examples/ipython/issue-516.ipynb b/examples/ipython/issue-516.ipynb
new file mode 100644
index 00000000..5360e912
--- /dev/null
+++ b/examples/ipython/issue-516.ipynb
@@ -0,0 +1,59 @@
+{
+ "cells": [
+  {
+   "cell_type": "code",
+   "execution_count": 1,
+   "id": "ee7c8882-7356-441a-920e-3603d081002b",
+   "metadata": {},
+   "outputs": [
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "[1]\n",
+      "[1]\n"
+     ]
+    }
+   ],
+   "source": [
+    "from galgebra.ga import Ga\n",
+    "import sympy as S\n",
+    "\n",
+    "ga = Ga('e', g=[1,1,1], coords=S.symbols(f\"0:3\", real=True), wedge=False)\n",
+    "x,y,z = ga.mv()\n",
+    "a = ga.mv('A','mv')\n",
+    "print((a.grade(0)*x).grades)\n",
+    "print((a.grade(0)^x).grades)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "8920d631-eb4f-4955-9558-9dc158b5171e",
+   "metadata": {},
+   "outputs": [],
+   "source": []
+  }
+ ],
+ "metadata": {
+  "kernelspec": {
+   "display_name": "Python 3 (ipykernel)",
+   "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.18"
+  }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 5
+}
diff --git a/galgebra/mv.py b/galgebra/mv.py
index ddf5933c..32c2dd2d 100644
--- a/galgebra/mv.py
+++ b/galgebra/mv.py
@@ -368,6 +368,7 @@ def __init__(self, *args, ga: 'Ga', recp=None, coords=None, **kwargs):
                 self.obj = x.obj
                 self.is_blade_rep = x.is_blade_rep
                 self.i_grade = x.i_grade
+                self.characterise_Mv()
             else:
                 if isinstance(x, Expr):  # copy constructor for obj expression
                     self.obj = x