We already have our own triangular solver kernels, extending them to a Gauß-Seidel kernel should be straightforward.