Implement some functions to compute Bulk modulus (imcomplete)
Trying to impose constraints algebraically.
Algebraic constraints implemented. Currently, only the constraints (self) are considered. This limitation may be fixed in a future release. (Inclusion of const (cross) for rotational invariance).
Merge branch 'feature/algebraic_constraint' into develop