We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
1 parent 4484536 commit 8d0b9ebCopy full SHA for 8d0b9eb
source/module_ri/module_exx_symmetry/symmetry_rotation.cpp
@@ -210,7 +210,7 @@ namespace ModuleSymmetry
210
// gmatc should be a rotation matrix, i.e. det(gmatc)=1
211
TCdouble Symmetry_rotation::get_euler_angle(const ModuleBase::Matrix3& gmatc) const
212
{
213
- double threshold = 1e-8;
+ double threshold = this->eps_;
214
double alpha, beta, gamma;
215
if (std::fabs(gmatc.e32) > threshold || std::fabs(gmatc.e31) > threshold) // sin(beta) is not zero
216
0 commit comments