210   double theta = 1.0e-2;  in TEST()  local
211   double axis_angle[3] = { theta, 0, 0 };  in TEST()
213   double expected[4] = { cos(theta/2), sin(theta/2.0), 0, 0 };  in TEST()
222   double theta = pow(numeric_limits<double>::min(), 0.75);  in TEST()  local
223   double axis_angle[3] = { theta, 0, 0 };  in TEST()
225   double expected[4] = { cos(theta/2), sin(theta/2.0), 0, 0 };  in TEST()
272   double theta = 1.0e-2;  in TEST()  local
273   double quaternion[4] = { cos(theta/2), sin(theta/2.0), 0, 0 };  in TEST()
275   double expected[3] = { theta, 0, 0 };  in TEST()
283   double theta = pow(numeric_limits<double>::min(), 0.75);  in TEST()  local
284   double quaternion[4] = { cos(theta/2), sin(theta/2.0), 0, 0 };  in TEST()
286   double expected[3] = { theta, 0, 0 };  in TEST()
326     double theta = kPi * 2 * RandDouble() - kPi;  in TEST()  local
328       axis_angle[i] = axis_angle[i] * theta / norm;  in TEST()
436     double theta = kPi - kMaxSmallAngle * RandDouble();  in TEST()  local
439       in_axis_angle[i] *= (theta / norm);  in TEST()
534     double theta = kPi * 2 * RandDouble() - kPi;  in TEST()  local
536       axis_angle[i] = axis_angle[i] * theta / norm;  in TEST()
695     double theta = pow(10.0, i);  in TEST()  local
696     J3 axis_angle[3] = { J3(theta, 0), J3(0, 1), J3(0, 2) };  in TEST()
699         MakeJ3(cos(theta/2), -sin(theta/2)/2, 0, 0),  in TEST()
700         MakeJ3(sin(theta/2), cos(theta/2)/2, 0, 0),  in TEST()
701         MakeJ3(0, 0, sin(theta/2)/theta, 0),  in TEST()
702         MakeJ3(0, 0, 0, sin(theta/2)/theta),  in TEST()
715     double theta = pow(10.0, i);  in TEST()  local
716     J3 axis_angle[3] = { J3(theta, 0), J3(0, 1), J3(0, 2) };  in TEST()
751     double theta = pow(10.0, i);  in TEST()  local
752     double s = sin(theta);  in TEST()
753     double c = cos(theta);  in TEST()
757         MakeJ4(s, -2*theta, 2*theta*c, 0, 0),  in TEST()
758         MakeJ4(0, 0, 0, 2*theta/s, 0),  in TEST()
759         MakeJ4(0, 0, 0, 0, 2*theta/s),  in TEST()
771     double theta = pow(10.0, i);  in TEST()  local
772     double s = sin(theta);  in TEST()
773     double c = cos(theta);  in TEST()
780         MakeJ4(theta, -2*theta, 2.0, 0, 0),  in TEST()
894     double theta = (2.0 * i * 0.0011 - 1.0) * kPi;  in TEST()  local
903       const double inv_norm = theta / sqrt(norm2);  in TEST()
943     double theta = (2.0 * i * 0.0001  - 1.0) * 1e-16;  in TEST()  local
944     const double inv_norm = theta / sqrt(norm2);  in TEST()