Overview

In this lecture we will first optimize our Stokes solver by implementing a blocked matrix assembly scheme. This will allow us to solve larger problems more efficiently. We will then apply our Stokes solver to simulate viscous folding instabilities in a 2D box setup.