Matrix3 Matrix3::Diag(double a11, double a22, double a33) { Matrix3 m; m.a11() = a11; m.a22() = a22; m.a33() = a33; return m; }
Matrix3 Matrix3::Eye() { Matrix3 m; m.a11() = m.a22() = m.a33() = 1; return m; }