diff --git a/dune/gfe/rotation.hh b/dune/gfe/rotation.hh index aa023750e5188f60422f3d30f516b97b443e5cec..e767e14776e49837220842138d81b1cdde1876fd 100644 --- a/dune/gfe/rotation.hh +++ b/dune/gfe/rotation.hh @@ -632,17 +632,17 @@ public: { const Quaternion<T>& q = (*this); - dd_dq[0][0][0] = 2*q[0]; dd_dq[0][0][1] = -2*q[1]; dd_dq[0][0][2] = -2*q[2]; dd_dq[0][0][3] = 2*q[3]; - dd_dq[0][1][0] = 2*q[1]; dd_dq[0][1][1] = 2*q[0]; dd_dq[0][1][2] = 2*q[3]; dd_dq[0][1][3] = 2*q[2]; - dd_dq[0][2][0] = 2*q[2]; dd_dq[0][2][1] = -2*q[3]; dd_dq[0][2][2] = 2*q[0]; dd_dq[0][2][3] = -2*q[1]; + dd_dq[0][0] = { 2*q[0], -2*q[1], -2*q[2], 2*q[3]}; + dd_dq[0][1] = { 2*q[1], 2*q[0], 2*q[3], 2*q[2]}; + dd_dq[0][2] = { 2*q[2], -2*q[3], 2*q[0], -2*q[1]}; - dd_dq[1][0][0] = 2*q[1]; dd_dq[1][0][1] = 2*q[0]; dd_dq[1][0][2] = -2*q[3]; dd_dq[1][0][3] = -2*q[2]; - dd_dq[1][1][0] = -2*q[0]; dd_dq[1][1][1] = 2*q[1]; dd_dq[1][1][2] = -2*q[2]; dd_dq[1][1][3] = 2*q[3]; - dd_dq[1][2][0] = 2*q[3]; dd_dq[1][2][1] = 2*q[2]; dd_dq[1][2][2] = 2*q[1]; dd_dq[1][2][3] = 2*q[0]; + dd_dq[1][0] = { 2*q[1], 2*q[0], -2*q[3], -2*q[2]}; + dd_dq[1][1] = {-2*q[0], 2*q[1], -2*q[2], 2*q[3]}; + dd_dq[1][2] = { 2*q[3], 2*q[2], 2*q[1], 2*q[0]}; - dd_dq[2][0][0] = 2*q[2]; dd_dq[2][0][1] = 2*q[3]; dd_dq[2][0][2] = 2*q[0]; dd_dq[2][0][3] = 2*q[1]; - dd_dq[2][1][0] = -2*q[3]; dd_dq[2][1][1] = 2*q[2]; dd_dq[2][1][2] = 2*q[1]; dd_dq[2][1][3] = -2*q[0]; - dd_dq[2][2][0] = -2*q[0]; dd_dq[2][2][1] = -2*q[1]; dd_dq[2][2][2] = 2*q[2]; dd_dq[2][2][3] = 2*q[3]; + dd_dq[2][0] = { 2*q[2], 2*q[3], 2*q[0], 2*q[1]}; + dd_dq[2][1] = {-2*q[3], 2*q[2], 2*q[1], -2*q[0]}; + dd_dq[2][2] = {-2*q[0], -2*q[1], 2*q[2], 2*q[3]}; }