I have seen this kind implementations in Maxima for general relativity. So there must be examples of its applications.