tjfulle opened this issue 3 years ago · 0 comments
symmetric3x3::isotropic_part has auto d = trace(a) / T(3.0) which results in a unitless quantity. It should be simply auto d = trace(a) / 3.0 or similar
symmetric3x3::isotropic_part
auto d = trace(a) / T(3.0)
auto d = trace(a) / 3.0