Add method to compute Bt*C*B faster while calculating element stiffness matrix