Thanks for the question. Yes the example code you linked to (involving the BondGate class) ought to work well for non-Hermitian operators. I think the Trotter expansion itself would still be justified in this case and the part of the code which does the exponentiation, which is inside the BondGate class itself, uses a general enough algorithm that it should handle exponentiating any type of operator. (Internally, BondGate sums a very high order Taylor series to exponentiate the local operator.)
Regarding part (2) of your question, there is a built-in function to exponentiate local Hermitian operators, which is called expHermitian, but only for that case. However, as mentioned above the code you linked to uses the BondGate class which does the exponentiation internally using a different (high-order Taylor) algorithm. So that may be all you need for your use case.
But we have not written tests for the non-Hermitian case, so I'd say please proceed with some case. Please test some exactly solvable or known limits of your problem such as small system sizes where you could do everything exactly with matrix methods. This is a good practice anyway even when using fully supported and tested cases such as Hermitian time evolution, because every new physical system presents its own challenges.
Please let us know if you don't get the results you expect when doing those initial checks.