Normal vector computation, local stiffness matrix computation and COO format of the stiffness matrix were implemented