Thanks for the replies. I wasn't expecting a response from what use to be usenet (last post here was 2021),
so I also posted an issue at github:
I was aware that David Lucas and David Joyner were both involved in grs_code.py and sent an email to both. I got
a reply from David Joyner, to post here at google groups. Although they added
GRSKeyEquationSyndromeDecoder
to Sage grs_code.py, I'm not sure who invented the algorithm or when. Doing a web search, it's mentioned
in other documents unrelated to Sage, such as:
The next reference dates back to 1992, but unless I buy the PDF, I can't be sure it is not a reference to
the BCH style encoding by a fixed generator polynomial, which these days use Berlekamp Massey or
Sugiyama decoders.
The Sage grs_code.py implements decoders for RS original encoding scheme, Berlekamp Welch (1986),
Gao (2002), and the syndrome decoder.
The other mystery is it didn't take much to add post correction code to handle a set of evaluation points that
includes zero, but I haven't seen this documented anywhere, so at this point it is a self-discovery.
Jeff