do_scf_relax : if .TRUE. the code makes a self-consistent relax calculation at
the equilibrium geometry to find the optimized atomic
coordinates, before computing the elastic constants.
This step is needed only for solids that have internal
degrees of freedom in the unstrained configuration.
If .FALSE. the coordinates of the input geometry are strained
uniformly to the equilibrium geometry.
Default: logical .FALSE.
Number of tasks for this option: Maximum between the number of tasks
needed by the what='mur_lc' option and the number of tasks
needed for the what='scf_elastic_constants' option.