Generation of non-Gaussian tensor-valued random fields using an ISDE-based algorithm