Complex fluids are challenging to simulate as they are characterized by stationary or evolving microstructure. We will discuss the discretization of integral equation formulations for complex fluids. In particular, we will discuss the formulation, numerical challenges and scalability of algorithms for volume integral equations and we will present a new open-source library for such problems. We will compare their performance to other state-of-the art codes.