Abstract: Feynman integrals beyond multiple polylogarithms have been of great interest recently. In this talk, I will discuss how to obtain the $\varepsilon$-form systematically for two non-planar triangle integrals at two loops, which have non-trivial sub-topology dependences. The methodology depends on the essential geometry information related to the Feynman integrals and can be generalized to more cases.