The paper introduces a new tree-level matrix element generator called COMIX, which is based on the color-dressed Berends-Giele recursive relations. The authors discuss two new algorithms for phase space integration, designed to handle large multiplicities and color sampling efficiently. The first algorithm is a general recursive approach for phase space generation, while the second is tailored for QCD and QCD-associated processes, utilizing the HAAG generator. The paper also covers the implementation details of COMIX, including multi-threading to reduce computation time and the choice of spinor basis. The effectiveness of the algorithms is demonstrated through comparisons with other matrix element generators and a detailed analysis of the phase space integration techniques.The paper introduces a new tree-level matrix element generator called COMIX, which is based on the color-dressed Berends-Giele recursive relations. The authors discuss two new algorithms for phase space integration, designed to handle large multiplicities and color sampling efficiently. The first algorithm is a general recursive approach for phase space generation, while the second is tailored for QCD and QCD-associated processes, utilizing the HAAG generator. The paper also covers the implementation details of COMIX, including multi-threading to reduce computation time and the choice of spinor basis. The effectiveness of the algorithms is demonstrated through comparisons with other matrix element generators and a detailed analysis of the phase space integration techniques.