214   Matrix<Scalar,4,4> coeffMatrix = Matrix<Scalar,4,4>::Zero();  in solveAuxiliaryEquation()  local
215   coeffMatrix.coeffRef(0,0) = A.coeff(0,0) + B.coeff(0,0);  in solveAuxiliaryEquation()
216   coeffMatrix.coeffRef(1,1) = A.coeff(0,0) + B.coeff(1,1);  in solveAuxiliaryEquation()
217   coeffMatrix.coeffRef(2,2) = A.coeff(1,1) + B.coeff(0,0);  in solveAuxiliaryEquation()
218   coeffMatrix.coeffRef(3,3) = A.coeff(1,1) + B.coeff(1,1);  in solveAuxiliaryEquation()
219   coeffMatrix.coeffRef(0,1) = B.coeff(1,0);  in solveAuxiliaryEquation()
220   coeffMatrix.coeffRef(0,2) = A.coeff(0,1);  in solveAuxiliaryEquation()
221   coeffMatrix.coeffRef(1,0) = B.coeff(0,1);  in solveAuxiliaryEquation()
222   coeffMatrix.coeffRef(1,3) = A.coeff(0,1);  in solveAuxiliaryEquation()
223   coeffMatrix.coeffRef(2,0) = A.coeff(1,0);  in solveAuxiliaryEquation()
224   coeffMatrix.coeffRef(2,3) = B.coeff(1,0);  in solveAuxiliaryEquation()
225   coeffMatrix.coeffRef(3,1) = A.coeff(1,0);  in solveAuxiliaryEquation()
226   coeffMatrix.coeffRef(3,2) = B.coeff(0,1);  in solveAuxiliaryEquation()
235   result = coeffMatrix.fullPivLu().solve(rhs);  in solveAuxiliaryEquation()