From: Michael Orlitzky <michael@orlitzky.com>
Date: Fri, 15 Nov 2019 13:15:53 +0000 (-0500)
Subject: eja: use peirce_decomposition() method to get the Peirce decomposition.
X-Git-Url: https://gitweb.michael.orlitzky.com/?a=commitdiff_plain;h=d4dda0de8204273c70b599bfa532e54563cf0af1;p=sage.d.git

eja: use peirce_decomposition() method to get the Peirce decomposition.
---

diff --git a/mjo/eja/eja_element.py b/mjo/eja/eja_element.py
index f44bee1..4d3b071 100644
--- a/mjo/eja/eja_element.py
+++ b/mjo/eja/eja_element.py
@@ -620,19 +620,10 @@ class FiniteDimensionalEuclideanJordanAlgebraElement(IndexedFreeModuleElement):
             True
 
         """
-        # TODO: when the Peirce decomposition is implemented for real,
-        # we can use that instead of finding this eigenspace manually.
-        #
-        # Trivial eigenspaces don't appear in the list, so we default to the
-        # trivial one and override it if there's a nontrivial space in the
-        # list.
         if not self.is_idempotent():
             return False
 
-        J1 = VectorSpace(self.parent().base_ring(), 0)
-        for (eigval, eigspace) in self.operator().matrix().left_eigenspaces():
-            if eigval == 1:
-                J1 = eigspace
+        (_,_,J1) = self.parent().peirce_decomposition(self)
         return (J1.dimension() == 1)