I was just tweaking a plasma simulation in Python, and I hit a performance wall—your thoughts on optimizing the numerical solver?