Microscopic dynamics in the cellular cytoskeleton must be carefully understood to predict structural and hydrodynamic behavior at the scale of the whole cell. The stochastic nature of sub-cellular processes coupled with the ever-changing topology of cytoskeletal networks make the microscopic dynamics of a cell difficult or impossible to interrogate through lab measurements or existing simulation techniques. In this talk I’ll discuss efficient methods to numerically simulate the types of flexible, fluctuating, inextensible microfilament gels that make up the cytoskeleton. I’ll introduce a method in which fibers are represented as curves on the unit sphere so that strict inextensibility is maintained as they evolve, and go on to describe how the fluctuating hydrodynamics of a fiber suspension can be captured with a carefully designed temporal integration scheme. I’ll briefly describe how fast linear algebra can scale simulations to thousands of fibers and demonstrate the efficacy of this approach through a few examples.