There doesn't seem to be a good way to accomplish this directly in linalg.generic. One possible solution is to build the structure of the output in numpy and use values from one of the inputs. Then perform an operation in linalg.generic, computing the value based on the other input using modulo indexing.