More robust atomic strain calculation (avoid singular matrix during matrix inverse by using dsylsd). Added test.